| 查看: 482 | 回复: 0 | ||
[求助]
udf程序调试
|
|
各位前辈： 大家好，我这里有一个程序，调试时显示以下信息： Error: C:\Users\user\Desktop\udf unsteady temperature w 1.c: line 27: parse error. Error: C:\Users\user\Desktop\udf unsteady temperature w 1.c: line 38: y: undeclared variable 程序如下： #include <udf.h> #include <math.h> #define ERF_PI 3.141592653589793 #define ERF_N 100 double erf(double x)////erf(x) = the cumulation of { 2/sqrt(pi)*exp(-z*z) } from 0 to x; //error function { // 2/sqrt(pi) * { ∑[(-1)^n / n! * x^(2n+1)/(2n+1)] + x } double res = x; double factorial = 1; //n! double x_pow = x; int one = 1, n; for( n=1; n<100; n++ ){ factorial *= n; one *= -1; x_pow *= x*x; res += one / factorial * x_pow / ( 2*n+1 ); } res *= 2 / sqrt(ERF_PI); return res; } double norm_cdf(double y)//cumulation distribution function of standard normal distribution { real time = CURRENT_TIME; y=(time-5.98023)/3.3487-3.3487/5.90205 return ( 1 + erf( y / sqrt(2) ) ) / 2; } DEFINE_PROFILE(unsteady_temperature, t, i) { real time = CURRENT_TIME; real temperature; face_t f; begin_f_loop(f, t) { temperature = -4.6376+6.9448*exp(1.1742-time/5.90205)*(erf(y/sqrt(2))+1)/2; F_PROFILE(f, t, i) = temperature; } end_f_loop(f, thread) } 望各位前辈给我指点 |
回复此楼» 猜你喜欢
求调剂 材料与工程 324分 专硕
已经有3人回复
-
化学工程与技术324调剂
已经有15人回复
-
一志愿矿大，材料工程专硕314分，0856可调都可以
已经有12人回复
-
环氧灌封胶 抗沉剂
已经有4人回复
-
268分085602化学工程调剂
已经有22人回复
-
考研调剂
已经有18人回复
-
还有化工二轮调剂的学校吗
已经有39人回复
-
277 数一104，学硕，求调剂
已经有8人回复
-
0854调剂
已经有5人回复
-
材料工程302分求调剂
已经有15人回复
