Fortran Coder

标题: 结构体传出来的数组会改变 [打印本页]

作者: zjk0112 时间: 2022-9-14 13:33
标题: 结构体传出来的数组会改变
本帖最后由 zjk0112 于 2022-9-14 13:35 编辑

我用MKl中CSR格式的结构体遇到了如下：在调用subroutine identity_matrix_sparse1时，打印出来的n为3，但export传出来三元组中的CSR_val值异常（图2）；若在identity_matrix_sparse1中重新赋值常量nn=3替换所有n，则export传出的三元组中的CSR_VAl的值正常（图1）。
我的问题：1.这里面哪里影响了n.2在subroutine中创建结构体，再将出来是否会破坏结构体。3.这里面涉及到C到Fortran的指针转换，会不会对常量有影响？

2.png (10.18 KB, 下载次数: 145)

1.png (9.2 KB, 下载次数: 147)

BadFile.f90

4.45 KB, 下载次数: 3

作者: zjk0112 时间: 2022-9-14 13:52
我又尝试了，如果用动态数组则替代固定数组则可以避免这样的问题。

作者: zjk0112 时间: 2022-9-20 16:58
经过了一周的研究和尝试，我终于找到了错误产生的原因：在subroutine identity_matrix_sparse1产生的结构体需要用到局部变量，而局部变量在默认会清空，导致结构体错误。所以，解决的方法也很简单，只需要在子例程subroutine identity_matrix_sparse1中加save语句即可。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)