/*
 * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 * Copyright (C) 2007 - INRIA - Bruno JOFRET
 *
 * This file must be used under the terms of the CeCILL.
 * This source file is licensed as described in the file COPYING, which
 * you should have received as part of this distribution.  The terms
 * are also available at
 * http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
 *
 */

#ifndef __LU_H__
#define __LU_H__

#define spINSIDE_SPARSE
#include "spConfig.h"
#include "spmatrix.h"
#include "spDefs.h"

#include "cerro.h"
#include "machine.h"
#include "localization.h"
#include "MALLOC.h"
#include "dynlib_sparse.h"

SPARSE_IMPEXP int addluptr (char *ptr); /* */
SPARSE_IMPEXP int getluptr (int sel, char **ptr);
SPARSE_IMPEXP int removeluptr (int sel);
SPARSE_IMPEXP void resetluptr (void); /* to be used to free the lu pointer table */

SPARSE_IMPEXP void C2F(lufact1)(double *val, int *lln, int *col, int *n, int *nel,
                                int *fmatindex, double *eps, double *releps, int *nrank, int *ierr);

SPARSE_IMPEXP void C2F(lusolve1)(int *fmatindex, double *b, double *x, int *ierr);

SPARSE_IMPEXP void C2F(ludel1)(int *fmatindex, int *ierr);

SPARSE_IMPEXP void C2F(lusiz1)(int *fmatindex, int* lsize, int* usize, int *ierr);

SPARSE_IMPEXP void C2F(luget1)(int *fmatindex, int *indP, double *P,
                               int *indl, double *l, int *indu, double *u,
                               int *indQ, double *Q, int *ierr);
#endif /* !__LU_H__ */