typedef	enum status {success, failure, error} status;

extern void	vectorcopy(int, double *, double *);
extern void	vectorfill(int, double *, double);
extern double	vectorvector(int, double *, double *);
extern void	scalarvector(int, double *, double, double *);
extern void	vectoradd(int, double *, double *, double *);
extern void	vectorsub(int, double *, double *, double *);

extern void	matrixvector(int, int, double *, double *, double *);
extern void	matrixrowvector(int, int, double *, double *, int);
extern void	matrixcolumnvector(int, int, double *, double *, int);

extern void	matrixcopy(int, int, double *, double *);
extern void	matrixunit(int, double *);
extern void	scalarmatrix(int, int, double *, double, double *);
extern void	matrixadd(int, int, double *, double *, double *);
extern void	matrixsub(int, int, double *, double *, double *);
extern void	matrixmatrix(int, int, int, double *, double *, double *);
extern void	matrixtranspose(int, int, double *, double *);

extern void	matrixrowscalar(int, int, double *, int, double);
extern void	matrixrowexchange(int, int, double *, int, int);
extern void	matrixrowadd(int, int, double *, int, double, int);
extern void	matrixrowaddtwo(int, int, int, double *, double *, int, double, int);
extern void	matrixsearchcolumnmaxabs(int, int, double *, int, int, int,
					 int *, double *);

extern void	vectorfprint(FILE *, int, double *);
extern void	vectorfscan(FILE *, int, double *);
extern void	matrixfprint(FILE *, int, int, double *);
extern void	matrixfscan(FILE *, int, int, double *);
extern void	matrixprintformat(char *);

extern void	matrixinverse(int, double *, double *, double);

extern int	nullspace(int, int, double *, double *, int *, double);
extern int	svd(int, int, double *, double *, int *, double);