(TIPsは、断片的に見つけたことの覚書です)
逆行列[A]-1を行っているクラスはここにありました。
matrix.cpp
matrix class のMatrix::Invert( )
この中では、やはり高速数理計算用のライブラリ LAPACKを使用していました。
//m行n列の行列[A]を行交換による部分枢軸選択を利用してLU分解
DGETRF(&n,&n,Aptr,&ldA,iPIV,&info);
//LU分解後の行列を利用して、逆行列[A]-1を計算
DGETRI(&n,Aptr,&ldA,iPIV,Wptr,&workSize,&info);
OpenSees のソースコード解析に挑戦してみる : 目次 |
0 件のコメント:
コメントを投稿