我想对syevd非连续数组调用 Lapack95 子程序，如下所示：
real :: mat(15000, 15000), vec(15000)
mat=1.d0
associate(eig_vects=>mat(:10000,:10000),eig_vals=>vec(:10000))
call syevd(eig_vects,eig_evals,'V')
end associate
这样做安全吗？之前，我在将不连续的数组切片传递给子例程时遇到了问题（无法分配该子例程创建的临时数组）。使用指针eig_vects作为参数时是否会出现类似的问题？
我正在努力用 C++ 代码填充字符串数组。以下是我尝试做的：
主要.f90：
program main
implicit none
character(len=100) arr(3)
call get_string_array(arr, 100, 3)
write(*,'(3(a,3x))') arr
end
cpp_函数.cpp：
#include <cstring>
#include <algorithm>
using namespace std;
void ConvertToFortran(char* fstring, size_t fstring_len, const char* cstring)
{
size_t inlen = strlen(cstring);
size_t cpylen = min(inlen, fstring_len);
fill(fstring, fstring + fstring_len, ' ');
copy(cstring, cstring + cpylen, fstring);
}
extern "C"
{
void get_string_array_(char** arr, int len, int n)
{
const char* strings[] = {"duma1", "duma2", "duma3"};
for (int i = 0; i < n; i++)
ConvertToFortran(arr[i], len, strings[i]);
}
}
编译这两个文件后
ifx -c main.f90
icx -c cpp_func.cpp
ifx main.o cpp_func.o
我收到分段错误。您能帮我如何arr在其中正确打印数组main.f90吗？
*** 编辑 ***
任何能够将字符串数组从 C++ 传递到 Fortran 的解决方案都将受到赞赏。