DAE_wrapper.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
// C++ Interface: DAE_wrapper.cc
//
// Description: Class used as a wrapper for all applications that will be
//                              solved using the DAE_solver.
//
// Author: Peter Dobson <pdobson@ualberta.ca>, (C) 2011
//              University of Alberta
//
// Copyright: See COPYING file that comes with this distribution
//
//---------------------------------------------------------------------------

#ifndef FUEL_CELL__DAE_WRAPPER__H
#define FUEL_CELL__DAE_WRAPPER__H

//------------------------------
// STD DECLARATIONS
//------------------------------
#include <iostream>


#include <boost/shared_ptr.hpp>
//------------------------------
// DEAL.II DECLARATIONS
//------------------------------
#include <deal.II/base/parameter_handler.h>

//------------------------------
// ALGLIB DECLARATIONS
//------------------------------
#include <integration.h>

//------------------------------
// FUEL CELL DECLARATIONS
//------------------------------
#include <contribs/DAE_solver.h>
#include <utils/logging.h>

// Must be very careful with this... it is now 'declared' as a global variable within FCST
// The 'definition' of the pointer in included in the .cc file.
// The implementation (assignment) of an object to the pointer should be within the derived class
#ifdef _OPENMP
        #include <omp.h>
#else
        #define omp_get_thread_num() 0
#endif
extern void *ptr_DAE_object[20];

using namespace dealii;
using namespace alglib;

namespace FuelCell
{
namespace ApplicationCore
{
        class DAEWrapper
        {
        public:

        DAEWrapper();

        ~DAEWrapper(){};

        double integrate(double lb, double ub, std::vector<double>& W, std::vector<double>& F);

        void get_quadrature_points (double lb, double ub, std::vector<double>& X, std::vector<double>& W, FuelCell::ApplicationCore::DAESolver* prob);

        virtual void setup_DAE_solver () = 0;

        virtual void fsub (double &, double [], double [], double []) = 0;


        virtual void dfsub (double &, double [], double [], double []) = 0;

        virtual void gsub (int &, double [], double &) = 0;

        virtual void dgsub (int &, double [], double []) = 0;

        virtual void guess (double &, double [], double [], double []) = 0;

        inline void verbosity(int i)
        {n_output = i;}

        void DAE_Error(int flag);

        void clear_memory()
        {
                delete [] mm;
                delete [] zeta;
                delete [] fixpnt;
                delete [] tol;
                delete [] ltol;
                delete prob;
                if (ptr_DAE_object[omp_get_thread_num()] != NULL)
                {
                        ptr_DAE_object[omp_get_thread_num()] = NULL;
                        //delete ptr_DAE_object;
                }
        }

        protected:

        int n_mesh;

        double *mesh;

        int n_colloc;

        int n_output;

        int n_comp;

        int n_y;

        int *mm;

        int m_star;

        double boundary_0;

        double boundary_1;

        double *zeta;

        FuelCell::ApplicationCore::DAESolver *prob;

        double *fixpnt;

        int *ltol;

        double *tol;

        double cm_to_m;

        double cm2_to_m2;

        double cm3_to_m3;
        };

}
}
#endif