35#include <lal/LALDatatypes.h>
36#include <lal/AVFactories.h>
37#include <lal/UserInput.h>
38#include <lal/LogPrintf.h>
39#include <lal/SFTfileIO.h>
40#include <lal/LALPulsarVCSInfo.h>
65main(
int argc,
char *argv[] )
88 XLALPrintInfo(
"OK. Maximal relative error %e is within tolerance of %e.\n",
maxd, uvar.relErrorMax );
131 REAL8 reldiff, reldiff_max, reldiff_avg;
133 if ( !ts1 || !ts2 || !ts1->
data || !ts2->
data ) {
139 if ( ts2->
length != numSteps ) {
148 for (
i = 0;
i < numSteps;
i ++ ) {
159 power = fabs( ts1->
data[
i] );
160 if ( power > maxpower ) {
164 total_power += power;
168 reldiff_max =
maxdiff / maxpower;
169 reldiff_avg = sumdiff / total_power;
172 LogPrintf(
LOG_DEBUG,
"%s: total difference = %g, total summed amplitude = %g ==> relative avg error %g\n",
__func__, sumdiff, total_power, reldiff_avg );
174 reldiff =
fmax( reldiff_max, reldiff_avg );
176 return (
REAL4 )reldiff;
197 if ( (
fp = fopen( fname,
"rb" ) ) == NULL ) {
202 if ( ( fread( &magic,
sizeof( magic ), 1,
fp ) != 1 ) || ( magic != 1234.5 ) ) {
208 if ( fread( &len,
sizeof( len ), 1,
fp ) != 1 ) {
220 if ( len != fread( vect->
data,
sizeof( *vect->
data ), len,
fp ) ) {
#define __func__
log an I/O error, i.e.
void LALCheckMemoryLeaks(void)
const LALVCSInfoList lalPulsarVCSInfoList
NULL-terminated list of VCS and build information for LALPulsar and its dependencies
int main(int argc, char *argv[])
REAL4Vector * XLALREAL4VectorFromFile(const CHAR *fname)
Load timeseries from binary input file into REAL4Vector.
REAL4 XLALcompareREAL4Vectors(REAL4Vector *ts1, REAL4Vector *ts2)
Compare two REAL8 vectors, returns a measure of the difference.
REAL4 getMaxErrSFT(const SFTtype *sft1, const SFTtype *sft2)
int vrbflg
defined in lal/lib/std/LALError.c
int initUserVars(UserVar *uvar)
REAL4 getMaxErrSFTVector(const SFTVector *sftvect1, const SFTVector *sftvect2)
#define XLAL_INIT_DECL(var,...)
void LogPrintf(LogLevel_t, const char *format,...) _LAL_GCC_PRINTF_FORMAT_(2
REAL4Vector * XLALCreateREAL4Vector(UINT4 length)
void XLALDestroyREAL4Vector(REAL4Vector *vector)
#define XLAL_ERROR_REAL8(...)
#define XLAL_ERROR_NULL(...)
int int int XLALPrintInfo(const char *fmt,...) _LAL_GCC_PRINTF_FORMAT_(1
#define XLAL_CHECK(assertion,...)
#define XLAL_CHECK_MAIN(assertion,...)
int XLALPrintError(const char *fmt,...) _LAL_GCC_PRINTF_FORMAT_(1
#define XLAL_IS_REAL8_FAIL_NAN(val)
A vector of COMPLEX8FrequencySeries.