includeLinks/MeshSet.hpp

Go to the documentation of this file.

/* ***************************************************************** 
    MESQUITE -- The Mesh Quality Improvement Toolkit

    Copyright 2004 Sandia Corporation and Argonne National
    Laboratory.  Under the terms of Contract DE-AC04-94AL85000 
    with Sandia Corporation, the U.S. Government retains certain 
    rights in this software.

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

    This library 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
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public License 
    (lgpl.txt) along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
    diachin2@llnl.gov, djmelan@sandia.gov, mbrewer@sandia.gov, 
    pknupp@sandia.gov, tleurent@mcs.anl.gov, tmunson@mcs.anl.gov      
   
  ***************************************************************** */
#ifndef Mesquite_MeshSet_hpp 
#define Mesquite_MeshSet_hpp

#ifdef MSQ_USE_OLD_C_HEADERS
#  include <stddef.h>
#else
#  include <cstddef>
#endif

#ifdef MSQ_USE_OLD_STD_HEADERS
#include <list.h>
#else
#include <list>
#endif

#include "Mesquite.hpp"
#include "PatchDataUser.hpp"

namespace Mesquite
{
  class MsqError;
  class Mesh;
  class MeshDomain;
  
  class MeshSet
  {
  public:
    
    MeshSet();
    ~MeshSet();
    

    void add_mesh(Mesquite::Mesh* mesh, MsqError &err);


    void set_domain_constraint(MeshDomain* domain, MsqError &err);

    Mesquite::MeshDomain* get_domain_constraint()
      { return mDomain; }
    
    void get_meshes(msq_std::list<Mesquite::Mesh*> &mesh_list) const
      { mesh_list = meshSet; }

      // We want to rename this to geometric_dimension() const
    int space_dim() const
      { return spaceDim; }
    

    bool get_next_patch(PatchData &pd,
                        PatchDataUser* pd_user,
                        MsqError &err)
      {
        return get_next_patch(pd, pd_user->get_all_parameters(), err);
      }
    
    bool get_next_patch(PatchData &pd,
                        PatchDataParameters &pd_params,
                        MsqError &err);
    
    void reset(MsqError &err);

    void update_mesh(const PatchData &pd, MsqError &err);

    void set_flag_to_cull(MsqVertex::FlagMask f_m)
      {
        cullFlag=f_m;
      }
    
    bool clear_all_soft_fixed_flags(MsqError &err);

    void write_vtk(const char* out_filebase, MsqError &err);
    
    void write_gnuplot(const char* out_filebase, MsqError &err);
    
    Mesh* get_current_mesh()
      { return meshSet.empty() ? 0 : *currentMesh; }
    
  private:
  
    bool get_next_elem_on_vert_patch( PatchData&, PatchDataParameters&, MsqError& );
    bool get_next_global_patch( PatchData&, PatchDataParameters&, MsqError& );
    
    msq_std::list<Mesquite::Mesh*> meshSet;
    msq_std::list<Mesquite::Mesh*>::iterator currentMesh;
    Mesquite::VertexIterator *vertexIterator;
    int spaceDim;
    
    Mesquite::EntityTopology elementType;

    msq_stdc::size_t *csrOffsets;
    //msq_stdc::size_t *csrData;
    //Mesh::VertexHandle *vertArray;
    //Mesh::ElementHandle *elemArray;
    Mesquite::EntityTopology *elemTopologies;
    bool *vertexOnBoundary;
    msq_stdc::size_t csrOffsetsSize;
    //msq_stdc::size_t csrDataSize;
    msq_stdc::size_t vertArraySize;
    msq_stdc::size_t elemArraySize;

      // This is the domain we snap everything back to.
    Mesquite::MeshDomain *mDomain;

      //Flag to tell the MeshSet which MSQ culling flag to use
      // to determine whether to build a local patch around a vertex
    Mesquite::MsqVertex::FlagMask cullFlag;

  };
  
    // -********** AOMD tmp TEST **********
  void test_aomd(void);
  
} //namespace

#endif