[交流] 【原创】MKL程序编译与连接:Lapack95篇

ifndef static
    EXT = so
else
    EXT = a
endif

MKLLIB = $(MKLROOT)/lib/32
MKLMOD = $(MKLROOT)/include/32
STATIC = $(MKLLIB)/libmkl_lapack95.a $(MKLLIB)/libmkl_intel.$(EXT)
ifeq ($(EXT),so)
SO = $(MKLLIB)/libmkl_lapack.so
endif
FCFLAGS = -Wl,--start-group \
      $(MKLLIB)/libmkl_intel_thread.$(EXT) \
      $(SO) \
      $(MKLLIB)/libmkl_core.$(EXT)\
      -Wl,--end-group \
      -liomp5 -lpthread
SRC = test_mkl.f90
TARGET = test

all:
>---ifort -xK -w -I$(MKLMOD) $(SRC) $(STATIC) $(FCFLAGS) -o $(TARGET

program test_mkl
      USE mkl95_LAPACK, ONLY: GETRF, GETRS

      IMPLICIT NONE
      REAL(KIND = 8), DIMENSION(3, 3) :: A
      REAL(KIND = 8), DIMENSION(3) :: B
      INTEGER, DIMENSION(3) :: IPIV
      INTEGER :: INFO, I

      DATA A/1,6,9,6,5,2,8,3,4/
      DATA B/2,4,7/

      PRINT '(3F5.1)', (A(I, :), I=1,3)  ! 这个才能 输出真实的A的样子
      PRINT '(F5.1)' , B ! 输出一个列向量

      CALL GETRF(A, IPIV, INFO)
      IF(INFO.NE.0) THEN
          PRINT *, "THERE IS SOMETHING WRONG WITH THE MATRIX: ", INFO
      ENDIF
      CALL GETRS(A, IPIV, B, 'N', INFO)
      IF(INFO.NE.0) THEN
          PRINT *, "THERE IS SOMETHING WRONG: ", INFO
      ENDIF

      PRINT *, ""
      PRINT '(3F6.2)', B