版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

jesseliuxl

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 3820.5
散金: 301
紅花: 1
帖子: 192
在線: 105.7小時
蟲號: 1174887
注冊: 2010-12-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

[求助] 100個金幣：udf改成并行的，求高人

本例是通過出口的濃度與進(jìn)口關(guān)聯(lián)，串行沒問題，并行就一直沒搞通，求高人

CODE:

#include "udf.h"

#include "mem.h" 

#define q 998

#define out 960

#define film-out 28

#define INT_YI 0.03

real AVG_YI;

DEFINE_ADJUST(outlet_avg_yi, d)

{

 real sum_yi_A=0.0;

 int j=0;

 face_t f; 

 Thread *thread_out=Lookup_Thread(d,5);

  begin_f_loop(f,thread_out)

 {

  sum_yi_A+=F_YI(f,thread_out,0); 

  j++;       

 }

 end_f_loop(f,thread_out)

  AVG_YI=sum_yi_A/j;

printf("%f\n",sum_yi_A);/* write out kount to data file */

printf("%f\n",AVG_YI); 

printf("%d\n",j);

} 

DEFINE_PROFILE(inlet,thread_in,index)

{

    real inlet_yi=0.0;

    face_t inlet_face;

    real time=RP_Get_Real("flow-time");   

    begin_f_loop(inlet_face,thread_in)

        {

          if(time<9)

        {

         inlet_yi=INT_YI;

          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;

        }

          else 

        {

          inlet_yi=AVG_YI+0.2*out/q;

          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;

        }

      }

        end_f_loop(inlet_face,thread_in)   

     printf("%f\n",inlet_yi);                      

     }

[ Last edited by xiegangmai on 2012-4-11 at 22:36 ]

回復(fù)此樓

» 猜你喜歡

一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有3人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有8人回復(fù)
材料學(xué)碩297已過四六級求調(diào)劑推薦已經(jīng)有4人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有14人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有8人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有5人回復(fù)
271材料工程求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有8人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

關(guān)于一個高級酸和高級醇的酯化反應(yīng)，各種求助，有金幣伺候！已經(jīng)有17人回復(fù)
調(diào)劑功能納米研究院怎么樣，望高人指點(diǎn)。。。已經(jīng)有24人回復(fù)
沒金幣求助，菌體在液體培養(yǎng)基中聚集成團(tuán)絮狀已經(jīng)有21人回復(fù)
求助：怎樣由鄰硝基甲苯合成2-硝基-4-溴苯甲酸采納答案贈金幣已經(jīng)有3人回復(fù)
【求教】有哪位高人做過長鏈伯胺與環(huán)氧基團(tuán)的固化反應(yīng)？已經(jīng)有23人回復(fù)
請教各位牛人：合成水性聚氨酯失敗原因（內(nèi)附配方）已經(jīng)有9人回復(fù)
請高人指點(diǎn)，如果CSC系統(tǒng)外語選不達(dá)標(biāo)可否已經(jīng)有14人回復(fù)
電腦悲催了在線等請求高人支招！！多謝多謝！已經(jīng)有22人回復(fù)
公司域用戶，無管理員權(quán)限，USB接口禁用、注冊表禁用如何解禁？如解決再送100金幣！已經(jīng)有7人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
y470電池模式改成最佳保養(yǎng)。一重啟就恢復(fù)成最長續(xù)航了�。。� 已經(jīng)有4人回復(fù)
【求助】單機(jī)多核UDF怎么并行計算已經(jīng)有15人回復(fù)
【求助】50金幣求幫寫一個并行的UDF 已經(jīng)有4人回復(fù)
【求助】物料屬性可以用UDF修改嗎已經(jīng)有8人回復(fù)
【求助】100金幣求助高人一個程序解方程已經(jīng)有15人回復(fù)

生命不息折騰不止

1樓 2012-04-01 09:06:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

