带内部 goto 的 Fortran do 循环
我有一个 Fortran77 代码片段,如下所示: DO 1301 N=NMLK-2,2,-1 Some code... IF(NB1(N).EQ.50) GOTO 1300 Some code... IF(BS(N).EQ.0.0) GOTO 130…
纯 Fortran 过程中的 I/O
我正在尝试将错误检查合并到我正在编写的纯过程中。我想要这样的东西: pure real function func1(output_unit,a) implicit none integer :: a, outpu…
从 COO 转换为 CSR 稀疏矩阵格式时对重复值求和
从 COO 格式转换为 CSR 时,如何有效地总结重复值。 fortran 子例程中是否存在类似于 scipy 实现(http://docs.scipy.org/doc/scipy-0.9.0/reference/…
通过 Python 代码使用 F2Py
如何通过Python使用F2Py而不是在命令行中调用它? 我正在使用 便携式 Python,但无法安装 Python。因此,如果我解压缩包并将所有内容复制到文件夹“si…
BLAS 2 的 C 与 Fortran
我有一个应用程序,需要执行大量范数、点积,最重要的是矩阵向量乘法。 矩阵和向量巨大。矩阵维度往往是 100000x100000 循环结构是: while(condition…
Python 和 Fortran 语言中的托马斯算法
Python 代码: import math import numpy n = input('Enter the dimension') print 'Matrix size',n dd = [ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,11,12,…
错误:(1) 处的数组引用排名不匹配 (2/1)
我恳请您对我不断得到的这段代码提供帮助 错误:(1) (2/1) 处的数组引用排名不匹配。我的目标是遍历立方体中的每个点(p = i+(j-1)*N + (k-1)NN)并计…
Fortran 变量类型转换
快问。 我有一个带有变量声明的 fortran77 子例程, DIMENSIONS HH(13, 1000) 我假设由于没有指定类型,因此该变量是一个整数数组。稍后在程序中,我…
用 Fortran 将复数矩阵写入文件
如何将 Fortran 中的复数 (n×n) 矩阵写入文件? 例如: DO I=1,N write(14,'(100g15.5)') ( M(i,j), j=1,n ) ENDDO 在此示例中,将 2n×n 个元素写入文…
FORTRAN 中子程序的指针
可能的重复: Fortran 中的函数指针数组 如何在 Fortran 中为函数名称添加别名 在 FORTRAN 中,如何我可以创建并使用一个指向子例程的指针吗? 此外…
使用静态链接对象中的模块 (Fortran)
亲爱的 stackoverflow 专家, 我想(使用 FORTRAN)创建一个共享对象文件,比如说 myso.so。在这个 .so 中,我有一个模块文件 modso.f 和一个“子例程…
FORTRAN 77 程序比 Fortran 90 程序更快吗?
今天,我正在阅读一些用 FORTRAN 77 编写的非常流行的数值库中的代码,例如 QUADPACK(最后更新于 1987 年),我想知道是否有任何理由不在 Fortran 90…
在 C++ 中传递多维数组对于 Fortran 程序员
我习惯用Fortran编写程序(这是f90),但我需要用C++编写程序。我对如何从 C++ 函数传递多维数组感到非常困惑。举个例子,我想在子程序中读取 XYZ 格…
GotoBLAS2 性能
我有一些代码,使用 LAPACK 例程 DPPTRF、DPPTRI 和 DSPMV 执行压缩对称矩阵求逆和乘法。 这里是一个较旧的主题,您可以在其中看到我用来调用 LAPACK …
MinGW:与 LAPACK 和 BLAS 链接会导致 C++未处理的异常
情况很简单,但是很奇怪。当我在没有 LinearAlgebra.o 源代码(需要链接到 LAPACK)的情况下编译程序时,会捕获并处理 C++ 异常。当我不包含该编译单…