求助MKL库求解三对角矩阵 - Fortran Coder

[Fortran] 纯文本查看 复制代码



PROGRAM MAIN

    IMPLICIT NONE

    INTEGER num

    PARAMETER (num = 5)

    INTEGER INFO

    REAL(8) A(1:3*num-2)

    REAL(8) B(1:num,1)    

    A(1) = 2.0000; A(2) = 1.0000;

    A(3) = 0.3571; A(4) = 2.0000; A(5) = 0.6429;

                   A(6) = 0.6000; A(7) = 2.0000; A( 8) = 0.4000;

                                  A(9) = 0.4286; A(10) = 2.0000; A(11) = 0.5714

                                                 A(12) = 1.0000; A(13) = 2.0000

    B(1,1) = 5.5200; B(2,1) = 4.3144; B(3,1) = 3.2661; B(4,1) = 2.4287; B(5,1) = 2.1150

    INFO = -1

    CALL DDTTRSB('T', num, 1, A(3:13:3), A(1:13:3), A(2:13:3), B, num, INFO)

    WRITE(*, *) INFO

END PROGRAM MAIN