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
|
/* CALTABLE - calibration table interface to PolConvert
Copyright (C) 2013 Ivan Marti-Vidal
Nordic Node of EU ALMA Regional Center (Onsala, Sweden)
Max-Planck-Institut fuer Radioastronomie (Bonn, Germany)
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 3 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, see <http://www.gnu.org/licenses/>
*/
#include <sys/types.h>
#include <iostream>
#include <complex>
/* Class to read a calibration table and interpolate
its gains in frequency and time. */
class CalTable {
public:
CalTable(int kind, double **R1,double **I1,double **R2,double **I2, double *freqs, double **times, int Na, long *Nt, long Nc, bool **flag, bool islinear, FILE *logF, bool verbose);
~CalTable();
int getNant();
long getNchan();
bool isBandpass();
void getTimeRange(double *JD);
void getFreqRange(double *Fr);
long getNEntries(int ant);
void getFrequencies(double *freqs);
void getTimes(int ant, double *times);
void getGains(int ant, long timeidx, double *gain[4]);
void setChanged(bool ch);
/* Prepares the instance for the frequency interpolation. User must provide
the array of frequencies to interpolate (*freqs) and the number
of channels in that array (mschan). */
void setMapping(long mschan, double *freqs);
/* Prepares the instance for the time interpolation. Returns False if the interpolation coefficients have not changed (so it would be a waste of resources to recompute everything). */
bool setInterpolationTime(double itime);
/* Interpolates the gains of an antenna (iant) at a given time, itime,
and applies them to the array gain[2] (elements of this array
are the different polarizations (X,Y)). User should have run "setMapping" before,
in case that the table channel frequencies do not match with the array of frequencies
where the interpolation is desired. Then,
If mode==0, the gains are just written to gain.
If mode==1, the gains are ADDED to the already-existing values in gain.
If mode==2, the gains are MULTIPLIED to the already-existing values in gain.
-------------------
*/
void applyInterpolation(int iant, int mode, std::complex<float>* gain[2]);
// Same as above, but one gain (for one channel) is returned:
bool getInterpolation(int iant, int ichan, std::complex<float> gain[2]);
private:
FILE *logFile;
char message[512];
void fillGaps(); // Fills flagged gains with interpolated values.
static const int Nmax = 256; // Maximum number of antennas.
std::string name;
int Nants;
long *Ntimes;
long Nchan;
bool SignFreq, success;
double ***GainAmp[2];
double ***GainPhase[2];
bool **flags, *firstTime;
double *BuffPhase[2];
double *BuffAmp[2];
double **Time;
double currTime;
double *Freqs;
double JDRange[2];
double *K0;
long *I0;
long *I1;
long MSChan;
double *preKt;
long *pret0, *pret1;
bool isDelay, gainChanged, isLinear, isDterm, isTsys, Verbose;
double deltaNu0, deltaNu;
std::complex<float>** bufferGain[2];
};
|