31#include <lal/AVFactories.h>
32#include <lal/PulsarDataTypes.h>
33#include <lal/ExtrapolatePulsarSpins.h>
34#include <lal/LALInitBarycenter.h>
35#include <lal/LogPrintf.h>
76 REAL8 minCoverFreq, maxCoverFreq;
80 REAL8 coveringFreqBand = maxCoverFreq - minCoverFreq;
81 fprintf(
fp,
"coveringFreqBand = %g;\n", coveringFreqBand );
94 REAL8 minFreq, maxFreq;
96 doppler.refTime = refTime;
97 doppler.Alpha = skypos_maxdoppler.
longitude;
98 doppler.Delta = skypos_maxdoppler.
latitude;
99 memcpy( doppler.fkdot, fkdot,
sizeof( fkdot ) );
100 doppler.asini = asini;
101 doppler.period = period;
107 REAL8 maxFreqBandWidth = maxFreq - minFreq;
109 fprintf(
fp,
"Alpha_max = %g; Delta_max = %g; minFreq = %g; maxFreq = %g;\nmaxFreqBandWidth = %g;\n",
110 doppler.Alpha, doppler.Delta, minFreq, maxFreq, maxFreqBandWidth );
112 "(maxFreqBandWidth = %g) > (coveringFreqBand = %g)\n",
113 maxFreqBandWidth, coveringFreqBand );
115 fprintf(
fp,
"%%%% Alpha Delta minFreq maxFreq FreqBandWidth\n" );
117 REAL8 time_XLALCWSignalBand = 0;
118 for (
UINT4 i = 0;
i < Ntrials;
i ++ ) {
120 doppler.Alpha =
LAL_TWOPI * ( 1.0 * rand() / ( RAND_MAX + 1.0 ) );
121 doppler.Delta =
LAL_PI_2 - acos( 1 - 2.0 * rand() / RAND_MAX );
125 REAL8 FreqBandWidth = maxFreq - minFreq;
127 "Alpha = %g, Delta = %g: (FreqBandWidth = %g) > (maxFreqBandWidth = %g)\n",
128 doppler.Alpha, doppler.Delta, FreqBandWidth, maxFreqBandWidth );
129 fprintf(
fp,
"%10g, %10g %10.8g %10.8g %10.2e;\n", doppler.Alpha, doppler.Delta, minFreq, maxFreq, FreqBandWidth );
132 time_XLALCWSignalBand /= Ntrials;
134 fprintf(
fp,
"time_XLALCWSignalCoveringBand = %.2e;\n", time_XLALCWSignalCoveringBand );
135 fprintf(
fp,
"time_XLALPrepareCWSignalBand = %.2e;\n", time_XLALPrepareCWSignalBand );
136 fprintf(
fp,
"time_XLALCWSignalBand = %.2e;\n", time_XLALCWSignalBand );
int test_SignalCoveringBand(UINT4 Ntrials)
void LALCheckMemoryLeaks(void)
void XLALDestroyDetectorStateSeries(DetectorStateSeries *detStates)
Get rid of a DetectorStateSeries.
EphemerisData * XLALInitBarycenter(const CHAR *earthEphemerisFile, const CHAR *sunEphemerisFile)
XLAL interface to reading ephemeris files 'earth' and 'sun', and return ephemeris-data in old backwar...
void XLALDestroyEphemerisData(EphemerisData *edat)
Destructor for EphemerisData struct, NULL robust.
#define XLAL_INIT_DECL(var,...)
REAL8 XLALGetCPUTime(void)
REAL8 PulsarSpins[PULSAR_MAX_SPINS]
Typedef for fixed-size array holding GW frequency and derivatives fk = d^k Freq/dt^k|(tau_ref)
const LALDetector * XLALGetSiteInfo(const CHAR *name)
Find the site geometry-information 'LALDetector' for given a detector name (or prefix).
#define XLAL_CHECK_MAIN(assertion,...)
LIGOTimeGPS * XLALGPSAdd(LIGOTimeGPS *epoch, REAL8 dt)
Timeseries of DetectorState's, representing the detector-info at different timestamps.
This structure contains all information about the center-of-mass positions of the Earth and Sun,...
Type containing the 'Doppler-parameters' affecting the time-evolution of the phase.
Contains a "spin-range", ie spins and corresponding bands at a given (SSB) reference GPS-time .