Commit 6b7088a1 authored by Karl Wette's avatar Karl Wette

SuperskyMetrics: improved equal-area patch bounds

- Divide reduced supersky parameter space (two unit
  disks) into roughly square patches of equal area
Original: 4e2c021995167de708d725c6eabfcdad231f78b6
parent 707224f3
This diff is collapsed.
......@@ -216,14 +216,12 @@ int XLALSetSuperskyPhysicalSkyBounds(
);
///
/// Divide the physical sky into \p patch_count equal-area patches, and return physical sky position
/// bounds \f$(\alpha, \delta)\f$ for the patch indexed by \p patch_index.
///
int XLALComputePhysicalSkyEqualAreaPatch(
double *alpha1, ///< [out] First bound on sky position right ascension
double *alpha2, ///< [out] Second bound on sky position right ascension
double *delta1, ///< [out] First bound on sky position declination
double *delta2, ///< [out] Second bound on sky position declination
/// Divide the reduced supersky parameter space into \p patch_count equal-area patches, and
/// set parameter-space bounds on the reduced supersky coordinates \f$(n_a,n_b)\f$ for the
/// patch indexed by \p patch_index.
///
int XLALSetSuperskyEqualAreaSkyBounds(
LatticeTiling *tiling, ///< [in] Lattice tiling
const UINT4 patch_count, ///< [in] Number of equal-area patches to divide sky into
const UINT4 patch_index ///< [in] Index of the patch for which to compute bounds
);
......
......@@ -670,9 +670,8 @@ static int SuperskyTests(
XLAL_CHECK( semi_tiling != NULL, XLAL_EFUNC );
// Add bounds
double alpha1 = 0, alpha2 = 0, delta1 = 0, delta2 = 0;
const double alpha1 = 0, alpha2 = LAL_PI, delta1 = -LAL_PI_2, delta2 = LAL_PI_2;
const double freq_min = freq_max - 5e-5, f1dot = -5e-9;
XLAL_CHECK( XLALComputePhysicalSkyEqualAreaPatch( &alpha1, &alpha2, &delta1, &delta2, 2, 0 ) == XLAL_SUCCESS, XLAL_EFUNC );
for ( size_t n = 0; n < metrics->num_segments; ++n ) {
XLAL_CHECK( XLALSetSuperskyPhysicalSkyBounds( coh_tiling[n], metrics->coh_rssky_metric[n], metrics->coh_rssky_transf[n], alpha1, alpha2, delta1, delta2 ) == XLAL_SUCCESS, XLAL_EFUNC );
XLAL_CHECK( XLALSetSuperskyPhysicalSpinBound( coh_tiling[n], metrics->coh_rssky_transf[n], 0, freq_min, freq_max ) == XLAL_SUCCESS, XLAL_EFUNC );
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment