RHS.h

Go to the documentation of this file.

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

#ifndef __RHS_H
#define __RHS_H

#include <limits>
#include <dolfin/constants.h>
#include <dolfin/Vector.h>
#include <dolfin/Event.h>

namespace dolfin {

  class ODE;
  class Solution;
  class Function;


  class RHS {
  public:

    RHS(ODE& ode, Solution& solution);

    RHS(ODE& ode, Function& function);

    ~RHS();

    unsigned int size() const;

    real operator() (unsigned int index, real t);

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

    // Compute derivative dfi/duj
    real dfdu(unsigned int i, unsigned int j, real t);

  private:

    // Update components that influence the current component at time t
    void update(unsigned int index, unsigned int node, real t);

    // Update when we use Solution
    void updateSolution(unsigned int index, unsigned int node, real t);

    // Update when we use Solution
    void updateFunction(unsigned int index, real t);

    // Check computed value
    inline real check(real value)
    {
      if ( value > -std::numeric_limits<real>::max() && 
           value <  std::numeric_limits<real>::max() )
        return value;

      illegal_number();
      return 0.0;
    }

    // Number of components
    unsigned int N;

    // The ODE
    ODE& ode;

    // Solution
    Solution* solution;

    // Function
    Function* function;
    
    // Solution vector
    Vector u;

    // Event for illegal value of right-hand side
    Event illegal_number;

  };
    
}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.