sphere.h
1.94 KB
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
#ifndef SPHERESTRUCT_H
#define SPHERESTRUCT_H
#include <ostream>
#include <sstream>
#include <vector>
#include <complex>
#include "fieldslice.h"
#include "dataTypes.h"
struct sphere
{
//sphere position
bsPoint p;
//sphere radius
ptype a;
//sphere material index
int iMaterial;
//GPU pointer to the scattered field produced by a plane wave
// this is a function of cos(theta) and |r| (distance from sphere center)
fieldslice Usp;
fieldslice Uip;
ptype d_min;
ptype d_max;
//sphere order
int Nl;
//refractive index for the current lambda
bsComplex n;
//external scattering coefficients
std::vector<bsComplex> B;
//internal scattering coefficients
std::vector<bsComplex> A;
sphere(ptype x = 0.0f, ptype y = 0.0f, ptype z = 0.0f, ptype a = 0.0f, int m = 0, int ang = 128)
{
this->p = bsPoint(x, y, z);
this->a = a;
this->iMaterial = m;
//surface = fieldslice(ang, ang/2);
}
//assignment operator
sphere & operator=(const sphere &rhs);
//copy constructor
sphere(const sphere &rhs);
std::string toStr()
{
std::stringstream ss;
ss<<p<<", "<<a<<", "<<iMaterial;
return ss.str();
}
//compute the order required to represent the scattered field
void calcNl(ptype lambda)
{
Nl = ceil( (2 * PI * a) / lambda + 4 * pow( (2 * PI * a) / lambda, 1.0/3.0) + 2);
}
//compute the scattering coefficients
void calcCoeff(ptype lambda, bsComplex n);
//compute the bessel function look-up tables
void calcLut(bsComplex* j, bsComplex* h, ptype lambda, bsComplex n, int aR, int rR);
void calcBesselLut(bsComplex* j, ptype k, bsComplex n, int aR);
void calcHankelLut(bsComplex* h, ptype k, int rR);
//calculate the scattering domain Us(theta, r)
void calcUp(ptype lambda, bsComplex n, rts::rtsQuad<ptype, 3> nfPlane, unsigned int R);
void scalarUsp(bsComplex* h, int rR, int thetaR);
void scalarUip(bsComplex* j, int aR, int thetaR);
};
#endif