mesh_loop.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: mesh_loop.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__mesh_loop_h
#define __deal2__appframe__mesh_loop_h

#include <lac/solver_control.h>
#include <appframe/base.h>
#include <numerics/vector_tools.h>

#include <set>

namespace AppFrame
{
  template <int dim> class DoFApplication;
  
  class MeshLoop :
      public ApplicationCopy
  {
    public:
      static const Event new_mesh;
      
      MeshLoop(ApplicationBase& app);
      virtual void declare_parameters(ParameterHandler& param);
      virtual void initialize (ParameterHandler& param);      
      virtual void solve(FEVector& start, const FEVectors& rhs);
      virtual void data_out(const std::string& filename, const FEVectors& src);
      void set_name(std::string new_name);      
    private:
      virtual void setup(FEVector& v);

      virtual void do_remesh(FEVector& v);

      virtual void do_evaluate(const FEVectors&);
      
      ReductionControl control;
      std::string base_name;
  };

  template <int dim>
  class DoFMeshLoop :
      public MeshLoop
  {
    public:
      
      DoFMeshLoop(ApplicationBase& app, DoFApplication<dim>& dof_app);
      
      ~DoFMeshLoop();
      virtual void declare_parameters(ParameterHandler& param);      
      virtual void initialize (ParameterHandler& param);
      virtual void data_out(const std::string& filename, const FEVectors& src);
      
    private:
      virtual void setup(FEVector&);
      virtual void do_remesh(FEVector&);
      virtual void do_evaluate(const FEVectors&);
      SmartPointer<DoFApplication<dim> > dof_app;
      
      bool transfer_solution;

      bool intrinsic_error;

      bool output_difference;
      
      FEVector* old_solution;      
  };
}

#endif