newton.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: newton.h 1348 2013-08-16 00:45:07Z secanell $
//
//    Copyright (C) 2006, 2007, 2008, 2009 by Guido Kanschat
//
//    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_h
#define __deal2__appframe__newton_h

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


namespace dealii
{
  class ParameterHandler;
  template <typename T> class BlockVector;
}

namespace AppFrame
{
  template <int dim> class DoFApplication;
  template <int dim> class LocalIntegrator;
  class Newton : public ApplicationCopy
  {
    public:
      static const Event bad_derivative;
      
      Newton (ApplicationBase& app);

      Newton (ApplicationBase& app, const Residual& res);
      
      virtual double residual(FEVector& start,
                              const FEVectors& rhs);
        
      virtual void solve(FEVector& start,
                         const FEVectors& rhs);
        
      virtual void declare_parameters (ParameterHandler& param);
        
      void _initialize (ParameterHandler& param);
        
      virtual void initialize (ParameterHandler& param);
        
      double threshold(double new_value);

      ReductionControl control;
    private:
      ApplicationResidual default_residual;
      
      SmartPointer<const Residual> compute_residual;
      
      bool assemble_now;

      unsigned int n_stepsize_iterations;
      
      double assemble_threshold;
      bool debug_solution;
      bool debug_update;
      bool debug_residual;
      unsigned int debug;       
  };

}




#endif