Fortran 编码建议答案

【问题标题】：Fortran Coding AdviceFortran 编码建议
【发布时间】：2015-10-25 17:28:47
【问题描述】：

我必须开发一个线性插值程序，但不断收到这些错误。以下是源代码：

!Interpolation program for exercise 1 of portfolio
PROGRAM interpolation
IMPLICIT NONE
!Specify table 1 for test of function linter
REAL, DIMENSION (5):: x,f
!Specify results for table 1 at intervals of 1
REAL, DIMENSION (10):: xd, fd
!Specify table 2 to gain linter results 
REAL, DIMENSION (9):: xx,ff
!Specify results for table 2 of at intervals of 0.25
REAL, DIMENSION (36):: xxd, ffd

INTEGER :: i, j

!Write values for table dimensions

!Enter x values for Table 1
x(1)=-4.0
x(2)=-2.0
x(3)=0.0
x(4)=2.0
x(5)=4.0
f(1)=28.0
f(2)=11.0
f(3)=2.0
f(4)=1.0
f(5)=8.0
xd(1)=-4.0
xd(2)=-3.0
xd(3)=-1.0
xd(4)=0.0
xd(5)=1.0
xd(6)=2.0
xd(7)=3.0
xd(9)=4.0

!Print Table 1 Array
PRINT *,"Entered Table Values are", x,f
PRINT *,"Interpolation Results for Table 1", xd, fd

END PROGRAM

SUBROUTINE interpol(x,f, xd,fd)
DO i=1, 5
    DO j=1, 5
        IF (x(j) <  xd(i) .AND. xd(i) <=  x(j+1)) THEN
        fd=linterp (x(j),x(j+1),f(j))
        END IF
    END DO 
END DO
END SUBROUTINE interpol


!Linear Interpolation function
FUNCTION linterp(x(i),x(i+1),f(i),f(i+1),x)
    linterp=f(i)+((x-x(i))/(x(i+1)-x(i)))*(f(i+1)-f(i))
END FUNCTION

它给出了这些错误；

lin.f90:55:18: Error: Expected formal argument list in function definition at (1)
lin.f90:56:19:

  linterp=f(i)+((x-x(i))/(x(i+1)-x(i)))*(f(i+1)-f(i))
                   1
Error: Expected a right parenthesis in expression at (1)
lin.f90:57:3:

 END FUNCTION
   1
Error: Expecting END PROGRAM statement at (1)

谁能指点我正确的方向？

【问题讨论】：

你有一个shortage of implicit none（正如 Alexander Vogt 的回答一样）：你可能会对 linterp 返回一个整数值感到惊讶。

标签： fortran fortran90 gfortran

【解决方案1】：

这正是编译器所抱怨的：你缺少一个右括号。

要么去掉多余的左边(：

linterp=f(i)+ ( x-x(i) ) / ( x(i+1)-x(i) )* ( f(i+1)-f(i) )

或添加另一个)

linterp=f(i)+ ( (x-x(i)) ) / ( x(i+1)-x(i) )* ( f(i+1)-f(i) )

请注意，我在中间部分删除了另一个错误放置的)。

除此之外，您的函数声明已损坏！声明中不能有x(i)！

试试：

real FUNCTION linterp(xI,xIp1,fI,fIp1,x)
    implicit none
    real, intent(in)    :: xI,xIp1,fI,fIp1,x
    linterp = fI + (x-xI)/(xIp1-xI)*(fIp1-fI)
END FUNCTION

或者，您可以提供整个数组（包括其长度N）和当前索引：

real FUNCTION linterp(x,f,N,i,xx)
    implicit none
    integer, intent(in) :: N
    real, intent(in)    :: x(N), f(N), xx
    integer, intent(in) :: i
    linterp = f(i) + (xx-x(i))/( x(i+1)-x(i) )*( f(i+1)-f(i) )
END FUNCTION

【讨论】：

太好了，我知道我哪里出错了，我已经编辑了函数来适应这个，谢谢！

【解决方案2】：

除了亚历山大所说的一切。您还需要确保函数声明中的输入数量与调用它时相同：

fd=linterp (x(j),x(j+1),f(j))

比你的函数声明少两个输入：

FUNCTION linterp(x(i),x(i+1),f(i),f(i+1),x)

另外，别忘了给 fd 添加一个索引，i 或 j：

fd(i)=linterp (x(j),x(j+1),f(j))

否则你每次都用 linterp 结果替换整个数组。

【讨论】：