weight_utils.h File Reference

#include <vector>
#include <memory>
#include "healpix_map.h"

Go to the source code of this file.

Functions
template<typename T >
void	apply_fullweights (Healpix_Map< T > &map, const std::vector< double > &wgt)
std::vector< double >	get_fullweights (int nside, int lmax, double epsilon, int itmax, double &epsilon_out)
std::vector< double >	get_ringweights (int nside, int lmax, double epsilon, int itmax, double &epsilon_out)

Detailed Description

Copyright (C) 2016-2019 Max-Planck-Society

Author

Martin Reinecke

Definition in file weight_utils.h.

Function Documentation

apply_fullweights()

template<typename T >

void apply_fullweights	(	Healpix_Map< T > &	map,
		const std::vector< double > &	wgt
	)

Applies the vector wgt containing compressed full weights to map.

get_fullweights()

std::vector<double> get_fullweights	(	int	nside,
		int	lmax,
		double	epsilon,
		int	itmax,
		double &	epsilon_out
	)

Computes full weights for a map of the given nside at a maximum multipole lmax. The solution is obtained via CGNE iteration, which stops once the norm of the residual falls below epsilon times the norm of the initial residual, or once itmax iterations have been performed.

Returns

a vector containing the compressed full weights.

Note

lmax must be even; odd l do not contribute to the weights.

Definition at line 340 of file weight_utils.cc.

get_ringweights()

std::vector<double> get_ringweights	(	int	nside,
		int	lmax,
		double	epsilon,
		int	itmax,
		double &	epsilon_out
	)

Computes ring weights for a map of the given nside at a maximum multipole lmax. The solution is obtained via CGNE iteration, which stops once the norm of the residual falls below epsilon times the norm of the initial residual, or once itmax iterations have been performed.

Returns

a vector of size 2*nside containing the ring weights.

Note

lmax must be even; odd l do not contribute to the weights.

Definition at line 380 of file weight_utils.cc.