24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1563  |  回复: 1

bmw1818

新虫 (初入文坛)

[交流] 咨询一个调用IMSL的问题

各位大侠好!小弟刚刚开始编程,在利用IMSL求矩阵逆时出现一个不理解的问题。
运行IMSL帮助里LINCG求逆程序是可以运行的,但是放到自己程序里面却报错,找不到这个程序。
error #6285: There is no matching specific subroutine for this generic subroutine call.   [CSSCAL]
想问下哪里出问题了?


附例子:
       include 'link_fnl_shared.h'
       USE LINCG_INT
       USE WRCRN_INT
       USE CSSCAL_INT
! Declare variables
       PARAMETER (LDA=3, LDAINV=3, N=3)
       REAL THIRD
       COMPLEX A(LDA,LDA), AINV(LDAINV,LDAINV),ccc(LDA,LDA)
! Set values for A
!
! A = ( 1.0+1.0i 2.0+3.0i 3.0+3.0i)
! ( 2.0+1.0i 5.0+3.0i 7.0+4.0i)
! ( -2.0+1.0i -4.0+4.0i -5.0+3.0i)
!
        DATA A/(1.0,1.0), (2.0,1.0), (-2.0,1.0), (2.0,3.0), (5.0,3.0),
     & (-4.0,4.0), (3.0,3.0), (7.0,4.0), (-5.0,3.0)/
!
! Scale A by dividing by three
        THIRD = 1.0/3.0
        DO 10 I=1, N
        CALL CSSCAL (N, THIRD, A(:,I), 1)
10     CONTINUE
! Calculate the inverse of A
        CALL LINCG (A, AINV)
! Print results
       ccc= matmul(rn,ra3)
       CALL WRCRN ('AINV', AINV)
        read(*,*)
        END
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bmw1818

新虫 (初入文坛)

自己解决了,是csscal函数的精度问题,我的数据是双精度的,函数默认是单精度。
2楼2017-01-23 20:28:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bmw1818 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见