aar_inv_gjd.c

一些矩阵运算的C代码实现：矩阵一般运算

/*         .h         freeware           http://groups.yahoo.com/group/mathc/ */

/* --------------------------------- INCLUDES ------------------------------- */
#include "v_ahfile.h"

/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
double pbA[4][4]=
{
1, 3, 2, 3,
9, 1, 1, 2,
2, 3, 8, 5,
3, 4, 8, 1,
};
mR mA={4,4,&pbA[0][0]};

double pbInvA[4][4];mR mInvA={4,4,&pbInvA[0][0]};
double pbID  [4][4];mR mID  ={4,4,&pbID  [0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
 printf("\n\nMatrix mA follow by mInvA \n M\n");
 p_mR(&mA);

 inv_gjd_mR(&mA,&mInvA);
      p_mR   (&mInvA);

 printf("\n\n mA * mInvA \n");
 mul_mR(&mA,&mInvA,&mID);
  p_mR(           &mID);

 printf("\n Press Return to stop\n");
 getchar();

 return 0;
}