我正在使用 R,但需要做大量的数字运算,我想在 Fortran 中做。我对 R 比较陌生,对 Fortran 来说是个新手...我已经有一个可用的 R 程序,我想对其进行优化。我创建了一个 Fortran 程序来求解 ODE 系统,并将其保存在子例程中。我还使用一个名为 aux.f90 的模块来存储参数和一个函数,该函数创建一个输入方程的信号。这按预期工作,数据保存在 .txt 文件中。
我现在想做的是创建一个 R 前端,向 Fortran 程序传递参数,例如模拟的长度或解决方案中使用的步骤数。然后 Fortran 完成繁重的工作,将结果保存在文件中,我可以使用 R 可视化文件中的数据。请参阅下面的 Fortran 代码:
! The auxiliary module contains all parameters
module aux
implicit none
integer,parameter :: n = 2**8 ! number of steps
real(kind=4) :: jout = 0.5 ! for normal metabolism
real(kind=4) :: alp = 4.0 ! factor for growth
real(kind=4) :: bet = 1.0 ! benefit value
real(kind=4) :: etay = 0.1 …Run Code Online (Sandbox Code Playgroud)