关于在DEFINE_DPM_OUTPUT中使用Cell Macros的问题

Estufa

请问坛友有没有在DEFINE_DPM_OUTPUT中使用Cell macros的经历？我在里面使用C_STRAIN_RATE_MAG(c,t）时遇到了问题。

我用TP_CELL(tp)取的c（只用这个宏时是没有问题的），t在DEFINE_DPM_OUTPUT的输入形参里我就直接用了。然后Console报错Process xxxx: Received signal SIGSEGV.并且fluent停止工作。在想是不是DEFINE_DPM_OUTPUT形参表里的Thread* t和C_STRAIN_RATE_MAG所要的Thread* t并不一样？

下附出现问题的代码。将与s有关的两行（第28行，第17行）注释掉后挂载使用是没有问题的。

1. #include "udf.h"
   
2. #include "dpm_types.h"
   
3. #include "mem.h"
   
4. /******************************************************************/
   
5. /* UDF that samples discrete phase size and velocity distributions*/
   
6. /* within the domain.            */
   
7. /******************************************************************/
   
8. #define REMOVE_PARTICLES FALSE
   
9. 
   
10. DEFINE_DPM_OUTPUT(hemocellOutput,header,fp,tp,t,plane)
    
11. {
    
12.          #if RP_2D
    
13.                 return;
    
14.          #else
    
15.                  
    
16.     cell_t c;
    
17.         real s;
    
18.     if(header)
    
19.     {
    
20.       par_fprintf_head(fp," ID, time, x-velocity[m/s], ");
    
21.       par_fprintf_head(fp,"y-velocity[m/s], z-velocity[m/s], ");
    
22.       par_fprintf_head(fp,"x-position[m], y-position[m], ");
    
23.       par_fprintf_head(fp,"z-position[m], strain rate \n");
    
24.     }
    
25.     if(NULLP(tp))
    
26.       return;
    
27.         c = TP_CELL(tp);
    
28.         s = C_STRAIN_RATE_MAG(c,t);
    
29.         par_fprintf(fp,"%d %" int64_fmt " %" int64_fmt ", %e, %f, %f, %f, %f, %f, %f,  \n",
    
30.    P_INJ_ID(TP_INJECTION(tp)), TP_ID(tp), TP_ID(tp), TP_TIME(tp),TP_VEL(tp)[0],
    
31.   TP_VEL(tp)[1],TP_VEL(tp)[2],TP_POS(tp)[0],TP_POS(tp)[1],TP_POS(tp)[2]);
    
32.          #endif
    
33. 
    
34.    #if REMOVE_PARTICLES
    
35.       MARK_TP(tp, P_FL_REMOVED);
    
36.    #endif
    
37. }

复制代码

Estufa

已解决，此tread非彼thread

1. #include "udf.h"
   
2. #include "dpm_types.h"
   
3. #include "mem.h"
   
4. /******************************************************************/
   
5. /* UDF that samples discrete phase size and velocity distributions*/
   
6. /* within the domain.            */
   
7. /******************************************************************/
   
8. #define REMOVE_PARTICLES FALSE
   
9. Domain *domain1;
   
10. 
    
11. DEFINE_DPM_OUTPUT(hemocellOutput,header,fp,tp,t,plane)
    
12. {
    
13.          #if RP_2D
    
14.                 return;
    
15.          #else
    
16.                  
    
17.     cell_t c;
    
18.         real s;
    
19. 
    
20.     int zone_ID = 3;
    
21.     Thread *thread_name = Lookup_Thread(domain1,zone_ID);
    
22.        
    
23.     if(header)
    
24.     {
    
25.       par_fprintf_head(fp," ID, time, x-velocity[m/s], ");
    
26.       par_fprintf_head(fp,"y-velocity[m/s], z-velocity[m/s], ");
    
27.       par_fprintf_head(fp,"x-position[m], y-position[m], ");
    
28.       par_fprintf_head(fp,"z-position[m], strain rate \n");
    
29.     }
    
30.     if(NULLP(tp))
    
31.       return;
    
32.         c = TP_CELL(tp);
    
33.         s = C_STRAIN_RATE_MAG(c,thread_name);
    
34.         par_fprintf(fp,"%d %" int64_fmt " %" int64_fmt ", %e, %f, %f, %f, %f, %f, %f, %f \n",
    
35.    P_INJ_ID(TP_INJECTION(tp)), TP_ID(tp), TP_ID(tp), TP_TIME(tp),TP_VEL(tp)[0],
    
36.   TP_VEL(tp)[1],TP_VEL(tp)[2],TP_POS(tp)[0],TP_POS(tp)[1],TP_POS(tp)[2],s);
    
37.          #endif
    
38. 
    
39.    #if REMOVE_PARTICLES
    
40.       MARK_TP(tp, P_FL_REMOVED);
    
41.    #endif
    
42. }
    
43. 
    
44. DEFINE_ON_DEMAND(my_udf)
    
45. {
    
46.     Domain *domain;          /* domain is declared as a variable  */
    
47.     domain = Get_Domain(1);  /* returns fluid domain pointer      */
    
48.     domain1 = domain;
    
49. }

复制代码