RecursiveTimeSlab.h

Go to the documentation of this file.

// Copyright (C) 2003 Johan Hoffman and Anders Logg.
// Licensed under the GNU GPL Version 2.
//
// Updates by Johan Jansson 2003

#ifndef __RECURSIVE_TIME_SLAB_H
#define __RECURSIVE_TIME_SLAB_H

#include <dolfin/NewArray.h>
#include <dolfin/TimeSlab.h>

namespace dolfin {

  class Element;
  class Adaptivity;
  class Partition;
  class RHS;
  class Solution;
  class FixedPointIteration;


  class RecursiveTimeSlab : public TimeSlab {
  public:

    RecursiveTimeSlab(real t0, real t1, Solution& u, RHS& f, 
                      Adaptivity& adaptivity, FixedPointIteration& fixpoint,
                      Partition& partition, int offset);
    
    ~RecursiveTimeSlab();

    bool accept(RHS& f, real TOL);
    
    void countElementGroups(unsigned int& size);
    
    void addElementGroups(NewArray<ElementGroup*>& groups, unsigned int& pos);
    
    void show(unsigned int depth = 0) const;

  private:
    
    // Create new time slab
    void create(Solution& u, RHS& f, Adaptivity& adaptivity,
                FixedPointIteration& fixpoint,
                Partition& partition, int offset);
    
    // Create list of time slabs within the time slab
    void createTimeSlabs(Solution& u, RHS& f, Adaptivity& adaptivity,
                         FixedPointIteration& fixpoint,
                         Partition& partition, int offset);

    // Create list of elements within the time slab
    void createElements(Solution& u, RHS& f, Adaptivity& adaptivity,
                        FixedPointIteration& fixpoint,
                        Partition& partition, int offset, int end);

    // Compute residuals and new time steps
    void computeResiduals(RHS& f, Adaptivity& adaptivity);

    //--- Time slab data ---

    // List of time slabs within this time slab
    NewArray<TimeSlab*> timeslabs;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.