86 REAL8 modXi, invModXi;
114 out->nFreqValid =
size->nFreqValid;
123 spinOrder =
par->spin.length;
127 timeDiff =
par->timeDiff;
129 spinF =
par->spin.data;
131 for (
i = 0;
i < spinOrder; ++
i ) {
132 timeDiffProd *= timeDiff;
133 vFactor += spinF[
i] * timeDiffProd;
137 xiX = vFactor * (
par->veloC.x );
138 xiY = vFactor * (
par->veloC.y );
139 xiZ = vFactor * (
par->veloC.z );
144 modXi = sqrt( xiX * xiX + xiY * xiY + xiZ * xiZ );
146 invModXi = 1. / modXi;
148 xiInit.
delta = asin( xiZ * invModXi );
152 xiInit.
alpha = atan2( xiY, xiX );
167 &( *out ).xi, &xiInit, &( *par ).skyPatch ),
status );
174 freqOffset = -modXi * sin(
delta );
175 offset =
out->offset = floor( 0.5 + freqOffset /
deltaF );
176 out->cosPhiMax0 = invModXi * ( offset + 0.5 ) *
deltaF;
177 out->cosPhiMin0 = invModXi * ( offset - 0.5 ) *
deltaF;
static double double delta
#define TRY(func, statusptr)
#define ATTATCHSTATUSPTR(statusptr)
#define ASSERT(assertion, statusptr, code, mesg)
#define DETATCHSTATUSPTR(statusptr)
#define INITSTATUS(statusptr)
#define RETURN(statusptr)
void LALRotatePolarU(LALStatus *status, REAL8UnitPolarCoor *out, REAL8UnitPolarCoor *in, REAL8UnitPolarCoor *par)
void LALNDHOUGHParamPLUT(LALStatus *status, HOUGHParamPLUT *out, HOUGHSizePar *size, HOUGHDemodPar *par)
Demodulation parameters needed for the Hough transform; all coordinates are assumed to be with respec...
Parameters needed to construct the partial look up table.
required for constructing patch
Polar coordinates of a unitary vector on the sphere.
REAL8 delta
In the interval [ ].