| 查看: 478 | 回复: 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) } 望各位前辈给我指点 |
回复此楼» 猜你喜欢
303求调剂
已经有4人回复
-
311求调剂
已经有3人回复
-
307求调剂
已经有3人回复
-
求调剂
已经有6人回复
-
课题组接收材料类调剂研究生
已经有4人回复
-
295复试调剂
已经有3人回复
-
一志愿中南大学理学化学
已经有4人回复
-
298求调剂
已经有12人回复
-
291分工科求调剂
已经有9人回复
-
284求调剂
已经有6人回复
