matlabBinaryGraphics.hh

Go to the documentation of this file.

 
#ifndef matlab_binary_hh
#define matlab_binary_hh
 
#include "matfile/matfileIO.hh"
 
#include "basics/outputMatlab.hh"
#include "basics/vectorsMatrices.hh"
#include "basics/vectorsMatricesForward.hh"
#include "denseMatrixGraphics.hh"
#include "matfile/matfileIO.hh"
 
namespace graphics {
 
  // ************************************************** MatlabBinaryGraphics **
 
  class MatlabBinaryGraphics: public concepts::MatfileIO,
                              public DenseMatrixGraphics {
  public:
 
    MatlabBinaryGraphics(const std::string filename);
 
    template<class G>
    MatlabBinaryGraphics(const concepts::Space<G>& spc,
                         const std::string filename);
 
    MatlabBinaryGraphics(concepts::Mesh& msh, const std::string filename,
                         const uint points = 5);
 
    template<class G, class F>
    MatlabBinaryGraphics(const concepts::Space<G>& spc,
                         const std::string filename,
                         const concepts::Vector<F>& sol,
                         const concepts::ElementFunction<F, G>* fun = 0);
 
    template<class G, class F>
    MatlabBinaryGraphics(const concepts::Space<G>& spc,
                         const std::string filename,
                         const concepts::Vector<F>& sol,
                         const std::string varName,
                         const concepts::ElementFunction<F, G>* fun = 0);
 
    template<class G, class F>
    MatlabBinaryGraphics(const concepts::Space<G>& spc,
                         const std::string filename,
                         std::vector<concepts::Vector<F> >& sol,
                         const concepts::ElementFunction<F, G>* fun = 0);
 
    template<class G, class F>
    MatlabBinaryGraphics(const concepts::Space<G>& spc,
                         const std::string filename,
                         const concepts::ElementFormula<F, G>& frm);
 
    virtual ~MatlabBinaryGraphics();
 
    template<class F, class G>
    void addSolution(const concepts::Space<G>& spc, const std::string name,
        const concepts::Vector<F>& sol,
        const concepts::ElementFunction<F, G>* fun = 0);
 
    template<class F, class G>
    void addFunction(const concepts::Space<G>& spc, const std::string name,
        const concepts::ElementFormula<F, G>& frm);
  protected:
    template<class G>
    void addSpace_(const concepts::Space<G>& spc);
 
    template<class G>
    void storeData_(DenseMatrixCollection<G>& dense_ptr_);
 
    uint data_;
  };
 
  template<class G, class F>
  MatlabBinaryGraphics::MatlabBinaryGraphics(const concepts::Space<G>& spc,
                                             const std::string filename,
                                             const concepts::Vector<F>& sol,
                                             const concepts::ElementFunction<F, G>* fun)
    : concepts::MatfileIO(filename), DenseMatrixGraphics()
  {
    this->addSpace_(spc);
    this->addSolution(spc, "u", sol, fun);
  }
 
 
  template<class G, class F>
  MatlabBinaryGraphics::MatlabBinaryGraphics(const concepts::Space<G>& spc,
                       const std::string filename,
                       const concepts::Vector<F>& sol,
                       const std::string varName,
                       const concepts::ElementFunction<F, G>* fun)
    :concepts::MatfileIO(filename), DenseMatrixGraphics()
  {
    this->addSpace_(spc);
    this->addSolution(spc, varName, sol, fun);
  }
 
  template<class G>
  MatlabBinaryGraphics::MatlabBinaryGraphics(const concepts::Space<G>& spc,
      const std::string filename) :
      concepts::MatfileIO(filename), DenseMatrixGraphics() {
 
    this->addSpace_(spc);
  }
 
  template<class G, class F>
  MatlabBinaryGraphics::MatlabBinaryGraphics(const concepts::Space<G>& spc,
      const std::string filename, const concepts::ElementFormula<F, G>& frm) :
      concepts::MatfileIO(filename) {
    this->addSpace_(spc);
    this->addFunction(spc, "u", frm);
  }
 
  template<class F, class G>
  void MatlabBinaryGraphics::addSolution(const concepts::Space<G>& spc,
                                         const std::string name,
                                         const concepts::Vector<F>& sol,
                                         const concepts::ElementFunction<F, G>* fun)
  {
    // get matrix of the solution
    concepts::RCP<concepts::ElementMatrix<F> > 
      function = this->getSolution(spc, sol, fun);
    // store matrix
    this->add(*function, name);
  }
 
  template<class F, class G>
  void MatlabBinaryGraphics::addFunction(const concepts::Space<G>& spc,
      const std::string name, const concepts::ElementFormula<F, G>& frm) {
    // get matrix for the formula
    concepts::RCP<concepts::ElementMatrix<typename concepts::Datatype<F>::type> > function =
        this->getFunction(spc, frm);
    // store matrix
    this->add(*function, name);
  }
 
  template<class G>
  void MatlabBinaryGraphics::addSpace_(const concepts::Space<G>& spc) {
    DenseMatrixCollection<G> dmc = this->getSpace(spc);
    //store matrix
    this->storeData_(dmc);
  }
 
  //Save the matrices that hold informations
  template<class G>
  void MatlabBinaryGraphics::storeData_(DenseMatrixCollection<G>& dense_ptr) {
    typename DenseMatrixCollection<G>::iterator i = dense_ptr.begin();
    for (; i != dense_ptr.end(); ++i)
      this->add(*i->second, i->first);
  } // storeData_
 
}
// namespace graphics
 
#endif // matlab_binary_hh