FlatMaster Namespace Reference

This modules contains routines to set up the utilities for performing Master Algorithm on flat sky using the J-function approximation to the Wigner 3-j symbols. More...

Classes
struct	PowSpecAtIntegers
	Interpolated Power At Integral l. More...
struct	PQMatrices
	The binning and unbinning matrices. More...
struct	FMUtils
Functions
subroutine	FMCreatePQfromBinFile (bin_def_file, power_of_l, lMax, PQM)
	Create Binning and Unibinning Matrices for power specta defined at integral multipoles (assumes cl starts at l=2).
subroutine	FMFreePQ (PQM)
subroutine	FMInterpolPowerOntoInteg (KPower, PowerAtInt)
	Given a 2-D power spectrum returns as power spectrum Interpolated onto integers.
subroutine	FMFreePowAtI (PowAtI)
subroutine	FMCreateMCMFromWindowPower (WPower, binningFileForInterpolation, Mllprime)
	Given the power spectrum of the window, generates the ModeCoupling Matrix M_l_l'.
subroutine	FMCreateBinnedMCM (Mllprime, binning_file, power_of_l, Mbbprime, doSVD)
subroutine	FMFreeMCM (Mll)
subroutine	FMMasterPowerSpectrum (Mbbprime, Kpower, BinnedMasteredPowSpec)
subroutine	GetCouplingKernel (l1, l2, lmax, wbar, CouplingKernel)
subroutine	ReturnJFunc (k1, k2, JFunc)
real(DP)	pythag (a, b)
subroutine	svbksb (u, w, v, b, x)
Variables
integer(i4b), parameter	MaxL = 7000
	The maximum l beyond which no computation will be done for mastering.

---

Detailed Description

This modules contains routines to set up the utilities for performing Master Algorithm on flat sky using the J-function approximation to the Wigner 3-j symbols.

Author:

Sudeep Das and Amir Hajian, Princeton University

Version:

Time-stamp: <2009-05-19 15:23:29 sudeep>

---

Function Documentation

subroutine FlatMaster::FMCreateBinnedMCM	(	type(FMUtils)	Mllprime,
		character(len = FILENAMELEN )	binning_file,
		real(dp)	power_of_l,
		type(FMUtils)	Mbbprime,
		logical(LGT),optional	doSVD
	)

Parameters:

Mllprime	The unbinned MCM
binning_file	A binnning file definig lower,upper bounds , and centers of desired bins
power_of_l	l^power_of_l will be multiplied when binning
Mbbprime	The binned MCM

subroutine FlatMaster::FMCreateMCMFromWindowPower	(	type(KMap)	WPower,
		character(len = FILENAMELEN)	binningFileForInterpolation,
		type(FMUtils)	Mllprime
	)

Given the power spectrum of the window, generates the ModeCoupling Matrix M_l_l'.

subroutine FlatMaster::FMCreatePQfromBinFile	(	character(len = FILENAMELEN)	Bin_def_file,
		real(dp),dimension(bin_upper(ibin)	power_of_l,
		integer(i4b)	lMax,
		type(PQMatrices)	PQM
	)

Create Binning and Unibinning Matrices for power specta defined at integral multipoles (assumes cl starts at l=2).

subroutine FlatMaster::FMFreeMCM

(

type(FMUtils)

Mll

)

subroutine FlatMaster::FMFreePowAtI

(

type(PowSpecAtIntegers)

PowAtI

)

subroutine FlatMaster::FMFreePQ

(

type(PQMatrices)

PQM

)

subroutine FlatMaster::FMInterpolPowerOntoInteg	(	type(KMap)	KPower,
		type(PowSpecAtIntegers)	PowerAtInt
	)

Given a 2-D power spectrum returns as power spectrum Interpolated onto integers.

subroutine FlatMaster::FMMasterPowerSpectrum	(	type(FMUtils)	MbbPrime,
		type(KMap)	KPower,
		type(BinnedKMap)	BinnedMasteredPowSpec
	)

subroutine FlatMaster::GetCouplingKernel	(	integer(i4b),dimension(abs(l1-l2)	l1,
		integer(i4b),dimension(abs(l1-l2)	l2,
		integer(i4b)	lmax,
		real(dp),dimension(0:lmax)	wbar,
		real(dp)	CouplingKernel
	)

real(DP) FlatMaster::pythag	(	real(DP),intent(in)	a,
		real(DP),intent(in)	b
	)

subroutine FlatMaster::ReturnJFunc	(	integer(i4b),dimension(abs(k1-k2)	k1,
		integer(i4b),dimension(abs(k1-k2)	k2,
		real(dp),dimension(abs(k1-k2)	JFunc
	)

subroutine FlatMaster::svbksb	(	real(DP),dimension(:,:),intent(in)	u,
		real(DP),dimension(:),intent(in)	w,
		real(DP),dimension(:,:),intent(in)	v,
		real(DP),dimension(:),intent(in)	b,
		real(DP),dimension(:),intent(out)	x
	)

---

Variable Documentation

integer(i4b),parameter FlatMaster::MaxL = 7000

The maximum l beyond which no computation will be done for mastering.

---