在Fortran中不匹配的接口90

我很难理解Fortran接口。在问题结尾处粘贴的测试代码代码在ifort < =版本15上编译得很好，但在版本16和更高版本上会导致灾难性编译器错误。同样gfortran（编译时出现seg fault）。在G95，然而，编译，出现以下错误在Fortran中不匹配的接口90

In file test.f90:79 

call rk4(y, dydx, x, h, yout1, der, pars,*10) 
           1 
Error: Interface of actual procedure does not match interface of dummy procedure at (1) 
In file test.f90:81 

call rk4(y, dydx, x, h/2.d0, yout2, der, pars,*10) 
            1 
Error: Interface of actual procedure does not match interface of dummy procedure at (1) 
In file test.f90:85 

call rk4(yout2, dydx2, x+h/2.d0, h/2.d0, yout2, der, pars,*10) 
               1 
Error: Interface of actual procedure does not match interface of dummy procedure at (1)

结束，但我不知道我做错了......

module RK 

implicit none 

contains 

SUBROUTINE rk4(y, dydx, x, h, yout, der, pars, *) 
IMPLICIT NONE 
integer, parameter :: dp=kind(1.d0) 
REAL(dp), DIMENSION(:), INTENT(IN) :: y,dydx,pars 
REAL(dp), INTENT(IN) :: x,h 
REAL(dp), DIMENSION(:), INTENT(OUT) :: yout 
REAL(dp) :: h6,hh,xh 
REAL(dp), DIMENSION(size(y)) :: dym,dyt,yt 

INTERFACE 
SUBROUTINE der(x,y,dydx,pars,*) 
IMPLICIT NONE 
integer, parameter :: dp=kind(1.d0) 
REAL(dp), INTENT(IN) :: x 
REAL(dp), DIMENSION(:), INTENT(IN) :: pars 
REAL(dp), DIMENSION(:), INTENT(IN) :: y 
REAL(dp), DIMENSION(:), INTENT(OUT) :: dydx 
END SUBROUTINE der 
END INTERFACE 


hh = h*0.5d0 
h6 = h/6.d0 
xh = x + hh 

yt(:)=y(:)+hh*dydx(:) 

call der(xh, yt, dyt, pars,*10) 

yt(:)=y(:)+hh*dyt(:) 

call der(xh, yt, dym, pars,*10) 

yt(:)=y(:)+h*dym(:) 

dym(:)=dyt(:)+dym(:) 

call der(x+h, yt, dyt, pars,*10) 

yout(:)=y(:)+h6*(dydx(:)+dyt(:)+2.d0*dym(:)) 

return 
10 return 1 


end subroutine 


subroutine adaptive_RK4(y,dydx,x,h,yout2,yerr,der,pars,*) 
implicit none 
integer, parameter :: dp=kind(1.d0) 
REAL (dp) :: h, x 
REAL (dp), DIMENSION (:) :: y, dydx, yout2, yerr, pars 
REAL (dp), DIMENSION (size(y)) :: dydx2, yout1 
intent(in) :: y, dydx, x, h, pars 
intent(out) :: yout2, yerr 
INTERFACE 
SUBROUTINE der(x,y,dydx,pars,*) 
IMPLICIT NONE 
integer, parameter :: dp=kind(1.d0) 
REAL(dp), INTENT(IN) :: x 
REAL(dp), DIMENSION(:), INTENT(IN) :: pars 
REAL(dp), DIMENSION(:), INTENT(IN) :: y 
REAL(dp), DIMENSION(:), INTENT(OUT) :: dydx 
END SUBROUTINE der 
END INTERFACE 


call rk4(y, dydx, x, h, yout1, der, pars,*10) 

call rk4(y, dydx, x, h/2.d0, yout2, der, pars,*10) 

call der(x+h/2.d0, yout2, dydx2, pars,*10) 

call rk4(yout2, dydx2, x+h/2.d0, h/2.d0, yout2, der, pars,*10) 

yerr(:)=yout2(:)-yout1(:) 

yout2(:)=(16.d0*yout2(:)-yout1(:))/15.d0 

return 

10 return 1 

end subroutine 

end module 

