Fortran Coder

查看: 151|回复: 4

[混编] Fortran调用C返回实数数组问题

[复制链接]

24

帖子

7

主题

0

精华

熟手

F 币
130 元
贡献
79 点
发表于 2017-5-7 22:26:35 | 显示全部楼层 |阅读模式
您好,请教一个Fortran调用C返回实数数组问题。之前群主在论坛里写的一个Fortran调用C返回一个数的例题,我想问问,如何返回一个数组呢?在群主之前的例题上修改了一点,结果还是有问题。我只是想实现在对数组中每个元素加1的结果再返回到Fortran中来,结果是只实现了第一个元素加1,别的不行
C程序:
[C] 纯文本查看 复制代码
#include <stdlib.h>
#include <stdio.h>
#include <string.h> 
double c_add(double a[], int n)
{
int i;
double s=0;
double *xx ;
xx = calloc(n,sizeof(double));
for(i=0;i<n;++i)
xx=a+1;
printf("n/", xx );
return *xx;
}

Fortran程序:
[Fortran] 纯文本查看 复制代码
module interfaces
interface
real(8) function add(a,n) Bind( C , Name = "c_add" )
  use ,Intrinsic::ISO_C_Binding
  type(C_PTR) , value :: a 
  integer , value :: n
end function add
end interface
end module interfaces

program main
  use interfaces
  use , intrinsic :: ISO_C_Binding
  implicit none
  integer::i;
  real(8) , target :: x(5) = [0.1,0.2,0.3,0.4,0.5]
  do i=1,size(x)
  write(*,*)  add( c_loc(x) , size(x) ) 
  end do
end program main
回复

使用道具 举报

1096

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂

F 币
847 元
贡献
746 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2017-5-7 22:43:22 | 显示全部楼层
混编时,请尽量遵循以下原则:
谁申请,谁释放;谁打开,谁关闭;谁构造,谁析构。

因此,你这样在C里面calloc一个数组,返回给fortran。先不管行不行得通,这样做非常容易导致内存泄漏(因为没法释放)

所以,最佳的方法是:由fortran定义(或分配)数组,传递给C语言,C语言对其赋值。然后fortran获得对应的赋值并继续使用。最终由定义者(或分配者)也就是fortran来释放。

334

帖子

2

主题

0

精华

大宗师

F 币
2306 元
贡献
1400 点

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

发表于 2017-5-8 09:07:40 | 显示全部楼层
本帖最后由 pasuka 于 2017-5-8 09:09 编辑

改成传递函数指针
参看,(没有C和C++基础的话,多花点功夫,两个月时间应该能搞明白)
The GNU Fortran Compiler: Working with Pointers
https://gcc.gnu.org/onlinedocs/g ... rking-with-Pointers

24

帖子

7

主题

0

精华

熟手

F 币
130 元
贡献
79 点
 楼主| 发表于 2017-5-8 09:22:33 | 显示全部楼层
fcode 发表于 2017-5-7 22:43
混编时,请尽量遵循以下原则:
谁申请,谁释放;谁打开,谁关闭;谁构造,谁析构。

好的,我试试,谢谢啦

24

帖子

7

主题

0

精华

熟手

F 币
130 元
贡献
79 点
 楼主| 发表于 2017-5-8 09:22:52 | 显示全部楼层
pasuka 发表于 2017-5-8 09:07
改成传递函数指针
参看,(没有C和C++基础的话,多花点功夫,两个月时间应该能搞明白)
The GNU Fortran Co ...

嗯嗯,我再看看
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|QQ群|Fcode

GMT+8, 2017-8-22 03:35

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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