aar_random1.c

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

/*         .h         freeware           http://groups.yahoo.com/group/mathc/ */

/* --------------------------------- INCLUDES ------------------------------- */
#include "v_ahfile.h"
/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void AplsB(int n,int p)
{
mR mA    ={n,n,init_mR(n,n)};
/*-------------------------------- PROGRAM ---------------------------------- */

  clrscrn();
  rnd_symetric_mR(&mA,9);
  printf("  A : rnd_symetric_mR(&mA,9); ");
  p_mR(&mA);
  getchar();

  clrscrn();
  rnd_Ltrig_mR(&mA,9);
  printf("  A : rnd_Ltrig_mR(&mA,9); ");
  p_mR(&mA);
  rnd_Utrig_mR(&mA,9);
  printf("  A : rnd_Utrig_mR(&mA,9); ");
  p_mR(&mA);
  getchar();

  clrscrn();
  rnd0_Ltrig_mR(&mA,9);
  printf("  A : rnd0_Ltrig_mR(&mA,9); ");
  p_mR(&mA);
  rnd0_Utrig_mR(&mA,9);
  printf("  A : rnd0_Utrig_mR(&mA,9); ");
  p_mR(&mA);
  getchar();

  clrscrn();
  rnd_r0_mR(&mA,9);
  printf("  A : rnd_r0_mR(&mA,9); ");
  p_mR(&mA);
  rnd_c0_mR(&mA,9);
  printf("  A : rnd_c0_mR(&mA,9); ");
  p_mR(&mA);
  getchar();

  clrscrn();
  rnd0_diag_mR(&mA,9);
  printf("  A : rnd0_diag_mR(&mA,9); ");
  p_mR(&mA);
  rnd_diag_mR(&mA,9);
  printf("  A : rnd_diag_mR(&mA,9); ");
  p_mR(&mA);
  rndp_diag_mR(&mA,9);
  printf("  A : rndp_diag_mR(&mA,9); ");
  p_mR(&mA);
  getchar();

  free(   mA.pb);
}


void hello(void)
{
  printf("\n Freeware by : xhunga bernard 06000 Nice France.\n");
  printf("\n                                                  ");
  printf("\n            random matrices                      ");
  printf("\n            ****************                    \n");
  printf("\n The size of mA, mB, mAB : 1x1 ... 6x6.         \n");
  printf("\n     ex : \"aplsbi.exe\"                        \n");
  printf("\n You can choose the size :                      \n");
  printf("\n     ex : \"aplsbf.exe 2 3\" -> 2x3               ");
  printf("\n          \"aplsbf.exe 1 5\" -> 1x5               ");
  printf("\n          \"aplsbf.exe 4 2\" -> 4x2               ");
  printf("\n          \"aplsbf.exe 5 5\" -> 5x5               ");
  printf("\n\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;
int    p;
/*-------------------------------- PROGRAM ---------------------------------- */
  srand(time(&t));
  hello();
do
{
if( argc == 3 && atoi(argv[1])&& atoi(argv[2]))
 {n = atoi(argv[1]); if( n > 6){n=6;}
  p = atoi(argv[2]); if( p > 6){p=6;}}
else
 {    n=rndp_I(6);       p=rndp_I(6); }
   if(n<1){n=1;}     if(p<1){p=1;}


  AplsB(n,p);


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

  return 0;
}