DAE_solver.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: DAE_solver.h 2605 2014-08-15 03:36:44Z secanell $
//
//    Copyright (C) 2011 by Jason Boisvert
//
//    This file is subject to QPL and may not be  distributed
//    without copyright and license information. Please refer
//    to the file deal.II/doc/license.html for the  text  and
//    further information on this license.
//
//---------------------------------------------------------------------------

#ifndef __deal2__appframe__DAE_solver_h
#define __deal2__appframe__DAE_solver_h

#include <stdio.h>
#include <stdexcept>

#ifdef _OPENMP
        #include <omp.h>
#else
        #define omp_get_thread_num() 0
#endif

//Load COLDAE subroutines
extern "C" {
        void coldae_(int &, int &, int [], double &, double &,
                double [], int [], int [], double [], double [], int [],
                double [], int &,
                void (*)(double &, double [], double [], double []),
                void (*)(double &, double [], double [], double []),
                void (*)(int &, double [], double &),
                void (*)(int &, double [], double []),
                void (*)(double &, double [], double [], double []));
        void appsln_ (double &, double [], double [], double [], int []);
}



namespace FuelCell
{
namespace ApplicationCore
{

        //Types
        typedef void (*fsub_ptr)(double &, double [], double [], double []);
        typedef void (*dfsub_ptr)(double &, double [], double [], double []);
        typedef void (*gsub_ptr)(int &, double [], double &);
        typedef void (*dgsub_ptr)(int &, double [], double []);
        typedef void (*guess_ptr)(double &, double [], double [], double []);

        //functions

        void DAE_dummy_guess (double &x, double z[], double y[], double df[]);
        void for_to_c_matrix(int rows, int cols, double *fmat, double **cmat);
        void c_to_for_matrix(int rows, int cols, double ** cmat, double *fmat);



        class DAESolver
        {
                public:

                DAESolver(int m_comp, int ny, int m[], double a, double b, double zeta[],
                void (*fsub)(double &, double [], double [], double []),
                void (*dfsub)(double &, double [], double [], double []),
                void (*gsub)(int &, double [], double &),
                void (*dgsub)(int &, double [], double []),
                void (*guess)(double &, double [], double [], double [])=NULL);

                ~DAESolver() {clear_mem();}

                void set_tolerance(int ltol_size = 0, int *ltol  = NULL, double *tol = NULL);
                void set_linear(void);
                void set_collocation_points(int pnts);
                void set_initial_mesh_size(int pnts);
                void set_integer_space(int ispace_size = 0);
                void set_float_space(int fspace_size= 0, double *fspace = NULL);
                void set_output(int level);
                void set_fixpnts(int fixpnt_size= 1, double *fixpnt = NULL);


                void set_solver_control(int control);


                void set_DAE_index(int index);



                void use_simple_cont(void);

                inline int get_ODE_order(void) {return *ODEs_Orders;}


                int DAE_solve(void);

                void DAE_solution(double x, double z[], double y[]);

                int *return_integer_space(void);

                double *return_float_space(void);

                inline int get_size_final_mesh(void) { return (this->ispace[0] +1); }

                void get_copy_final_mesh (double *mesh);

//              void operator delete (void *p);

                inline void use_sol_as_guess() { this->ispace[8] = 3; }

                protected:
                void set_fsub(void (*fsub)(double &, double [], double [], double []));
                void set_dfsub(void (*dfsub)(double &, double [], double [], double []));
                void set_gsub(void (*)(int &, double [], double &));
                void set_dgsub(void (*)(int &, double [], double []));
                void set_guess(void (*)(double &, double [], double [], double []));


                void (*DAE_fsub) (double &, double [], double [], double []);
                void (*DAE_dfsub) (double &, double [], double [], double []);
                void (*DAE_gsub) (int &, double [], double &);
                void (*DAE_dgsub) (int &, double [], double []);
                void (*DAE_guess) (double &, double [], double [], double []);

                void set_prob_size(int num_ODEs, int num_Alg_Const, int *ODEs_Orders);
                void set_boundary_points(int a, int b);
                void set_side_conditions(int zeta_size, double *zeta );
                void set_ipar(void);

                void clear_mem(void);


                // user-supllied subroutines flags
                bool have_fsub;
                bool have_dfsub;
                bool have_gsub;
                bool have_dgsub;
                bool have_guess;


                double *zeta;
                int zeta_size;
                bool set_zeta;


                int *ltol;
                bool set_ltol;
                int ltol_size;

                double *tol;
                bool set_tol;

                int *ipar;

                bool linear;

                int collpnts;
                bool set_collpnts;

                int intialmeshsize;
                bool set_intialmeshsize;

                int *ispace;
                int ispace_size;
                bool set_ispace;

                double *fspace;
                int fspace_size;
                bool set_fspace;
                bool use_old_fspace;

                int output_level;

                double *fixpnt;
                bool set_fixpnt;
                int fixpnt_size;

                bool set_solvercontrol;
                int solvercontrol;

                int DAE_index;

                int num_ODEs;
                int num_Alg_Const;
                int *ODEs_Orders;
                double a;
                double b;
                bool use_cont;
        };


}

}

#endif