aar_qr_s.c

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

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

/* --------------------------------- INCLUDES ------------------------------- */
#include "v_ahfile.h"
/* --------------------------------- MAIN ----------------------------------- */
int main(void)
{
double pba[5][5]=
{
     3,          7,          3,          3,          3,
     1,          1,          3,          3,          1,
    -2,          1,         -2,         -2,         -2,
    -2,          3,          3,          1,          1,
     3,          3,         -2,          1,          3,
};
mR A={5,5,&pba[0][0]};

double pbQ [MXR][MXC];mR  Q={MXR,MXC,&pbQ [0][0]};
double pbR [MXR][MXC];mR  R={MXR,MXC,&pbR [0][0]};
double pbQR[MXR][MXC];mR QR={MXR,MXC,&pbQR[0][0]};

     Q.rows = 5;            Q.cols = 5;
     R.rows = 5;            R.cols = 5;
    QR.rows = 5;           QR.cols = 5;
/*-------------------------------- PROGRAM ---------------------------------- */
   clrscrn();
   QR_mR(&A,&Q,&R);
   printf("                                          A : ");
   p_mR(&A);
   printf("                                          Q : ");
   p_mR(&Q);
   printf("                                          R : ");
   p_mR(&R);
   getchar();

   clrscrn();
   printf("                                          A : ");
   p_mR(&A);
   printf("                                     A = QR : ");
   mul_mR(&Q,&R,&QR);
   p_mR(&QR);

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

 return 0;
}