jesseliuxl

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 3820.5
散金: 301
紅花: 1
帖子: 192
在線: 105.7小時
蟲號: 1174887
注冊: 2010-12-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

★ ★
jesseliuxl: 回帖置頂 2012-04-06 09:45:48
xiegangmai: 金幣+2, 鼓勵分享解決經(jīng)驗(yàn)，金幣獎勵！ 2012-04-06 10:11:20
xiegangmai: 編輯內(nèi)容 2012-04-06 10:11

并行的udf

CODE:

#include "udf.h"

#include "mem.h"  

#define q 989.57

#define film_out 72.59

#define cir_out 916.98

#define INT_YI 0.022

real AVG_YI=0.0;

DEFINE_ADJUST(out, d)

     {

      real sum_yi_A=0.0;

      real sum_flux_A=0.0;

      int j=0;

   #if !RP_HOST

      face_t f; 

      Thread *thread_out;

   #endif

   #if !RP_HOST

      thread_out=Lookup_Thread(d,19);

      begin_f_loop(f,thread_out)

       if (PRINCIPAL_FACE_P(f,thread_out))

       {

         sum_yi_A+=F_YI(f,thread_out,0)*F_FLUX(f,thread_out); 

         sum_flux_A+=F_FLUX(f,thread_out);       

       }

      end_f_loop(f,thread_out)

   #endif  

     #if RP_NODE

      sum_yi_A=PRF_GRSUM1(sum_yi_A);

      sum_flux_A=PRF_GRSUM1(sum_flux_A);

     #endif

      node_to_host_real_2(sum_yi_A,sum_flux_A);  

      AVG_YI=sum_yi_A/sum_flux_A;

}

DEFINE_PROFILE(in,thread_in,index)

   {

    real inlet_yi=0.0;

    face_t inlet_face;

    #if !RP_HOST 

    real time=RP_Get_Real("flow-time");

    begin_f_loop(inlet_face,thread_in)

     if (PRINCIPAL_FACE_P(inlet_face,thread_in))  

       {

          if(time<9)

           {

             inlet_yi=INT_YI;

             F_PROFILE(inlet_face,thread_in,index)=inlet_yi;

           }

          else 

           {

             inlet_yi=AVG_YI*cir_out/q+INT_YI;

             F_PROFILE(inlet_face,thread_in,index)=inlet_yi;

           }

        }

    end_f_loop(inlet_face,thread_in)   

    #endif

    }

[ Last edited by xiegangmai on 2012-4-6 at 10:11 ]

贊一下(7人)

回復(fù)此樓

生命不息折騰不止

13樓2012-04-06 09:45:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

5超7

至尊木蟲 (著名寫手)

菜鳥

仿真EPI: 1
應(yīng)助: 488 (碩士)
金幣: 25690.5
散金: 41
紅花: 41
帖子: 2465
在線: 880.2小時
蟲號: 1371942
注冊: 2011-08-18
性別: GG
專業(yè): 能源化工

并行這東西在流體網(wǎng)上面有所討論的把！浪子大哥好像開過帖子

贊一下(2人)

回復(fù)此樓

我要奮斗到底

8樓2012-04-02 08:08:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yp1990

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 219.2
散金: 17
紅花: 2
帖子: 44
在線: 47.8小時
蟲號: 2007436
注冊: 2012-09-17
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

13樓: Originally posted by jesseliuxl at 2012-04-06 09:45:42
并行的udf

#include "udf.h"
#include "mem.h"
#define q 989.57
#define film_out 72.59
#define cir_out 916.98

#define INT_YI 0.022
real AVG_YI=0.0;

