TobiR's questions

我有以下示例：

LIST_JOB = job1 job2 job3

define macro
$(info "JOB: $(1)")
job_$(1): job.f90
   ifort -c job.f90
   $(shell touch job.f90)
endef

all: rule

rule: $(foreach job,$(LIST_JOB), job_$(job))

$(foreach job,$(LIST_JOB), \
   $(eval $(call macro,$(job))) \
)

首先打印$(info "JOB: $(1)")行，然后打印编译命令：

"JOB: job1"
"JOB: job2"
"JOB: job3"
ifort -c job.f90
ifort -c job.f90
ifort -c job.f90

有没有办法在编译之前打印标题行(当然，同时保持macro循环foreach)？这意味着以下顺序：

"JOB: job1"
ifort -c job.f90
"JOB: job2"
ifort -c job.f90
"JOB: job3"
ifort -c job.f90

我有以下 Makefile 示例：

BUILD = build
PATH_POT3 = src/pot/pot3/
MPIFC = mpiifort
FFLAGS = -O3

MOL_LIST = $(shell find $(PATH_POT3) -mindepth 1 -maxdepth 1 -type d)

$(foreach MOL,$(MOL_LIST), \
  $(eval LIST_MOL := $(shell find $(MOL) -name "*.f90")) \
  $(info Fortran files: $(LIST_MOL)) \
  $(foreach SRC,$(LIST_MOL), \
    $(eval OBJ := $(BUILD)/$(notdir $(SRC:.f90=.o))) \
    $(OBJ): $(SRC) ; \
    $(info Compiling: $(SRC)) \
    $(info Output: $(OBJ)) \
    $(MPIFC) $(FFLAGS) -c $(SRC) -o $(OBJ) \
  ) \
)

all: $(foreach MOL,$(MOL_LIST),\
    $(foreach SRC,$(shell find $(MOL) -name "*.f90"),\
        $(BUILD)/$(notdir $(SRC:.f90=.o))\
    )\
)

这里我想PATH_POT3以顺序的方式（一个接一个地）编译每个子目录的（递归）内容。遗憾的是，这种构造并没有提供预期的结果，因为它只编译列表中的第一个目标文件，而跳过了所有其他目标文件。Makefile 为我的测试用例打印了以下信息：

Fortran files: src/pot/pot3/ar3/pot_ar3.f90 src/pot/pot3/ar3/ar43/pot_ar43.f90
Compiling: mpiifort -O3 -c src/pot/pot3/ar3/pot_ar3.f90 -o build/pot_ar3.o
Output: build/pot_ar3.o
Compiling: mpiifort -O3 -c src/pot/pot3/ar3/ar43/pot_ar43.f90 -o build/pot_ar43.o
Output: build/pot_ar43.o
Fortran files: src/pot/pot3/ar33/pot_ar33.f90 src/pot/pot3/ar33/pot_ar333.f90 src/pot/pot3/ar33/ar53/pot_ar53.f90
Compiling: mpiifort -O3 -c src/pot/pot3/ar33/pot_ar33.f90 -o build/pot_ar33.o
Output: build/pot_ar33.o
Compiling: mpiifort -O3 -c src/pot/pot3/ar33/pot_ar333.f90 -o build/pot_ar333.o
Output: build/pot_ar333.o
Compiling: mpiifort -O3 -c src/pot/pot3/ar33/ar53/pot_ar53.f90 -o build/pot_ar53.o
Output: build/pot_ar53.o
make: 'build/pot_ar3.o' is up to date.

一切似乎都设置正确，但编译仅针对执行build/pot_ar3.o。

我应该如何更改此文件以获得所需的行为？

我想对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 的解决方案都将受到赞赏。