SpectralMeshTool.hpp

Go to the documentation of this file.

#ifndef SPECTRALMESHTOOL_HPP
#define SPECTRALMESHTOOL_HPP

#include "moab/Interface.hpp" // needs to be here to support inline query_interface
#include "moab/Error.hpp" // needs to be here to support inline query_interface
#include <vector>

namespace moab {

class SpectralMeshTool
{
public:

  SpectralMeshTool(Interface *impl, int order = 0);

  ~SpectralMeshTool();

  Tag spectral_vertices_tag(const bool create_if_missing = false);
  
  Tag spectral_order_tag(const bool create_if_missing = false);
  
  ErrorCode convert_to_fine(EntityHandle spectral_set);
  
  ErrorCode convert_to_coarse(int order = 0, EntityHandle spectral_set = 0);

  template <class T>
  ErrorCode create_spectral_elems(const T *conn, int num_fine_elems, int dim, 
                                  Range &output_range, int start_idx = 0, Range *local_gids = NULL);
  
  void spectral_order(int order) {spectralOrder = order; spectralOrderp1 = order+1;}
  
  int spectral_order() {return spectralOrder;}
/*
  struct ConnMap 
  {
    const short a[16];
  };
  */
  static const short int permute_array[];
  
  static const short int lin_permute_array[];
  
private:

  Interface* mbImpl;

  Error *mError;
  
  Tag svTag;
  
  Tag soTag;

  int spectralOrder;

  int spectralOrderp1;
};

inline SpectralMeshTool::SpectralMeshTool(Interface *impl, int order) 
        : mbImpl(impl), svTag(0), soTag(0), spectralOrder(order), spectralOrderp1(order+1)
{
  impl->query_interface(mError);
}
    
inline SpectralMeshTool::~SpectralMeshTool() 
{}

} // namespace moab 

#endif