Fortran Coder

查看: 191|回复: 5

[数学库] 指针作为型参数传入函数以及接口的编写问题

[复制链接]

14

帖子

3

主题

0

精华

入门

F 币
66 元
贡献
34 点
发表于 2024-6-27 23:41:34 | 显示全部楼层 |阅读模式
最近在调用IMSL函数库的曲面拟合问题,其代码中用到了指针变量作为结构体的一个成员变量随着结构体一起定义作为形参在函数中调用,我写了一个类似的代码,有幸QQ群里的高人指点了下,有几个地方还不明白:1、Type Api中,Real(Kind=8),Pointer::Bus(:)为什么要这样定义?能否定义为Real(Kind=8),Pointer::Bus(10)?
2、下述代码运行是OK的,但是逐步跟踪时,在函数Hello外部,Tiger%Bus是有值的,但是到了函数内部,Tiger%Bus 就变成了 Undefined Pointer/Array,现实没有进行内存分配,这是什么原因,但是计算结果又是对的。

      
[Fortran] 纯文本查看 复制代码
Module Pointer_Mod
          Implicit None

          Type Api
          Integer::Id
          Real(Kind=8),Pointer::Bus(:)
          End Type Api

          Contains

          Subroutine Hello(Moon,Tiger,Fish)
          Integer,Intent(In)::Moon
          Type(Api),Intent(In)::Tiger
          Real(Kind=8),Intent(Out)::Fish

          Integer::I
          Fish=0.d0
          If(Moon==1)Then
          Do I=1,10,1
          Fish=Fish+Tiger%Bus(I)
          End Do 
          Else if(Moon==2)Then
          Do I=1,10,1
          Fish=Fish-Tiger%Bus(I)
          End Do 
          End If   !!
          End Subroutine        
End Module    

    program Test_Pointer
    Use Pointer_Mod
    implicit none

    Integer::Moon,I
    Type(Api)::Tiger
    Real(Kind=8)::Fish

    Moon=1
    Tiger%Id=10
    Allocate(Tiger%Bus(10))

    Do I=1,10,1
    Tiger%Bus(I)=I+1.d0
    End Do    
    Call Hello(Moon,Tiger,Fish)
    Write(*,*)Fish
    Pause
    end program Test_Pointer

Ps:
3、还有一个问题,如果该帖子有幸被堂主见到,在前面的贴中,关于IMSL的Surface_Fitting和Surface_Values两个函数,堂主都进行了接口的重新编写,可以实现IMSL库中对应例子的调用,但是有如下的扩展问题,对于Surface_Values的接口,当时堂主编写如下:

[Fortran] 纯文本查看 复制代码
interface SURFACE_VALUES    
    function DV_SURFACE_VALUES(DERIVATIVE, VARIABLESX, VARIABLESY, KNOTSX, KNOTSY, COEFFS)
      import
      implicit none
      integer :: DERIVATIVE(2)
      real(kind(1d0)) :: VARIABLESX(:), VARIABLESY(:)
      type (d_spline_knots) :: knotsx, knotsy
      real(kind(1d0)) :: COEFFS(:,:)
      real(kind(1d0)) :: DV_SURFACE_VALUES(size(VARIABLESX),size(VARIABLESY))
    end function DV_SURFACE_VALUES
  end interface SURFACE_VALUES


请问,上面的函数名字为什么要叫DV_SURFACE_VALUES,而改成D_SURFACE_VALUES时,调用就会报错?其次,VARIABLESX, VARIABLESY以及函数本身的返回值一定要是一维数组类型吗?能否使单纯的Real(Kind(1d0))类型?如果可以的话,是否可以做一个函数重载,再写一个function DV_SURFACE_VALUES(DERIVATIVE, VARIABLESX, VARIABLESY, KNOTSX, KNOTSY, COEFFS),其中对应的返回值与型参的是浮点型呢?
4、本人愚钝,堂主能否伸出您那充满光芒的双手再将另外一个IMSL库函数 Surface_Constraints也一并编写一个接口,我写了半天,调用老是不成功,不胜感激!跪谢!


1995

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1481 元
贡献
636 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2024-6-28 08:47:25 | 显示全部楼层
1. 指针的定义必须大小为假定的 (:) 。而不能是固定大小。
   不管是不是在一个Type中。
