Fortran Coder

查看: 24609|回复: 11
打印 上一主题 下一主题

[混编] C++与fortran混编 subroutine和function在module里

[复制链接]

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
跳转到指定楼层
楼主
发表于 2017-7-12 21:03:44 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
6F 币
先贴代码:
[Fortran] 纯文本查看 复制代码
module hydromodelUtil
  implicit none
  double precision, parameter :: pi=3.141593d0
  type singleinfo
    integer :: controlno
    integer :: unittype
    integer :: solid
    double precision :: coearea
    double precision :: volume
    double precision, dimension(3) :: center
    integer :: pointcount, trianglecount
    double precision, dimension(:,:), allocatable :: pointchain
    double precision, dimension(:,:,:), allocatable :: trianglechain
  end type singleinfo
contains
  subroutine cylindermeshing(r, h, cylinderinfo)
    implicit none
    double precision :: r, h
    type(singleinfo) :: cylinderinfo
    integer :: polyedge=60
    double precision, dimension(:,:), allocatable :: point
    double precision :: dtheta
    integer :: meshi

    dtheta=2*pi/real(polyedge)

    cylinderinfo%volume=(0.5d0*polyedge*sin(dtheta)*r**2)*h
    cylinderinfo%center=(/ 0.0d0, 0.0d0, 0.5d0*h /)

    cylinderinfo%pointcount=2*polyedge
    cylinderinfo%trianglecount=4*polyedge
    allocate(point(2,polyedge+1) )

    do meshi=1, polyedge, 1
      point(1,meshi)=r*cos( (real(meshi)-1.0d0)*dtheta)
      point(2,meshi)=r*sin( (real(meshi)-1.0d0)*dtheta)
    end do
    point(:, polyedge+1)=point(:, 1)

    do meshi=1, polyedge, 1
      cylinderinfo%pointchain(1:2,meshi)=point(1:2,meshi)
      cylinderinfo%pointchain(3,meshi)=0.0d0
      cylinderinfo%pointchain(1:2,polyedge+meshi)=point(1:2,meshi)
      cylinderinfo%pointchain(3,polyedge+meshi)=h
    end do

  end subroutine cylindermeshing

  function thirddet(deta)
    double precision, dimension(3,3) :: deta
    double precision :: thirddet
    thirddet=deta(1,1)*(deta(2,2)*deta(3,3)-deta(3,2)*deta(2,3)) &
      & -deta(2,1)*(deta(1,2)*deta(3,3)-deta(3,2)*deta(1,3)) &
      & +deta(3,1)*(deta(1,2)*deta(2,3)-deta(2,2)*deta(1,3))
  end function thirddet

end module hydromodelUtil

请问怎么才能让c++调用生成的dll 里面的子程序和函数啊   我自己改的程序老是出现各种莫名奇妙的错误,拜求大神给看看,多谢了

最佳答案

查看完整内容

1. 我的代码是通用代码。不分动态库,静态库。 2. 用静态库的话,fortran编译得到obj或lib文件,加入C++的项目即可。 3. 用动态库的话,对于 Intel Fortran,需要加一句 !DEC$ ATTRIBUTES DLLEXPORT :: cylindermeshing 4. 动态库可以有隐式调用和显式调用两种用法。这是 VC++ 的事情。我这里的代码适合隐式调用。
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

259

帖子

0

主题

0

精华

版主

World Analyser

F 币
717 元
贡献
510 点

新人勋章美女勋章元老勋章热心勋章规矩勋章管理勋章

QQ
沙发
发表于 2017-7-12 21:03:45 | 只看该作者
1. 我的代码是通用代码。不分动态库,静态库。
2. 用静态库的话,fortran编译得到obj或lib文件,加入C++的项目即可。
3. 用动态库的话,对于 Intel Fortran,需要加一句 !DEC$ ATTRIBUTES  DLLEXPORT  ::   cylindermeshing
4. 动态库可以有隐式调用和显式调用两种用法。这是 VC++ 的事情。我这里的代码适合隐式调用。
回复

使用道具 举报

259

帖子

0

主题

0

精华

版主

World Analyser

F 币
717 元
贡献
510 点

新人勋章美女勋章元老勋章热心勋章规矩勋章管理勋章

QQ
板凳
发表于 2017-7-13 07:51:07 | 只看该作者
module中的函数,和普通的外部函数一样的。直接写 Bind( C , Name ) 就可以了。
你目前的设计,最难的地方在于,你要在 type singleinfo 结构中包含 allocatable 的数组 pointchain ,这是比较困难的。
这一实现,由于fortran和C均没有标准的底层实现。所以标准的混编也就无法实现。

一般的做法是,把allocatable数组单独放在外面当做另外的参数进行传递。fortran不必关心它在C++里是静态的还是动态分配的。

排除这几个数组之外,我写了一个范例给你:

[C++] 纯文本查看 复制代码
#include <iostream>

struct singleinfo{
  int controlno;
  int unittype;
  int solid;
  double coearea;
  double volume;
  double center[3];
  int pointcount , trianglecount;
};

extern "C" {
  void cylindermeshing( double & , double & , struct singleinfo *  );
  double thirddet( double * );
}

using namespace std;

