skyviewer/boundary.cpp at main · nasa-lambda/skyviewer · 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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
/* ============================================================================
'boundary.cpp' defines a boundary class.

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
============================================================================ */
/*
			Fetch header files.
*/
#include <math.h>
#include "boundary.h"
#include "map_exception.h"

using namespace std;
/* ----------------------------------------------------------------------------
'Boundary' is the class constructor.

Arguments:
	None.

Returned:
	Nothing.

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
---------------------------------------------------------------------------- */
Boundary::Boundary(void)
{
	a = b = 0.;
	face = 0;
 	eq = true;
	return;
}
/* ----------------------------------------------------------------------------
'set_eq' defines an equatorial boundary.

Arguments:
	z         -
	phi       -
	direction -

Returned:
	Nothing.

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
---------------------------------------------------------------------------- */
void Boundary::set_eq(const double z, const double phi, const int direction)
{
	b = direction*0.375*M_PI;
	a = phi - b*z;
 	eq = true;
	return;
}
/* ----------------------------------------------------------------------------
'set_np' defines a north pole boundary.

Arguments:
	z     -
	phi   -
	face_ -

Returned:
	Nothing.

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
---------------------------------------------------------------------------- */
void Boundary::set_np(const double z, const double phi, const int face_)
{
 	eq = false;
	face = face_;
	a = sqrt(1-z)*(0.5*M_PI*(face+1) - phi);
	return;
}
/* ----------------------------------------------------------------------------
'set_sp' defines a south pole boundary.

Arguments:
	z     -
	phi   -
	face_ -

Returned:
	Nothing.

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
---------------------------------------------------------------------------- */
void Boundary::set_sp(const double z, const double phi, const int face_)
{
 	eq = false;
	face = face_;
	a = sqrt(1+z)*(phi - 0.5*M_PI*(face-8));
	return;
}
/* ----------------------------------------------------------------------------
'operator()' computes phi for a given z.

Arguments:
	z   -

Returned:
	phi -

Written by Nicholas Phillips.
QT4 adaption by Michael R. Greason, ADNET, 27 August 2007
---------------------------------------------------------------------------- */
double Boundary::operator()(const double z) const
{
	double phi;
	if( eq )
		phi = a + b*z;
	else if( face < 4) {
		phi = 0.5*M_PI*(face+1);
		if( z < 1 )
			phi += -a/sqrt(1-z);
	}
	else if( face > 7) {
		phi = 0.5*M_PI*(face-8);
		if( z > -1 )
			phi += a/sqrt(1+z);
	}
	else {
		phi = 0;
		throw MapException(MapException::Other, 0,
			"Boundary::operator():Bad call");
	}
	return phi;
}