meshImport2D.hh

Go to the documentation of this file.

 
#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);
 
    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