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

dongtian99

您好，请教一个Fortran调用C返回实数数组问题。之前群主在论坛里写的一个Fortran调用C返回一个数的例题，我想问问，如何返回一个数组呢？在群主之前的例题上修改了一点，结果还是有问题。我只是想实现在对数组中每个元素加1的结果再返回到Fortran中来，结果是只实现了第一个元素加1，别的不行
C程序：

[C] 纯文本查看 复制代码

查看代码

打印?

01	#include <stdlib.h>

02	#include <stdio.h>

03	#include <string.h>

04	double c_add(double a[], int n)

05

{

06

int i;

07	double s=0;

08	double *xx ;

09	xx = calloc(n,sizeof(double));

10	for(i=0;i<n;++i)

11

xx=a+1;

12	printf("n/", xx );

13	return *xx;

14

}

Fortran程序：

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	module interfaces

02

interface

03	real(8) function add(a,n) Bind( C , Name = "c_add" )

04	use ,Intrinsic::ISO_C_Binding

05	type(C_PTR) , value :: a

06	integer , value :: n

07	end function add

08	end interface

09	end module interfaces

10

11	program main

12	use interfaces

13	use , intrinsic :: ISO_C_Binding

14	implicit none

15	integer::i;

16	real(8) , target :: x(5) = [0.1,0.2,0.3,0.4,0.5]

17	do i=1,size(x)

18	write(*,*)  add( c_loc(x) , size(x) )

19

end do

20	end program main

dongtian99

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

嗯嗯，我再看看

dongtian99

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

好的，我试试，谢谢啦

pasuka

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

fcode

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

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

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