Fortran Coder

标题: Fortran调用C返回实数数组问题 [打印本页]

作者: dongtian99    时间: 2017-5-7 22:26
标题: Fortran调用C返回实数数组问题
您好,请教一个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

作者: fcode    时间: 2017-5-7 22:43
混编时,请尽量遵循以下原则:
谁申请,谁释放;谁打开,谁关闭;谁构造,谁析构。

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

所以,最佳的方法是:由fortran定义(或分配)数组,传递给C语言,C语言对其赋值。然后fortran获得对应的赋值并继续使用。最终由定义者(或分配者)也就是fortran来释放。
作者: pasuka    时间: 2017-5-8 09:07
本帖最后由 pasuka 于 2017-5-8 09:09 编辑

改成传递函数指针
参看,(没有C和C++基础的话,多花点功夫,两个月时间应该能搞明白)
The GNU Fortran Compiler: Working with Pointers
https://gcc.gnu.org/onlinedocs/g ... rking-with-Pointers
作者: dongtian99    时间: 2017-5-8 09:22
fcode 发表于 2017-5-7 22:43
混编时,请尽量遵循以下原则:
谁申请,谁释放;谁打开,谁关闭;谁构造,谁析构。

好的,我试试,谢谢啦
作者: dongtian99    时间: 2017-5-8 09:22
pasuka 发表于 2017-5-8 09:07
改成传递函数指针
参看,(没有C和C++基础的话,多花点功夫,两个月时间应该能搞明白)
The GNU Fortran Co ...

嗯嗯,我再看看




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2