newton_w_3ppC.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: newton_w_3ppC.h 1354 2013-08-17 00:01:22Z secanell $
//
//    Copyright (C) 2009 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__newton_w_3ppC_h
#define __deal2__appframe__newton_w_3ppC_h

#include <appframe/base.h>
#include <appframe/residual.h>
#include <lac/solver_control.h>
#include <newton_base.h>


namespace AppFrame
{
class Newton3ppC : public newtonBase
  {
    public:

        static const AppFrame::Event bad_derivative; 

        Newton3ppC(ApplicationBase& app);

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

  private:

        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;
        int step;
       };
}

#endif