AraSim_old/RaySolver.h at master · ara-software/AraSim_old · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#ifndef RAYSOLVER_H
#define RAYSOLVER_H


#include <vector>
#include <iostream>
#include "RayTrace.h"
#include "RayTrace_IceModels.h"
#include "Vector.h"
#include <fstream>
#include <iomanip>
#include "Settings.h"
#ifndef __CINT__
#include <boost/lexical_cast.hpp>
#include <boost/program_options.hpp>
#endif

class Position;
class IceModel;

//--------------------------------------------------
// struct underline{
// 	static char esc;
// 	const std::string& str;
// 	underline(const std::string& s):str(s){}
// 	friend std::ostream& operator<<(std::ostream& os, const underline& u);
// };
//
// char underline::esc=0x1B;
//
// std::ostream& operator<<(std::ostream& os, const underline& u){
// 	return(os << underline::esc << "[4m" << u.str << underline::esc << "[0m");
// }
//--------------------------------------------------

template<typename positionType>

class pathPrinter{
public:
	pathPrinter(){}
	void operator()(const positionType& p, RayTrace::RKStepType stepType){
		std::cout << p.x << ' ' << p.z << '\n';
	}
};

class RaySolver {

    private:
        void Earth_to_Flat_same_depth(Position &source, Position &target, IceModel* antarctica);
        void Earth_to_Flat_same_angle(Position &source, Position &target, IceModel* antarctica);

    public:
        RaySolver();
//        RaySolver(int argc, char* argv[]);

        void test1();
        void Solve_Ray_org(Position &source, Position &target, std::vector < std::vector <double> > &outputs, Settings *settings1);

        //void Solve_Ray_org (double source_x, double source_y, double source_z, double target_x, double target_y, double target_z, double &travel_time, double &travel_dist, int &no_sol, Settings *settings1);
        void Solve_Ray_org (double source_x, double source_y, double source_z, double target_x, double target_y, double target_z, bool print_path=false );

        void Solve_Ray(Position &source, Position &target, IceModel *antarctica, std::vector < std::vector <double> > &outputs, Settings *settings1);

        int source_over_surface;
        int solution_toggle;    // no solution : 0  solution exist : 1
};


#endif //RAYSOLVER_H