Solution.h

Go to the documentation of this file.

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

#ifndef __SOLUTION_H
#define __SOLUTION_H

#include <fstream>
#include <dolfin/constants.h>
#include <dolfin/Element.h>
#include <dolfin/NewArray.h>
#include <dolfin/Variable.h>

namespace dolfin {

  class ODE;
  class RHS;
  class Function;
  class ElementData;


  class Solution : public Variable
  {
  public:
    
    Solution(ODE& ode, Function& u);

    ~Solution();

    Element* createElement(Element::Type type, unsigned int q, unsigned int index, real t0, real t1);

    Element* element(unsigned int i, real t);

    Element* first(unsigned int i);

    Element* last(unsigned int i);

    real operator() (unsigned int i, unsigned int node, real t);

    real u(unsigned int i, real t);

    real u(unsigned int i, unsigned int node, real t);

    real k(unsigned int i, real t);

    real r(unsigned int i, real t, RHS& f);

    unsigned int size() const;

    Element::Type method(unsigned int i);

    unsigned int order(unsigned int i);

    void setlast(unsigned int i, real value);

    void shift(real t0);

    void save(real t0);

    void reset();

    enum Action { create = 0, update };
    void debug(Element& element, Action action);

  private:

    // The ODE
    ODE& ode;

    // Element data
    ElementData& elmdata;

    // Initial values (propagated values)
    NewArray<real> u0;

    // Time where current initial values are specified
    real t0;
    
    // Save debug info to file 'timesteps.debug'
    bool _debug;
    std::ofstream file;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.