aaz_sub.c

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

/*         .h         freeware           http://groups.yahoo.com/group/mathc/ */
/* --------------------------------- INCLUDES ------------------------------- */
#include "y_ahfile.h"
/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
double pbA[3][3*C2]=
{
  1, 2,   3, 4,   5, 6,
  1, 1,   3, 4,   5, 6,
  1, 1,   1, 1,   1, 6,
};mZ mA={3,3*C2,&pbA[0][0]};

double pbB[3][3*C2]=
{
  9, 6,   -9, -6,  9, -6,
  8, 5,   -8, -5,  9, -6,
  7, 4,   -7, -4,  9, -6,
};mZ mB={3,3*C2,&pbB[0][0]};

double pbAmnsB[3][3*C2];mZ mAmnsB={3,3*C2,&pbAmnsB[0][0]};
/*-------------------------------- PROGRAM ---------------------------------- */
  clrscrn();

  printf("\n A");
  p_mZ(&mA);
  printf("\n B");
  p_mZ(&mB);

  sub_mZ(&mA,&mB,&mAmnsB);
  printf("\n A-B");
  p_mZ(&mAmnsB);

 printf("\nPress return to continue.\n");
 getchar();

 return 0;
}