Solver.h

Go to the documentation of this file.

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

#ifndef __SOLVER_H
#define __SOLVER_H

#include <dolfin.h>

namespace dolfin {

  class Solver {
  public:

    Solver();

    Solver(Mesh& mesh_);

    Solver(ODE& ode_);
    
    virtual ~Solver();

    virtual const char* description() = 0;
    
    virtual void solve() = 0;
    
  protected:
    
    Mesh& mesh;
    ODE& ode;

  private:
    
    // Every solver needs to have both a mesh and an ODE (since the same solver
    // interface is used for different equations). We create dummy objects to
    // be able to initialize the references if they are not used.

    class DummyODE : public ODE {
    public:
      DummyODE() : ODE(1) {}
      real u0(unsigned int i) { return 0.0; }
      real f(const Vector& u, real t, unsigned int i) { return 0.0; }
    };

    Mesh     dummy_mesh;
    DummyODE dummy_ode; 
    
  };
  
}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.