| 查看: 535 | 回复: 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) } 望各位前辈给我指点 |
回复此楼» 猜你喜欢
申博发邮件
已经有8人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有9人回复
-
国社科系统bug了，是不是要放榜了？
已经有7人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有8人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复