【求助】e元件批处理命令
我做的非线性抛物方程,这几天的程序老是要吗没初值,要吗迭代一次就退出循环,今天看了看e元件批处理命令, `:o}:x}3Vl1。觉得不是书上教的那样?有点不理解? 7a&b2^ XO
tarta disp0 id0 nv disp0 bfd time0 time coor0 unod npa.io elem0 sys diag YdJ Qyud
if exist stop del stop
:1 \7P5Lr s1P@Q0T-y3v"Y
bft disp0 time coor0 nv bfdt stop
if exist end del end %P$`$eeir)d-MC
:2 #y2kW@)M
enpa disp0 sys time nv coor0 bfd diag elem0 f displace end 1sX sF8\lS&`@U
if not exist end goto 2 c?M3g,zw
call post.bat
if not exist stop goto 1 eh7uK#r'Aj0I$F2~
2。我的目录下没有生成end文件,为什么迭代后就退出,开始执行call post.bat? &cI4vpgz"|L^
3。我读了读enpa.for,发现程序中并没有用到enpa批处理命令最后两个文件displace ,end。为什么? ]0q%IR#kL,k/Y$K
附件是我的enpa.for !HL;j:?xq;owe
请斑竹指点!!学生我非常感激!!
回复: 【求助】e元件批处理命令
疑点可能找到,1。 4^B{"mo?uY2l
if (err .lt. 1.0e-8 .or. err .lt. 1.0e-8*ul) then h k)sN2M!e
open(11,file='end',status='unknown')
close(11)
我觉得能不能写成
if (err .lt. 1.0e-8) then。。。。 HU"Wa,Sy;Xt L!q
2。 ]M/yPs,F
每次生成的end文件在下一次时步开始就被delete 7BCB;x9n,d3o){z
if exist end del exist -N8m@7~8hGX]|
所以我找不到这个文件
{8jGjrk`DYX
请教斑竹上面对不对?谢谢
回复: 【求助】e元件批处理命令
但是即使我修改了程序,执行时还是按照原来的,怎么样才能执行修改后的程序?谢谢回复: 【求助】e元件批处理命令
我重新按照 !g@&v~/X5ffl32 /c enpa.for Buv$~'^|Mg#[2u7MJ
编译链接时,出现许多warning,主要是参数类型不匹配, ._*L)T{'@'e"Oy
是不是这个引起的?
回复: 【求助】e元件批处理命令
'qFIuqIif (err .lt. 1.0e-8 .or. err .lt. 1.0e-8*ul) then M"i@Vr5tkR,c
open(11,file='end',status='unknown') S4c:`)cv:n
close(11)
是收敛准则,可以修改 +RDw*e#TD X:Y7d
#H*RfOu8w
算完后应该有end文件
页:
[1]
