aaz_inv_maple.c

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

/* gjid03.c                      freeware  http://groups.yahoo.com/group/mathc/
*/

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


/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void DeterA(int n)
{
complexZ z_det;

mZ mA   ={n,n*C2,init_mZ(n,n*C2)};
mZ minvA={n,n*C2,init_mZ(n,n*C2)};
mZ mID  ={n,n*C2,init_mZ(n,n*C2)};
/*-------------------------------- PROGRAM ---------------------------------- */
  do
  {
   if(n==4) rnd0_mZ(&mA,2);
   if(n==3) rnd0_mZ(&mA,3);
   if(n==2) rnd0_mZ(&mA,4);
   if(n==1) rnd0_mZ(&mA,5);

   z_det = det_Z(&mA);

  }while( (!z_det.r) && (!z_det.i));

  clrscrn();
  printf(" Copy/Past into the maple windows \n\n");
  printf(" restart:\n");
  printf(" with(linalg):\n\n");
  p_mZ_maple(&mA,"a");
  printf("\n");
  printf(" inverse(a);\n");
  getchar();

  clrscrn();
  printf("     A\n");
  p_mZ(&mA);

  printf(" inv(A)\n");
  inverse_mZ(&mA,&minvA);
  p_mZ(&minvA);

  mul_mZ(&mA,&minvA,&mID);
  p_mZ(&mID);

 free(   mA.pb);
 free(minvA.pb);
 free(  mID.pb);
}

void hello(void)
{
  printf("\n Freeware by : xhunga bernard 06000 Nice France.\n");
  printf("\n                                                  ");
  printf("\n                 inverse(mA)                      ");
  printf("\n                 ***********                    \n");
  printf("\n The size of mA : 2x2 ... 5x5.                  \n");
  printf("\n           ex : \"inverse.exe\"                 \n");
  printf("\n You can choose the size :                      \n");
  printf("\n           ex : \"inverse.exe 2\" -> mA 2x2       ");
  printf("\n                \"inverse.exe 3\" -> mA 3x3       ");
  printf("\n\n                                                ");
  
  printf("\n                                                  ");
  printf("\n Press return to continue                          ");
  printf("\n                                                  ");

  getchar();
}

/* --------------------------------- MAIN ----------------------------------- */
int main(int argc, char *argv[])
{
time_t t;
int    n;

/*-------------------------------- PROGRAM ---------------------------------- */
  srand(time(&t));
  hello();
do
{
if( argc == 2 && atoi(argv[1]))
 {  n = atoi(argv[1]); if( n > 3){n=3;}}
else
 {  n = rndp_I(3);}
if( n < 1){n=1;}


  DeterA(n);


printf("\n\n Press 'X' Return to stop\n");
} while(toupper(getchar())!='X');

  return 0;
}