36 #include <sys/types.h>
42 #define DLL_EXPORT extern
353 ({ __typeof (a) _a = (a); \
354 __typeof (b) _b = (b); \
355 _a < _b ? _a : _b; })
360 ({ __typeof (a) _a = (a); \
361 __typeof (b) _b = (b); \
362 _a > _b ? _a : _b; })
372 #define hz2rad(hz) (2.0*PI*hz)
377 #define sinpsin(s1, s2) (2.0*sin((asin(s1)+asin(s2))/2.0)*cos((asin(s1)-asin(s2))/2.0))
382 #define sinmsin(s1, s2) (2.0*cos((asin(s1)+asin(s2))/2.0)*sin((asin(s1)-asin(s2))/2.0))
387 #define cospcos(c1, c2) (2.0*cos((acos(s1)+acos(s2))/2.0)*cos((acos(s1)-acos(s2))/2.0))
392 #define cosmcos(c1, c2) (2.0*sin((acos(s1)+acos(s2))/2.0)*sin((acos(s1)-acos(s2))/2.0))
397 #define sinxsin(s1, s2) ((cos(asin(s1)-asin(s2))-cos(asin(s1)+asin(s2)))/2.0)
402 #define cosxcos(c1, c2) ((cos(acos(s1)+acos(s2))+cos(acos(s1)-acos(s2)))/2.0)
407 #define sinxcos(s, c) ((sin(asin(s)+acos(c))+sin(asin(s)-acos(c)))/2.0)
412 #define sin2cos(s) cos(asin(s))
417 #define cos2sin(c) sin(acos(c))
420 #ifndef VLBI_VERSION_STRING
422 #define VLBI_VERSION_STRING "1.23.0"
425 #ifndef VLBI_CATALOG_PATH
427 #define VLBI_CATALOG_PATH "/usr/share/OpenVLBI/cat/index.txt"
432 #define CIRCLE_DEG 360.0
436 #define CIRCLE_AM (CIRCLE_DEG * 60.0)
440 #define CIRCLE_AS (CIRCLE_AM * 60.0)
444 #define RAD_AS (CIRCLE_AS/(PI*2.0))
446 #ifndef ONE_SECOND_TICKS
448 #define ONE_SECOND_TICKS 100000000
450 #ifndef ONE_MILLISECOND_TICKS
452 #define ONE_MILLISECOND_TICKS 100000
454 #ifndef ONE_MICROSECOND_TICKS
456 #define ONE_MICROSECOND_TICKS 100
460 #define SOLAR_DAY 86400
464 #define SIDEREAL_DAY 86164.0905
466 #ifndef TRACKRATE_SIDEREAL
468 #define TRACKRATE_SIDEREAL (CIRCLE_AS / SIDEREAL_DAY)
470 #ifndef TRACKRATE_SOLAR
472 #define TRACKRATE_SOLAR (CIRCLE_AS / SOLAR_DAY)
474 #ifndef TRACKRATE_LUNAR
476 #define TRACKRATE_LUNAR 14.511415
478 #ifndef EARTHRADIUSEQUATORIAL
480 #define EARTHRADIUSEQUATORIAL 6378137.0
482 #ifndef EARTHRADIUSPOLAR
484 #define EARTHRADIUSPOLAR 6356752.0
486 #ifndef EARTHRADIUSMEAN
488 #define EARTHRADIUSMEAN 6372797.0
492 #define AVOGADRO 6.02214076E+23
496 #define EULER 2.71828182845904523536028747135266249775724709369995
500 #define PLANK 6.62607015E-34
504 #define BOLTSMANN 1.380649E-23
506 #ifndef STEPHAN_BOLTSMANN
508 #define STEPHAN_BOLTSMANN (2.0*pow(PI, 5)*pow(BOLTSMANN, 2)/(pow(LIGHTSPEED, 2)*15*pow(PLANK, 3)))
512 #define GAS_R (BOLTSMANN * AVOGADRO)
516 #define ROOT2 1.41421356237309504880168872420969807856967187537694
520 #define PI 3.14159265358979323846
528 #define LIGHTSPEED 299792458.0
532 #define J2000 2451545.0
536 #define GAMMAJ2000 18.6971378528
540 #define ELECTRON 1.602176634E-19
546 #ifndef ASTRONOMICALUNIT
548 #define ASTRONOMICALUNIT 1.495978707E+11
552 #define PARSEC (ASTRONOMICALUNIT/sin(PI*2.0/CIRCLE_AS))
556 #define LY (LIGHTSPEED * SIDEREAL_DAY * 365.0)
564 #define AU2M(au) (au * ASTRONOMICALUNIT)
572 #define PARSEC2M(parsec) (parsec * PARSEC)
580 #define LY2M(ly) (ly * LY)
588 #define M2AU(m) (m / ASTRONOMICALUNIT)
596 #define M2PARSEC(m) (m / PARSEC)
604 #define M2LY(m) (m / LY)
612 #define RAD2AS(rad) (rad * RAD_AS)
620 #define AS2RAD(as) (as / RAD_AS);
643 return sqrt(pow(x, 2)+pow(y, 2));
655 double mag = sqrt(pow(x, 2)+pow(y, 2));
658 rad = acos (y / (mag > 0.0 ? mag : 1.0));
659 if(x < 0 && rad != 0)
851 DLL_EXPORT
void vlbi_get_offsets(
vlbi_context ctx,
double J2000Time,
const char* node1,
const char* node2,
double Ra,
double Dec,
double *offset1,
double *offset2);
873 DLL_EXPORT
void vlbi_get_uv_plot(
void *ctx,
const char *name,
int u,
int v,
double *target,
double freq,
double sr,
int nodelay,
int moving_baseline,
vlbi_func2_t delegate,
int *interrupt);
DLL_EXPORT double vlbi_astro_estimate_temperature(double wavelength, double flux)
Estimate the physical temperature with a given flux.
DLL_EXPORT dsp_stream_p vlbi_astro_create_reference_catalog(dsp_stream_p *catalog, int catalog_size)
Create a dsp_stream_p containing all the spectral lines of all elements of a catalog.
DLL_EXPORT double vlbi_astro_diff_spectra(dsp_stream_p spectrum0, dsp_stream_p spectrum, double decades)
Compare a spectrum to a reference one.
DLL_EXPORT double vlbi_astro_estimate_distance_ratio(double luminosity_ratio, double flux_ratio)
Returns the distance ratio of two celestial object.
DLL_EXPORT dsp_align_info vlbi_astro_align_spectra(dsp_stream_p spectrum, dsp_stream_p catalog, int max_lines, double decimals, double min_score)
Align a spectrum to a reference catalog.
DLL_EXPORT double vlbi_astro_estimate_flux(double wavelength, double temperature)
Estimate the flux from a physical temperature.
DLL_EXPORT double vlbi_astro_estimate_distance_parallax(double rad, double baseline)
Returns the distance of an object from its parallax.
DLL_EXPORT dsp_stream_p vlbi_astro_load_spectrum(char *filename)
Load a spectrum file.
DLL_EXPORT int vlbi_astro_load_spectra_catalog(char *path, dsp_stream_p **catalog, int *catalog_size)
Load a spectrum file catalog.
DLL_EXPORT double vlbi_astro_estimate_brightness_temperature(double wavelength, double flux)
Estimate the brightness temperature from a flux.
DLL_EXPORT double vlbi_astro_flux_ratio(double flux0, double flux, double delta_spectrum)
Returns the flux ratio of two objects.
DLL_EXPORT void vlbi_astro_alt_az_from_ra_dec(double J2000time, double Ra, double Dec, double Lat, double Long, double *Alt, double *Az)
Obtain the altitude and azimuth coordinate of a celestial coordinate at a specific time.
DLL_EXPORT void vlbi_astro_scan_spectrum(dsp_stream_p stream, int sample_size)
Scan a dsp_stream_p and detect the relevant spectral lines.
DLL_EXPORT void vlbi_astro_save_spectrum(dsp_stream_p stream, char *filename)
Save a spectrum stream into a file.
DLL_EXPORT void vlbi_astro_get_alt_az_coordinates(double hour_angle, double dec, double latitude, double *alt, double *az)
Returns alt-azimuth coordinates of an object.
DLL_EXPORT double vlbi_astro_estimate_luminosity_ratio(double size_ratio, double flux_ratio)
Returns the luminosity ratio between two celestial objects.
DLL_EXPORT double vlbi_astro_redshift_adjust(double distance, double redshift)
Returns the distance of a far object adjusted with its measured redshift.
DLL_EXPORT double vlbi_astro_estimate_temperature_ratio(double rad_ratio, double flux_ratio)
Returns the temperature ratio of two objects.
DLL_EXPORT double vlbi_astro_estimate_redshift(double wavelength0, int wavelength)
Returns the redshift of an object.
DLL_EXPORT dsp_stream_p vlbi_file_read_fits(char *filename)
Returns the distance of a far object adjusted with its measured redshift.
DLL_EXPORT double vlbi_astro_estimate_size_transient(double transient_object_velocity, double transit_time)
Returns the size of an object by the observation of a transient body.
DLL_EXPORT double vlbi_astro_estimate_size_ratio(double luminosity_ratio, double temperature_ratio)
Returns the size ratio of two objects.
DLL_EXPORT double vlbi_astro_get_local_hour_angle(double local_sideral_time, double ra)
Returns local hour angle of an object.
DLL_EXPORT double vlbi_astro_mean_speed(double speed)
Obtain or set the reference constant speed of the radiation to measure.
DLL_EXPORT dsp_stream_p * vlbi_file_read_sdfits(char *filename, long *n)
Returns the distance of a far object adjusted with its measured redshift.
DLL_EXPORT dsp_stream_p vlbi_get_baseline_stream(void *ctx, const char *node1, const char *node2)
Obtain the baseline dsp_stream structure containing the complex visibility data.
DLL_EXPORT void vlbi_unlock_baseline(void *ctx, const char *node1, const char *node2)
Unlock the baseline and get visibility from its nodes correlations.
DLL_EXPORT void vlbi_get_offsets(vlbi_context ctx, double J2000Time, const char *node1, const char *node2, double Ra, double Dec, double *offset1, double *offset2)
Get the offsets of a single baseline nodes to the farest node to the target.
DLL_EXPORT void vlbi_set_baseline_buffer(void *ctx, const char *node1, const char *node2, fftw_complex *buffer, int len)
Fill the buffer of a single baseline with complex visibility data. This function locks this baeline a...
DLL_EXPORT int vlbi_get_baselines(void *ctx, vlbi_baseline **baselines)
List all baselines of the current OpenVLBI context.
DLL_EXPORT void vlbi_set_baseline_stream(void *ctx, const char *node1, const char *node2, dsp_stream_p stream)
Set the baseline dsp_stream structure containing the complex visibility data. This function locks thi...
DLL_EXPORT void vlbi_set_location(void *ctx, double lat, double lon, double el)
Set the location of the reference station.
#define sinxsin(s1, s2)
Multiply a sine to a sine.
Definition: vlbi.h:397
double vlbi_default_delegate(double x, double y)
A placeholder delegate that simply multiplies the values received from vlbi_get_uv_plot.
Definition: vlbi.h:630
#define PI
Our PI constant.
Definition: vlbi.h:520
double vlbi_magnitude_correlator_delegate(double x, double y)
A magnitude correlator delegate for vlbi_get_uv_plot.
Definition: vlbi.h:672
double vlbi_phase_correlator_delegate(double x, double y)
A phase correlator delegate for vlbi_get_uv_plot.
Definition: vlbi.h:683
double vlbi_phase_delegate(double x, double y)
A phase calculator delegate for vlbi_get_uv_plot.
Definition: vlbi.h:654
double vlbi_magnitude_delegate(double x, double y)
A magnitude calculator delegate for vlbi_get_uv_plot.
Definition: vlbi.h:642
DLL_EXPORT void vlbi_exit(vlbi_context ctx)
Close a OpenVLBI instance.
DLL_EXPORT const char * vlbi_get_version(void)
Print the current version of OpenVLBI.
DLL_EXPORT unsigned long int vlbi_max_threads(unsigned long value)
get/set the maximum number of threads allowed
DLL_EXPORT vlbi_context vlbi_init(void)
Initialize a OpenVLBI instance.
DLL_EXPORT double * vlbi_matrix_calc_location(double *loc)
Convert geographic location into xyz location.
DLL_EXPORT double * vlbi_matrix_calc_3d_projection(double alt, double az, double *baseline)
Return The 3d projection of the current observation.
DLL_EXPORT double vlbi_matrix_estimate_geocentric_elevation(double latitude, double sea_level_elevation)
Returns an estimation of the actual geocentric elevation.
DLL_EXPORT double vlbi_matrix_estimate_resolution_zero(double frequency)
Estimate the angular resolution of a 1 meter baseline at a given frequency.
DLL_EXPORT double * vlbi_matrix_calc_baseline_center(double *loc1, double *loc2)
Return The baseline center in geographic coordinates.
DLL_EXPORT double * vlbi_matrix_calc_uv_coordinates(double *proj, double wavelength)
Return The UV coordinates of the current observation.
DLL_EXPORT double vlbi_matrix_estimate_resolution(double resolution0, double baseline)
Estimate Signal to noise ratio after a given integration time.
DLL_EXPORT void vlbi_add_model_from_fits(void *ctx, char *filename, const char *name)
Add a model from a fits file.
DLL_EXPORT void vlbi_get_model_to_fits(void *ctx, char *filename, const char *name)
Write a model to a fits file.
DLL_EXPORT void vlbi_add_model(vlbi_context ctx, dsp_stream_p Stream, const char *name)
Add a model into the current OpenVLBI context.
DLL_EXPORT int vlbi_get_models(void *ctx, dsp_stream_p **models)
List all models of the current OpenVLBI context.
DLL_EXPORT void vlbi_get_model_to_jpeg(void *ctx, char *filename, const char *name)
Write a model to a jpeg file.
DLL_EXPORT void vlbi_stack_models(vlbi_context ctx, const char *name, const char *model1, const char *model2)
Stack two models into a new one.
DLL_EXPORT void vlbi_diff_models(vlbi_context ctx, const char *name, const char *model1, const char *model2)
Diff two models into a new one.
DLL_EXPORT void vlbi_add_model_from_png(void *ctx, char *filename, const char *name)
Add a model from a png file.
DLL_EXPORT void vlbi_shift(vlbi_context ctx, const char *name)
Shift a model by its dimensions.
DLL_EXPORT void vlbi_get_ifft(vlbi_context ctx, const char *name, const char *magnitude, const char *phase)
Save into name an inverse fourier transform of the uv plot using its current magnitude and phase comp...
DLL_EXPORT void vlbi_get_fft(vlbi_context ctx, const char *model, const char *magnitude, const char *phase)
Get the fourier transform of the given model and save its magnitude and phase components into two new...
DLL_EXPORT dsp_stream_p vlbi_get_model(void *ctx, const char *name)
Get a single model from the current OpenVLBI context.
DLL_EXPORT void vlbi_del_model(vlbi_context ctx, const char *name)
Remove a model from the current OpenVLBI context.
DLL_EXPORT void vlbi_add_model_from_jpeg(void *ctx, char *filename, const char *name)
Add a model from a jpeg file.
DLL_EXPORT void vlbi_apply_mask(vlbi_context ctx, const char *name, const char *model, const char *mask)
Mask the stream with the content of the mask stream, by multiplication of each element.
DLL_EXPORT void vlbi_get_model_to_png(void *ctx, char *filename, const char *name)
Write a model to a png file.
DLL_EXPORT void vlbi_get_uv_plot(void *ctx, const char *name, int u, int v, double *target, double freq, double sr, int nodelay, int moving_baseline, vlbi_func2_t delegate, int *interrupt)
Fill a fourier plane with an aperture synthesis projection of the baselines during the integration ti...
DLL_EXPORT dsp_stream_p vlbi_get_node(void *ctx, const char *name)
Get a stream from the current OpenVLBI context.
DLL_EXPORT void vlbi_add_node(vlbi_context ctx, dsp_stream_p Stream, const char *name, int geographic_coordinates)
Add a stream into the current OpenVLBI context.
DLL_EXPORT void vlbi_add_node_from_fits(void *ctx, char *filename, const char *name, int geo)
Add a node from a 2d image fits file.
DLL_EXPORT void vlbi_add_nodes_from_sdfits(void *ctx, char *filename, const char *name, int geo)
Add nodes from each row of a single dish fits -SDFITS- file.
DLL_EXPORT void vlbi_del_node(vlbi_context ctx, const char *name)
Remove a stream from the current OpenVLBI context.
DLL_EXPORT int vlbi_get_nodes(void *ctx, vlbi_node **nodes)
List all nodes of the current OpenVLBI context.
DLL_EXPORT timespec_t vlbi_time_string_to_timespec(const char *time)
Obtain a timespec struct containing the date and time specified by a time string.
DLL_EXPORT timespec_t vlbi_time_J2000time_to_timespec(double secs_since_J2000)
Obtain a timespec struct containing the date and time specified by a J2000 time.
DLL_EXPORT timespec_t vlbi_time_mktimespec(int year, int month, int dom, int hour, int minute, int second, long nanosecond)
Obtain a timespec struct containing the date and time specified.
DLL_EXPORT double vlbi_time_timespec_to_J2000time(timespec_t tp)
Convert a timespec into J2000 time.
DLL_EXPORT double vlbi_time_J2000time_to_lst(double secs_since_J2000, double Long)
Obtain the local sidereal time at an arbitrary moment and location.
double(* vlbi_func2_t)(double, double)
The delegate function type to pass to vlbi_plot_uv_plane.
Definition: vlbi.h:337
void * vlbi_context
the OpenVLBI context object type
Definition: vlbi.h:340
struct timespec timespec_t
Definition of the timespec_t in a C type, just for convenience.
Definition: vlbi.h:343
Alignment informations needed.
Definition: dsp.h:269
Contains a set of informations and data relative to a buffer and how to use it.
Definition: dsp.h:362
A single baseline, returned by vlbi_get_baselines into an array.
Definition: vlbi.h:294
int relative
Whether this baseline's nodes use coordinates relative to the context reference station.
Definition: vlbi.h:296
char * Name
Baseline's name.
Definition: vlbi.h:322
double delay
Current delay in seconds.
Definition: vlbi.h:312
double * baseline
The baseline's 3d sizes in meters.
Definition: vlbi.h:306
double u
Current u perspective location in pixels.
Definition: vlbi.h:308
vlbi_node Node2
Latter node.
Definition: vlbi.h:320
int locked
Whether this baseline's nodes use coordinates relative to the context reference station.
Definition: vlbi.h:298
double WaveLength
Wavelength observed in meters.
Definition: vlbi.h:314
dsp_stream_p Stream
Baseline's DSP stream.
Definition: vlbi.h:324
double * Target
The baseline's current celestial target array.
Definition: vlbi.h:300
vlbi_node Node1
Earlier node.
Definition: vlbi.h:318
double Ra
The baseline's current celestial target right ascension.
Definition: vlbi.h:302
double SampleRate
Samples per second.
Definition: vlbi.h:316
double Dec
The baseline's current celestial target declination.
Definition: vlbi.h:304
double v
Current v perspective location in pixels.
Definition: vlbi.h:310
A single node, returned by vlbi_get_nodes into an array.
Definition: vlbi.h:278
int Geo
Whether this node uses geographic coordinates.
Definition: vlbi.h:284
double * GeographicLocation
Geographic coordinates array (latitude, longitude, elevation)
Definition: vlbi.h:280
dsp_stream_p Stream
Node's DSP stream.
Definition: vlbi.h:286
int Index
Node's index - oldest zero.
Definition: vlbi.h:290
double * Location
Current location array in meters.
Definition: vlbi.h:282
char * Name
Node's name.
Definition: vlbi.h:288