[求助]WENO编程问题！！！

endly111

文字文字最近自己在看有关WENO的一些东西，想自己编个程序，可是程序调试出现了两个问题，一个就是按照论文中所说的用那个WENO加权得到的两组数据不知道怎么用？还有就是时间离散的时候那个L(U)到底应该是个什么表达式？比如对于f=u^2/2的情形？希望各位高手指导一下，不胜感激啊！

wangaifeng

weno加权你可以直接选择拉格朗日插值中的数据，选择均匀网格，最常见的就是weno5阶，可以在写程序的时候写成一个子程序，循环调用；再一个lu=ui+1-ui-1一般是这样的

taibaisha

其实很简单的。先要选择好模板，按照论文上给的；然后你可以将守恒变量或者原始变量投影到特征变量空间（也可以不走这一步）；其次将这些投影好的变量进行插值，中间涉及到计算一些系数（比如光滑因子的计算），这就需要那两组数据了；插值完成以后就进行反变换，返回到守恒变量或者原始变量空间，这样，插值就算完成了。

yeclerk

我回复这个贴子吧，我把自己写的一段程序给粘上去，大家看看！
double weno(double data[5])//插(1/2)处的通量
{
double s1,s2,s3;
double err=1.0e-6;
double w1,w2,w3,w;
double value;
s1=(13.0/12.0)*pow((data[0]-2.0*data[1]+data[2]),2)+(1.0/4.0)*pow((data[0]-4.0*data[1]+3.0*data[2]),2);
s2=(13.0/12.0)*pow((data[1]-2.0*data[2]+data[3]),2)+(1.0/4.0)*pow((data[3]-data[1]),2);
s3=(13.0/12.0)*pow((data[2]-2.0*data[3]+data[4]),2)+(1.0/4.0)*pow((3.0*data[2]-4.0*data[3]+data[4]),2);
w1=1.0*pow((err+s2),2)*pow((err+s3),2);
w2=6.0*pow((err+s1),2)*pow((err+s3),2);
w3=3.0*pow((err+s1),2)*pow((err+s2),2);
/*w1=0.7;
w2=0;
w3=0.3;*/
/*
w1=1.0/(pow((err+s1),2));//没有除10,因为10*w1/(10*w1+10*w2+10*w3)=w1/(w1+w2+w3)
w2=6.0/(pow((err+s2),2));
w3=3.0/(pow((err+s3),2));*/
w=w1+w2+w3;
value=(w1*(2.0*data[0]-7.0*data[1]+11.0*data[2])+w2*((-1.0)*data[1]+5.0*data[2]+2.0*data[3])+w3*(2.0*data[2]+5.0*data[3]-data[4]))/(6.0*w);//数值通量f[i+1/2]
return value;
}[br][br][以下内容由 yeclerk 在 2007年10月20日 10:49am 时添加] [br]
望对大家有所帮助[br][br][以下内容由 yeclerk 在 2007年10月20日 10:49am 时添加] [br]
望对大家有所帮助