fmi-rs 0.3.0

#![allow(non_camel_case_types, non_snake_case)]

use std::ffi::c_void;

use crate::sundials::{sundials_nonlinearsolver::SUNNonlinearSolver, sundials_nvector::N_Vector, sundials_types::{SUNContext, sunrealtype}};

// /* -----------------------------------------------------------------
//  * Programmer(s): Allan G. Taylor, Alan C. Hindmarsh, Radu Serban,
//  *                and Aaron Collier @ LLNL
//  * -----------------------------------------------------------------
//  * SUNDIALS Copyright Start
//  * Copyright (c) 2025-2026, Lawrence Livermore National Security,
//  * University of Maryland Baltimore County, and the SUNDIALS contributors.
//  * Copyright (c) 2013-2025, Lawrence Livermore National Security
//  * and Southern Methodist University.
//  * Copyright (c) 2002-2013, Lawrence Livermore National Security.
//  * All rights reserved.
//  *
//  * See the top-level LICENSE and NOTICE files for details.
//  *
//  * SPDX-License-Identifier: BSD-3-Clause
//  * SUNDIALS Copyright End
//  * -----------------------------------------------------------------
//  * This is the header file for the main IDA solver.
//  * -----------------------------------------------------------------*/

// #ifndef _IDA_H
// #define _IDA_H

// #include <ida/ida_ls.h>
// #include <stdio.h>
// #include <sundials/sundials_core.h>

// #ifdef __cplusplus /* wrapper to enable C++ usage */
// extern "C" {
// #endif

// /* -----------------
//  * IDA Constants
//  * ----------------- */

// /* itask */
// #define IDA_NORMAL   1
pub const IDA_NORMAL: i32 = 1;
// #define IDA_ONE_STEP 2

// /* icopt */
// #define IDA_YA_YDP_INIT 1
// #define IDA_Y_INIT      2

// /* return values */

// #define IDA_SUCCESS      0
pub const IDA_SUCCESS: i32 = 0;
// #define IDA_TSTOP_RETURN 1
pub const IDA_TSTOP_RETURN: i32 = 1;
// #define IDA_ROOT_RETURN  2
pub const IDA_ROOT_RETURN: i32 = 2;

// #define IDA_WARNING 99

// #define IDA_TOO_MUCH_WORK -1
// #define IDA_TOO_MUCH_ACC  -2
// #define IDA_ERR_FAIL      -3
// #define IDA_CONV_FAIL     -4

// #define IDA_LINIT_FAIL  -5
// #define IDA_LSETUP_FAIL -6
// #define IDA_LSOLVE_FAIL -7
// #define IDA_RES_FAIL    -8
// #define IDA_REP_RES_ERR -9
// #define IDA_RTFUNC_FAIL -10
// #define IDA_CONSTR_FAIL -11

// #define IDA_FIRST_RES_FAIL  -12
// #define IDA_LINESEARCH_FAIL -13
// #define IDA_NO_RECOVERY     -14
// #define IDA_NLS_INIT_FAIL   -15
// #define IDA_NLS_SETUP_FAIL  -16
// #define IDA_NLS_FAIL        -17

// #define IDA_MEM_NULL     -20
// #define IDA_MEM_FAIL     -21
// #define IDA_ILL_INPUT    -22
// #define IDA_NO_MALLOC    -23
// #define IDA_BAD_EWT      -24
// #define IDA_BAD_K        -25
// #define IDA_BAD_T        -26
// #define IDA_BAD_DKY      -27
// #define IDA_VECTOROP_ERR -28

// #define IDA_CONTEXT_ERR -29

// #define IDA_TOO_CLOSE -60

// #define IDA_UNRECOGNIZED_ERROR -99

// /* ------------------------------
//  * User-Supplied Function Types
//  * ------------------------------ */

// typedef int (*IDAResFn)(sunrealtype tt, N_Vector yy, N_Vector yp, N_Vector rr,
//                         void* user_data);
pub type IDAResFn = unsafe extern "C" fn(tt: sunrealtype, yy: N_Vector, yp: N_Vector, rr: N_Vector,
    user_data: *mut c_void) -> i32;

