92 lines
2.1 KiB
C
92 lines
2.1 KiB
C
/*
|
||
*-----------------------------------------------------------------------------
|
||
* file: demo.c
|
||
* desc: demostrate how to use Patrick's matrix toolbox
|
||
* by: ko shu pui, patrick
|
||
* date: 24 nov 91 v0.1
|
||
* revi: 14 may 92 v0.2
|
||
* 24 may 92 v0.4
|
||
* ref:
|
||
* [1] Mary L.Boas, "Mathematical Methods in the Physical Sciene,"
|
||
* John Wiley & Sons, 2nd Ed., 1983. Chap 3.
|
||
*
|
||
* [2] Kendall E.Atkinson, "An Introduction to Numberical Analysis,"
|
||
* John Wiley & Sons, 1978.
|
||
*
|
||
* [3] Alfred V.Aho, John E.Hopcroft, Jeffrey D.Ullman, "The Design
|
||
* and Analysis of Computer Algorithms," 1974.
|
||
*
|
||
*-----------------------------------------------------------------------------
|
||
*/
|
||
#include <stdio.h>
|
||
#include <time.h>
|
||
#include "matrix.h"
|
||
|
||
int main()
|
||
{
|
||
MATRIX A, B, X, M;
|
||
FILE *fp;
|
||
double result;
|
||
time_t t1, t2;
|
||
int tinv, tdet, tmul;
|
||
|
||
A = mat_creat( 10, 10, UNDEFINED );
|
||
B = mat_creat( 10, 1, UNDEFINED );
|
||
|
||
if ((fp = fopen( "demo.dat", "r" )) == NULL)
|
||
{
|
||
fprintf( stderr, "file cannot be opened\n" );
|
||
exit (0);
|
||
}
|
||
|
||
fgetmat( A, fp );
|
||
printf( "|- Matrix A -|\n");
|
||
mat_dumpf( A, "%+06.1f " );
|
||
|
||
t1 = time(&t1);
|
||
result = mat_det(A);
|
||
t2 = time(&t2);
|
||
tdet = t2 - t1;
|
||
printf( "\n\nDet(A) = %f\n", result );
|
||
|
||
printf( "|- Inv A -|\n");
|
||
t1 = time(&t1);
|
||
X = mat_inv( A );
|
||
t2 = time(&t2);
|
||
tinv = t2 - t1;
|
||
|
||
if (X == NULL)
|
||
printf( "A is a singular matrix\n" );
|
||
else
|
||
{
|
||
mat_dumpf(X, "%+06.1f ");
|
||
|
||
printf( "|- A x Inv A -|\n");
|
||
t1 = time(&t1);
|
||
M = mat_mul( X, A );
|
||
t2 = time(&t2);
|
||
tmul = t2 - t1;
|
||
mat_dumpf( M, "%+06.1f " );
|
||
|
||
mat_free(M);
|
||
mat_free(X);
|
||
}
|
||
|
||
fgetmat( B, fp );
|
||
printf( "|- Matrix B -|\n");
|
||
mat_dumpf( B, "%+06.1f " );
|
||
|
||
printf( "|- A x B -|\n");
|
||
mat_free(mat_dumpf(mat_mul(A, B), "%+06.1f "));
|
||
|
||
printf( "time for finding 10 x 10 matrix inversion is less than %d secs\n", tinv );
|
||
printf( "time for finding 10 x 10 matrix determinant is less than %d secs\n", tdet );
|
||
printf( "time for finding 10 x 10 matrix multiplication is less than %d secs\n", tmul );
|
||
|
||
mat_free( A );
|
||
mat_free( B );
|
||
|
||
fclose(fp);
|
||
|
||
}
|
||
|