external_packages/hydro_from_external_file/SurfaceFinder.h

0001 // Copyright Chun Shen @ 2015
0002 #ifndef SRC_SurfaceFinder_H_
0003 #define SRC_SurfaceFinder_H_
0004
0005 #ifdef USE_HDF5
0006 #include "./Hydroinfo_h5.h"
0007 #endif
0008
0009 #include "./Hydroinfo_MUSIC.h"
0010 #include "./ParameterReader.h"
0011
0012 #include <fstream>
0013
0014 // using namespace std;
0015 using std::string;
0016 using std::vector;
0017 using std::ofstream;
0018 using std::ifstream;
0019 using std::abs;
0020
0021 class SurfaceFinder {
0022  private:
0023     int hydro_type;
0024 #ifdef USE_HDF5
0025     HydroinfoH5 *hydroinfo_ptr;
0026 #endif
0027     Hydroinfo_MUSIC *hydroinfo_MUSIC_ptr;
0028     ParameterReader *paraRdr;
0029     double T_cut;
0030
0031  public:
0032     SurfaceFinder(void* hydroinfo_ptr_in, ParameterReader* paraRdr_in);
0033     SurfaceFinder(void* hydroinfo_ptr_in, ParameterReader* paraRdr_in,
0034                   double T_cut_in);
0035     ~SurfaceFinder();
0036
0037     bool check_intersect(double T_cut, double tau, double x, double y,
0038                          double dt, double dx, double dy, double ***cube);
0039     int Find_full_hypersurface();
0040 };
0041
0042 #endif  // SRC_SurfaceFinder_H_