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功能的實現是否正確？

回答

相關問題