ceo 0.1.0

CUDA Engined Optics
Documentation 
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
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224

#ifndef __GMTMIRRORS_H_
#define __GMTMIRRORS_H_

#ifndef __RAYTRACING_H__
#include "rayTracing.h"
#endif

#ifndef __SOURCE_H__
#include "source.h"
#endif

__device__ rtd bending_modes_surface(vector *v, const int N_mode, rtd *b,
                                     double *BM,
                                     const int BM_N_SAMPLE, const double BM_radius);
__device__ rtd partial_x_bending_modes_surface(vector *v, const int N_mode, rtd *b,
                                               double *BM,
                                               const int BM_N_SAMPLE, const double BM_radius);
__device__ rtd partial_y_bending_modes_surface(vector *v, const int N_mode, rtd *b,
                                               double *BM,
                                               const int BM_N_SAMPLE, const double BM_radius);
__device__ rtd partial_z_bending_modes_surface(void);



struct bending_modes {
  double *d__x_BM, *d__y_BM;
  double *d__BM;
  double *d__BMS;
  double BM_radius;
  int BM_N_SAMPLE;
  double *d__BM_buffer;
  int n_mode;
  rtd *b, *d__b;
  int N;  
  void setup(int n_mode);
  void setup(int n_mode, int N_SET);
  void setupKL(int n_mode);
  void setupKL(int n_mode, int N_SET);
  void setupPolish(int n_mode, int N_SET);
  void cleanup(void);
  void load(void);
  void load_reg(void);
  void load_KL(void);
  void load_polish(void);
  void nearest_neighbor(rtd *d__BMi,
                   rtd *d__partial_x_BMi,
                   rtd *d__partial_y_BMi,
                   int NI, rtd di, int k_mode);
  void bilinear(rtd *d__BMi,
           rtd *d__partial_x_BMi,
           rtd *d__partial_y_BMi,
           int NI, rtd di, int k_mode);
  void update(rtd *b);
};
struct modes {
  double *d__x_BM, *d__y_BM;
  double *d__BM;
  double *d__BMS;
  double BM_radius;
  int BM_N_SAMPLE;
  double *d__BM_buffer;
  int n_mode;
  rtd *b, *d__b;
  int N;  
  char filename[256];
  int N_SET;
  int N_MODE;
  int *d__s2b;  
  void setup(char *filename, int N, int n_mode);
  void setup(int _BM_N_SAMPLE_, double _BM_radius_,
             int _N_SET_, int _N_MODE_, int *s2b,
             double *buffer, int _N_, int _n_mode_);
  void cleanup(void);
  void load(void);
  void load(int _BM_N_SAMPLE_, double _BM_radius_,
            int _N_SET_, int _N_MODE_, int *s2b,
            double *buffer);
  void reset_modes(double *buffer);
  void update(rtd *b);
};
struct gmt_m1 {
  int M_ID;
  rtd D_assembly;
  rtd D_clear;
  rtd D_full;
  rtd ri;
  rtd beta;
  rtd L;
    //rtd D;
    //int D_px;
  rtd area0, area_fraction;
  float area0_px;
  rtd area;
  int N;
  rtd depth;
  coordinate_system aperture_CS;
  coordinate_system conic_CS;
  vector conic_origin[7], *d__conic_origin;
  rtd conic_c, conic_k;
  rtd *d__conic_c, *d__conic_k;
  coordinate_system rigid_body_CS;
  coordinate_system motion_CS;
  rtd height;
  mask *V;
  int idx_offset;
  zernikeS *ZS;
  int *d__piston_mask;
  coordinate_system TT_CS;
  double *d__C, *d__D;
  cublasHandle_t handle;
  char *d__valid_segments;
    //bending_modes *BS;
  modes BS;
  float *d__segment_reflectivity;  
  
  void setup(void);
  //void setup(modes *_BS_);
  void setup(char *_filename_, int _N_, int _n_mode_);
  void setup(zernikeS *ZS);
  void cleanup(void);
  void update(vector _origin_, vector _euler_angles_,int idx);
  void reset(void);
  void trace(bundle *rays);
  void traceall(bundle *rays);
  void blocking(bundle *rays);
  void global_tiptilt(float tip, float tilt);
  void test_ray_tracing(void);
  void track(float *x, float *y, float *z, int N, int idx);
  void locate(float *x, float *y, float *z, int N, int idx);
  void remove(int *seg_ID, int N_ID);
  void keep(int *seg_ID, int N_ID);
  void update_conic_c(rtd *_conic_c_);
  void update_conic_k(rtd *_conic_k_);
  void set_reflectivity(float *reflectivity);
  void preset(bundle *rays, rtd margin);
  void edge_sensors_data(void);
};
struct gmt_m2 {
  int M_ID;
  rtd D_assembly;
  rtd D_clear;
  rtd D_full;
  rtd ri;
  rtd beta;
  rtd L;
    //rtd D;
    //int D_px;
  rtd area0, area_fraction;
  float area0_px;
  rtd area;
  int N;
  rtd depth;
  coordinate_system aperture_CS;
  coordinate_system conic_CS;
  vector conic_origin[7], *d__conic_origin;
  rtd conic_c, conic_k;
  rtd *d__conic_c, *d__conic_k;
  coordinate_system rigid_body_CS;
  coordinate_system motion_CS;
  rtd height;
  mask *V;
  int idx_offset;
  zernikeS *ZS;
  int *d__piston_mask;
  coordinate_system TT_CS;
  double *d__C, *d__D;
  cublasHandle_t handle;
  char *d__valid_segments;
    //bending_modes *BS;
  modes BS;
  float *d__segment_reflectivity;  
  void setup(void);
  //void setup(modes *_BS_);
  void setup(char *_filename_, int _N_, int _n_mode_);
  void setup(zernikeS *ZS);
  void cleanup(void);
  void update(vector _origin_, vector _euler_angles_,int idx);
  void reset(void);
  void trace(bundle *rays);
  void traceall(bundle *rays);
  void blocking(bundle *rays);
  void global_tiptilt(float tip, float tilt);
  void test_ray_tracing(void);
  void track(float *x, float *y, float *z, int N, int idx);
  void locate(float *x, float *y, float *z, int N, int idx);
  void remove(int *seg_ID, int N_ID);
  void keep(int *seg_ID, int N_ID);
  void update_conic_c(rtd *_conic_c_);
  void update_conic_k(rtd *_conic_k_);
  void set_reflectivity(float *reflectivity);
  void pointing_neutral(float tip, float tilt);
  void coma_neutral(float tip, float tilt);
};
struct stereoscopic_edge_sensors {
int N, N_DATA;
int n_byte, n_byte_DATA;
rtd alpha;
gmt_m1 *mirror;
vector *v0, *v, *d__v0, *d__v;
vector *dv0, *d__dv0;
vector *dv, *d__dv;
void setup(gmt_m1 *_mirror_);
void cleanup(void);
void data(void);
};
struct lateral_edge_sensors {
int N, N_DATA, N_HEIGHT;
rtd alpha;
rtd height;
gmt_m1 *mirror;
int n_byte_DATA;
rtd *x, *y, *d__x, *d__y, *d, *d__d;
vector *k_cam, *d__k_cam;
vector *k_laser, *d__k_laser;
vector *uP, *d__uP;
vector *vP, *d__vP;
vector *A, *A0, *B, *B0, *d__A, *d__A0, *d__B, *d__B0;
void setup(gmt_m1 *_mirror_);
void setup(gmt_m1 *_mirror_, rtd _height_);
void cleanup(void);
void data(void);
};

#endif // __GMTMIRRORS_H_