#include <iostream>
#include "announce.h"
#include "paramfile.h"
#include "healpix_map_fitsio.h"
#include "lsconstants.h"
#include "arr.h"
#include "fitshandle.h"
#include "vec3.h"
#include "string_utils.h"
#include "alm.h"
#include "alm_healpix_tools.h"
#include "planck_rng.h"
#include "healpix_map.h"

Functions
void	get_step (const PolarizationHolder &ph, vec3 &loc, vec3 &dir, double theta)

void	runge_kutta_step (vec3 &loc, vec3 &dir, const PolarizationHolder &ph, double theta)

void	runge_kutta_2 (const vec3 &location, const PolarizationHolder &ph, double theta, arr< vec3 > &locs)

void	make_kernel (arr< double > &kernel)

void	convolve (const arr< double > &kernel, const arr< double > &raw, arr< double > &convolution)

int	lic_function (Healpix_Map< double > &hitcount, Healpix_Map< double > &texture, const PolarizationHolder &ph, const Healpix_Map< double > &th, int steps, int kernel_steps, double step_radian)

void	lic_main (const Healpix_Map< double > &Q, const Healpix_Map< double > &U, const Healpix_Map< double > &th, Healpix_Map< double > &hit, Healpix_Map< double > &tex, Healpix_Map< double > &mag, int steps, int kernel_steps, double step_radian, double polmin, double polmax)

Detailed Description

Author: David Larson; Martin Reinecke

Definition in file alice3.cc.

Function Documentation

◆ get_step()

void get_step	(	const PolarizationHolder &	ph,
		vec3 &	loc,
		vec3 &	dir,
		double	theta
	)

Steps from loc in direction dir for an angle theta and updates loc and dir.

Definition at line 84 of file alice3.cc.

◆ runge_kutta_step()

void runge_kutta_step	(	vec3 &	loc,
		vec3 &	dir,
		const PolarizationHolder &	ph,
		double	theta
	)

Performs one Runge-Kutta second order step. Updates loc and dir.

Definition at line 92 of file alice3.cc.

◆ runge_kutta_2()

void runge_kutta_2	(	const vec3 &	location,
		const PolarizationHolder &	ph,
		double	theta,
		arr< vec3 > &	locs
	)

Second order Runge-Kutta integration on the sphere. Given a starting location, a qu map of the sky, and a step size theta, this subroutine returns an array of vectors extending in both directions from the starting location.

Definition at line 107 of file alice3.cc.

◆ make_kernel()

void make_kernel ( arr< double > & kernel )

Create a sinusoidal kernel.

Definition at line 132 of file alice3.cc.

◆ convolve()

void convolve	(	const arr< double > &	kernel,
		const arr< double > &	raw,
		arr< double > &	convolution
	)

Convolve an array with a kernel.

Definition at line 142 of file alice3.cc.

◆ lic_function()

int lic_function	(	Healpix_Map< double > &	hitcount,
		Healpix_Map< double > &	texture,
		const PolarizationHolder &	ph,
		const Healpix_Map< double > &	th,
		int	steps,
		int	kernel_steps,
		double	step_radian
	)

Perform line integral convolution on sphere.

Definition at line 155 of file alice3.cc.

◆ lic_main()

void lic_main	(	const Healpix_Map< double > &	Q,
		const Healpix_Map< double > &	U,
		const Healpix_Map< double > &	th,
		Healpix_Map< double > &	hit,
		Healpix_Map< double > &	tex,
		Healpix_Map< double > &	mag,
		int	steps,
		int	kernel_steps,
		double	step_radian,
		double	polmin,
		double	polmax
	)

Expose line integral convolution for external use.

Definition at line 188 of file alice3.cc.

Functions