aaz_lu_r.c

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

/*         .h         freeware           http://groups.yahoo.com/group/mathc/ */
/* --------------------------------- INCLUDES ------------------------------- */
#include "y_ahfile.h"
/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void DeterA(int n)
{
complexZ z_det;

double pbA [n][n*C2];mZ mA ={n,n*C2,&pbA [0][0]};
double pbU [n][n*C2];mZ mU ={n,n*C2,&pbU [0][0]};
double pbL [n][n*C2];mZ mL ={n,n*C2,&pbL [0][0]};
double pbLU[n][n*C2];mZ mLU={n,n*C2,&pbLU[0][0]};

                     mA.rows=n;                   mA.cols=n*C2;
                     mU.rows=n;                   mU.cols=n*C2;
                     mL.rows=n;                   mL.cols=n*C2;
                    mLU.rows=n;                  mLU.cols=n*C2;
/*-------------------------------- PROGRAM ---------------------------------- */
  do
  {
   if(n==4) rnd_mZ(&mA,2);
   if(n==3) rnd_mZ(&mA,2);
   if(n==2) rnd_mZ(&mA,4);

   z_det = det_Z(&mA);

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

  clrscrn();

  printf("  A :");
  p_mZ(&mA);

  LU_mZ( &mA,&mL, &mU);
  printf("  L :");
  p_mZ(&mL);

  printf("  U :");
  p_mZ(&mU);

  mul_mZ(&mL,&mU,&mLU);
  printf(" LU :");
  p_mZ(&mLU);
}

void hello(void)
{
  printf("\n Freeware by : xhunga bernard 06000 Nice France.\n");
  printf("\n                                                  ");
  printf("\n            LU decomposition.                     ");
  printf("\n            ****************                    \n");
  printf("\n The size of mA : 2x2 ... 4x4.                  \n");
  printf("\n           ex : \"LUFR.exe\"                    \n");
  printf("\n You can choose the size :                      \n");
  printf("\n           ex : \"LUFR.exe 2\" -> mA 2x2          ");
  printf("\n                \"LUFR.exe 3\" -> mA 3x3          ");
  printf("\n                \"LUFR.exe 4\" -> mA 4x4          ");
  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 > 4){n=4;}}
else
 {  n = rndp_I(4);}
if( n < 2){n=2;}


  DeterA(n);


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

  return 0;
}