// typedef int (*IDARootFn)(sunrealtype t, N_Vector y, N_Vector yp,
//                          sunrealtype* gout_1d, void* user_data);
pub type IDARootFn = unsafe extern "C" fn(t: sunrealtype, y: N_Vector, yp: N_Vector,
    gout_1d: *mut sunrealtype, user_data: *mut c_void) -> i32;

// typedef int (*IDAEwtFn)(N_Vector y, N_Vector ewt, void* user_data);

// /* -------------------
//  * Exported Functions
//  * ------------------- */

unsafe extern "C" {
// /* Initialization functions */

// SUNDIALS_EXPORT void* IDACreate(SUNContext sunctx);
pub unsafe fn IDACreate(sunctx: SUNContext) -> *mut c_void;

// SUNDIALS_EXPORT int IDAInit(void* ida_mem, IDAResFn res, sunrealtype t0,
//                             N_Vector yy0, N_Vector yp0);
pub unsafe fn IDAInit(ida_mem: *mut c_void, res: IDAResFn, t0: sunrealtype,
    yy0: N_Vector, yp0: N_Vector) -> i32;

// SUNDIALS_EXPORT int IDAReInit(void* ida_mem, sunrealtype t0, N_Vector yy0,
//                               N_Vector yp0);
pub unsafe fn IDAReInit(ida_mem: *mut c_void, t0: sunrealtype, yy0: N_Vector, 
    yp0: N_Vector) -> i32;

// /* Tolerance input functions */

// SUNDIALS_EXPORT int IDASStolerances(void* ida_mem, sunrealtype reltol,
//                                     sunrealtype abstol);
// SUNDIALS_EXPORT int IDASVtolerances(void* ida_mem, sunrealtype reltol,
//                                     N_Vector abstol);
pub unsafe fn IDASVtolerances(ida_mem: *mut c_void, reltol: sunrealtype,
     abstol: N_Vector) -> i32;

// SUNDIALS_EXPORT int IDAWFtolerances(void* ida_mem, IDAEwtFn efun);

// /* Initial condition calculation function */
// SUNDIALS_EXPORT int IDACalcIC(void* ida_mem, int icopt, sunrealtype tout1);

// /* Command-line control over scalar-valued set routines */
// SUNDIALS_EXPORT int IDASetOptions(void* ida_mem, const char* idaid,
//                                   const char* file_name, int argc, char* argv[]);

// /* Initial condition calculation optional input functions */

// SUNDIALS_EXPORT int IDASetNonlinConvCoefIC(void* ida_mem, sunrealtype epiccon);
// SUNDIALS_EXPORT int IDASetMaxNumStepsIC(void* ida_mem, int maxnh);
// SUNDIALS_EXPORT int IDASetMaxNumJacsIC(void* ida_mem, int maxnj);
// SUNDIALS_EXPORT int IDASetMaxNumItersIC(void* ida_mem, int maxnit);
// SUNDIALS_EXPORT int IDASetLineSearchOffIC(void* ida_mem, sunbooleantype lsoff);
// SUNDIALS_EXPORT int IDASetStepToleranceIC(void* ida_mem, sunrealtype steptol);
// SUNDIALS_EXPORT int IDASetMaxBacksIC(void* ida_mem, int maxbacks);

// /* Optional input functions */

// SUNDIALS_EXPORT int IDASetDeltaCjLSetup(void* ida_max, sunrealtype dcj);
// SUNDIALS_EXPORT int IDASetUserData(void* ida_mem, void* user_data);
pub unsafe fn IDASetUserData(ida_mem: *mut c_void, user_data: *mut c_void) -> i32;

// SUNDIALS_EXPORT int IDASetMaxOrd(void* ida_mem, int maxord);
// SUNDIALS_EXPORT int IDASetMaxNumSteps(void* ida_mem, long int mxsteps);
// SUNDIALS_EXPORT int IDASetInitStep(void* ida_mem, sunrealtype hin);
// SUNDIALS_EXPORT int IDASetMaxStep(void* ida_mem, sunrealtype hmax);
// SUNDIALS_EXPORT int IDASetMinStep(void* ida_mem, sunrealtype hmin);
// SUNDIALS_EXPORT int IDASetStopTime(void* ida_mem, sunrealtype tstop);
// SUNDIALS_EXPORT int IDAClearStopTime(void* ida_mem);
// SUNDIALS_EXPORT int IDASetMaxErrTestFails(void* ida_mem, int maxnef);
// SUNDIALS_EXPORT int IDASetSuppressAlg(void* ida_mem, sunbooleantype suppressalg);
// SUNDIALS_EXPORT int IDASetId(void* ida_mem, N_Vector id);
// SUNDIALS_EXPORT int IDASetConstraints(void* ida_mem, N_Vector constraints);
// SUNDIALS_EXPORT int IDASetMaxNumConstraintFails(void* ida_mem, int max_fails);

// /* Optional step adaptivity input functions */

// SUNDIALS_EXPORT
// int IDASetEtaFixedStepBounds(void* ida_mem, sunrealtype eta_min_fx,
//                              sunrealtype eta_max_fx);
// SUNDIALS_EXPORT
// int IDASetEtaMin(void* ida_mem, sunrealtype eta_min);
// SUNDIALS_EXPORT
// int IDASetEtaMax(void* ida_mem, sunrealtype eta_max);
// SUNDIALS_EXPORT
// int IDASetEtaLow(void* ida_mem, sunrealtype eta_low);
// SUNDIALS_EXPORT
// int IDASetEtaMinErrFail(void* ida_mem, sunrealtype eta_min_ef);
// SUNDIALS_EXPORT
// int IDASetEtaConvFail(void* ida_mem, sunrealtype eta_cf);

// /* Nonlinear solve input functions */

// SUNDIALS_EXPORT int IDASetMaxConvFails(void* ida_mem, int maxncf);
// SUNDIALS_EXPORT int IDASetMaxNonlinIters(void* ida_mem, int maxcor);
// SUNDIALS_EXPORT int IDASetNlsResFn(void* IDA_mem, IDAResFn res);
// SUNDIALS_EXPORT int IDASetNonlinConvCoef(void* ida_mem, sunrealtype epcon);
// SUNDIALS_EXPORT int IDASetNonlinearSolver(void* ida_mem, SUNNonlinearSolver NLS);
pub unsafe fn IDASetNonlinearSolver(ida_mem: *mut c_void, NLS: SUNNonlinearSolver) -> i32;

// /* Rootfinding initialization function */

// SUNDIALS_EXPORT int IDARootInit(void* ida_mem, int nrtfn, IDARootFn g);
pub unsafe fn IDARootInit(ida_mem: *mut c_void, nrtfn: i32, g: IDARootFn) -> i32;

// /* Rootfinding optional input functions */

// SUNDIALS_EXPORT int IDASetRootDirection(void* ida_mem, int* rootdir_1d);
// SUNDIALS_EXPORT int IDASetNoInactiveRootWarn(void* ida_mem);

// /* Solver function */
// SUNDIALS_EXPORT int IDASolve(void* ida_mem, sunrealtype tout, sunrealtype* tret,
//                              N_Vector yret, N_Vector ypret, int itask);
pub unsafe fn IDASolve(ida_mem: *mut c_void, tout: sunrealtype, tret: *mut sunrealtype,
        yret: N_Vector, ypret: N_Vector, itask: i32) -> i32;

// /* Utility functions to update/compute y and yp based on ycor */

// SUNDIALS_EXPORT int IDAComputeY(void* ida_mem, N_Vector ycor, N_Vector y);
// SUNDIALS_EXPORT int IDAComputeYp(void* ida_mem, N_Vector ycor, N_Vector yp);

// /* Dense output function */
// SUNDIALS_EXPORT int IDAGetDky(void* ida_mem, sunrealtype t, int k, N_Vector dky);

// /* Optional output functions */

// SUNDIALS_DEPRECATED_EXPORT_MSG(
//   "Work space functions will be removed in version 8.0.0")
// int IDAGetWorkSpace(void* ida_mem, long int* lenrw, long int* leniw);
// SUNDIALS_EXPORT int IDAGetNumSteps(void* ida_mem, long int* nsteps);
// SUNDIALS_EXPORT int IDAGetNumResEvals(void* ida_mem, long int* nrevals);
// SUNDIALS_EXPORT int IDAGetNumLinSolvSetups(void* ida_mem, long int* nlinsetups);
// SUNDIALS_EXPORT int IDAGetNumErrTestFails(void* ida_mem, long int* netfails);
// SUNDIALS_EXPORT int IDAGetNumBacktrackOps(void* ida_mem, long int* nbacktr);
// SUNDIALS_EXPORT int IDAGetConsistentIC(void* ida_mem, N_Vector yy0_mod,
//                                        N_Vector yp0_mod);
// SUNDIALS_EXPORT int IDAGetLastOrder(void* ida_mem, int* klast);
// SUNDIALS_EXPORT int IDAGetCurrentOrder(void* ida_mem, int* kcur);
// SUNDIALS_EXPORT int IDAGetCurrentCj(void* ida_mem, sunrealtype* cj);
// SUNDIALS_EXPORT int IDAGetCurrentY(void* ida_mem,
//                                    N_Vector* ycur); // nb::rv_policy::reference
// SUNDIALS_EXPORT int IDAGetCurrentYp(void* ida_mem,
//                                     N_Vector* ypcur); // nb::rv_policy::reference
// SUNDIALS_EXPORT int IDAGetActualInitStep(void* ida_mem, sunrealtype* hinused);
// SUNDIALS_EXPORT int IDAGetLastStep(void* ida_mem, sunrealtype* hlast);
// SUNDIALS_EXPORT int IDAGetCurrentStep(void* ida_mem, sunrealtype* hcur);
// SUNDIALS_EXPORT int IDAGetCurrentTime(void* ida_mem, sunrealtype* tcur);
// SUNDIALS_EXPORT int IDAGetTolScaleFactor(void* ida_mem, sunrealtype* tolsfact);
// SUNDIALS_EXPORT int IDAGetErrWeights(void* ida_mem, N_Vector eweight);
// SUNDIALS_EXPORT int IDAGetEstLocalErrors(void* ida_mem, N_Vector ele);
// SUNDIALS_EXPORT int IDAGetNumConstraintFails(void* ida_mem,
//                                              long int* num_fails_out);
// SUNDIALS_EXPORT int IDAGetNumConstraintCorrections(void* ida_mem,
//                                                    long int* num_corrections_out);
// SUNDIALS_EXPORT int IDAGetNumGEvals(void* ida_mem, long int* ngevals);
// SUNDIALS_EXPORT int IDAGetRootInfo(void* ida_mem, int* rootsfound_1d);
pub unsafe fn IDAGetRootInfo(ida_mem: *mut c_void, rootsfound_1d: *mut i32) -> i32;

// SUNDIALS_EXPORT int IDAGetIntegratorStats(void* ida_mem, long int* nsteps,
//                                           long int* nrevals, long int* nlinsetups,
//                                           long int* netfails, int* qlast,
//                                           int* qcur, sunrealtype* hinused,
//                                           sunrealtype* hlast, sunrealtype* hcur,
//                                           sunrealtype* tcur);
// SUNDIALS_EXPORT int IDAGetNonlinearSystemData(
//   void* ida_mem, sunrealtype* tcur, N_Vector* yypred, N_Vector* yppred,
//   N_Vector* yyn, N_Vector* ypn, N_Vector* res, sunrealtype* cj,
//   void** user_data); // nb::rv_policy::reference
// SUNDIALS_EXPORT int IDAGetNumNonlinSolvIters(void* ida_mem, long int* nniters);
// SUNDIALS_EXPORT int IDAGetNumNonlinSolvConvFails(void* ida_mem,
//                                                  long int* nnfails);
// SUNDIALS_EXPORT int IDAGetNonlinSolvStats(void* ida_mem, long int* nniters,
//                                           long int* nnfails);
// SUNDIALS_EXPORT int IDAGetNumStepSolveFails(void* ida_mem, long int* nncfails);
// SUNDIALS_EXPORT int IDAGetUserData(void* ida_mem, void** user_data);
// SUNDIALS_EXPORT int IDAPrintAllStats(void* ida_mem, FILE* outfile,
//                                      SUNOutputFormat fmt);
// SUNDIALS_EXPORT char* IDAGetReturnFlagName(long int flag);

// /* Free function */
// SUNDIALS_EXPORT void IDAFree(void** ida_mem);
pub unsafe fn IDAFree(ida_mem: *mut *mut c_void);

// /* IDALS interface function that depends on IDAResFn */
// SUNDIALS_EXPORT int IDASetJacTimesResFn(void* ida_mem, IDAResFn jtimesResFn);
}

// #ifdef __cplusplus
// }
// #endif

// #endif