fps/fptest/linpack.c

45 static int LINSUB(REAL * residn, REAL * resid,
46     REAL * eps, REAL * x11, REAL * xn1, int fps_verbose_msg);
47 static int MATGEN(REAL a[], int lda, int n, REAL b[], REAL * norma);
48 static REAL EPSLON(REAL x);
49 static void MXPY(int n1, REAL y[], int n2, int ldm, REAL x[], REAL m[]);
71 	REAL eps;
72 	REAL resid;
73 	REAL residn;
74 	REAL x11;
75 	REAL xn1;
76 	REAL EPS;
77 	REAL RESID;
78 	REAL RESIDN;
79 	REAL X11;
80 	REAL XN1;
202  * LINSUB(REAL *residn, REAL *resid, REAL *eps,
203  * REAL *x11, REAL *xn1, int fps_verbose_msg)begins
207 LINSUB(REAL *residn, REAL *resid,
208 	REAL *eps, REAL *x11, REAL *xn1,
215 	REAL *a;
216 	REAL abs;
217 	REAL *b;
218 	REAL norma;
219 	REAL normx;
220 	REAL *x;
233 	a = (REAL *) malloc((MAT_SIZE + 8) * (MAT_SIZE + 1) *
234 	    (size_t)sizeof (REAL));
235 	b = (REAL *) malloc(MAT_SIZE * (size_t)sizeof (REAL));
236 	x = (REAL *) malloc(MAT_SIZE * (size_t)sizeof (REAL));
276 	(void) MATGEN((REAL *) a, lda, n, b, &norma);
282 	MXPY(n, b, n, lda, x, (REAL *) a);
289 		abs = (REAL)fabs((double)b[i]);
291 		abs = (REAL)fabs((double)x[i]);
297 	*eps = EPSLON((REAL) LP_ONE);
364  * MATGEN(REAL a[], int lda, int n, REAL b[], REAL *norma)
370 MATGEN(REAL a[], int lda, int n, REAL b[], REAL *norma)
375 	REAL	value;
391 			value = (REAL) dran(init) - 0.5;
415  * EPSLON(REAL x) estimates unit roundoff in
418 static REAL
419 EPSLON(REAL x)
421 	REAL a;
422 	REAL abs;
423 	REAL b;
424 	REAL c;
425 	REAL eps;
433 		eps = (REAL)fabs((double)(c - LP_ONE));
436 	abs = (REAL)fabs((double)x);
442  * MXPY(int n1, REAL y[], int n2, int ldm, REAL x[], REAL m[])
447 MXPY(int n1, REAL y[], int n2, int ldm, REAL x[], REAL m[])