Fortran Coder

查看: 151|回复: 0

[数学库] 求助MKL库求解三对角矩阵

[复制链接]

1

帖子

1

主题

0

精华

新人

F 币
9 元
贡献
3 点
发表于 2018-5-9 15:23:06 | 显示全部楼层 |阅读模式
我在调用库函数DDTTRSB求解三对角矩阵的时候(满足对角占优),发现求解结果不对,可能是我对该函数的参数理解错误(尤其DL,D,DU的形式),求帮忙解答下,谢谢。
代码如下:
[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


回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2018-9-26 05:02

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表