int main()
{
  struct singleinfo cylinderinfo;
  double r = 3.0;
  double h = 4.0;
  cylinderinfo.controlno  = 1;
  cylinderinfo.unittype   = 2;
  cylinderinfo.solid      = 3;
  cylinderinfo.coearea    = 4;
  cylinderinfo.volume     = 5;
  cylinderinfo.unittype   = 6;
  cylinderinfo.pointcount = 7;
  cylinderinfo.trianglecount = 8;
  cylinderinfo.center[0] = 11;
  cylinderinfo.center[1] = 11;
  cylinderinfo.center[2] = 11;
  cylindermeshing(r, h, &cylinderinfo);
  double deta[3][3];
  for (int i = 0; i < 3; i++){
    for (int j = 0; j < 3; j++){
      deta[j][i] = j*4 + i;
    }
  }
  cout << thirddet( &deta[0][0] );
}

[Fortran] 纯文本查看 复制代码
module hydromodelUtil
  implicit none
  double precision, parameter :: pi=3.141593d0
  type , Bind( C ) :: singleinfo
    integer :: controlno
    integer :: unittype
    integer :: solid
    double precision :: coearea
    double precision :: volume
    double precision, dimension(3) :: center
    integer :: pointcount, trianglecount
  end type singleinfo
contains
  subroutine cylindermeshing(r, h, cylinderinfo) Bind(C,Name="cylindermeshing")
    implicit none
    double precision :: r, h
    type(singleinfo) :: cylinderinfo
    write(*,*) r , h , cylinderinfo
  end subroutine cylindermeshing

  function thirddet(deta) Bind( C , Name = "thirddet" )
    double precision, dimension(3,3) :: deta
    double precision :: thirddet
    thirddet=deta(1,1)*(deta(2,2)*deta(3,3)-deta(3,2)*deta(2,3)) &
      & -deta(2,1)*(deta(1,2)*deta(3,3)-deta(3,2)*deta(1,3)) &
      & +deta(3,1)*(deta(1,2)*deta(2,3)-deta(2,2)*deta(1,3))
  end function thirddet

end module hydromodelUtil
回复

使用道具 举报

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
地板
 楼主| 发表于 2017-7-13 14:14:06 | 只看该作者
本帖最后由 cfortran 于 2017-7-13 14:17 编辑

我看你的代码里没有加载动态链接库的代码,现在这样做是不需要做成动态链接库 而采用直接链接的方式混编吗   我用的是ifc的平台   看网上直接混编有个实例是用gfc平台做的。在VS里面具体该怎么做啊这是网上那个gfc平台做的
回复

使用道具 举报

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
5#
 楼主| 发表于 2017-7-13 14:24:34 | 只看该作者
kyra 发表于 2017-7-13 07:51
module中的函数,和普通的外部函数一样的。直接写 Bind( C , Name ) 就可以了。
你目前的设计,最难的地方 ...

我看你的代码里没有加载动态链接库的代码,现在这样做是不需要做成动态链接库 而采用直接链接的方式混编吗   我用的是ifc的平台   看网上直接混编有个实例是用gfc平台做的。在VS里面具体该怎么做啊 ?  这是网上那个gfc平台做的

捕获.PNG (25.36 KB, 下载次数: 1206)

捕获.PNG
回复

使用道具 举报

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
6#
 楼主| 发表于 2017-7-13 15:21:13 | 只看该作者
kyra 发表于 2017-7-12 21:03
1. 我的代码是通用代码。不分动态库,静态库。
2. 用静态库的话,fortran编译得到obj或lib文件,加入C++的 ...

嗯  多谢了    我再试试
回复

使用道具 举报

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
7#
 楼主| 发表于 2017-7-13 15:58:32 | 只看该作者
kyra 发表于 2017-7-12 21:03
1. 我的代码是通用代码。不分动态库,静态库。
2. 用静态库的话,fortran编译得到obj或lib文件,加入C++的 ...

链接时出现错误error LNK1104:无法打开文件‘ifconsol.lib’
我在附加依赖项中添加ifconsol.lib后还是不行。请问这个怎么解决   

ps:我已经把fortran编译得到的source1.obj文件添加到了C++的项目里了
回复

使用道具 举报

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

水王勋章元老勋章热心勋章

8#
发表于 2017-7-13 16:02:20 | 只看该作者
回复

使用道具 举报

14

帖子

4

主题

0

精华

入门

F 币
11 元
贡献
42 点
QQ
9#
 楼主| 发表于 2017-7-13 16:19:30 | 只看该作者
pasuka 发表于 2017-7-13 16:02
为啥不能先看看ivf帮助文档和范例呢?
Intel® Fortran Compiler XE 12.1 User and Reference GuidesMixed  ...

嗯  也在看modern fortran explained 这本书上的介绍   无奈英语渣的不行   只能一个词一个词的查。。。。。。  
回复

使用道具 举报

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

水王勋章元老勋章热心勋章

10#
发表于 2017-7-13 16:22:47 | 只看该作者
cfortran 发表于 2017-7-13 16:19
嗯  也在看modern fortran explained 这本书上的介绍   无奈英语渣的不行   只能一个词一个词的查。。。 ...

ivf帮助文档啃不动,ivf自带的sample 代码总归看得懂吧?!
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-11-15 01:26

Powered by Tencent X3.4

© 2013-2024 Tencent

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