如何在fortran中读取.mat
在Fortran中读取.mat文件，可以使用MATLAB的引擎库来实现。MATLAB引擎库提供了Fortran的接口，使得Fortran程序可以直接读取和处理MATLAB的.mat文件。
以下是在Fortran中读取.mat文件的步骤：
- 首先，确保已经安装了MATLAB软件，并且配置了MATLAB的编译器和库文件路径。
- 在Fortran程序中，使用
use, intrinsic :: iso_c_binding语句导入ISO_C_BINDING模块，以便使用C语言的数据类型和函数。 - 声明MATLAB引擎库的函数和数据类型，例如：
! 声明MATLAB引擎库的函数和数据类型
interface
subroutine engOpen(arg1) bind(C, name="engOpen")
use iso_c_binding
character(kind=c_char), dimension(*) :: arg1
end subroutine engOpen
subroutine engClose(ep) bind(C, name="engClose")
use iso_c_binding
type(c_ptr), value :: ep
end subroutine engClose
subroutine engGetVariable(ep, name, array) bind(C, name="engGetVariable")
use iso_c_binding
type(c_ptr), value :: ep
character(kind=c_char), dimension(*) :: name
type(c_ptr), value :: array
end subroutine engGetVariable
end interface- 在主程序中，打开MATLAB引擎并读取.mat文件中的变量。以下是一个示例代码：
program read_mat_file
use iso_c_binding
implicit none
character(kind=c_char), dimension(1:100) :: mat_file = "example.mat"
character(kind=c_char), dimension(1:100) :: var_name = "my_variable"
type(c_ptr) :: ep
type(c_ptr) :: array
! 打开MATLAB引擎
call engOpen(mat_file)
! 读取.mat文件中的变量
call engGetVariable(ep, var_name, array)
! 在这里可以对读取到的变量进行处理
! 关闭MATLAB引擎
call engClose(ep)
end program read_mat_file在上述示例代码中，mat_file变量指定了.mat文件的路径，var_name变量指定了要读取的变量名。engOpen函数用于打开MATLAB引擎，engGetVariable函数用于读取变量，engClose函数用于关闭MATLAB引擎。
需要注意的是，上述代码只是一个简单的示例，实际应用中可能需要根据具体情况进行适当的修改和扩展。
推荐的腾讯云相关产品：腾讯云弹性计算（Elastic Compute）产品，提供了弹性的计算资源，适用于各种计算任务和应用场景。您可以通过以下链接了解更多信息：腾讯云弹性计算产品介绍。
相关·内容
谁能告诉我怎么用fortran来读.mat？我需要使用fortran来处理来自matlab的二维矩阵。我使用
open(unit = 9006, file = 'zmatrix.mat',status = 'unknown', form = 'unformatted')
do i=1,296
do j=1,278
read(9006,*) zmatrix(i,j)
end do
end do
但据说fortran运行时错误:格式存在于未格式化的数据传输中。感谢您的任何建议！
浏览 10提问于2017-07-07得票数 1
我有一个很长(大约1000行)的复杂的Matlab函数(比如mat_test.m)，它可以很好地工作。另一方面，我有一组复杂的Fortran代码。我正在尝试将Matlab函数插入其中一段Fortran代码中。我对Fortran没有太多的能力，所以我不想用Fortran重写整个Matlab函数。有没有一种方法可以在Fortran中定义变量，并远程调用此Matlab函数并使用它生成的输出？
浏览 0提问于2015-08-06得票数 1
我正在使用带有Visual Studio 2008 SP1的英特尔Fortran。
我的主要问题是:我想从Matlab的.mat文件中读取二维数组到fortran中。另外，将Fortran 2D矩阵的输出保存为.mat文件，因为目前我可以使用以下命令将其保存为文本文件：
write(unit = #, <linelength>F22.8>),matrixname
这条线行得通，但我不确定我的双精度是否松动了。如果我不放松精确度，我可以坚持下去，否则我将需要帮助。我只需要一种方法来从Matlab文件中读取到英特尔-fortran，并保持精度。这些数组中没有字符，它们是数值。
我
浏览 1提问于2018-01-23得票数 0
0回答
我在下面有一个简单的Fortran代码，其中我从一个文本文件中读取一个矩阵，然后将其传递给子例程。但是当子程序读取它时，它会给出一个无用的值。我似乎找不出代码出了什么问题？
program main
implicit none
real, dimension(:,:), allocatable :: mat
open (unit=99, file='inp.txt', status='old', action='read')
allocate(mat(2,2))
read(99, *), mat
CALL fun(mat)
end
subrou
浏览 4提问于2016-07-14得票数 1
回答已采纳
2回答
在Fortran中，可以对数组进行操作，但是如何将派生类型的索引也视为数组的一部分呢？代码可以解释我最想做的事情：
type mytype
integer :: b(3,3)
real :: c(4)
endtype
integer :: a(3,3)
real :: d(2,4)
type(mytype) :: mat(2)
!do stuff so that 'mat' gets values
....
!usually one does this
a = matmul(mat(1)%b, transpose(mat(2)%b))
!multiply
浏览 0提问于2012-01-11得票数 2
回答已采纳
我使用VSCode与现代Fortran+fortls一起用PETSc库编写Fortran代码。使用PETSc关键词的声明似乎根本没有得到承认。我尝试使用PETSc等方法将.fortls库的路径和头文件包含在.fortls文件中，但没有解决这个问题。有人能给我一个提示吗，以便PETSc中的所有关键字和函数声明都可以由fortls处理？非常感谢
我在.fortls文件中使用了以下内容
{
"source_dirs": ["./**","./src/**",
浏览 10提问于2022-11-01得票数 0
回答已采纳
1回答
我尝试在Android4.4中使用.mat加载一个JMatIO文件。
当.mat文件在外部存储中时，它工作得很好，但是当文件位于应用程序的raw文件夹中时，它就不能工作了。目前，应用程序找不到任何*.mat文件。日志只打印消息“未能加载矩阵文件”。
下面是我使用的代码：
try {
mfr = new MatFileReader("android.resource://" + getPackageName() + "/raw/rm.mat");
} catch (IOException e) {
e.printStackTrace
浏览 2提问于2015-07-13得票数 0
回答已采纳
我在一个文件(number.txt)中存储了一个矩阵，如下所示：
12323456
54254311
76534522
我如何在Fortran中读取这样的矩阵，因此结果将是：
1 2 3 2 3 4 5 6
5 4 2 5 4 3 1 1
7 6 5 3 4 5 2 2
使用awk分隔这些列并在Fortran中读取它是非常容易的。但是，我想知道我是否可以只使用Fortran来完成所有这些工作。在我完成之后，我将需要将这个矩阵乘以它的转置。
浏览 13提问于2019-04-03得票数 2
我用ndarray.tobytes(order='F')在Python中保存了一个矩阵。
f_mat = open('mat.dat', 'wb')
f_mat.write(matrix.tobytes('F'))
f_mat.close()
在Fortran中是否有一种快速阅读它的方法，或者我应该做一个for循环并计算元素的顺序吗？
浏览 0提问于2017-06-22得票数 0
回答已采纳
现在，我有一个N*N矩阵，mat。我也有
maskmat=(a(1),a(2),...,a(i),...a(N))。
a(i)等于0或1。
如果为a(i)==0，则不做任何更改。因此，我们可以根据规则maskmat得到一个子矩阵。
如何在Fortran中实现？
浏览 1提问于2014-04-21得票数 0
3回答
Fortran 2003有用于数组连接的方括号语法，Intel fortran编译器也支持它。我在这里为矩阵连接编写了一个简单的代码：
program matrix
implicit none
real,dimension (3,3) :: mat1,mat2
real,dimension(3,6):: mat3
integer i
mat1=reshape( (/1,2,3,4,5,6,7,8,9/),(/3,3/))
mat2=reshape( (/1,2,3,4,5,6,7,8,9/),(/3,3/))
mat3=[mat1,mat2]
!display
do i=1,3,1
wri
浏览 4提问于2015-06-19得票数 3
回答已采纳
我正在使用Fortran的一个旧的有限元程序。这是大学里的一个研究项目。
我有一个文件.raw，它表示3D中的图像。
.raw文件中的数据存储在uint16或uint8中，并且知道整数的总数。
如何将此图像读取到Fortran程序中的整数数组？
类似于
allocate(imgarray(total_int))
call raw2array(filename,imgarray)
我目前正在使用python读取图像，并将其转换为整数向量。此向量用文本文件在Fortran中读取。
In python
imgarray = np.fromfile(fid, dtype=np.uint16,count
浏览 3提问于2018-04-26得票数 1
回答已采纳
我是个Fortran新手。我试图为矩阵运算编写fortran代码，但在设计我的函数的原型(接口)时却陷入了困境。
我的C/C++编程经验告诉我，要编写适用于所有维数的矩阵的代码，需要将数组和维数分别传递给函数。示例(可能不是最好的)：
int * matx_op(int *mat_a, int arows, int acols, int *mat_b, int brows, int bcols);
但在fortran中，matmul()函数似乎是自动执行此操作的。想要了解它是如何做到的。
print *, 'Enter 16 elements of matrix A&#
浏览 2提问于2015-04-09得票数 1
1回答
我正在通过Ubuntu用Fortran编写代码，并使用make文件来运行编译器-当涉及到编译器时，这有点像新手，但我认为这就是问题所在。我做了一些研究，我的问题和这里问的是一样的：。然而，在我的例子中，-C (检查边界)标志一直处于打开状态，并且我最近才遇到这个问题。我一点也没有改变编译器，所以很困惑。举个例子，我写道：
REAL, ALLOCATABLE :: mat_b(:,:)
INTEGER :: i
ALLOCATE(mat_b(5,5))
mat_b = 0.0
DO i = -1, 10
WRITE(*,*) ma
浏览 1提问于2013-05-05得票数 2
1回答
最近，我用python编写了一个耗时的程序，并决定用fortran重写最耗时的部分。
然而，用f2py包装的fortran代码比python代码慢，谁能告诉我这里发生了什么？
作为参考，下面是python函数：
def iterative_method(alpha0, beta0, epsilon0, epsilons0, omega, smearing=0.01, precision=0.01, max_step=20, flag=0):
# alpha0, beta0, epsilon0, epsilons0 are numpy arrays
m, n = np.shape
浏览 3提问于2016-03-12得票数 0
回答已采纳
我正在用Fortran编写一个简单的矩阵向量乘法程序。在这个程序中，我使用了cuBLAS的Fortran绑定例程，解释了。当我的数组包含少于90个元素时，程序就像一个魅力。我读过其他几篇建议使用ulimit -s unlimited或-heap-array [size]的文章，但这些在我的例子中行不通。还有别的方法可以解决这个我不知道的问题吗？我正在用ifort和gfortran编译，使用nvcc编译cublas例程。
ifort
$ifort -heap-arrays 1 -names uppercase -assume nounderscore -c -o cublasf.o main.f
浏览 3提问于2014-03-28得票数 0
回答已采纳
3回答
我用f2py包装了一些fortran代码。以下是fortran代码：
MODULE iteration
implicit none
contains
SUBROUTINE iterate(alpha, beta, e, es, rank, omega, smearing, prec, max_step)
REAL(kind=8), INTENT(in) :: omega, smearing, prec
INTEGER :: max_step, step, rank, cnt
COMPLEX(kind=16) :: alpha(rank,rank), beta(rank
浏览 8提问于2016-03-11得票数 3
回答已采纳
2回答
我在编辑表格时遇到了问题。单击“编辑”按钮无法获得表单值。只有这些值在控制台中可见。更新按钮也不起作用。
prepaid.component.html
<div><button style="margin:10px 0 0 50px;"mat-raised-button color="warn" routerLink="/addPlan"><mat-icon >add_circle_outline</mat-icon>Add Plan</button>
</div>
<
浏览 5提问于2022-03-17得票数 0
回答已采纳
我试图读取这一行中的整数：
# 14 14 10
在Fortran，2008年。
我试图使用以下代码：
read(21, "(A, I,I,I)") garbage, a, b, c
但这不符合标准。英特尔Fortran发出警告"Fortran 2008 does not allow this edit descriptor. [I]“和其他问题解释这个问题：
如何使用Fortran 2008正确读取未知宽度的整数？我不能简单地指定I2，因为我事先不知道整数的宽度。
浏览 3提问于2017-10-04得票数 4
回答已采纳
1回答
当我编译频谱解算器时，我正在尝试使用Cmake来管理Fortran和C代码的依赖关系。但是似乎Cmake找不到编译器，尽管我已经下载了Fortran和GCC。这就是我得到的错误。我希望有人能告诉我如何解决这个问题。 the error screenshot `CMake Error at /usr/share/cmake-3.16/Modules/CMakeDetermineCompiler.cmake:58 (set_property):
set_property could not find CACHE variable CMAKE_Fortran_COMPILER. Perhap
浏览 383提问于2020-12-02得票数 0
我正在尝试使用openMP编写一个fortran应用程序。我在MPI方面有丰富的经验，但仍然在与OMP同步进行斗争。请考虑以下例子：
module ps_mod
contains
subroutine cg(A,b,N)
use iso_fortran_env
implicit none
real(real64) :: A(:,:), b(:)
integer(int32) :: N
integer :: i
real(
浏览 7提问于2021-03-08得票数 0
回答已采纳
我是MPI和Fortran 77。我有一个fortran 77代码FKRPRO.f，我想使用OpenMPI并行化它。代码需要许多参数，这些参数是在运行时从一个单独的文件输入到它的。编译和运行如下所示
gfortran -o FKRPRO FKRPRO.f
./FKRPRO < Modelfile.txt
代码中的等效行(不是我的代码)是
PARAMETER(LIN=5)
INTEGER ERROR
LOGICAL PRNTA
PRNTA=.FALSE.
READ(LIN,'(L3)') PRNTA
R
浏览 1提问于2013-10-03得票数 0
回答已采纳
2回答
我的问题是在Matlab中对数组进行整形。
我在Matlab中读取来自Fortran的"diegm.MAT“文件。这个数组的大小是12x3，我需要一个4x3x3。
我尝试了整形功能，但不起作用。
这是我正在读取的数组：
5 2 5
2 1 2
4 3 2
5 3 3
5 2 4
4 2 3
1 1 3
4 5 1
3 3 1
2 1 4
2 3 1
4 2 4
这是我需要的
浏览 2提问于2018-05-31得票数 0
回答已采纳
我正在尝试从转向fortran和c++。我主要处理的是巨大的图像，坐标方向是x:从左到右，y:自上而下，像素被存储在行中，这是图像格式中常见的。
Numpy说，它也按行存储矩阵。就好而言，它与fortran/c++保存图像是一样的。但是，请参阅下面的示例，灰色从黑色到白色依次排列。让我们创建第一行为123等的图像3x3。
mat = np.array([1,2,3,4,5,6,7,8,9], 'i1').reshape(3,3)
接下来，我使用下面的代码获取有关矩阵的信息。例程打印:在位置x=1, y=0的值，在内存中的位置，看数组是否被复制，大步和迭代器。
def pr(a
浏览 4提问于2022-01-17得票数 1
最近，我在Fortran代码中收到以下错误
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0 0x2AD9B0F8FE08
#1 0x2AD9B0F8EF90
#2 0x2AD9B12D44AF
#3 0x401A3E in MAIN__ at tstreadin.f90:?
我的代码如下
Program www
implicit none
integer ::i,j,rows,cols,row
real(ki
浏览 1提问于2017-12-18得票数 1
我有3个垫标签与垫滑块切换。在第三个mat选项卡上，使用ngIf删除mat-滑块切换。然而，这似乎重置了垫上滑块上的检查值。在前两个选项卡之间来回移动，滑块值保持良好，但第三个选项卡始终重置垫块滑块切换。如果在mat选项卡1上选择了滑块，则移动到mat选项卡3，然后回到mat选项卡1，滑块将移回未选中状态。
选项卡component.html
<div class="container">
<mat-tab-group
(selectedTabChange)="focusChange($event)"
[selectedIndex]=
浏览 21提问于2018-07-21得票数 1
回答已采纳
我必须阅读我在这篇文章末尾写的简单的文本文件(它只是一个结构化的网格)。在fortran中很容易做到这一点，你只需要做：
read(fileunit,*)
read(fileunit,*) mc,nc
do j = 1, nc
read (fileunit, *) dummy, dummy, (xcor(j,i), i=1,mc)
enddo
在matlab中有没有一个等价的函数，可以像fortran一样逐个读取元素，并在换行符之后继续读取？我找不到它，所有的函数，如fscanf，textscan等逐行读取，然后我必须解析每一行。这是文件。谢
浏览 1提问于2013-02-06得票数 2
回答已采纳
1回答
如何在openGL中沿Z轴画线
、、、、
我被困在Z轴上画了条线。我已经检查了的相关话题，但是即使我改变了相机的位置，我仍然看不到我的线；同时，我可以看到XZ平面上的正方形。
这是我的代码：
int main()
{
glfwInit();
glfwWindowHint(GLFW_CONTEXT_VERSION_MAJOR, 2);
glfwWindowHint(GLFW_CONTEXT_VERSION_MINOR, 1);
GLFWwindow* window = glfwCreateWindow(1024, 768, "some program", NULL, NULL);
浏览 3提问于2022-07-02得票数 0
回答已采纳
2回答
我正在尝试使用CodeBlocks作为我的集成开发环境，GNU Fortran compiler作为我的编译器来运行fengdemo.f。它运行include fine，但它无法构建对matlab引擎的任何方法调用
C:\Program Files\MATLAB\R2012b\extern\examples\eng_mat\fengdemo.F|51|undefined reference to `engopen_'|
C:\Program Files\MATLAB\R2012b\extern\examples\eng_mat\fengdemo.F|58|undefined refe
浏览 8提问于2013-03-19得票数 0
回答已采纳
3回答
我有一个C函数来规范化对数空间中数组的行(这可以防止数值下溢)。
我的C函数的原型如下所示：
void normalize_logspace_matrix(size_t nrow, size_t ncol, double* mat);
您可以看到，它接受一个指向数组的指针，并在适当的位置修改它。当然，C代码假定数据保存为C连续数组，即行连续数组。
我使用Cython将函数包装如下(imports和cdef extern from被省略)：
def normalize_logspace(np.ndarray[np.double_t, ndim=2] mat):
cdef Py_ssize
浏览 1提问于2010-12-12得票数 10
回答已采纳
当将数组的部分传递给Fortran中的子例程时，例如f(a, b, c(2:5,4:6)) (都是2D数组)，程序是否首先创建c的临时副本，然后将其传递给子例程(作为引用、指针等)，还是整个程序都是动态处理的？
我正在尝试将一些Fortran代码转换为C++，并且我看到对具有数组部分作为参数传递的子程序的调用。据我所知，C++不允许这样做，所以我试图在C++ (mat2d = std::vector<std::vector<T>>)中绕过这个问题：
f(mat2d &a, mat2d &b, mat2d *a, int rows, int rows, i
浏览 1提问于2020-03-04得票数 3
回答已采纳
1回答
我是Fortran新手，我需要编写一段代码，允许我读取一个文本文件(.dat)，其中包含一个数据数组(1876年列x3126行)，其中有不规则的空格/不规则的分隔符，并将其写入格式一致的新文件。简而言之，我只是希望数据的列彼此对齐。
例如：
<4 spaces> 0.2341 15.0769 <2 spaces> -10.0672 1.6278 <3 spaces> 66.9851 -14.7123 <3 spaces> -0.4468 -2.5673
上面的例子基本上是数据列的样子，每个值之前对应的空格数如下所示。这些值的格式是“f10
浏览 1提问于2013-03-21得票数 1
回答已采纳
我尝试用gfortran编译基于matlab "fintrf.h“的fortran代码，以便在fortran程序中导入MAT文件。下面是编译命令： gfortran -mcmodel=medium -Ofast $1.F90 -I/mnt/f/MATLAB/extern/include/ -L/mnt/f/MATLAB/bin/win64/ -cpp -lmat -lmx -Wl,-rpath /mnt/f/MATLAB/bin/win64/ -o $1.exe 我得到了这些错误： /usr/bin/ld: cannot find -lmat
/usr/bin/ld: cannot
浏览 50提问于2020-06-17得票数 0
1回答
多点测井方法
、、、、
我正在构建一个带有Java前端和性能非常关键的Fortran库的桌面应用程序。我们最近开始将更多的代码从Java中推到Fortran中(参见下面)
不久，我们可能会在SCALA中有模块(使用akka.io或JPPF的分布式计算系统)。
我想知道我能做些什么来确保正常的、合并的、格式正确的日志。现在，我们相当广泛的java前端正在使用java.util.logging (7月)。正如我确信你们中的许多人也已经处理过，引导我们的记录器配置已经是一个比我想要的更复杂的过程，并且没有自动化的测试。
我们的Fortran二进制文件，根据Intel Fortran库的默认行为，偶尔会低声说出标准的内容，这
浏览 0提问于2015-07-24得票数 3
我想知道如何使用RcppArmadillo从QR分解返回“紧凑”表单。使用RcppArmadillo在矩阵上运行QR分解相对简单，如下所示：
library(inline)
src <- '
using namespace Rcpp;
using namespace arma;
mat X = as<mat>(A);
mat Q, R;
qr(Q,R,X);
return Rcpp::List::create(Rcpp::Named("Q")=Q,Rcpp::Named("R")=R);
'
fun <- cxxf
浏览 4提问于2015-03-02得票数 2
我试图用flex & bison自动生成一个fortran和解析器，但是在扫描以下fortran语言时遇到了一个错误："if(i.le.20.and.j.le.10)“
我发现的原因是“.and”是fortran语法中的逻辑运算符，而浮点数可以写为“20”。这样我的雷克萨斯就能认出“20”作为一个浮点数，基于“匹配最长的可能字符串”的规则，而左字符串"and.j.le.10“没有匹配的规则。那我该怎么解决这个问题呢？
浏览 0提问于2019-01-08得票数 1
1回答
我不知道我是碰到了编译器的错误还是遗漏了什么。我试图运行以下代码：
program test
implicit none
integer, parameter :: dp=kind(1.d0)
integer, allocatable :: int_mat(:,:)
integer, allocatable :: int_mat_2(:,:)
real(dp), allocatable :: real_mat(:,:)
allocate(int_mat(2,2))
int_mat = 0
int_mat_2 = int_mat
real_mat =
浏览 5提问于2017-09-04得票数 4
回答已采纳
我试图在我的角度测试用例中模拟viewchild。
这是我的html文件
<mat-form-field appearance="fill">
<mat-label>Date</mat-label>
<input matInput #dateInput [matDatepicker]="picker"
(dateChange)="dateChange($event)" readonly>
<mat-datepicker-toggle matSuffix [for]="pick
浏览 5提问于2022-04-06得票数 0
在作业作业中，我的任务是将C#中的一个循环转换为Fortran 95。
outerLoop:
for(row = 0; row < numRows; rows++){
for(col = 0; col < numCols; col++){
if(mat[row][col] == 0)
continue outerLoop;
sum += mat[row][col];
}
}
正如一些人所看到的，这看起来类似于Java和Perl中使用的exit语句标签规范，据我所知，这些规范用于打破具有嵌套循环或“if”语句的循环，而不
浏览 5提问于2012-04-21得票数 2
回答已采纳
有人能帮助我阅读从Visual生成的.mat文件吗？您可以从链接下载.mat文件：
您可以在链接中的zip中获得一个.mat文件，而该文件正是我所请求的帮助。好像是个ASCII文件。我不知道如何读取文件中的数据。我试图用scipy.io.loadmat()加载.mat文件中的数据，但是出现了一个错误，如下所示：
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
raise ValueError('Unknown mat file type, versio
浏览 7提问于2020-07-09得票数 0
回答已采纳
我有一个Fortran程序，我使用f2c将其转换为C，因为我需要在iOS上运行它，所以这将是我唯一的选择。该程序依赖于文件的读取和写入，特别是它生成一个未格式化的二进制文件，它随后读取这些文件以生成人类可读的结果文件。
当我在Windows上运行C代码时，一切都按预期工作，这意味着C程序生成一个未格式化的二进制文件，然后使用该文件生成一些结果。但是，当我将相同的C代码作为iOS应用程序在iPad上运行时，会生成未格式化的二进制文件，但文件大小略有不同，而且当对二进制文件使用各种比较工具时，文件内容显然是不可读的。问题是，iOS生成的文件导致C代码与：Thread 1:EXC_BAD_ACCES
浏览 14提问于2022-08-19得票数 -1
回答已采纳
2回答
我有一个21x21矩阵，在大多数情况下，它包含一个矩阵。某些元素需要用离散数字初始化。现有的解决方案是两步初始化:首先将所有元素设置为1，然后将选定的元素重新分配到离散数字。
这是在遗留的FORTRAN77代码中实现的，请参见下面的代码，它使用gfortran-4.9 (即使没有警告)进行很好的编译。
但是，我需要在运行gfortran-4.4的旧服务器上编译它。这与f2c或g95编译器等替代工具一起失败，原因是“重叠初始化”。
现在，我如何重新编码这个问题，它避免了这种重叠在特殊的块数据程序单元？我也认为这是不推荐的，虽然可能的英特尔Fortran编译器，因为初始化的顺序不能保证。
浏览 5提问于2015-12-02得票数 1
回答已采纳
我已经工作了相当一段时间的科学(学术)软件的机械分析。现在，代码从八十年代开始(Fortran 77)，以混合/混合形式的Fortran 90/95到达我这里。然而，由于需要添加像MKL这样的强大工具，我决定从旧的11.072 (使用VS2008)转移到“最近”14.0 (在ComposerXE中)。F77核心没有问题地编译，但是我在处理变量延迟定义的子程序中遇到了一些接口问题。而不是烦人的庞大的例程，我做了一个MWE能够重复问题。
迷你程序复制如下，因此您可以修补它：
program main
implicit none
print *, 'Start of the p
浏览 2提问于2016-01-17得票数 1
回答已采纳
我刚开始尝试寻找NxN矩阵的行列式。我使用了来自dualm.wordpress.com/2012/01/06/computing-determinant-in-fortran/的代码
program lapackdet
implicit none
integer :: N
complex, allocatable, dimension(:,:) :: mat
complex :: det
integer :: i, info
integer, allocatable :: ipiv(:)
real :: sgn
N=2
allocate(ipiv(N)
浏览 2提问于2013-08-16得票数 2
2回答
我正在尝试优化一些python代码(以加速一些矩阵运算)，我的代码类似于这个(我的真实数据集也类似于'gps')，
import numpy as np
gps = [np.random.rand(50,50) for i in xrange(1000)]
ips = np.zeros( (len(gps),len(gps)), dtype='float32')
for i in xrange(len(gps)):
for j in xrange(0,i+1):
ips[i,j]= f.innerProd(gps[i],gps[j])
ip
浏览 1提问于2011-05-01得票数 3
回答已采纳
1回答
我有以下Fortran代码：
module class_MAT
implicit none
private
public :: load_coo
type Mat_CSR
real, dimension(:), allocatable :: val, colInd, rowPtr
end type Mat_CSR
contains
subroutine load_coo(filename, this)
implicit none
type(mat_csr) :: this
character(len=50) , intent(i
浏览 1提问于2016-04-22得票数 0
回答已采纳
2回答
我正在编写一个非常基本的Fortran代码来创建Ising模型。然而，我被困在最后一步-重复计算，直到达到最稳定的条件。
do
!Calculation (omitted here)
!Making a decision
if (E1 /= -64) then !not yet stable
if(dE > 0.0) then
call seed
call random_number(prob) ! random generate a number which is 0 <= Prob <= 1
浏览 0提问于2015-12-07得票数 4
3回答
我正在尝试用Dymola的.mat结果进行数据处理。我的计划是使用MATLAB。我有几个关于.mat文件的问题：
如果我直接将.mat文件加载到MATLAB中，数据结构非常奇怪，我必须使用Dymola附带的MATLAB脚本来加载.mat文件。有关于数据如何存储在.mat文件中的解释吗？
当用结果绘制图表时，我想改变这个单元，但我不知道如何使Dymola输出数据与我想要使用的单元。当Dymola将数据输出到.mat文件中时，是否有允许我更改单元的设置？
浏览 8提问于2020-08-01得票数 2
回答已采纳
我有一个问题，如何将非常大的矩阵导入fortran 6.5。例如，我的.for文件矩阵非常大，为120MB，矩阵的形式为
A(1,1)=691627.009959795*x(1)**2 - 3428.41117804251*x(1)*x(2) -
# 156.91392014999*x(2)**2 + 428391.21915942*x(1)*x(3) - 953.4687530
#88931*x(2)*x(3) + 66317.7918091961*x(3)**2 + 159956.172489745*x(1)
#*x(4) - 427.609460294073*x(2)*x(4)
浏览 0提问于2012-05-08得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
