- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
如何访问 Fortran 中的二维数组,该数组作为 C 函数的指针传递。下面是我的代码
program linkFwithC
use iso_c_binding
implicit none
interface
subroutine my_routine(p,r) bind(c,name='print2')
import :: c_ptr
import :: c_int
type(c_ptr), value :: p
integer(c_int), value :: r
end subroutine
end interface
REAL, POINTER :: PTR2
integer i,j
integer,parameter ::n=3
real (c_double), allocatable, target :: xyz(:,:)
real (c_double), target :: abc(3,3)
type(c_ptr) :: cptr
allocate(xyz(n,n))
cptr = c_loc(xyz(1,1))
!Inputing array valyes
xyz(1,1)= 1
xyz(1,2)= 2
xyz(1,3)= 3
xyz(2,1)= 4
xyz(2,2)= 5
xyz(2,3)= 6
xyz(3,1)= 7
xyz(3,2)= 8
xyz(3,3)= 9
call my_routine(cptr,n)
do j=1,n
do i=1,n
print*,"xyz(i,j)",i,j,xyz(j,i)
enddo
enddo
deallocate(xyz)
! pause
end program linkFwithC
下面是我的C代码
#include <stdio.h>
void print2(double *p, int n)
{
printf("Array from C is \n");
double *dptr;
int i,j;
dptr = (double *)p;
for ( i = 0; i < n; i++)
{
for ( j = 0; j<n; j++)
{
printf("%.6g \t",dptr[i*n+j]);
*(p+i*n+j)=1;
printf("\n");
}
}
}
下面是输出
Array from C is
1 4 7
1 5 8
1 6 9
xyz(i,j) 1 1 1.00000000000000
xyz(i,j) 2 1 1.00000000000000
xyz(i,j) 3 1 1.00000000000000
xyz(i,j) 1 2 4.00000000000000
xyz(i,j) 2 2 5.00000000000000
xyz(i,j) 3 2 6.00000000000000
xyz(i,j) 1 3 7.00000000000000
xyz(i,j) 2 3 8.00000000000000
xyz(i,j) 3 3 9.00000000000000
*** glibc detected *** ./main.exe: free(): invalid next size (normal): 0x000000000093c290 ***
======= Backtrace: =========
/lib64/libc.so.6[0x320ea75f4e]
/lib64/libc.so.6[0x320ea78cf0]
./main.exe[0x408786]
./main.exe[0x4240bc]
./main.exe[0x429f54]
./main.exe[0x402e63]
/lib64/libc.so.6(__libc_start_main+0xfd)[0x320ea1ed5d]
./main.exe[0x402d69]
======= Memory map: ========
我无法在 Fortran 中打印修改后的值。谁能提出工 Assets 出的可能原因是什么?
最佳答案
首先,你的代码不是FORTRAN,而是Fortran; C 语言互操作性设施的出现要晚得多。二、如下代码
#include <stdio.h>
void print2(double *p, int n)
{
double *dptr;
int i,j;
/* Why are you making a copy here?
* Did you intent to pass by value or reference?
*/
dptr = p;
printf("Array from C is \n");
for(i = 0; i < n; i++){
for(j = 0; j < n; j++){
printf("%.6g \t", dptr[i*n+j]);
*(p+i*n+j)=1;
printf("\n");
}
}
}
module mymod
use ISO_C_binding, only: &
c_ptr, &
c_int
! Explicit typing only
implicit none
interface
subroutine my_routine(p, r) bind(c, name='print2')
import :: c_ptr, c_int
type(c_ptr), value :: p
integer(c_int), value :: r
end subroutine my_routine
end interface
end module mymod
program main
use ISO_Fortran_env, only: &
compiler_version, &
compiler_options
use ISO_C_binding, only: &
c_double, &
c_ptr, &
c_loc
use mymod, only: &
my_routine
! Explicit typing only
implicit none
integer i, j
integer, parameter :: N = 3
real(c_double), allocatable, target :: xyz(:,:)
type(c_ptr) :: cptr
! Allocate memory
allocate(xyz(N, N))
! Get C-language address
cptr = c_loc(xyz(1,1))
! Inputting array values
xyz(1,1)= 1
xyz(1,2)= 2
xyz(1,3)= 3
xyz(2,1)= 4
xyz(2,2)= 5
xyz(2,3)= 6
xyz(3,1)= 7
xyz(3,2)= 8
xyz(3,3)= 9
call my_routine(cptr, N)
do j=1, N
do i=1, N
print *, "xyz(i,j)", i, j, xyz(j, i)
end do
end do
! Release memory
deallocate(xyz)
print '(/4a/)', &
' This file was compiled using ', compiler_version(), &
' using the options ', compiler_options()
end program main
产量
gfortran -Wall -o main.exe print2.c mymod.f90 main.f90
./main.exe
Array from C is
1
4
7
2
5
8
3
6
9
xyz(i,j) 1 1 1.0000000000000000
xyz(i,j) 2 1 1.0000000000000000
xyz(i,j) 3 1 1.0000000000000000
xyz(i,j) 1 2 1.0000000000000000
xyz(i,j) 2 2 1.0000000000000000
xyz(i,j) 3 2 1.0000000000000000
xyz(i,j) 1 3 1.0000000000000000
xyz(i,j) 2 3 1.0000000000000000
xyz(i,j) 3 3 1.0000000000000000
This file was compiled using GCC version 6.1.1 20160802 using the options -mtune=generic -march=x86-64 -Wall
关于c - 在 Fortran 中访问从 C 作为指针传递的二维数组,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/39799557/
我正在尝试创建一个包含 int[][] 项的数组 即 int version0Indexes[][4] = { {1,2,3,4}, {5,6,7,8} }; int version1Indexes[
我有一个整数数组: private int array[]; 如果我还有一个名为 add 的方法,那么以下有什么区别: public void add(int value) { array[va
当您尝试在 JavaScript 中将一个数组添加到另一个数组时,它会将其转换为一个字符串。通常,当以另一种语言执行此操作时,列表会合并。 JavaScript [1, 2] + [3, 4] = "
根据我正在阅读的教程,如果您想创建一个包含 5 列和 3 行的表格来表示这样的数据... 45 4 34 99 56 3 23 99 43 2 1 1 0 43 67 ...它说你可以使用下
我通常使用 python 编写脚本/程序,但最近开始使用 JavaScript 进行编程,并且在使用数组时遇到了一些问题。 在 python 中,当我创建一个数组并使用 for x in y 时,我得
我有一个这样的数组: temp = [ 'data1', ['data1_a','data1_b'], ['data2_a','data2_b','data2_c'] ]; // 我想使用 toStr
rent_property (table name) id fullName propertyName 1 A House Name1 2 B
这个问题在这里已经有了答案: 关闭13年前。 Possible Duplicate: In C arrays why is this true? a[5] == 5[a] array[index] 和
使用 Excel 2013。经过多年的寻找和适应,我的第一篇文章。 我正在尝试将当前 App 用户(即“John Smith”)与他的电子邮件地址“jsmith@work.com”进行匹配。 使用两个
当仅在一个边距上操作时,apply 似乎不会重新组装 3D 数组。考虑: arr 1),但对我来说仍然很奇怪,如果一个函数返回一个具有尺寸的对象,那么它们基本上会被忽略。 最佳答案 这是一个不太理
我有一个包含 GPS 坐标的 MySQL 数据库。这是我检索坐标的部分 PHP 代码; $sql = "SELECT lat, lon FROM gps_data"; $stmt=$db->query
我需要找到一种方法来执行这个操作,我有一个形状数组 [批量大小, 150, 1] 代表 batch_size 整数序列,每个序列有 150 个元素长,但在每个序列中都有很多添加的零,以使所有序列具有相
我必须通过 url 中的 json 获取文本。 层次结构如下: 对象>数组>对象>数组>对象。 我想用这段代码获取文本。但是我收到错误 :org.json.JSONException: No valu
enter code here- (void)viewDidLoad { NSMutableArray *imageViewArray= [[NSMutableArray alloc] init];
知道如何对二维字符串数组执行修剪操作,例如使用 Java 流 API 进行 3x3 并将其收集回相同维度的 3x3 数组? 重点是避免使用显式的 for 循环。 当前的解决方案只是简单地执行一个 fo
已关闭。此问题需要 debugging details 。目前不接受答案。 编辑问题以包含 desired behavior, a specific problem or error, and the
我有来自 ASP.NET Web 服务的以下 XML 输出: 1710 1711 1712 1713
如果我有一个对象todo作为您状态的一部分,并且该对象包含数组列表,则列表内部有对象,在这些对象内部还有另一个数组listItems。如何更新数组 listItems 中 id 为“poi098”的对
我想将最大长度为 8 的 bool 数组打包成一个字节,通过网络发送它,然后将其解压回 bool 数组。已经在这里尝试了一些解决方案,但没有用。我正在使用单声道。 我制作了 BitArray,然后尝试
我们的数据库中有这个字段指示一周中的每一天的真/假标志,如下所示:'1111110' 我需要将此值转换为 boolean 数组。 为此,我编写了以下代码: char[] freqs = weekday
我是一名优秀的程序员,十分优秀!