/* Steven Andrews, 9/13/98 */
/* See documentation called RnLU doc */
/* Copyright 2003 by Steven Andrews.  Permission is granted
   for non-commercial use of and modifications to the code. */

#ifndef __RnLU_h
#define __RnLU_h

float LUdecomp(float *a,int n,int **indx);
void LUsolveV(float *a,int *indx,float *b,int n);
void LUsolveM(float *a,int *indx,float *b,int n,int nb);
float LUlndetM(float *a,int *indx,int n,int *sgn);
float LUimprV(float *a,float *alu,int *indx,float *x,float *b,int n);
float LUimprM(float *a,float *alu,int *indx,float *x,float *b,int n,int nb);

#endif