大家好，现在为了模拟一个流体力学现象，需要实现以下过程
某一空间内，参数 a 的当前值 是上一时间步内 参数 b 和 a 的函数 ，用数学表达式表示就是：a（i+1）= a（i）+b （i）
代码如下,大家帮忙看看有没有什么问题
DEFINE_PROFILE(function,t,i)
{
real x[ND_ND];
real b;
real a=100; /*a的初始值*/
cell_t c;
begin_c_loop(c,t)
{
C_CENTROID(x,c,t);
b=C_V_M1(c,t); /*前一时间步的y 方向速度*/
a=10.*b+a;
F_PROFILE(c,t,i) = a;
}
end_c_loop(c,t)
}
不甚感激，金币奉上 |
回复此楼
点击这里搜索更多相关资源