ReducedModel.h

Go to the documentation of this file.

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

#ifndef __REDUCED_MODEL_H
#define __REDUCED_MODEL_H

#include <dolfin/constants.h>
#include <dolfin/NewArray.h>
#include <dolfin/ODE.h>

namespace dolfin
{
  class Vector;
  class Function;

  
  class ReducedModel : public ODE
  {
  public:
    
    ReducedModel(ODE& ode);
    
    ~ReducedModel();

    real f(const Vector& u, real t, unsigned int i);

    real u0(unsigned int i);

    Element::Type method(unsigned int i);

    unsigned int order(unsigned int i);

    real timestep(unsigned int i);
    
    void update(RHS& f, Function& u, real t);

    void update(Solution& u, Adaptivity& adaptivity, real t);

    void save(Sample& sample);

  private:

    // The component-specific model
    class Model
    {
    public:
      
      Model();

      ~Model();

      real operator() () const;

      bool active() const;

      void inactivate();

      void computeModel(Vector& ubar, Vector& fbar, unsigned int i, 
                        real tau, ODE& ode);

    private:

      real g;
      bool _active;

    };

    // Compute averages
    void computeAverages(RHS& f, Function& u, Vector& fbar, Vector& ubar);
    
    // The given model
    ODE& ode;

    // The reduced model
    NewArray<Model> g;

    // True if model has been reduced
    bool reduced;

    // Length of running average
    real tau;

    // Number of samples for running average
    unsigned int samples;

    // Tolerance for active components
    real tol;
    
    // Averages of u and f
    Vector ubar;
    Vector fbar;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.