在Fedora Linux上使用gfortran与LAPACK链接

Question

我需要有关在Fedora机器上安装LAPACK（用于Fortran 95）的一些帮助。

我在尝试过在此论坛中针对相关查询发布的所有答案后，正在撰写此帖子。

我尝试了两种不同的选择：

方法1

我直接从此处下载了lapack95.a库（ http://www.netlib.org/lapack95/ ，文件名为lapack95_linux_redhat.tgz。解压缩后，我可以获得文件lapack95.a。

在找到有关如何链接和使用库的答案之后，我编写了一个试用代码并将其保存在文件try.f90中，然后执行以下操作

gfortran -c try.f90

gfortran -o try try.o -L/(name of directory where lapack95.a
                          was present. I put it in the same directory
                          as try.f90)/ -llapack95

（ fortran文件try.f90和lapack95.a位于同一目录中 ）

我收到以下错误：

/usr/bin/ld: cannot find -llapack95
collect2: error: ld returned 1 exit status

我在互联网上进行搜索，发现此处（ https://forums.fedoraforum.org/archive/index.php/t-248227.html ）需要安装lapack-devel。 我做到了，结果也一样。

方法二

我也尝试运行上述命令（在方法1中），但在/和-llapack95之间没有空格，即：

gfortran -o try try.o -L/(name of directory where lapack95.a 
                          was present I put it in the same directory
                         as try.f90)/-llapack95

并收到以下错误：

try.o: In function `MAIN__':
try.f90:(.text+0xdb): undefined reference to `sgesv_'
collect2: error: ld returned 1 exit status

方法3

在没有成功之后，我认为直接使用lapack95.a并不是一个好主意。 我应该从原始文件中构建它。

我意识到可能是我需要通过正确的方法来安装lapack95，然后从同一链接http://www.netlib.org/lapack95

我下载了文件lapack95.tgz并按照自述文件中的说明进行操作。（该文件只是为了在SRC中运行make文件）

我做了一些更改，如：

FC1      = f95 -fixed  to 

FC1      = f95 

和

OPTS0    = -u -V -dcfuns -dusty -ieee=full 

至

OPTS0    = -u -V 

（如https://www.reddit.com/r/linux4noobs/comments/7tgmb9/fortran_library_netlib_lapack95_installation_error/上的建议）

我期望该库将以LAPACK3的形式在/ usr / lib /或usr / local / lib中创建，但找不到。

实际上有一个错误（因此不能有任何.a文件）

f95 -free -c -u -V--I./../lapack95_modules f77_lapack_double.f90 f95：错误：当输入来自标准输入时需要-E或-x make：*** [../make.inc： 45：f77_lapack_double.o]错误1

我不确定我的方法是否还可以。 我希望方法1可以正常工作，因为我正在本地链接库，但是在所有3种情况下都失败了。

请提示。 谢谢

Answer 1

Lapack库是用于快速执行矩阵运算的非常有用的工具。 在linux机器上，使用内置工具dnf / yum或apt-get [2]安装lapack非常简单。 在我的Fedora桌面上，安装如下：

sudo dnf install lapack64-devel.x86_64

这是一个有关如何使用lapack lib的小示例。 源代码是将线性方程求解为linear_eqs.f [3]。

    Program LinearEquations
c solving the matrix equation A*x=b using LAPACK
    Implicit none
c declarations, notice single precision
    Real*4 A(3,3), b(3)
    integer i, j, pivot(3), ok
c define matrix A
    A(1,1)=3.1
    A(1,2)=1.3
    A(1,3)=-5.7
    A(2,1)=1.0
    A(2,2)=-6.9
    A(2,3)=5.8
    A(3,1)=3.4
    A(3,2)=7.2
    A(3,3)=-8.8
c define vector b, make b a matrix and you can solve multiple
c equations with the same A but different b
    b(1)=-1.3
    b(2)=-0.1
    b(3)=1.8
c find the solution using the LAPACK routine SGESV
    call SGESV(3, 1, A, 3, pivot, b, 3, ok)
c
c parameters in the order as they appear in the function call
c    order of matrix A, number of right hand sides (b), matrix A,
c    leading dimension of A, array that records pivoting, 
c    result vector b on entry, x on exit, leading dimension of b
c    return value   
c   
c print the vector x
    do i=1, 3
       write(*,*) b(i)
    end do
    end

在此代码中，使用了一个lapack函数'sgesv（）'。 要执行此文件，请运行以下命令：

gfortran linear_eqs.f -llapack -o lineqs.out
./lineqs.out

输出将是：

1.00000024    
1.00000036    
1.00000036 

参考

[1] https://github.com/fxcoudert/gfortran-for-macOS/releases

[2] BLAS，ATLAS，LAPACK共享库的最小示例

[3] http://physics.oregonstate.edu/~landaur/nacphy/lapack/codes/linear-f.html