2. 调试器是人类设计最复杂的软件之一,所以经常有bug,仅作为参考。有时候要结合 write 输出到屏幕来判断。
3. 为什么要叫 DV_SURFACE_VALUES,我不知道。得问 IMSL 厂家为什么这么设定。有可能是历史遗留问题。我也是通过分析IMSL的*.mod和*.lib文件才知道的。
4. IMSL是闭源软件,你不能在他的基础上实现重载(除非你有独立书写SURFACE_VALUES函数的能力。)
闭源函数库,你必须严格遵循它本身的函数接口,故而VARIABLESX, VARIABLESY以及函数本身的返回值一定要是一维数组类型。否则调用就会失败。
5. 看到你的需求了,这两天有点忙,我看周末抽个时间帮你弄弄。

14

帖子

3

主题

0

精华

入门

F 币
66 元
贡献
34 点
 楼主| 发表于 2024-6-28 10:06:20 | 显示全部楼层
十分感谢版主大人的耐心细致的讲解,醍醐灌顶,受益匪浅!
对版主大人的无私付出,表示由衷的感谢和致敬!

14

帖子

3

主题

0

精华

入门

F 币
66 元
贡献
34 点
 楼主| 发表于 2024-6-30 09:43:28 | 显示全部楼层
fcode 发表于 2024-6-28 08:47
1. 指针的定义必须大小为假定的 (:) 。而不能是固定大小。
   不管是不是在一个Type中。
2. 调试器是人类设 ...

Dear 堂主!
打扰您周末休息了,麻烦您看看我的需求,您有时间的话,麻烦帮我看看Surface_Constraints的接口文件的重写封装,谢谢您了!

1995

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1481 元
贡献
636 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2024-7-3 09:14:02 | 显示全部楼层
[Fortran] 纯文本查看 复制代码
Module surface_fitting_int
  USE norm_int
  implicit none
  interface surface_fitting
    function D_SURFACE_FITTING(DATA,KNOTSX,KNOTSY,CONSTRAINTS,COVARIANCE,IOPT,epack)
      import
      implicit none
      type (d_spline_knots) :: knotsx, knotsy
      real(kind(1d0)) :: DATA(:,:)
      type (D_SURFACE_CONSTRAINTS) , optional :: CONSTRAINTS(:)
      real(kind(1d0)) , optional :: COVARIANCE(:,:)
      type (d_options) , optional :: iopt(:)
      type (d_error) , optional :: epack
      real(kind(1d0)) :: D_SURFACE_FITTING(size(KNOTSX%D_KNOTS)-KNOTSX%SPLINE_DEGREE-1,size(KNOTSY%D_KNOTS)-KNOTSY%SPLINE_DEGREE-1)
    end function D_SURFACE_FITTING    
  end interface surface_fitting

  interface SURFACE_VALUES    
    function DV_SURFACE_VALUES(DERIVATIVE, VARIABLESX, VARIABLESY, KNOTSX, KNOTSY, COEFFS)
      import
      implicit none
      integer :: DERIVATIVE(2)
      real(kind(1d0)) :: VARIABLESX(:), VARIABLESY(:)
      type (d_spline_knots) :: knotsx, knotsy
      real(kind(1d0)) :: COEFFS(:,:)
      real(kind(1d0)) :: DV_SURFACE_VALUES(size(VARIABLESX),size(VARIABLESY))
    end function DV_SURFACE_VALUES
  end interface SURFACE_VALUES
  
  interface SURFACE_CONSTRAINTS
    function DSURFACE_CONSTRAINTS(derivative,point,value,type,periodic)
      import
      type (d_surface_constraints) :: DSURFACE_CONSTRAINTS
      integer , optional :: derivative(:)
      real(kind(1d0)) :: point(:) 
      real(kind(1d0)) , optional :: periodic(:)
      real(kind(1d0)) , optional :: value
      character(len=*):: type      
    end function DSURFACE_CONSTRAINTS
  end interface SURFACE_CONSTRAINTS

End Module surface_fitting_int

14

帖子

3

主题

0

精华

入门

F 币
66 元
贡献
34 点
 楼主| 发表于 2024-7-3 14:55:25 | 显示全部楼层
fcode 发表于 2024-7-3 09:14
[mw_shl_code=fortran,true]Module surface_fitting_int
  USE norm_int
  implicit none

感谢堂主,您辛苦了,您费心了!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-7-23 20:45

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表