VBMicrolensing 5.1.1__cp310-cp310-macosx_11_0_arm64.whl → 5.3.1__cp310-cp310-macosx_11_0_arm64.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Potentially problematic release.
This version of VBMicrolensing might be problematic. Click here for more details.
- VBMicrolensing/VBMicrolensing.so +0 -0
- VBMicrolensing/__init__.py +1 -0
- VBMicrolensing/data/SunEphemeris.txt +22130 -0
- VBMicrolensing/lib/VBMicrolensingLibrary.cpp +594 -178
- VBMicrolensing/lib/VBMicrolensingLibrary.h +21 -7
- VBMicrolensing/lib/python_bindings.cpp +339 -2
- {vbmicrolensing-5.1.1.dist-info → vbmicrolensing-5.3.1.dist-info}/METADATA +4 -3
- vbmicrolensing-5.3.1.dist-info/RECORD +14 -0
- {vbmicrolensing-5.1.1.dist-info → vbmicrolensing-5.3.1.dist-info}/WHEEL +1 -1
- vbmicrolensing-5.1.1.dist-info/RECORD +0 -13
- {vbmicrolensing-5.1.1.dist-info → vbmicrolensing-5.3.1.dist-info}/licenses/LICENSE +0 -0
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
// VBMicrolensing v5.
|
|
1
|
+
// VBMicrolensing v5.3 (2025)
|
|
2
2
|
//
|
|
3
3
|
// This code has been developed by Valerio Bozza (University of Salerno) and collaborators.
|
|
4
4
|
// Check the repository at https://github.com/valboz/VBMicrolensing
|
|
@@ -75,6 +75,8 @@ class VBMicrolensing
|
|
|
75
75
|
|
|
76
76
|
int* ndatasat;
|
|
77
77
|
double** tsat, *** possat;
|
|
78
|
+
double** posEar, startEar, stepEar;
|
|
79
|
+
int ndataEar;
|
|
78
80
|
double Mag0;
|
|
79
81
|
double* dist_mp, * q;
|
|
80
82
|
int nim0, n, n2, nnm1, nroots, nrootsmp, * nrootsmp_mp;
|
|
@@ -95,9 +97,10 @@ class VBMicrolensing
|
|
|
95
97
|
int iastro;
|
|
96
98
|
double Obj[3], rad[3], tang[3], t0old;
|
|
97
99
|
double Eq2000[3], Quad2000[3], North2000[3];
|
|
98
|
-
double Et0[2], vt0[2], Et[2], Ehel[2];
|
|
100
|
+
double Et0[2], vt0[2], Et[2], Ehel[2], lighttravel, lighttravel0;
|
|
99
101
|
double ESPLout[__rsize_ESPL][__zsize_ESPL], ESPLin[__rsize_ESPL][__zsize_ESPL], ESPLoutastro[__rsize_ESPL][__zsize_ESPL], ESPLinastro[__rsize_ESPL][__zsize_ESPL];
|
|
100
|
-
bool
|
|
102
|
+
bool coordinates_set;
|
|
103
|
+
bool multidark;
|
|
101
104
|
double* LDtab, * rCLDtab, * CLDtab;
|
|
102
105
|
double scr2, sscr2;
|
|
103
106
|
int npLD;
|
|
@@ -146,12 +149,18 @@ public:
|
|
|
146
149
|
bool astrometry;
|
|
147
150
|
bool turn_off_secondary_source;
|
|
148
151
|
bool turn_off_secondary_lens;
|
|
152
|
+
bool ESPLoff;
|
|
153
|
+
bool t_in_HJD;
|
|
149
154
|
|
|
150
155
|
static char ESPLtablefile[1024];
|
|
151
156
|
static void SetESPLtablefile(char* instring) { strcpy(ESPLtablefile, instring); }
|
|
152
|
-
|
|
157
|
+
static char Suntablefile[1024];
|
|
158
|
+
static void SetSuntablefile(char* instring) { strcpy(Suntablefile, instring); }
|
|
159
|
+
double Tol, RelTol, a1, a2,corrquad, corrquad2, safedist;
|
|
153
160
|
double mass_radius_exponent, mass_luminosity_exponent, lens_mass_luminosity_exponent;
|
|
154
|
-
int satellite, parallaxsystem, t0_par_fixed, nsat;
|
|
161
|
+
int satellite, parallaxsystem, t0_par, t0_par_fixed, nsat;
|
|
162
|
+
bool suntable, parallaxephemeris;
|
|
163
|
+
int parallaxextrapolation;
|
|
155
164
|
int minannuli, maxannuli, nannuli, NPS, NPcrit;
|
|
156
165
|
int newtonstep;
|
|
157
166
|
double y_1, y_2, av, therr, astrox1, astrox2;
|
|
@@ -163,6 +172,7 @@ public:
|
|
|
163
172
|
// Initialization for parallax calculation
|
|
164
173
|
void SetObjectCoordinates(char* Coordinates_file, char* Directory_for_satellite_tables);
|
|
165
174
|
void SetObjectCoordinates(char* CoordinateString);
|
|
175
|
+
bool AreCoordinatesSet();
|
|
166
176
|
// Skowron & Gould root calculation
|
|
167
177
|
void cmplx_roots_gen(complex*, complex*, int, bool, bool);
|
|
168
178
|
void cmplx_roots_multigen(complex*, complex**, int, bool, bool);
|
|
@@ -202,7 +212,7 @@ public:
|
|
|
202
212
|
void SetMethod(Method);
|
|
203
213
|
|
|
204
214
|
//ESPL functions
|
|
205
|
-
void LoadESPLTable(char* tablefilename);
|
|
215
|
+
void LoadESPLTable(const char* tablefilename);
|
|
206
216
|
double ESPLMag(double u, double rho);
|
|
207
217
|
double ESPLMag2(double u, double rho);
|
|
208
218
|
double ESPLMagDark(double u, double rho);
|
|
@@ -211,6 +221,7 @@ public:
|
|
|
211
221
|
|
|
212
222
|
// New (v2) light curve functions, operating on arrays
|
|
213
223
|
|
|
224
|
+
void LoadSunTable(char* tablefilename);
|
|
214
225
|
void PSPLLightCurve(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
215
226
|
void PSPLLightCurveParallax(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
216
227
|
void ESPLLightCurve(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
@@ -229,6 +240,7 @@ public:
|
|
|
229
240
|
void BinSourceExtLightCurveXallarap(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, double* y1_array2, double* y2_array2, int np);
|
|
230
241
|
void BinSourceSingleLensXallarap(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, double* y1_array2, double* y2_array2, int np);
|
|
231
242
|
void BinSourceBinLensXallarap(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
243
|
+
void BinSourceBinLensLightCurve(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, double* y1_array2, double* y2_array2, double* seps_array, int np);
|
|
232
244
|
|
|
233
245
|
void TripleLightCurve(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
234
246
|
void TripleLightCurveParallax(double* parameters, double* t_array, double* mag_array, double* y1_array, double* y2_array, int np);
|
|
@@ -253,7 +265,8 @@ public:
|
|
|
253
265
|
double BinSourceExtLightCurveXallarap(double* parameters, double t);
|
|
254
266
|
double BinSourceBinLensXallarap(double* parameters, double t);
|
|
255
267
|
double BinSourceSingleLensXallarap(double* parameters, double t);
|
|
256
|
-
double
|
|
268
|
+
double BinSourceBinLensLightCurve(double* parameters, double t);
|
|
269
|
+
|
|
257
270
|
|
|
258
271
|
double TripleLightCurve(double* parameters, double t);
|
|
259
272
|
double TripleLightCurveParallax(double* parameters, double t);
|
|
@@ -266,6 +279,7 @@ public:
|
|
|
266
279
|
void BinaryAstroLightCurveOrbital(double* parameters, double* t_array, double* mag_array, double* centroid_s1_array, double* centroid_s2_array, double* centroid_l1_array, double* centroid_l2_array, double* y1_array, double* y2_array, double* seps_array, int np);
|
|
267
280
|
void BinaryAstroLightCurveKepler(double* parameters, double* t_array, double* mag_array, double* centroid_s1_array, double* centroid_s2_array, double* centroid_l1_array, double* centroid_l2_array, double* y1_array, double* y2_array, double* seps_array, int np);
|
|
268
281
|
void BinSourceAstroLightCurveXallarap(double* parameters, double* t_array, double* mag_array, double* centroid_s1_array, double* centroid_s2_array, double* centroid_l1_array, double* centroid_l2_array, double* y1_array, double* y2_array, double* y1_array2, double* y2_array2, int np);
|
|
282
|
+
void BinSourceBinLensAstroLightCurve(double* parameters, double* t_array, double* mag_array, double* centroid_s1_array, double* centroid_s2_array, double* centroid_l1_array, double* centroid_l2_array, double* y1_array, double* y2_array, double* y1_array2, double* y2_array2, double* seps_array, int np);
|
|
269
283
|
void TripleAstroLightCurve(double* parameters, double* t_array, double* mag_array, double* centroid_s1_array, double* centroid_s2_array, double* centroid_l1_array, double* centroid_l2_array, double* y1_array, double* y2_array, int np);
|
|
270
284
|
|
|
271
285
|
// Constructor and destructor
|