main.cpp
1.15 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
#include "sphere.h"
#include "rts/optics/material.h"
#include <vector>
//#include "rts/complex.h"
#include "nearfield.h"
//nearfieldStruct* NF;
#include "microscope.h"
microscopeStruct* SCOPE;
#include "fieldslice.h"
#include "fileout.h"
#include "options.h"
#include "montecarlo.h"
#include "rts/math/point.h"
#include "rts/math/spherical_bessel.h"
#include "rts/math/matrix.h"
#include "rts/math/quaternion.h"
#include "rts/envi/envi.h"
#include "warnings.h"
fileoutStruct gFileOut;
using namespace std;
int cbessjyva(double v,complex<double> z,double &vm,complex<double>*cjv,
complex<double>*cyv,complex<double>*cjvp,complex<double>*cyvp);
int main(int argc, char *argv[])
{
//test Envi loading and saving
//EnviFile envi("testenvi", "w");
//float* data = (float*)malloc(sizeof(float) * 100 * 100);
//envi.addBand(data, 100, 100, 100);
//envi.close();
//return 0;
SCOPE = new microscopeStruct();
cout<<SCOPE->nf.Uf.R[0]<<endl;
LoadParameters(argc, argv);
//TestSimulation(NF, SCOPE, &gFileOut);
//initialize GPU memory for fields
SCOPE->init();
OutputOptions();
gFileOut.Save(SCOPE);
//NF->destroy();
SCOPE->destroy();
}