C函式当参数回传到Fortran Dll

Radiosan · 发表于 2016-12-5 23:43:41

各位好! 小弟我想在C++下使用Fortran套件 ODRPACK95。我的作业环境为 Win10，使用VS2013，Fortran编译器为Intel Visual Fortran。
首先下面是我的Fortran程序，用来将C++下的变量值传递到ODRPACK做运算。

[Fortran] 纯文本查看 复制代码

subroutine wrapper_ODR(FCN,N,M,NP,NQ,BETA,Y,X,&
        DELTA,WE,WD,IFIXB,IFIXX,JOB,NDIGIT,TAUFAC,&
        SSTOL,PARTOL,MAXIT,IPRINT,LUNERR,LUNRPT,&
     STPB,STPD,SCLB,SCLD,WORK,IWORK,INFO,LOWER,UPPER) bind(C, name='wrapper_ODR')
!DEC$ ATTRIBUTES DLLEXPORT :: wrapper_ODR    
    use iso_c_binding
    use ODRPACK95
    implicit none
    
!自訂意凾式介面宣告    
    interface 
        subroutine FCN(N,M,NP,NQ,LDN,LDM,LDNP,BETA,XPLUSD,IFIXB,IFIXX,LDIFX,&
        IDEVAL,F,FJACB,FJACD,ISTOP) bind(C)
        
        use, intrinsic :: iso_c_binding
        implicit none
        
        integer(c_int) :: IDEVAL,ISTOP,LDIFX,LDM,LDN,LDNP,M,N,NP,NQ
        
        real (c_double) :: BETA(1:NP),F(1:LDN,1:NQ),FJACB(1:LDN,1:LDNP,1:NQ), &
        FJACD(1:LDN,1:LDM,1:NQ),XPLUSD(1:LDN,1:M)
        
        integer(c_int) :: IFIXB(1:NP),IFIXX(1:LDIFX,1:M)
    
        end subroutine
    end interface
  
    
    integer(c_int),value :: N,M,NP,NQ
    real(c_double) :: BETA(1:NP),Y(1:N,1:NQ),X(1:N,1:M)
    
!!!!!Optional variable
    
    integer(c_int), intent(in), optional :: IFIXB(:),IFIXX(:,:),JOB,NDIGIT,MAXIT&
    ,IPRINT,LUNERR,LUNRPT,IWORK(:),INFO
    
    real(c_double), intent(in), optional :: DELTA(:,:),&
        WE(:,:,:),WD(:,:,:),TAUFAC,SSTOL,PARTOL,&
        STPB(:),STPD(:,:),SCLB(:),SCLD(:,:),&
         WORK(:),LOWER(:),UPPER(:)
    
!!!!!Call ODR    
    call ODR(inter_func,N,M,NP,NQ,BETA,Y,X)
    
contains 

!!!!!Local subroutine
subroutine inter_func(N,M,NP,NQ,LDN,LDM,LDNP,BETA,XPLUSD,IFIXB,IFIXX,LDIFX,&
   IDEVAL,F,FJACB,FJACD,ISTOP)
    
        use REAL_PRECISION

        integer :: IDEVAL,ISTOP,LDIFX,LDM,LDN,LDNP,M,N,NP,NQ
       
        real (kind=R8) :: BETA(1:NP),F(1:LDN,1:NQ),FJACB(1:LDN,1:LDNP,1:NQ), &
        FJACD(1:LDN,1:LDM,1:NQ),XPLUSD(1:LDN,1:M)
       
        integer :: IFIXB(1:NP),IFIXX(1:LDIFX,1:M)
       
        integer(c_int) :: inter_IDEVAL,inter_ISTOP,inter_LDIFX,inter_LDM,&
        inter_LDN,inter_LDNP,inter_M,inter_N,inter_NP,inter_NQ
        
        real (c_double) :: inter_BETA(1:NP),inter_F(1:LDN,1:NQ),&
        inter_FJACB(1:LDN,1:LDNP,1:NQ),&
        inter_FJACD(1:LDN,1:LDM,1:NQ),inter_XPLUSD(1:LDN,1:M)
        
        integer(c_int) :: inter_IFIXB(1:NP),inter_IFIXX(1:LDIFX,1:M)
     
        inter_IDEVAL =  IDEVAL
        inter_ISTOP = ISTOP
        inter_LDIFX = LDIFX
        inter_LDM = LDM
        inter_LDN = LDN
        inter_LDNP = LDNP
        inter_M = M
        inter_N = N
        inter_NP = NP
        inter_NQ = NQ
       
        !!!!REAL array
        inter_BETA = BETA
        inter_F = F
        inter_FJACB = FJACB
        inter_FJACD = FJACD
        inter_XPLUSD = XPLUSD
        !!!!INTEGER array
        inter_IFIXB = IFIXB
        inter_IFIXX = IFIXX
          
   end subroutine inter_func