DEFINE_ADJUST(out, d)
   {
   ...

你好，我按照你的模型也修改了我的udf，但是單機(jī)多核并行的時候還是出錯，不知道怎么回事？

贊一下(1人)

回復(fù)此樓

愛拼才會贏

20樓2015-05-03 21:51:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

zxbok1979

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 1858.1
散金: 25
紅花: 2
沙發(fā): 3
帖子: 172
在線: 56.7小時
蟲號: 1489252
注冊: 2011-11-13
性別: GG
專業(yè): 流體力學(xué)

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
jesseliuxl: 金幣+2, ★有幫助, 要是能這樣，還折騰個啥啊 2012-04-01 13:45:51
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-02 15:41:08

在單核求解器中編譯UDF，保存dat和case，然后再用并行求解器打開就可以了。

贊一下

回復(fù)此樓

SILENCEISALSOAWAYTOLEARN

2樓2012-04-01 12:02:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17462
散金: 883
紅花: 18
帖子: 6207
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
jesseliuxl: 金幣+10, ★有幫助, 能編譯通過，但是編譯后出現(xiàn)錯誤啦 2012-04-01 18:41:20
臭水溝: 金幣+2, 謝謝交流~~ 2012-04-02 15:41:20

#include "udf.h"
#include "mem.h"
#define q 998
#define out 960
#define film-out 28
#define INT_YI 0.03
real AVG_YI=0.0;

DEFINE_ADJUST(outlet_avg_yi, d)
{#if RP_NODE

real sum_yi_A=0.0;
int j=0;
face_t f;
Thread *thread_out=Lookup_Thread(d,5);
begin_f_loop(f,thread_out)
{
sum_yi_A+=F_YI(f,thread_out,0);
  j++;
}
end_f_loop(f,thread_out)
sum_yi_A = PRF_GRSUM1(sum_yi_A);
j = PRF_GRSUM1(j);
  AVG_YI=sum_yi_A/j;

Message0("%f\n",sum_yi_A);
Message0("%f\n",AVG_YI);
Message0("%d\n",j);
#endif
}

DEFINE_PROFILE(inlet,thread_in,index)
{#if RP_NODE
   real inlet_yi=0.0;
   face_t inlet_face;

   real time=RP_Get_Real("flow-time" ) ;
   begin_f_loop(inlet_face,thread_in)
      {
         if(time<9)
      {
      inlet_yi=INT_YI;
         F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
      }

         else
      {
         inlet_yi=AVG_YI+0.2*out/q;
         F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
      }
   }
      end_f_loop(inlet_face,thread_in)

   Message0("%f\n",inlet_yi);
   #endif
   }

贊一下

回復(fù)此樓

3樓2012-04-01 16:20:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jesseliuxl

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 3820.5
散金: 301
紅花: 1
帖子: 192
在線: 105.7小時
蟲號: 1174887
注冊: 2010-12-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

3樓: Originally posted by kylafree at 2012-04-01 16:20:13:
#include "udf.h"
#include "mem.h"
#define q 998
#define out 960
#define film-out 28
#define INT_YI 0.03
real AVG_YI=0.0;

DEFINE_ADJUST(outlet_avg_yi, d)
{#if R ...

編譯能通過，不過之后出現(xiàn)這問題，不知是什么原因

你并行里試過沒有�。�

贊一下

回復(fù)此樓

生命不息折騰不止

4樓2012-04-01 18:44:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jesseliuxl

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 3820.5
散金: 301
紅花: 1
帖子: 192
在線: 105.7小時
蟲號: 1174887
注冊: 2010-12-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

并行的udf只能編譯？

贊一下

回復(fù)此樓

生命不息折騰不止

5樓2012-04-01 22:32:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17462
散金: 883
紅花: 18
帖子: 6207
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

★
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-02 15:41:40

#include "udf.h"
#include "mem.h"
#define q 998
#define out 960
#define film_out 28
#define INT_YI 0.03
real AVG_YI=0.0;

DEFINE_ADJUST(outlet_avg_yi, d)
{#if RP_NODE

real sum_yi_A=0.0;
int j=0;
face_t f;
Thread *thread_out=Lookup_Thread(d,5);
begin_f_loop(f,thread_out)
{
sum_yi_A+=F_YI(f,thread_out,0);
  j++;
}
end_f_loop(f,thread_out)
sum_yi_A = PRF_GRSUM1(sum_yi_A);
j = PRF_GRSUM1(j);
AVG_YI=sum_yi_A/j;

Message0("%f\n",sum_yi_A);
Message0("%f\n",AVG_YI);
Message0("%d\n",j);
#endif
}

DEFINE_PROFILE(inlet,thread_in,index)
{#if RP_NODE
   real inlet_yi=0.0;
   face_t inlet_face;

   real time=RP_Get_Real("flow-time" ) ;
   begin_f_loop(inlet_face,thread_in)
      {
         if(time<9)
      {
         inlet_yi=INT_YI;
         F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
      }

         else
      {
         inlet_yi=AVG_YI+0.2*out/q;
         F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
      }
      }
      end_f_loop(inlet_face,thread_in)

   Message0("%f\n",inlet_yi);
   #endif
      }

贊一下

回復(fù)此樓

6樓2012-04-01 22:55:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kylafree

至尊木蟲 (知名作家)

應(yīng)助: 166 (高中生)
金幣: 17462
散金: 883
紅花: 18
帖子: 6207
在線: 1161.7小時
蟲號: 1061079
注冊: 2010-07-20
性別: GG
專業(yè): 流體力學(xué)

你連你的東西是干什么的都沒說你讓別人怎么試

贊一下

回復(fù)此樓

7樓2012-04-01 22:57:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jesseliuxl

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 3820.5
散金: 301
紅花: 1
帖子: 192
在線: 105.7小時
蟲號: 1174887
注冊: 2010-12-22
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

jesseliuxl: 回帖置頂 2012-04-02 08:59:39

引用回帖:

7樓: Originally posted by kylafree at 2012-04-01 22:57:27:
你連你的東西是干什么的都沒說你讓別人怎么試

組分輸送，把進(jìn)口組分濃度與出口關(guān)聯(lián)，你這高手應(yīng)該一看就知道啊

贊一下

回復(fù)此樓

生命不息折騰不止

9樓2012-04-02 08:32:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dhw0075

木蟲 (知名作家)

至尊木蟲

應(yīng)助: 58 (初中生)
金幣: 6561.1
散金: 354
紅花: 16
帖子: 5576
在線: 284.4小時
蟲號: 1107824
注冊: 2010-09-25
性別: GG
專業(yè): 防災(zāi)工程

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-02 18:03:36
jesseliuxl: 金幣+2 2012-04-02 19:57:20

在單核求解器中編譯UDF，保存dat和case，然后再用并行求解器打開就可以了。

祝福~

贊一下

回復(fù)此樓

學(xué)無止境。

10樓2012-04-02 16:40:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 jesseliuxl 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料080500調(diào)劑求收留 +5	一顆meteor 2026-03-13	5/250	2026-03-20 09:14 by xingguangj
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +6	Doleres 2026-03-19	6/300	2026-03-20 09:03 by xingguangj
[考研] 279分求調(diào)劑一志愿211 +7	chaojifeixia 2026-03-19	7/350	2026-03-20 09:00 by sunny300
[考研] 085700資源與環(huán)境308求調(diào)劑 +8	墨墨漠 2026-03-18	8/400	2026-03-20 08:59 by lbsjt
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	7/350	2026-03-20 08:53 by 申子申申
[考研] 281求調(diào)劑（0805） +11	煙汐憶海 2026-03-16	21/1050	2026-03-20 08:45 by 無際的草原
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 288求調(diào)劑 +15	于海海海海 2026-03-19	15/750	2026-03-19 22:41 by 學(xué)員8dgXkO
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	10/500	2026-03-19 16:05 by 余麟余
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu

24小時熱門版塊排行榜

[求助] 100個金幣：udf改成并行的，求高人

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 100個金幣：udf改成并行的，求高人

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: