timestepping.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: timestepping.h 1348 2013-08-16 00:45:07Z secanell $
//
//    Copyright (C) 2006, 2007, 2008 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__timestepping_h
#define __deal2__appframe__timestepping_h

#include <appframe/base.h>
#include <appframe/residual.h>
#include <appframe/timestep_control.h>

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

namespace AppFrame
{
  class ThetaTimestepping : public ApplicationCopy
  {
    public:
      static const Event new_time;
      
      static const Event new_step_size;
      
      ThetaTimestepping (ApplicationBase& app);
      
      ThetaTimestepping (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);
      virtual void initialize (ParameterHandler& param);
      
      void _initialize (ParameterHandler& param);

    protected:
      TimestepControl control;
      
    private:
      double theta;
      bool adaptive;

      ApplicationResidual default_residual;
      
      SmartPointer<const Residual> compute_residual;      
  };


  class ThetaResidual : public Residual
  {
    public:
      ThetaResidual(ApplicationBase& application);
      
      virtual double operator() (FEVector& dst, const FEVectors& src) const;
      
    private:
      SmartPointer<ApplicationBase> app;
  };



}


#endif