vlbi.h

/*  OpenVLBI - Open Source Very Long Baseline Interferometry
*   Copyright © 2017-2023  Ilia Platone
*
*   This program is free software; you can redistribute it and/or modify
*   it under the terms of the GNU General Public License as published by
*   the Free Software Foundation; either version 2 of the License, or
*   (at your option) any later version.
*
*   This program is distributed in the hope that it will be useful,
*   but WITHOUT ANY WARRANTY; without even the implied warranty of
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*   GNU General Public License for more details.
*
*   You should have received a copy of the GNU General Public License along
*   with this program; if not, write to the Free Software Foundation, Inc.,
*   51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*/
 
#ifndef _VLBI_H
#define _VLBI_H
 
#ifdef  __cplusplus
extern "C" {
#endif
 
#ifndef DLL_EXPORT
#ifdef _WIN32
#define DLL_EXPORT __declspec(dllexport)
#else
#define DLL_EXPORT extern
#endif
#endif
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <errno.h>
#include <unistd.h>
#include <string.h>
#include <math.h>
#include <float.h>
#include <libgen.h>
#include <sys/types.h>
#include <assert.h>
#include <pthread.h>
#include <dsp.h>
 
typedef struct {
    double *GeographicLocation;
    double *Location;
    int Geo;
    dsp_stream_p Stream;
    char *Name;
    int Index;
} vlbi_node;
 
typedef struct {
    int relative;
    int locked;
    double *Target;
    double Ra;
    double Dec;
    double *baseline;
    double u;
    double v;
    double delay;
    double WaveLength;
    double SampleRate;
    int nodes_count;
    vlbi_node *Nodes;
    char *Name;
    dsp_stream_p Stream;
} vlbi_baseline;
 
typedef double(* vlbi_func2_t)(double, double);
 
typedef void* vlbi_context;
 
typedef struct timespec timespec_t;
#define VLBI_BASELINE_NAME_SIZE DSP_NAME_SIZE*2+1
 
#ifndef Min
#define Min(a,b) \
   ({ __typeof (a) _a = (a); \
       __typeof (b) _b = (b); \
     _a < _b ? _a : _b; })
#endif
#ifndef Max
#define Max(a,b) \
   ({ __typeof (a) _a = (a); \
       __typeof (b) _b = (b); \
     _a > _b ? _a : _b; })
#endif
#ifndef Log
#define Log(a,b) \
( log(a) / log(b) )
#endif
 
#ifndef hz2rad
#define hz2rad(hz) (2.0*PI*hz)
#endif
 
#ifndef sinpsin
#define sinpsin(s1, s2) (2.0*sin((asin(s1)+asin(s2))/2.0)*cos((asin(s1)-asin(s2))/2.0))
#endif
 
#ifndef sinmsin
#define sinmsin(s1, s2) (2.0*cos((asin(s1)+asin(s2))/2.0)*sin((asin(s1)-asin(s2))/2.0))
#endif
 
#ifndef cospcos
#define cospcos(c1, c2) (2.0*cos((acos(s1)+acos(s2))/2.0)*cos((acos(s1)-acos(s2))/2.0))
#endif
 
#ifndef cosmcos
#define cosmcos(c1, c2) (2.0*sin((acos(s1)+acos(s2))/2.0)*sin((acos(s1)-acos(s2))/2.0))
#endif
 
#ifndef sinxsin
#define sinxsin(s1, s2) ((cos(asin(s1)-asin(s2))-cos(asin(s1)+asin(s2)))/2.0)
#endif
 
#ifndef cosxcos
#define cosxcos(c1, c2) ((cos(acos(s1)+acos(s2))+cos(acos(s1)-acos(s2)))/2.0)
#endif
 
#ifndef sinxcos
#define sinxcos(s, c) ((sin(asin(s)+acos(c))+sin(asin(s)-acos(c)))/2.0)
#endif
 
#ifndef sin2cos
#define sin2cos(s) cos(asin(s))
#endif
 
#ifndef cos2sin
#define cos2sin(c) sin(acos(c))
#endif
 
#ifndef VLBI_VERSION_STRING
#define VLBI_VERSION_STRING "3.0.3"
#endif
 
#ifndef VLBI_CATALOG_PATH
#define VLBI_CATALOG_PATH "/usr/share/OpenVLBI/cat/index.txt"
#endif
 
#ifndef CIRCLE_DEG
#define CIRCLE_DEG 360.0
#endif
#ifndef CIRCLE_AM
#define CIRCLE_AM (CIRCLE_DEG * 60.0)
#endif
#ifndef CIRCLE_AS
#define CIRCLE_AS (CIRCLE_AM * 60.0)
#endif
#ifndef RAD_AS
#define RAD_AS (CIRCLE_AS/(PI*2.0))
#endif
#ifndef ONE_SECOND_TICKS
#define ONE_SECOND_TICKS 100000000
#endif
#ifndef ONE_MILLISECOND_TICKS
#define ONE_MILLISECOND_TICKS 100000
#endif
#ifndef ONE_MICROSECOND_TICKS
#define ONE_MICROSECOND_TICKS 100
#endif
#ifndef SOLAR_DAY
#define SOLAR_DAY 86400
#endif
#ifndef SIDEREAL_DAY
#define SIDEREAL_DAY 86164.0905
#endif
#ifndef TRACKRATE_SIDEREAL
#define TRACKRATE_SIDEREAL (CIRCLE_AS / SIDEREAL_DAY)
#endif
#ifndef TRACKRATE_SOLAR
#define TRACKRATE_SOLAR (CIRCLE_AS / SOLAR_DAY)
#endif
#ifndef TRACKRATE_LUNAR
#define TRACKRATE_LUNAR 14.511415
#endif
#ifndef EARTHRADIUSEQUATORIAL
#define EARTHRADIUSEQUATORIAL 6378137.0
#endif
#ifndef EARTHRADIUSPOLAR
#define EARTHRADIUSPOLAR 6356752.0
#endif
#ifndef EARTHRADIUSMEAN
#define EARTHRADIUSMEAN 6372797.0
#endif
#ifndef AVOGADRO
#define AVOGADRO 6.02214076E+23
#endif
#ifndef EULER
#define EULER 2.71828182845904523536028747135266249775724709369995
#endif
#ifndef PLANK
#define PLANK 6.62607015E-34
#endif
#ifndef BOLTSMANN
#define BOLTSMANN 1.380649E-23
#endif
#ifndef STEPHAN_BOLTSMANN
#define STEPHAN_BOLTSMANN (2.0*pow(PI, 5)*pow(BOLTSMANN, 2)/(pow(LIGHTSPEED, 2)*15*pow(PLANK, 3)))
#endif
#ifndef GAS_R
#define GAS_R (BOLTSMANN * AVOGADRO)
#endif
#ifndef ROOT2
#define ROOT2 1.41421356237309504880168872420969807856967187537694
#endif
#ifndef PI
#define PI 3.14159265358979323846
#endif
#ifndef AIRY
#define AIRY 1.21966
#endif
#ifndef LIGHTSPEED
#define LIGHTSPEED 299792458.0
#endif
#ifndef J2000
#define J2000 2451545.0
#endif
#ifndef GAMMAJ2000
#define GAMMAJ2000 18.6971378528
#endif
#ifndef ELECTRON
#define ELECTRON 1.602176634E-19
#endif
#ifndef PROTON
#define PROTON 1.67262192E-24
#endif
#ifndef CANDLE
#define CANDLE 0.683
#endif
#ifndef ASTRONOMICALUNIT
#define ASTRONOMICALUNIT 1.495978707E+11
#endif
#ifndef PARSEC
#define PARSEC (ASTRONOMICALUNIT/sin(PI*2.0/CIRCLE_AS))
#endif
#ifndef LY
#define LY (LIGHTSPEED * SIDEREAL_DAY * 365.0)
#endif
#ifndef AU2M
#define AU2M(au) (au * ASTRONOMICALUNIT)
#endif
#ifndef PARSEC2M
#define PARSEC2M(parsec) (parsec * PARSEC)
#endif
#ifndef LY2M
#define LY2M(ly) (ly * LY)
#endif
#ifndef M2AU
#define M2AU(m) (m / ASTRONOMICALUNIT)
#endif
#ifndef M2PARSEC
#define M2PARSEC(m) (m / PARSEC)
#endif
#ifndef M2LY
#define M2LY(m) (m / LY)
#endif
#ifndef RAD2AS
#define RAD2AS(rad) (rad * RAD_AS)
#endif
#ifndef AS2RAD
#define AS2RAD(as) (as / RAD_AS);
#endif
 
inline double vlbi_default_delegate(double x, double y) {
    return x*y;
}
 
inline double vlbi_magnitude_delegate(double x, double y) {
    return sqrt(pow(x, 2)+pow(y, 2));
}
 
inline double vlbi_phase_delegate(double x, double y) {
    double mag = sqrt(pow(x, 2)+pow(y, 2));
    double rad = 0.0;
    if(mag > 0.0) {
        rad = acos (y / (mag > 0.0 ? mag : 1.0));
        if(x < 0 && rad != 0)
            rad = PI*2-rad;
    }
    return rad;
}
 
inline double vlbi_magnitude_correlator_delegate(double x, double y) {
    return x*y;
}
 
inline double vlbi_phase_correlator_delegate(double x, double y) {
    return fmod(x+y, M_PI*2);
}
 
DLL_EXPORT unsigned long int vlbi_max_threads(unsigned long value);
 
DLL_EXPORT const char *vlbi_get_version(void);
 
DLL_EXPORT vlbi_context vlbi_init(void);
 
DLL_EXPORT void vlbi_exit(vlbi_context ctx);
 
DLL_EXPORT void vlbi_add_node(vlbi_context ctx, dsp_stream_p Stream, const char *name, int geographic_coordinates);
 
DLL_EXPORT void vlbi_copy_node(void *ctx, const char *name, const char *node);
 
DLL_EXPORT dsp_stream_p vlbi_get_node(void *ctx, const char *name);
 
DLL_EXPORT int vlbi_has_node(void *ctx, const char *name);
 
DLL_EXPORT void vlbi_del_node(vlbi_context ctx, const char *name);
 
DLL_EXPORT int vlbi_get_nodes(void *ctx, vlbi_node** nodes);
 
DLL_EXPORT void vlbi_add_node_from_fits(void *ctx, char *filename, const char *name, int geo);
 
DLL_EXPORT void vlbi_add_nodes_from_sdfits(void *ctx, char *filename, const char *name, int geo);
 
DLL_EXPORT void vlbi_filter_lp_node(void *ctx, const char *name, const char *node, double radians);
 
DLL_EXPORT void vlbi_filter_hp_node(void *ctx, const char *name, const char *node, double radians);
 
DLL_EXPORT void vlbi_filter_bp_node(void *ctx, const char *name, const char *node, double lo_radians, double hi_radians);
 
DLL_EXPORT void vlbi_filter_br_node(void *ctx, const char *name, const char *node, double lo_radians, double hi_radians);
 
DLL_EXPORT void vlbi_set_correlation_order(void *ctx, int order);
 
DLL_EXPORT int vlbi_get_baselines(void *ctx, vlbi_baseline** baselines);
 
DLL_EXPORT void vlbi_free_baselines(vlbi_baseline** baselines, int num_baselines);
 
DLL_EXPORT void vlbi_set_baseline_buffer(void *ctx, const char **nodes, complex_t *buffer, int len);
 
DLL_EXPORT dsp_stream_p vlbi_get_baseline_stream(void *ctx, const char**nodes);
 
DLL_EXPORT void vlbi_unlock_baseline(void *ctx, const char**nodes);
 
DLL_EXPORT void vlbi_set_baseline_stream(void *ctx, const char**nodes, dsp_stream_p stream);
 
DLL_EXPORT void vlbi_set_location(void *ctx, double lat, double lon, double el);
 
DLL_EXPORT double vlbi_get_offset(vlbi_context ctx, double J2000Time, const char* node, double Ra, double Dec, double Distance);
 
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 void vlbi_add_model(vlbi_context ctx, dsp_stream_p Stream, const char *name);
 
DLL_EXPORT void vlbi_copy_model(void *ctx, const char *name, const char *node);
 
DLL_EXPORT void vlbi_del_model(vlbi_context ctx, const char *name);
 
DLL_EXPORT int vlbi_get_models(void *ctx, dsp_stream_p** models);
 
DLL_EXPORT dsp_stream_p vlbi_get_model(void *ctx, const char *name);
 
DLL_EXPORT int vlbi_has_model(void *ctx, const char *name);
 
DLL_EXPORT void vlbi_get_ifft(vlbi_context ctx, const char *name, const char *magnitude, const char *phase);
 
DLL_EXPORT void vlbi_get_fft(vlbi_context ctx, const char *model, const char *magnitude, const char *phase);
 
DLL_EXPORT void vlbi_apply_mask(vlbi_context ctx, const char *name, const char *model, const char *mask);
 
DLL_EXPORT void vlbi_apply_convolution_matrix(vlbi_context ctx, const char *name, const char *model, const char *matrix);
 
DLL_EXPORT void vlbi_stack_models(vlbi_context ctx, const char *name, const char *model1, const char *model2);
 
DLL_EXPORT void vlbi_diff_models(vlbi_context ctx, const char *name, const char *model1, const char *model2);
 
DLL_EXPORT double vlbi_compare_models(vlbi_context ctx, const char *model1, const char *model2);
 
DLL_EXPORT void vlbi_shift(vlbi_context ctx, const char *name);
 
DLL_EXPORT void vlbi_add_model_from_png(void *ctx, char *filename, const char *name);
 
DLL_EXPORT void vlbi_add_model_from_jpeg(void *ctx, char *filename, const char *name);
 
DLL_EXPORT void vlbi_add_model_from_fits(void *ctx, char *filename, const char *name);
 
DLL_EXPORT void vlbi_get_model_to_png(void *ctx, char *filename, const char *name);
 
DLL_EXPORT void vlbi_get_model_to_jpeg(void *ctx, char *filename, const char *name);
 
DLL_EXPORT void vlbi_get_model_to_fits(void *ctx, char *filename, const char *name);
 
DLL_EXPORT double* vlbi_matrix_calc_baseline_center(double *loc1, double *loc2);
 
DLL_EXPORT double* vlbi_matrix_calc_3d_projection(double alt, double az, double *baseline);
 
DLL_EXPORT double* vlbi_matrix_calc_parametric_projection(double *target, double *baseline);
 
DLL_EXPORT double* vlbi_matrix_calc_uv_coordinates(double *proj, double wavelength);
 
DLL_EXPORT double* vlbi_matrix_calc_location(double *loc);
 
DLL_EXPORT double vlbi_matrix_estimate_geocentric_elevation(double latitude, double sea_level_elevation);
 
DLL_EXPORT double vlbi_matrix_estimate_resolution_zero(double frequency);
 
DLL_EXPORT double vlbi_matrix_estimate_resolution(double resolution0, double baseline);
 
DLL_EXPORT timespec_t vlbi_time_mktimespec(int year, int month, int dom, int hour, int minute, int second, long nanosecond);
 
DLL_EXPORT double vlbi_time_timespec_to_J2000time(timespec_t tp);
 
DLL_EXPORT double vlbi_time_J2000time_to_lst(double secs_since_J2000, double Long);
 
DLL_EXPORT timespec_t vlbi_time_string_to_timespec(const char *time);
 
DLL_EXPORT timespec_t vlbi_time_J2000time_to_timespec(double secs_since_J2000);
 
DLL_EXPORT double vlbi_astro_mean_speed(double speed);
 
DLL_EXPORT void vlbi_astro_alt_az_from_ra_dec(double J2000time, double Ra, double Dec, double Lat, double Long, double* Alt, double *Az);
 
DLL_EXPORT double vlbi_astro_get_local_hour_angle(double local_sideral_time, double ra);
 
DLL_EXPORT void vlbi_astro_get_alt_az_coordinates(double hour_angle, double dec, double latitude, double* alt, double *az);
 
DLL_EXPORT dsp_stream_p vlbi_astro_load_spectrum(char *filename);
 
DLL_EXPORT int vlbi_astro_load_spectra_catalog(char *path, dsp_stream_p **catalog, int *catalog_size);
 
DLL_EXPORT dsp_stream_p vlbi_astro_create_reference_catalog(dsp_stream_p *catalog, int catalog_size);
 
DLL_EXPORT void vlbi_astro_save_spectrum(dsp_stream_p stream, char *filename);
 
DLL_EXPORT void vlbi_astro_scan_spectrum(dsp_stream_p stream, int sample_size);
 
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);
 
