不能在FORTRAN子程序中使用SUNDIALS-KINSOL？

Question

感谢您查看此问题。

问题： SEG。从包含KINSOL求解过程的f90子程序返回时出现错误， 计算结果已生成。当主程序中有相同的求解过程时没问题。

环境： linux下， 海湾合作委员会， 日su静态库

如何解决问题：

获取附加的REDUCED测试代码

module moduleNonlinearSolve
  integer,save::nEq
contains
  subroutine solveNonlinear(u)
    double precision::u(*)
    integer iout(15),ier
    double precision rout(2),koefScal(nEq)
    koefScal(:)=1d0
    call fnvinits(3,nEq,ier)
    call fkinmalloc(iout,rout,ier)
    call fkinspgmr(50,10,ier)
    call fkinsol(u,1,koefScal,koefScal,ier)
    call fkinfree()
    do i=1,nEq
      write(*,*),i,u(i)
    end do
  end subroutine
end module

subroutine fkfun(u,fval,ier)
  use moduleNonlinearSolve
  double precision::u(*)
  double precision::fval(*)
  integer::ier
  forall(i=2:nEq-1)
    fval(i)=-u(i-1)+2d0*u(i)-u(i+1)-1d0
  end forall
  fval(1)=u(1)+2d0*u(1)-u(2)-1d0
  fval(nEq)=-u(nEq-1)+2d0*u(nEq)+u(nEq)-1d0
  ier=0
end subroutine

program test
  use moduleNonLinearSolve
  double precision u(10)
  nEq=size(u)
  u(:)=10d0
  call solveNonlinear(u)
end program``

编译

$ gfortran -c -Wall -g test.f90
$ gfortran -Wall -g -o test test.o -lsundials_fkinsol -lsundials_fnvecserial -lsundials_kinsol -lsundials_nvecserial -llapack -lblas

运行

$ ./test

注意：如果将所有SUNDIALS程序放在主程序中，它将完美无缺。

非常感谢您的任何意见。

绵纸

Answer 1

根据KINSOL文档，fkinmalloc的第一个参数必须与C类型long int具有相同的整数类型。在您的情况下，long int长度为8个字节，但您传入的是一个4字节整数的数组。这将导致fkinmalloc尝试超出数组边界并进入其他内存。这通常会导致内存损坏，其症状就像您正在观察的那样：稍后随机崩溃，例如从函数返回时。您应该能够通过valgrind运行程序来确认这一点，这可能会报告大小为8的无效写入。无论如何，替换

integer :: iout(15)

与

integer*8 :: iout(15)

应该解决问题。