end subroutine wrapper_ODR

其中subroutine inter_func是将上面函式subroutine FCN里，将使用iso_cbinding的变量值指派给inter_func里的变量用
C++程序代码如下:

[C++] 纯文本查看 复制代码

#include "stdafx.h"
#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;
//宣告ODR呼叫接口
extern "C" {

        void wrapper_ODR(void(*)(int*, int*, int*, int*, int*, int*, int*, \
                double [2] ,double [1][4],int [2], int [1][4],int* ,\
                int*, double [1][4],double [1][2][4],double [1][1][4] ,int* ),\

                int N,int M,int NP,int NQ,double BETA[],double Y[][4],double X[][4],\
                double DELTA[][4],double WE[],double WD[],int IFIXB[],int IFIXX[],\
                int *JOB,int *NDIGIT,double *TAUFAC,double *SSTOL, double *PARTOL,\
                int *MAXIT, int *IPRINT, int *LUNERR, int *LUNRPT,double STPB[],\
                double STPD[], double SCLB[], double SCLD[], double WORK[], double IWORK[],\
                int *INFO, double LOWER[], double UPPER[]);

}
//自定义函数
void FCN(int *N, int *M, int *NP, int *NQ, int *LDN, int *LDM, int *LDNP,\
        double BETA[2], double XPLUSD[1][4], int IXIFB[], int IFIXX[1][4],int *LDIFX,\
        int *IDEVAL, double F[1][4], double FJACB[1][2][4], double FJACD[1][1][4],int *ISTOP){

         *ISTOP = 2;
        
        if (fmod(*IDEVAL,10)!=0) {
                for (int i = 0; i < 4; i++) {
                        F[0][i] = BETA[0] * XPLUSD[0][i] + BETA[1]; //Fitting model BETA[0]*x+BETA[1]
                }

        }

}


int main(){

        int NP = 2, N = 4, M = 1, NQ = 1;
        double BETA[2] = { 2.0, 0.5 };
        double X[1][4] = { 0.0, 1.0, 2.0, 3.0 };
        double Y[1][4] = { 2.0, 5.0, 8.0, 11.0 };

        wrapper_ODR(FCN, N, M, NP, NQ, BETA, Y, X, NULL, NULL, NULL, NULL, NULL, NULL, NULL, NULL, \
                NULL, NULL, NULL, NULL, NULL, NULL, NULL,NULL, NULL, NULL, NULL, NULL, NULL, \
                NULL, NULL);

        system("pause");
        return 0;

}

根据ODRPACK文件定义数组大小如下:

BETA[NP],XPLUSD[M][LDN], IXIFB[NP], IFIXX[M][LDIFX], F[NQ][LDN], FJACB[NQ][LDNP][LDN], FJACD[NQ][LDM][LDN]
NP = 2, N = 4,M = 1, NQ = 1, LDN = 4, LDM = 1, LDNP = 2,LDIFX = 4

执行结果

NP, N, M, NQ, BETA, X, Y皆成功传到ODRPACK去，但它并没有算出正确的BETA值，应该说它好像没有执行，且红线部分，
结果显示ISTOP值却是0，但我在C下定义ISTOP值应为2，似乎我在C下面的自定义义函数的变量值没有传到Fortran去。
不知该怎么解决这问题?需要各位帮忙指正一下，谢谢!

Radiosan · 发表于 2016-12-6 12:37:00

pasuka 发表于 2016-12-6 08:32
为啥不考虑Python，直接吃现成饭
Orthogonal distance regression (scipy.odr) — SciPy v0.18.1 Reference ...

您好! 在一年前我老师丢给我ODR要我能让他在VC++下可以用，在这期间我发现可以用Python，我也进入Python的大世界，用的很愉快，但是我却把给当初给老师交代的是丢在后头了。到了最近，发觉自己几乎都是自己搞自己的事，说来也惭愧，老师交代的事没一件事完成的。所以希望能把ODR这事情给完成，当然，虽然有ISO_C_BINDING，但是做接口这事情，仍然很复杂。从一个月前开始摸Fortran到会ISO_C_BINDING让C去呼叫，一路上跌跌撞撞，几乎都是上网发问，慢慢查资料，一步一步整理出头绪来。但到现在已经想不到该怎么办了?偶然之下，发现这个论坛，因此在这发文希望有人能给个指点。

[混编] C函式当参数回传到Fortran Dll