Fortran代码：PHI功能的实现是否正确？

REAL FUNCTION PHI(D) 
    COMMON FAC,DELTA,ER,T 
    DOW=2.*D/DELTA 
    TOW=4.*T/DELTA 
    DOWM=DOW-1. 
    DOWP=DOW+1. 
    IF(D.EQ.0.) GO TO 1 
    PHI=DOWM*ALOG(DOWM)-DOWP*ALOG(DOWP) 
1-0.5*DOWM*ALOG(DOWM**2+TOW**2) 
2+0.5*DOWP*ALOG(DOWP**2+TOW**2) 
3-TOW*(ATAN(DOWM/TOW)-ATAN(DOWP/TOW)) 
    PHI=PHI*FAC/2. 
    GO TO 2 
    PHI=FAC*(0.5*ALOG(1.+TOW*TOW)+TOW*ATAN(1./TOW)) 
    CONTINUE 
    RETURN 
    END

我的作业是将上面的代码转换成Matlab .m文件。Fortran代码：PHI功能的实现是否正确？

但我不能明白，我从来没有见过这样的实现：

.... 
    PHI=DOWM*ALOG(DOWM)-DOWP*ALOG(DOWP) 
1-0.5*DOWM*ALOG(DOWM**2+TOW**2) 
2+0.5*DOWP*ALOG(DOWP**2+TOW**2) 
3-TOW*(ATAN(DOWM/TOW)-ATAN(DOWP/TOW)) 
...

没有任何续行标志！这很混乱。我应该如何将这些行转换为Matlab代码？我知道Matlab，但我不知道Fortran。

来源

2017-05-27 valgrind

它看起来像这个功能是从[这里]复制（http://read.pudn.com/downloads110/sourcecode/book/453805/Crosstalk%20Fortran%20Programs/MSTRP.FOR__.htm），并在遭受过程（例如，线标签'1'和'2'丢失）。 –

我认为这是以固定形式提供的源代码，所以请注意代码的缩进（这里的编辑会吞下前四个空格）。如果是这样，那么确实存在行续标记：您引用的第6列中的“1”，“2”和“3”。 – francescalus

避免所有的帽子，在Fortrsn看起来很糟糕，但它是一个旧代码。但不要在新代码中使用它。还要在=运算符和括号之间使用空格，以便将运算组合在一起。 –

的PHI计算是在四大行中，第二线标有，等等...

这些都是的延续，标志着你正在寻找（ 1，2，3）

你的计算公式为：

PHI=DOWM*ALOG(DOWM)-DOWP*ALOG(DOWP)-0.5*DOWM*ALOG(DOWM**2+TOW**2)+0.5*DOWP*ALOG(DOWP**2+TOW**2)-TOW*(ATAN(DOWM/TOW)-ATAN(DOWP/TOW))

但请记住，原来的Fortran无法理解第73列及以上（仅1至72）的字符，这就是为什么在第6列中使用连续标记（通常是星号，但数字也可以使用）的原因。

另外代码只能在第7到第72列之间，因为第1到第5列是为标签保留的。

我的猜测是，标签1是在PHI第二个定义的前面，并在该CONTINUE前面的标签2，这是应该如何看待用Fortran

GO TO 2 
1 PHI=FAC*(0.5*ALOG(1.+TOW*TOW)+TOW*ATAN(1./TOW)) 
2 CONTINUE

中的Fortran代码应该是这样的：

 REAL FUNCTION PHI(D) 
     COMMON FAC,DELTA,ER,T 
     DOW=2.*D/DELTA 
     TOW=4.*T/DELTA 
     DOWM=DOW-1. 
     DOWP=DOW+1. 
     IF(D.EQ.0.) GO TO 1 
     PHI=DOWM*ALOG(DOWM)-DOWP*ALOG(DOWP) 
    1-0.5*DOWM*ALOG(DOWM**2+TOW**2) 
    2+0.5*DOWP*ALOG(DOWP**2+TOW**2) 
    3-TOW*(ATAN(DOWM/TOW)-ATAN(DOWP/TOW)) 
C 
C Previous four (4) lines have this formula (C mark for comment) 
C 
C PHI=DOWM*ALOG(DOWM)-DOWP*ALOG(DOWP)-0.5*DOWM*ALOG(DOWM**2+TOW**2)+0.5*DOWP*ALOG(DOWP**2+TOW**2)-TOW*(ATAN(DOWM/TOW)-ATAN(DOWP/TOW)) 
C 
     PHI=PHI*FAC/2. 
     GO TO 2 
1  PHI=FAC*(0.5*ALOG(1.+TOW*TOW)+TOW*ATAN(1./TOW)) 
2  CONTINUE 
     RETURN 
     END

来源

2017-05-27 16:17:36 TrustworthySystems

令人困惑的部分是跳转**转到1 **和**转到2 **，因为您的代码中没有显示（标签1和2应该在第1列的某处）。你也有一个** CONTINUE **子句，但没有循环或块，是无害的，但显然这里有一些遗漏的代码。 – TrustworthySystems

如果可以，请提供原始的Fortran代码以帮助您更好地进行操作 – TrustworthySystems

A.continue不需要执行循环。 –

Fortran代码：PHI功能的实现是否正确？

回答

相关问题