DLL_EXPORT double vlbi_astro_diff_spectra(dsp_stream_p spectrum0, dsp_stream_p spectrum, double decades);
 
DLL_EXPORT double vlbi_astro_flux_ratio(double flux0, double flux, double delta_spectrum);
 
DLL_EXPORT double vlbi_astro_estimate_brightness_temperature(double wavelength, double flux);
 
DLL_EXPORT double vlbi_astro_estimate_temperature(double wavelength, double flux);
 
DLL_EXPORT double vlbi_astro_estimate_flux(double wavelength, double temperature);
 
DLL_EXPORT double vlbi_astro_estimate_temperature_ratio(double rad_ratio, double flux_ratio);
 
DLL_EXPORT double vlbi_astro_estimate_size_ratio(double luminosity_ratio, double temperature_ratio);
 
DLL_EXPORT double vlbi_astro_estimate_luminosity_ratio(double size_ratio, double flux_ratio);
 
DLL_EXPORT double vlbi_astro_estimate_distance_ratio(double luminosity_ratio, double flux_ratio);
 
DLL_EXPORT double vlbi_astro_estimate_distance_parallax(double rad, double baseline);
 
DLL_EXPORT double vlbi_astro_estimate_redshift(double wavelength0, int wavelength);
 
DLL_EXPORT double vlbi_astro_estimate_size_transient(double transient_object_velocity, double transit_time);
 
DLL_EXPORT double vlbi_astro_redshift_adjust(double distance, double redshift);
 
DLL_EXPORT dsp_stream_p vlbi_file_read_fits(char *filename);
 
DLL_EXPORT dsp_stream_p *vlbi_file_read_sdfits(char * filename, long *n);
 
#ifdef __cplusplus
}
#endif
 
#endif //_VLBI_H