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

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

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

double pb[5][5*TWOCOL]=
{
5,1, 4,1, 1,1, 0,1, 5,1,
2,1, 1,1, 7,1, 4,1, 0,1,
3,1, 2,1, 6,1, 1,1, 1,1,
2,1, 3,1, 2,1, 2,1, 2,1,
7,1, 5,1, 0,1, 9,1, 1,1,
};
matrix m={5,5*TWOCOL,&pb[0][0]};

double pbAT [5][5*TWOCOL];  matrix mAT  = {5,5*TWOCOL,&pbAT [0][0]};
double pbT  [5][5*TWOCOL];  matrix mT   = {5,5*TWOCOL,&pbT  [0][0]};

double pbID [5][5*TWOCOL];  matrix mID  = {5,5*TWOCOL,&pbID [0][0]};
double pbIDT[5][5*TWOCOL];  matrix mIDT = {5,5*TWOCOL,&pbIDT[0][0]};
double pbID0[5][5*TWOCOL];  matrix mID0 = {5,5*TWOCOL,&pbID0[0][0]};
double pbID1[5][5*TWOCOL];  matrix mID1 = {5,5*TWOCOL,&pbID1[0][0]};
double pbID2[5][5*TWOCOL];  matrix mID2 = {5,5*TWOCOL,&pbID2[0][0]};
double pbID3[5][5*TWOCOL];  matrix mID3 = {5,5*TWOCOL,&pbID3[0][0]};
double pbID4[5][5*TWOCOL];  matrix mID4 = {5,5*TWOCOL,&pbID4[0][0]};

/*-------------------------------- PROGRAM ---------------------------------- */
 printf(" The matrix mA \n");
 printmF(&m);
 copymF(&m,&mAT);

idF3(&mID,&m,0);copymF(&mID,&mID0);printf("mID0");printmF(&mID0);getchar();
idF3(&mID,&m,1);copymF(&mID,&mID1);printf("mID1");printmF(&mID1);getchar();
idF3(&mID,&m,2);copymF(&mID,&mID2);printf("mID2");printmF(&mID2);getchar();
idF3(&mID,&m,3);copymF(&mID,&mID3);printf("mID3");printmF(&mID3);getchar();
idF3(&mID,&m,4);copymF(&mID,&mID4);printf("mID4");printmF(&mID4);getchar();

 printf("\n mID4 * mID3 : \n");
 multmF(&mID4, &mID3, &mID);
 printmF(&mID );
 getchar();

 printf("\n (mID4 * mID3) * mID2 : \n");
 multmF(&mID,  &mID2, &mIDT);
 printmF(&mIDT );
 getchar();

 printf("\n mID = ((mID4 * mID3) * mID2) * mID1 : \n");
 multmF(&mIDT, &mID1, &mID);
 printmF(&mID );
 getchar();

 printf("\n mID = (((mID4 * mID3) * mID2) * mID1) * mDI0 : \n");
 multmF(&mID, &mID0, &mIDT);
 printmF(&mIDT );
 getchar();

 copymF(&mAT,&m);
 printf("\n mID * mA : \n");
 multmF(&mIDT,&m,&mT);
 printmF(        &mT);

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

 return 0;
}