app-kersten/Models.hh

Go to the documentation of this file.

#ifndef hp2DModels_hh
#define hp2DModels_hh

#include "basics/exceptions.hh"
#include "function/vector.hh"
#include "graphics/matlab.hh"
#include "hp2D/space.hh"
#include "toolbox/inputOutput.hh"
#include "toolbox/inputParam.hh"
#include "formula/piecewiseFormula.hh"

#include "basics/debug.hh"

#define Models_formulaOut_D 0

namespace hp2D {

  class InputModels;

  // ********************************************************* SpaceNotBuilt **

  class ModelNotSolved : public concepts::ExceptionBase {};

  // ********************************************************* SpaceNotBuilt **

  class MatrixNotBuilt : public concepts::ExceptionBase {};

  // **************************************************************** Models **
  
  template<class F>
  class Models : public concepts::OutputOperator {
  public:
    Models(const std::string problemName, const uint geomRefAttr = 100);

    virtual ~Models() {}
    void solve(const uint l = 0, const uint p = 1, const uint g = 0,
         const std::string vectorFile = "",
         const std::string matrixFile = "");
    void solve(InputModels& input);
    void matrixBuild(const uint l = 0, const uint p = 1, const uint g = 0,
         const std::string matrixFile = "");
    void matrixBuild(InputModels& input);
    void spaceBuild(const uint l = 0, const uint p = 1, const uint g = 0);
    virtual void loadSolution(const std::string vectorFile);
    virtual void loadSolution(const uint l = 0, const uint p = 1,
            const uint g = 0, const std::string path = "");
    virtual void graphicsOut(const uint numPoints,
           const std::string path = "") = 0;
    virtual void dataOut(concepts::InOutParameters& parameters,
       const bool store = false,
       const std::string path = "");
    void basisFunctions(const std::string path = "");
    const hp2D::Space* space() const { return spc_.get(); }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual concepts::Mesh2& mesh_() = 0;
    virtual const std::string mshAbbr_() = 0;
    const std::string problemName_;
    const uint geomRefAttr_;
    std::auto_ptr<hp2D::Space> spc_;
    std::auto_ptr<concepts::BoundaryConditions> bc_;
    std::auto_ptr<concepts::Vector<F> > sol_;
    std::string matrixFile_;
    std::string vectorFile_;
    virtual void solve_() = 0;
    virtual void matrices_() = 0;
    virtual void dataOut_() = 0;
    std::string filename_(const std::string varname,
        const std::string subpath, const bool p = false);
    template<class G>
    void formulaOut_(const concepts::PiecewiseFormulaBase<G>& frm,
                     const std::string frmstr, const std::string path = "");
    concepts::InOutParameters parameters_;
  private:
    void refine_(uint attrib, int pMax[2]);
    void newfile_(std::string& name, const std::string varname,
      const std::string ending);
    uint nextCounter_(const std::string path = "") const;
    std::string nodenameLetter_() const;
  };

  template<class F>
  template<class G>
  void Models<F>::formulaOut_(const concepts::PiecewiseFormulaBase<G>& frm,
                              const std::string frmstr,
                              const std::string path) {
    if (!spc_.get()) spaceBuild();
    std::string dname =
      filename_(frmstr, path + "/" + problemName_ + "/Grafik");
    DEBUGL(Models_formulaOut_D,
     "Writing formula " << frmstr << " into " << dname << ".m");
    graphics::MatlabGraphics(*spc_, dname.c_str(), frm);

    parameters_.addString((std::string("grafik_") + frmstr).c_str(),
        (dname + ".m").c_str());
  }

  // *********************************************************** InputModels **

  class InputModels : public concepts::InputParameter {
  public:
    InputModels(concepts::InOutParameters& input);
    virtual std::ostream& letters(std::ostream& os) const;
    virtual std::ostream& arguments(std::ostream& os) const;
    virtual std::ostream& description(std::ostream& os) const;
    virtual int input(int opt, const char* optarg);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };

} // namespace hp2D

#endif // hp2DModels_hh