如何在Visual Studio Code中为gfortran编译器安装LAPACK？

一、确认基础环境 #

• 确保gfortran已安装：打开终端执行gfortran --version，能显示版本号即为正常
• VS Code中已安装Fortran相关插件（如Fortran IntelliSense、Fortran Breakpoint Support）

二、按操作系统安装LAPACK #

1. Linux（Ubuntu/Debian为例） #

• 终端执行命令安装LAPACK及依赖的BLAS：

  sudo apt update
  sudo apt install liblapack-dev libblas-dev
  

• 验证安装：执行find /usr/lib -name liblapack.so，能找到对应文件即为成功

2. macOS #

• 通过Homebrew安装：

  brew install lapack
  

• 验证安装：Intel芯片执行find /usr/local/lib -name liblapack.dylib，M系列芯片执行find /opt/homebrew/lib -name liblapack.dylib，找到文件即为成功

3. Windows（MinGW-w64为例） #

• 方法一：MinGW包管理器安装
  打开MinGW Installation Manager，找到lapack和liblapack-dev类包，标记为安装后点击Installation > Apply Changes完成安装
• 方法二：手动编译安装
  1. 解压LAPACK源码包到本地
  2. 打开MinGW终端进入源码目录，执行：

     cmake -G "MinGW Makefiles" .
     mingw32-make
     mingw32-make install
     

  3. 将LAPACK库路径（如C:\MinGW\lib）添加到系统环境变量PATH

三、VS Code中配置编译任务 #

1. 打开Fortran项目文件夹，按Ctrl+Shift+B（Windows/Linux）或Cmd+Shift+B（macOS），选择创建tasks.json文件并指定gfortran
2. 修改tasks.json，添加LAPACK链接参数：
   • Linux/macOS版：

     {
       "version": "2.0.0",
       "tasks": [
         {
           "label": "gfortran build",
           "type": "shell",
           "command": "gfortran -o ${fileBasenameNoExtension} ${file} -llapack -lblas",
           "group": {
             "kind": "build",
             "isDefault": true
           },
           "problemMatcher": ["$gcc"]
         }
       ]
     }
     

   • Windows版：

     {
       "version": "2.0.0",
       "tasks": [
         {
           "label": "gfortran build",
           "type": "shell",
           "command": "gfortran -o ${fileBasenameNoExtension}.exe ${file} -llapack -lblas",
           "group": {
             "kind": "build",
             "isDefault": true
           },
           "problemMatcher": ["$gcc"]
         }
       ]
     }
     

3. 保存文件后，按Ctrl+Shift+B即可编译调用LAPACK的代码

四、测试配置有效性 #

编写测试代码验证：

program test_lapack
  implicit none
  real :: A(2,2), A_inv(2,2)
  integer :: ipiv(2), info

  A = reshape([1.0, 2.0, 3.0, 4.0], [2,2])
  
  call dgetrf(2, 2, A, 2, ipiv, info)
  call dgetri(2, A, 2, ipiv, A_inv, 4, info)

  if (info == 0) then
    print *, "矩阵A的逆矩阵："
    print *, A_inv(1,:)
    print *, A_inv(2,:)
  else
    print *, "求逆失败，错误码：", info
  end if

end program test_lapack

编译运行后，若能输出正确逆矩阵，说明配置完成

内容的提问来源于stack exchange，提问作者Dilpreet sandhu