
double ddot0(int n, double dx[], int incx, double dy[], int incy);
void dgesl0(double **a, int n, int * ipvt, double b[], int job);
void dgefa0(double ** a, int n, int * ipvt, int * info);
void daxpy0(int n, double da, double dx[], int incx, double dy[], int incy);
int idamax0(int n, double dx[], int incx);
void dscal0(int n, double da, double dx[], int incx);
double vmnorm0(int n, double *v, double *w);
double fnorm0(int n, double **a, double *w);
#if 0
static double vmnorm0(int n, double *v, double *w)

/*
   This function routine computes the weighted max-norm
   of the vector of length n contained in the array v, with weights
   contained in the array w of length n.

   vmnorm = max( i = 1, ..., n ) fabs( v[i] ) * w[i].
*/

{
	int             i;
	double          vm;

	vm = 0.;
	for (i = 1; i <= n; i++)
		vm = fmax(vm, fabs(v[i]) * w[i]);
	return vm;

}
#endif
