MeshReader.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stéphane Del Pino

//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2, or (at your option)
//  any later version.

//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.

//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  

//  $Id: MeshReader.hpp,v 1.18 2007/12/15 18:33:46 delpinux Exp $

#ifndef MESH_READER_HPP
#define MESH_READER_HPP

#include <MeshGenerator.hpp>

#include <ReferenceCounting.hpp>
#include <Vector.hpp>

#include <VerticesSet.hpp>

#include <Vertex.hpp>
#include <Hexahedron.hpp>
#include <Tetrahedron.hpp>
#include <Triangle.hpp>

#include <Stringify.hpp>
#include <ErrorHandler.hpp>

#include <string>
#include <cstdio>

class MeshReader
  : public MeshGenerator
{
protected:
  FILE * __ifh;                 
  std::string __fileName;       
  ReferenceCounting<VerticesSet> __vertices;

  ReferenceCounting<Vector<Hexahedron> > __hexahedra;

  ReferenceCounting<Vector<Quadrangle> > __quadrilaterals;

  ReferenceCounting<Vector<Tetrahedron> > __tetrahedra;

  ReferenceCounting<Vector<Triangle> > __triangles;

  MeshReader(const MeshReader& M);

  inline real_t __getReal()
  {
    // a double is used here instead of a real_t to allow scanf to
    // behave correctly
    register double r;
    int retval = fscanf(__ifh,"%le",&r);
    if (retval <= 0) {
      throw ErrorHandler(__FILE__,__LINE__,
                         "reading mesh '"+__fileName+"': expected real not found",
                         ErrorHandler::normal);
      return 0;
    }
    return r;
  }

  inline int __getInteger()
  {
    int i;
    int retval = fscanf(__ifh,"%d",&i);
    if (retval <= 0) {
      throw ErrorHandler(__FILE__,__LINE__,
                         "reading mesh '"+__fileName+"': expected integer not found",
                         ErrorHandler::normal);
      return 0;
    }
    return i;
  }

  virtual void __createMesh();

public:

  class Error
  { 
  private:
    const std::string __message;
  public:
    Error(const std::string & message) : __message(message) { }
    const std::string & message() const { return __message; }
  };

  ReferenceCounting<Mesh> mesh()
  {
    return __mesh;
  }

  MeshReader(const std::string & s);

  virtual ~MeshReader();
};

#endif // MESH_READER_HPP

---