112 lines
3.8 KiB
C++
112 lines
3.8 KiB
C++
/*
|
|
* Defines some prototypes for UMFPACK functions
|
|
*/
|
|
|
|
/*
|
|
* Copyright © 2013-2023 Dynare Team
|
|
*
|
|
* This file is part of Dynare.
|
|
*
|
|
* Dynare is free software: you can redistribute it and/or modify
|
|
* it under the terms of the GNU General Public License as published by
|
|
* the Free Software Foundation, either version 3 of the License, or
|
|
* (at your option) any later version.
|
|
*
|
|
* Dynare is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
* GNU General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU General Public License
|
|
* along with Dynare. If not, see <https://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#ifndef DYNUMFPACK_H
|
|
#define DYNUMFPACK_H
|
|
|
|
#ifdef OCTAVE_MEX_FILE
|
|
# ifdef HAVE_SUITESPARSE_UMFPACK_H
|
|
# include <suitesparse/umfpack.h>
|
|
# endif
|
|
# ifdef HAVE_UMFPACK_H
|
|
# include <umfpack.h>
|
|
# endif
|
|
#else
|
|
|
|
/* Under MATLAB, we have to provide our own header file for functions in
|
|
libmwumfpack, since there is no associated header */
|
|
|
|
# ifdef __cplusplus
|
|
extern "C"
|
|
{
|
|
# endif
|
|
|
|
/* -------------------------------------------------------------------------- */
|
|
/* size of Info and Control arrays */
|
|
/* -------------------------------------------------------------------------- */
|
|
|
|
/* These might be larger in future versions, since there are only 3 unused
|
|
* entries in Info, and no unused entries in Control. */
|
|
|
|
# define UMFPACK_INFO 90
|
|
# define UMFPACK_CONTROL 20
|
|
/* used in all UMFPACK_report_* routines: */
|
|
# define UMFPACK_PRL 0 /* print level */
|
|
/* returned by all routines that use Info: */
|
|
# define UMFPACK_OK (0)
|
|
# define UMFPACK_STATUS 0 /* UMFPACK_OK, or other result */
|
|
|
|
# ifdef _WIN64
|
|
# ifdef __cplusplus
|
|
# include <cstdint>
|
|
# else
|
|
# include <stdint.h>
|
|
# endif
|
|
typedef int64_t SuiteSparse_long;
|
|
# else
|
|
typedef long SuiteSparse_long;
|
|
# endif
|
|
|
|
void umfpack_dl_defaults(double Control[UMFPACK_CONTROL]);
|
|
|
|
SuiteSparse_long umfpack_dl_symbolic(SuiteSparse_long n_row, SuiteSparse_long n_col,
|
|
const SuiteSparse_long Ap[], const SuiteSparse_long Ai[],
|
|
const double Ax[], void** Symbolic,
|
|
const double Control[UMFPACK_CONTROL],
|
|
double Info[UMFPACK_INFO]);
|
|
|
|
SuiteSparse_long umfpack_dl_numeric(const SuiteSparse_long Ap[], const SuiteSparse_long Ai[],
|
|
const double Ax[], void* Symbolic, void** Numeric,
|
|
const double Control[UMFPACK_CONTROL],
|
|
double Info[UMFPACK_INFO]);
|
|
|
|
SuiteSparse_long umfpack_dl_solve(SuiteSparse_long sys, const SuiteSparse_long Ap[],
|
|
const SuiteSparse_long Ai[], const double Ax[], double X[],
|
|
const double B[], void* Numeric,
|
|
const double Control[UMFPACK_CONTROL],
|
|
double Info[UMFPACK_INFO]);
|
|
|
|
void umfpack_dl_report_info(const double Control[UMFPACK_CONTROL],
|
|
const double Info[UMFPACK_INFO]);
|
|
|
|
void umfpack_dl_report_status(const double Control[UMFPACK_CONTROL], SuiteSparse_long status);
|
|
|
|
void umfpack_dl_free_symbolic(void** Symbolic);
|
|
|
|
void umfpack_dl_free_numeric(void** Numeric);
|
|
|
|
SuiteSparse_long umfpack_dl_load_symbolic(void** Symbolic, char* filename);
|
|
|
|
SuiteSparse_long umfpack_dl_load_numeric(void** Numeric, char* filename);
|
|
|
|
SuiteSparse_long umfpack_dl_save_symbolic(void* Symbolic, char* filename);
|
|
|
|
SuiteSparse_long umfpack_dl_save_numeric(void* Numeric, char* filename);
|
|
|
|
# ifdef __cplusplus
|
|
} /* extern "C" */
|
|
# endif
|
|
|
|
#endif // OCTAVE_MEX_FILE
|
|
#endif
|