4tcmul0.c

常用数值算法程序集

  #include "stdio.h"
  #include "4tcmul.c"
  main()
  {  int i,j;
     static double cr[3][4],ci[3][4];
     static double ar[3][4]={ {1.0,2.0,3.0,-2.0},
            {1.0,5.0,1.0,3.0},{0.0,4.0,2.0,-1.0}};
     static double ai[3][4]={ {1.0,-1.0,2.0,1.0},
            {-1.0,-1.0,2.0,0.0},{-3.0,-1.0,2.0,2.0}};
     static double br[4][4]={ {1.0,4.0,5.0,-2.0},
            {3.0,0.0,2.0,-1.0},{6.0,3.0,1.0,2.0},
            {2.0,-3.0,-2.0,1.0}};
     static double bi[4][4]={ {-1.0,-1.0,1.0,1.0},
            {2.0,1.0,0.0,5.0},{-3.0,2.0,1.0,-1.0},
            {-1.0,-2.0,1.0,-2.0}};
     tcmul(ar,ai,br,bi,3,4,4,cr,ci);
     for (i=0; i<=2; i++)
       { for (j=0; j<=3; j++)
           printf("%13.7e ",cr[i][j]);
         printf("\n");
       }
     printf("\n");
     for (i=0; i<=2; i++)
       { for (j=0; j<=3; j++)
           printf("%13.7e ",ci[i][j]);
         printf("\n");
       }
  }