newton_w_3pp.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: newton_w_3pp.h
//    - Description: Base class for all variants of the Newton-Raphson solver
//    - Developers: Jason Boisvert and M. Secanell
//    - $Id: newton_w_3pp.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef __deal2__appframe__newton_w_3pp_h
#define __deal2__appframe__newton_w_3pp_h

#include <newton_base.h>

namespace FuelCell
{
namespace ApplicationCore
{
    class Newton3pp : public newtonBase
    {
    public:

        static const FuelCell::ApplicationCore::Event bad_derivative;

        Newton3pp(ApplicationBase& app);

        void declare_parameters(ParameterHandler& param);

        void initialize(ParameterHandler& param);

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

    private:
        inline void set_predefined_parameters()
        {
            // Depends in which refinement
            const double* ref = get_data()->scalar("Refinement");
            //FcstUtilities::log << "Ref: "<< *ref <<"\n";
            if (*ref == 0)
            {
                max_reduction_factor = 0.90;
                min_reduction_factor = 0.01;
                alpha=0.01;
            }
            else
            {
                max_reduction_factor = 0.90;
                min_reduction_factor = 0.50;
                alpha=0.001;
            }
        }

        double sum_of_squares(double residual_norm);

        double three_point_step (int line_search_iterations,double lambda_current, double lambda_previous, double mf_original, double mf_current, double mf_previous);

        bool armijo(double step_length,double residual_norm,double new_residual_norm);

        double max_reduction_factor;
        double min_reduction_factor;
        double alpha;

        bool with_prediction;

        bool use_predefined;
    };
}
}

#endif