mulmfr.c

[[ Complex Matrices : Language c]]

/* .c                         freeware                 [e-mail]
*/

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

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void AB(int n, int p, int q)
{
double pbA [MXR][MXC*C4];mzF mA ={MXR,MXC*C4,&pbA [0][0]};
double pbB [MXR][MXC*C4];mzF mB ={MXR,MXC*C4,&pbB [0][0]};
double pbAB[MXR][MXC*C4];mzF mAB={MXR,MXC*C4,&pbAB[0][0]};

                     mA.rows=n;                mA.cols=p*C4;
                     mB.rows=p;                mB.cols=q*C4;
                    mAB.rows=n;               mAB.cols=q*C4;
/*-------------------------------- PROGRAM ---------------------------------- */
  rnd0_mzF(&mA,3,3);
  rnd0_mzF(&mB,3,3);

  clrscrn();
  printf("                                                       A :");
  p_mizF(&mA);
  printf("                                                       B :");
  p_mizF(&mB);
  getchar();

  mul_mzF(&mA,&mB,&mAB);

  printf("                                                       A * B :");
  p_mizF(&mAB);
}


void hello(void)
{
  printf("\n Freeware by : xhunga bernard 06000 Nice France.\n");
  printf("\n                                                  ");
 printf("\n            Mutiply two complex matrices           ");
  printf("\n            ***************************         \n");
  printf("\n The size of mA, mB, mAB : 1x1 ... 6x6.         \n");
  printf("\n     ex : \"mulmFR.exe\"                        \n");
  printf("\n You can choose the size :                      \n");
  printf("\n     ex : \"mulmFR.exe 2 3 5\" -> 2x3 3x4 = 2x4   ");
  printf("\n          \"mulmFR.exe 1 5 2\" -> 1x4 4x2 = 1x2   ");
  printf("\n          \"mulmFR.exe 4 2 1\" -> 4x2 2x1 = 4x1   ");
  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;
int    q;

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


  AB(n,p,q);


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

  return 0;
}