Adaptivity.h

Go to the documentation of this file.

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

#ifndef __ADAPTIVITY_H
#define __ADAPTIVITY_H

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

namespace dolfin {

  class ODE;
  class TimeSlab;
  class Solution;
  class RHS;


  class Adaptivity {
  public:

    Adaptivity(ODE& ode);

    ~Adaptivity();

    Regulator& regulator(unsigned int i);

    const Regulator& regulator(unsigned int i) const;

    real tolerance() const;

    real maxstep() const;

    real minstep() const;
    
    void stabilize(real k, unsigned int m);
    
    bool fixed() const;

    real threshold() const;

    unsigned int size() const;

    void shift(Solution& u, RHS& f);

    bool accept(TimeSlab& timeslab, RHS& f);

    void adjustMaximumTimeStep(real kmax);

  private:

    // Regulators, one for each component
    NewArray<Regulator> regulators;

    // Tolerance
    real TOL;

    // Maximum and minimum allowed time step
    real kmax;

    // Current maximum time step
    real kmax_current;

    // Flag for fixed time steps
    bool kfixed;
    
    // Threshold for reaching end of interval
    real beta;

    // Remaining number of small time steps
    unsigned int m;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.