Fortran调用C返回实数数组问题

dongtian99 · 发表于 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

fcode · 发表于 2017-5-7 22:43:22

混编时，请尽量遵循以下原则：
谁申请，谁释放；谁打开，谁关闭；谁构造，谁析构。

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

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

pasuka · 发表于 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

dongtian99 · 发表于 2017-5-8 09:22:33

fcode 发表于 2017-5-7 22:43
混编时，请尽量遵循以下原则：
谁申请，谁释放；谁打开，谁关闭；谁构造，谁析构。

好的，我试试，谢谢啦

dongtian99 · 发表于 2017-5-8 09:22:52

pasuka 发表于 2017-5-8 09:07
改成传递函数指针
参看，（没有C和C++基础的话，多花点功夫，两个月时间应该能搞明白）
The GNU Fortran Co ...

嗯嗯，我再看看

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

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章