/* gjid03.c                         freeware                   [[Email Removed]]
*/

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

/* --------------------------------- MAIN ----------------------------------- */
 main()
{

double pb[3][3*TWOCOL]=
{
 15,2,  13,3,   0,1,
  4,1,   0,1, -11,2,
-13,2,   1,2, -15,3,
};
matrix m={3,3*TWOCOL,&pb[0][0]};

double pbmAID [3][6*TWOCOL]; matrix mAID  ={3,6*TWOCOL,&pbmAID [0][0]};
double pbmAIDT[3][6*TWOCOL]; matrix mAIDT ={3,6*TWOCOL,&pbmAIDT[0][0]};
double pbinv  [3][3*TWOCOL]; matrix minv  ={3,3*TWOCOL,&pbinv  [0][0]};
double pbT    [3][3*TWOCOL]; matrix mT    ={3,3*TWOCOL,&pbT    [0][0]};
double pbID   [3][3*TWOCOL]; matrix mID   ={3,3*TWOCOL,&pbID   [0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
 copymF  (&m,&mT);
 copymAID(&m,&mAID);

 clrscrn();
 printf(" The matrix  mA \n");
 printmF2(&m);
 getchar();

 printf("\n");
 printf(" The matrix  mA_ID \n");
 printmF2(&mAID);
 getchar();

 clrscrn();
 printf(" Find mID (below the pivots) \n");
 idF4(&mID,&m);
 printmF2(&mID);
 getchar();

 printf("  mID * mA_ID \n");
 multmF(&mID,&mAID,&mAIDT);
 printmF2(&mAIDT);
 getchar();

 multmF(&mID,&mT,&m);

 clrscrn();
 printf(" Find mID (above the pivots) \n");
 idF4a(&mID,&m);
 printmF2(&mID);
 getchar();

 printf("  mID * mA_ID  \n");
 multmF(&mID,&mAIDT,&mAID);
 printmF2(&mAID);
  getchar();

 clrscrn();
 printmF2(&mAID);
 printf("\n");
 printf("The inverse matrix is : \n");
 copyinv(&mAID, &minv);
 printmF(&minv);
 getchar();

 copymF(&mT,&m);

 printf("\n mA * minv give \n");
 multmF(&m, &minv,&mID);
 printmF(        &mID);

 printf("\n Press return to continue");
 getchar();

 return 0;
}