[問題] 如何用GNU scientific libray(GSL)求複數反矩陣
我常常要算複數矩陣的反矩陣
就是一個矩陣內每一個單元都是呈現a+bj的形式
我目前掛載了GNU scientific libray(GSL)
我用他求實數矩陣的反矩陣答案正確
但將指令改為複數指令結果是錯的
p = gsl_permutation_alloc (4);
int sin=1;
gsl_linalg_LU_decomp(&AA.matrix , p,&sin );
gsl_linalg_LU_invert (&AA.matrix , p, inverseAA ) ;
以上是實數的寫法
p = gsl_permutation_alloc (4);
int sin=1;
gsl_linalg_LU_complex_decomp(&AA.matrix , p,&sin );
gsl_linalg_LU_complex_invert (&AA.matrix , p, inverseAA ) ;
以上是我直接改為複數指令的方法但結果是錯誤的
我要解的是4x4的方陣
希望有高手能幫我解答 如何修改才能求出複數反矩陣ꨊ
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 211.74.215.161
C_and_CPP 近期熱門文章
PTT數位生活區 即時熱門文章
