aar_maxima_adjt.c

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

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

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

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void adjoint(int n)
{
int      row;
int      col;
double f;

double pbA    [MXR][MXC];mR mA    ={MXR,MXC,&pbA    [0][0]};
double pbCofac[MXR][MXC];mR mCofac={MXR,MXC,&pbCofac[0][0]};
double pbAdjnt[MXR][MXC];mR mAdjnt={MXR,MXC,&pbAdjnt[0][0]};

                     mA.rows=n;                   mA.cols= n;
                 mCofac.rows=n;               mCofac.cols= n;
                 mAdjnt.rows=n;               mAdjnt.cols= n;
/*-------------------------------- PROGRAM ---------------------------------- */
  rnd_mR(&mA,9);

  clrscrn();
  printf(" Copy/Past into the maxima windows \n\n\n\n");
  printf(" kill(all); \n\n");
  p_mR_maxima(&mA,"a");
  printf(" b : adjoint(a); \n\n");
  getchar();

  clrscrn();
  printf("                                                     A :\n");
  p_mR(&mA);
  printf("                                                mAdjnt :\n");
  cofactor_mR(&mA,&mCofac);
   adjoint_mR(    &mCofac,&mAdjnt);
        p_mR(            &mAdjnt);
}

void hello(void)
{
  printf("\n Freeware by : xhunga bernard 06000 Nice France.\n");
  printf("\n                                                  ");
  printf("\n            Find the adjoint matrix of mA        ");
  printf("\n            ******************************      \n");
  printf("\n The size of mA : 2x2 ... 4x4.                  \n");
  printf("\n           ex : \"adjointf.exe\"                \n");
  printf("\n You can choose the size :                      \n");
  printf("\n           ex : \"adjointf.exe 2\" -> mA 2x2      ");
  printf("\n                \"adjointf.exe 3\" -> mA 3x3      ");
  printf("\n                \"adjointf.exe 4\" -> mA 4x4      ");
  printf("\n                                                  ");
  printf("\n Distributed  by :  http://www.simtel.net/        ");
  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;}


  adjoint(n);


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

  return 0;
}