newton_base.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: newton_base.h 1354 2013-08-17 00:01:22Z secanell $
//
//    Copyright (C) 2007 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_base_h
#define __deal2__appframe__newton_base_h

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

using namespace dealii;
using namespace AppFrame;

namespace AppFrame
{
        class newtonBase:public ApplicationCopy
                {
                public:
                        static const AppFrame::Event bad_derivative; 
                        newtonBase(ApplicationBase& app);
                        
                        virtual void assemble();
                        
                        virtual double residual(AppFrame::FEVector& dst,
                                                                        const AppFrame::FEVectors& rhs);
                        
                        virtual void solve(AppFrame::FEVector& u,
                                                           const AppFrame::FEVectors& in_vectors);
                        virtual void declare_parameters (ParameterHandler& param);
                        
                        void _initialize (ParameterHandler& param);
                        
                        virtual void initialize (ParameterHandler& param);
                        
                        double threshold(double new_value);
                        
                        void initialize_initial_guess(BlockVector<double>& dst);

                        ReductionControl control;
                        
                        // Output data
                        //int numIterations;
                        double numIter;
                protected:
                        ApplicationResidual default_residual;
                        
                        SmartPointer<const Residual> compute_residual;
                        bool assemble_now;
                        
                        double assemble_threshold;
                        bool debug_solution;
                        bool debug_update;
                        bool debug_residual;
                        unsigned int debug;     
                };
}

#endif