innerprod6.c

来自「基本矩阵运算 : + - *, power, transpose, trace,」· C语言 代码 · 共 48 行

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

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

/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
fraction InnerProduct;

double pbu[6][1*C2]=
{
  1,5,
  3,2,
  5,3,
  1,2,
  3,1,
  1,3,
};
mF mu={6,1*C2,&pbu[0][0]};

double pbv[1][6*C2]=
{
  5,9,  2,8,  8,7, 1,2,  3,4,  5,1,
};
mF mv={1,6*C2,&pbv[0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
  InnerProduct = innerprod_F(&mu,&mv);

  printf("\n\n                                    u : \n");
  p_mF3(&mu);

  printf("\n\n                                    u : \n");
  p_mF3(&mv);

  printf("\n\n Inner Product =  %6.0f/%-6.0f ",
            InnerProduct.n,
            InnerProduct.d);

 printf("\n\n Press Return to continue");
 getchar();

 return 0;
}