geometry/meshImport2D.hh

Go to the documentation of this file.

/* Meshes
 * Combining cells to a mesh to represent a computational domain
 * Importing 2D meshes from text files
 */

#ifndef meshImport2D_hh
#define meshImport2D_hh

#include <vector>
#include <map>
#include "mesh.hh"
#include "meshImport.hh"
#include "toolbox/array.hh"
#include "toolbox/multiArray.hh"

namespace concepts {

  // ****************************************************** Import2dMeshBase **

  class Import2dMeshBase : public Mesh2, public ImportMesh {
  public:
    Import2dMeshBase(const std::string coord, const std::string elms,
         const uint idxStart = 1);
    virtual ~Import2dMeshBase();
    virtual uint ncell() const { return cell_.size(); }
    virtual Scan2* scan();
  protected:
    virtual std::ostream& info(std::ostream& os) const;

    class S;

    std::vector<Cell2*> cell_;

    const std::string coord_, elms_;

    std::vector<Real3d> vertices_;
    uint dim_;

    std::vector<Connector2*> cntr2_;
    MultiArray<3, Triangle*> Tri_;
    MultiArray<4, Quad*> Quad_;

    virtual void createEntity_(const MultiIndex<1>& idx) 
    { ImportMesh::createEntity_(idx); }
    virtual void createEntity_(const MultiIndex<2>& idx) 
    { ImportMesh::createEntity_(idx); }
    virtual void createEntity_(const MultiIndex<3>& idx) 
      throw(concepts::MissingFeature);
    virtual void createEntity_(const MultiIndex<4>& idx)
      throw(concepts::MissingFeature);

    virtual void createCell_(const MultiIndex<3>& idx)
      throw(concepts::MissingFeature);
    virtual void createCell_(const MultiIndex<4>& idx)
      throw(concepts::MissingFeature);

    virtual void readVertices_();
    uint readVertexLine_(const std::string& i, uint* idx, Real3d* coord) const;

    void attributes_(const std::string boundary);
    void import_();
  };

  class Import2dMeshBase::S : public Scan<Cell2> {
  public:
    inline S(std::vector<Cell2*>& cell) :
      idx_(cell.begin()), cell_(cell) {}
    inline S(const S& scan) : idx_(scan.idx_), cell_(scan.cell_) {}
    inline bool eos() const { return idx_ == cell_.end(); }
    inline Cell2& operator++(int) { return *(*idx_++); }
    inline Scan2* clone() const { return new S(*this); }
  private:
    std::vector<Cell2*>::iterator idx_;
    std::vector<Cell2*>& cell_;
  };

  // ********************************************************** Import2dMesh **

  class Import2dMesh : public Import2dMeshBase {
  public:
    Import2dMesh(const std::string coord, const std::string elms,
     const std::string boundary, const uint idxStart = 1);

//my_new:
    Import2dMesh(const std::string coord, const std::string elms,
     const std::string boundary, const std::string edgecorr,
     const uint idxStart = 1);
    Import2dMesh(const std::string coord, const std::string elms,
     const uint idxStart = 1);
    virtual ~Import2dMesh() {}
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };

} // namespace concepts

#endif // meshImport2D_hh