Fortran Coder

查看: 8990|回复: 3
打印 上一主题 下一主题

[混编] 请教关于Fortran调C后返回值的问题

[复制链接]

24

帖子

7

主题

0

精华

熟手

F 币
132 元
贡献
80 点
跳转到指定楼层
楼主
发表于 2017-2-27 16:19:11 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
请教一个关于Fortran调用C后返回值的问题。主程序是Fortran,求解时是调用的C接口函数(lib),现在已经利用C求解器把结果计算出来了,但是却没能把C接口计算的结果返回到Fortran中。程序有点大,没法粘贴上。举个小例子说明一下自己的问题,如下面这个C语言sum函数,a,n是两个参数,s是返回值。利用main函数调用sum函数时,y1=sum(x,5),此时y1就是计算结果,这里的y1相当于sum函数中的s。
[C] 纯文本查看 复制代码
int sum(int a[], int n)
{
int i,s=0;
for(i=0;i<n;i++)
s+=a;
return s;
}

int main(void)
{
int sum(int a[], int n);
int x[5]={1,2,3,4,5},y1,y2;
y1=sum(x,5);
  return 0;
}

通过这样类比,目前,我在Fortran中写的call sum(a, n), 把结果输出到txt是正确的(求解是C完成的,输出语句也是在C中完成的),Fortran中不能写成y1=sum(x,5),所以不知道怎么把C计算的结果返回到Fortran中。而现在我是想在Fortran中想写这样一个语句call sum(a, n,s),a,n是两个参数,s是返回值。而在Fortran与C混合编程过程中要保持共用函数一致,那么C的接口函数也是这种形式的话,就变成 int sum(int a[], int n, int s ), 好像又不合理。请问大家,能帮忙解答一下吗?
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

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

沙发
发表于 2017-2-27 18:04:38 | 只看该作者
C里面数组等价于指针,int a[] 等价于 int *a
具体请参考gfortran的相关说明
https://gcc.gnu.org/onlinedocs/g ... tines-and-Functions

2015

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1574 元
贡献
676 点

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

板凳
发表于 2017-2-27 18:35:41 | 只看该作者
[Fortran] 纯文本查看 复制代码
module interfaces
interface
integer function sum(a,n) Bind( C , Name = "sum" )
  use ,Intrinsic::ISO_C_Binding
  type(C_PTR) , value :: a 
  integer , value :: n
end function sum
end interface
end module interfaces

program main
  use interfaces
  use , intrinsic :: ISO_C_Binding
  implicit none
  integer , target :: x(5) = [1,2,3,4,5]
  write(*,*) "Sum=" , sum( c_loc(x) , size(x) )
end program main

[C] 纯文本查看 复制代码
int sum(int a[], int n)
{
int i,s=0;
for(i=0;i<n;i++)
s+=a[i];
return s;
}




24

帖子

7

主题

0

精华

熟手

F 币
132 元
贡献
80 点
地板
 楼主| 发表于 2017-2-27 18:46:19 | 只看该作者
fcode 发表于 2017-2-27 18:35
[mw_shl_code=fortran,true]module interfaces
interface
integer function sum(a,n) Bind( C , Name = "su ...

谢谢您啦,我尝试一下
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-11-15 06:58

Powered by Tencent X3.4

© 2013-2024 Tencent

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