picard_base.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2009-13 by Energy Systems Design Laboratory, University of Alberta
//
//    This software is distributed under the MIT License.
//    For more information, see the README file in /doc/LICENSE
//
//    - Class: picard_base.h
//    - Description: Base class for all variants of the Picard solver
//    - Developers: Mayank Sabharwal
//
//---------------------------------------------------------------------------

#ifndef __deal2__appframe__picard_base_h
#define __deal2__appframe__Picard_base_h

#include <deal.II/lac/solver_control.h>

#include <application_core/application_wrapper.h>

using namespace dealii;
using namespace FuelCell::ApplicationCore;

namespace FuelCell
{
namespace ApplicationCore
{
    class PicardBase:public ApplicationWrapper
    {
    public:

        PicardBase(ApplicationBase& app);

        virtual void declare_parameters (ParameterHandler& param);

        virtual void initialize (ParameterHandler& param);

        void initialize_initial_guess(BlockVector<double>& dst){};

        virtual void assemble();

        virtual double residual(FuelCell::ApplicationCore::FEVector& dst,
                                const FuelCell::ApplicationCore::FEVectors& rhs);

        virtual void solve(FuelCell::ApplicationCore::FEVector& u,
                           const FuelCell::ApplicationCore::FEVectors& in_vectors) = 0;

        


    protected:
        void debug_output(const FEVector& sol,
                          const FEVector& update,
                          const FEVector& residual) const;

        bool assemble_now;

        double abs_tolerance;
        
        double rel_tolerance;
        
        int maxsteps;
        bool debug_solution;

        bool debug_residual;
        unsigned int debug;
        unsigned int step;
        
        double numIter;

        std::vector<unsigned int> blocks;

        unsigned int n_blocks;
    };
}
}

#endif