dirichlet.hh

Go to the documentation of this file.

 
#ifndef dirichlet_hh
#define dirichlet_hh
 
#include "basics/outputOperator.hh"
#include "basics/typedefs.hh"
#include "space/space.hh"
#include "space/element.hh"
#include "geometry/mesh.hh"
 
 
namespace concepts {
 
  // ************************************************************* Dirichlet **
 
  template<class F>
  class Dirichlet : public virtual OutputOperator {
  public:
    Dirichlet(Space<Real>& spc) : coeff_(0)  {}
    std::map<int,F>* operator[](const Element<F>* elm) { 
      return (coeff_) ? (*coeff_)[elm] : 0; }
    virtual ~Dirichlet();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual std::ostream& datainfo(std::ostream& os) const;
    // --> Idee:
    // Schlaufe über alle Formfunktionen, Eintrag für konkrete Funktion
    // nichtbenutzte Indizes bleiben Null
    // Alternative: Array, Schlaufe über benutzte Indizes,
    // Problematisch Zuordnung Index zur Formfunktion
 
    // nicht schlimm, laufe ueber alle Eintraege in Data
    //    Joiner<Array<Data>*, 1>* coeff_;
 
//      Joiner<SparseArray<F>*, 1>* coeff_;
 
    // Map okay, aber keine Ausgabe (info), geschrieben!
    //Joiner<std::map<int,F>*, 1>* coeff_;
 
 
    std::map<const Element<F>*, std::map<int,F>*>* coeff_;
    // schafft es Sparse'Matrix' auch?
    // es gibt Funktionen zum Multiplizieren mit Vektoren
    // das waere gut zum Ausrechnen der Funktionswerte
    // jedoch ist ein Raum (spc) nötig
    // Joiner<SparseMatrix<Real>*, 1>* coeff_;
  };
 
 
 
} // namespace concepts
 
#endif // dirichlet_hh