Serial Spectral Extraction¶

The serial spectral extraction tools in HARP are useful debugging new data plugins or testing algorithms on small toy problems. If you are performing production spectral extraction, you should use the parallel extraction tools.

class harp::spec_slice_region

Public Functions

spec_slice_region()

~spec_slice_region()

spec_slice_region(spec_slice_region const & orig)

spec_slice_region & operator=(spec_slice_region const & rhs)

Public Members

size_t overlap_spec

size_t overlap_lambda

size_t first_spec

size_t first_lambda

size_t first_good_spec

size_t first_good_lambda

size_t n_spec

size_t n_lambda

size_t n_good_spec

size_t n_good_lambda

Friends

friend class boost::serialization::access

class harp::spec_slice

Public Functions

spec_slice()

spec_slice(size_t nworker, size_t nspec, size_t nlambda, size_t chunk_nspec, size_t chunk_nlambda, size_t overlap_spec, size_t overlap_lambda)

~spec_slice()

spec_slice_region const & full_region() const

std::vector< spec_slice_region > regions(size_t const & worker) const

Friends

friend class boost::serialization::access

void sub_spec(spec_slice_region const & full_region, spec_slice_region const & sub_region, vector_double const & full_data, bool use_good_sub, vector_double & sub_data)

void accum_spec(spec_slice_region const & sub_region, spec_slice_region const & full_region, vector_double const & sub_data, bool use_good_sub, vector_double & full_data)

void noise_weighted_spec(matrix_double_sparse const & AT, vector_double const & invnoise, vector_mask const & mask, vector_double const & img, vector_double & z)

void inverse_covariance(matrix_double_sparse const & AT, vector_double const & invnoise, vector_mask const & mask, matrix_double & invC)

void resolution(vector_double const & D, matrix_double const & W, vector_double & S, matrix_double & R)

void accum_diag_resolution(spec_slice_region const & region, matrix_double const & R, matrix_double_sparse & diag)

void extract(vector_double const & D, matrix_double const & W, vector_double const & S, vector_double const & z, vector_double & Rf, vector_double & f)

void extract_slices(spec_slice_p slice, psf_p design, vector_double const & img, vector_double const & img_inv_var, vector_double const & truth, vector_double & Rf, vector_double & f, vector_double & err, vector_double & Rtruth, std::map< std::string, double > & profile, bool region_threads, bool lambda_mask, std::string const & status_prefix = "")