module derivative 
implicit none 
contains 
    SUBROUTINE derr(x,y,dydx,pars,*) 
    IMPLICIT NONE 
integer, parameter :: dp=kind(1.d0) 
    REAL(dp), INTENT(IN) :: x 
    REAL(dp), DIMENSION(:), INTENT(IN) :: pars 
    REAL(dp), DIMENSION(:), INTENT(IN) :: y 
    REAL(dp), DIMENSION(:), INTENT(OUT) :: dydx 


dydx=(pars+y)*x 
return 

END SUBROUTINE derr 

end module 

program test 
use rk 
use derivative 
implicit none 
integer, parameter :: dp=kind(1.d0) 
real(dp), dimension(2) :: y,pars,dydx, yout, yerr 
real(dp) :: x 

y=(/0.1d0,2.d0/) 
pars=(/0.7d0,3.d0/) 
x=2.1d0 

call derr(x,y,dydx,pars,*10) 

write(*,*) dydx 
call adaptive_RK4(y,dydx,x,0.0001d0,yout,yerr,derr,pars,*10) 

stop 
10 write(*,*) "some error" 

end program

来源

2017-10-09 Enrico B

欢迎光临。务必采取欢迎[游览]。针对所有Fortran问题使用标签[tag：fortran]。您可以为特定于版本的问题添加特定的版本标签，但请注意，Fortran 90是旧的和过时的版本，通常您不希望受限于此。 –

请使用一些缩进（每行开始处的空格）。你的代码很难阅读。我们看不到结构。 –

而OMG停止使用备用返回功能（https：// stackoverflow。com/questions/7835270/asterisks-in-fortran-syntax-in-argument-list-at-1）现在，这是一种下地狱，相信我们的方式。 –

的代码运行在gfortran版本6和7

罚款

Gfortran版本4.8崩溃，但这是编译器的问题。即使代码不好，编译器也不应该崩溃。

对于英特尔Fortran也是如此。如果在编译过程中崩溃，则是编译器的错误。

现在，代码真的不是很好，并且使用编译器可能很难的功能。当涉及到对标准的东西，gfortran 6抱怨的：

interface2.f90:81:9: 

call rk4(yout2, dydx2, x+h/2.d0, h/2.d0, yout2, der, pars,*10) 
     1 
Warning: Same actual argument associated with INTENT(IN) argument ‘y’ and INTENT(OUT) argument ‘yout’ at (1)

它不是用Fortran允许通过相同的变量中的参数以这种方式。编译器将假定intent(in)参数不会更改，但在修改yout时确实会更改。这可能会导致严重的问题。

对于G95编译时错误，我认为编译器是错误的，错误信息是假的。 rk4和adaptive_rk4中的接口块完全相同。

我强烈建议您清理代码，摆脱备用退出功能，编译器编译代码会更容易。我预计这是由于您过度使用未经太多测试的备用返回功能。

你的下一个问题，肯定读[MCVE]非常谨慎。有必要使您的代码短得多为您的未来的问题。这是实际MCVE你应该询问，并可以作为一个bug报告基地：

module m 

    implicit none 

contains 

    subroutine a(d) 
    interface 
     subroutine d(*) 
     end subroutine d 
    end interface 

    10 return 1 
    end subroutine 


    subroutine b(d) 
    interface 
     subroutine d(*) 
     end subroutine d 
    end interface 

    call a(d) 
    end subroutine 
end module

你给了我们很长的代码进行比较。如果没有编译器崩溃消失（ifort 17和gcc 4.8），这不会真的被减少。这个问题在虚拟过程中很明显（尽管它实际上并没有被调用），它使用了备用返回特性，编译器无法正确处理它。

来源

2017-10-09 11:54:01

谢谢！我没有注意到intent（out）属性的问题，后者在代码编译时确实发生了一些奇怪的事情。对于过时的功能/长代码抱歉，我将来会记得 –

我从DERR/DER实例中替换了备用返回，然后编译在英特尔Fortran（XE 2017）中运行。我同意，建议（1）清理代码并（2）避免使用替代退货。

来源

2017-10-09 13:28:22

在Fortran中不匹配的接口90

回答

相关问题