Eliteruirui 2022-12-03 11:18 采纳率: 0%
浏览 30
已结题

freefem++ gnuplot文件生成有问题

问题遇到的现象和发生背景

freefem++怎么能让这个程序成功编译呢?

用代码块功能插入代码,请勿粘贴截图
verbosity=0;
real Dx=.2, R=10.;
border C(t=0.,2.*pi){x=R*cos(t);y=R*sin(t);label=1;};
mesh Th=buildmesh(C(floor(2.*pi*R/Dx)));
fespace Vh(Th,P2);
real dt=0.01,Tf=10., lambda=-1., p=3.,alpha=0.;
Vh<complex> uh, vh, uh0=exp(-x^2-y^2-5.*1i*x), uhk=uh0, TNL, B;
varf a(u,v) = int2d(Th)(u*v*1i/dt + u*v*1i*alpha/2. + (dx(u)*dx(v) + dy(u)*dy(v))/2.) + on(1,u=0);
matrix<complex> A = a(Vh,Vh);
varf b(u,v) = int2d(Th)(uh0*v*1i/dt - uh0*v*1i*alpha/2. - (dx(uh0)*dx(v) + dy(uh0)*dy(v))/2. - lambda*TNL*(uhk+uh0)*v/2.) + on(1,u=0);
Vh ABSU;
int kk=0;
real[int] NORML2(floor(Tf/dt)+1);
for (real t=0.;t<=Tf;t+=dt){
    TNL=abs(uh0)^(p-1);
     for (int i=0;i<2;i++){
        B[] = b(0,Vh);
        set(A,solver=sparsesolver);
        uhk[] = A^-1*B[];
    }
    uh0[]=uhk[];
    ABSU=abs(uh0);
    NORML2[kk]=sqrt(int2d(Th)(abs(uh0)^2));        
    if ( !(kk % 10)){
        plot(ABSU,cmm="t="+t+"  ;||u||_L^2="+NORML2[kk], fill=true,value=true,dim=2);
    {
        ofstream gnufile("||u||_L2.gnu");
         for (int i=0;i<=kk;i++)
            gnufile<<i*dt<<" "<<NORML2(i)<<endl;
    }
exec("echo 'plot \"||u||_L2.gnu\" w lp \pause 5 \quit' | gnuplot");}
    kk++;
}
```c++


```

运行结果及报错内容

Error opening file ||u||_L2.gnu
current line = 32
Exec error : Error opening file
-- number :1
Exec error : Error opening file
-- number :1
err code 8 , mpirank 0

我的解答思路和尝试过的方法
我想要达到的结果
  • 写回答

2条回答 默认 最新

  • Jackyin0720 2022-12-03 12:18
    关注

    你报错信息看不到哦,

    评论

报告相同问题?

问题事件

  • 系统已结题 12月11日
  • 修改了问题 12月3日
  • 创建了问题 12月3日

悬赏问题

  • ¥15 安装quartus II18.1时弹出此error,怎么解决?
  • ¥15 keil官网下载psn序列号在哪
  • ¥15 想用adb命令做一个通话软件,播放录音
  • ¥30 Pytorch深度学习服务器跑不通问题解决?
  • ¥15 部分客户订单定位有误的问题
  • ¥15 如何在maya程序中利用python编写领子和褶裥的模型的方法
  • ¥15 Bug traq 数据包 大概什么价
  • ¥15 在anaconda上pytorch和paddle paddle下载报错
  • ¥25 自动填写QQ腾讯文档收集表
  • ¥15 DbVisualizer Pro 12.0.7 sql commander光标错位 显示位置与实际不符