xginvgjd.h

来自「基本矩阵运算 : + - *, power, transpose, trace,」· C头文件 代码 · 共 21 行

/* xginvgjd.h                    freeware                  [e-mail]  */

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void inv_gjd_mF(
PmF mA,
PmF minv)
{
double pAID [MXR][MXC*2*C2];mF mAID ={MXR,MXC*2*C2,&pAID [0][0]};

                  mAID.rows=mA->rows;         mAID.cols=2*mA->rows*C2;

     createAid_mF(mA,&mAID);
   gaussjordan_mF(   &mAID);
       copyinv_mF(   &mAID,minv);
}