mcnpmit.hpp

Go to the documentation of this file.

#include "moab/Core.hpp"
#include "moab/Range.hpp"
#include <iostream>
#define MCNP mc_instance()
#define BOXMIN_TAG "BOXMIN_TAG"
#define BOXMAX_TAG "BOXMAX_TAG"
#define TALLY_TAG "TALLY_TAG"
#define ERROR_TAG "ERROR_TAG"
 
#define MBI mb_instance()
 
enum MCNPError
{
 MCNP_SUCCESS,
 MCNP_FAILURE,
 DONE
};
enum
{
 NOSYS,
 CARTESIAN,
 CYLINDRICAL,
 SPHERICAL
};
 
class McnpData
{
 
 public:
 // Constructor and Destructor
 McnpData();
 ~McnpData();
 
 // Coordinate system and rotation matrix
 int coord_system;
 double rotation_matrix[16];
 
 // Vertices and elements
 std::vector< moab::EntityHandle > MCNP_vertices;
 std::vector< moab::EntityHandle > MCNP_elems;
 moab::Range vert_handles;
 moab::Range elem_handles;
 
 // Tally data
 moab::Tag box_min_tag, box_max_tag;
 moab::Tag tally_tag;
 moab::Tag relerr_tag;
 
 // MCNP Meshtal file name
 std::string MCNP_filename;
 
 // Setting and retrieving coordinate sysem
 MCNPError set_coord_system( int );
 int get_coord_system();
 
 // Setting and retrieving roation matrix
 MCNPError set_rotation_matrix( double[16] );
 double* get_rotation_matrix();
 
 // Set the filename
 MCNPError set_filename( std::string );
 std::string get_filename();
 
 // MCNP reading routines
 MCNPError read_mcnpfile( bool );
 MCNPError read_coord_system( std::string );
 MCNPError read_rotation_matrix( std::string, int );
 MCNPError make_elements( std::vector< double >[3], int* );
 MCNPError make_adjacencies( int* );
 MCNPError initialize_tags();
 MCNPError extract_tally_data( std::string, moab::EntityHandle );
 
 // Transformation routine
 MCNPError transform_point( double*, double*, int, double* );
};