#include "output.h"

Include dependency graph for output.c:

Functions
int	output_init (struct background pba, struct thermodynamics pth, struct perturbations ppt, struct primordial ppm, struct transfer ptr, struct harmonic phr, struct fourier pfo, struct lensing ple, struct distortions psd, struct output pop)

int	output_cl (struct background pba, struct perturbations ppt, struct harmonic phr, struct lensing ple, struct output *pop)

int	output_pk (struct background pba, struct perturbations ppt, struct fourier pfo, struct output pop, enum pk_outputs pk_output)

int	output_tk (struct background pba, struct perturbations ppt, struct output *pop)

int	output_heating (struct injection pin, struct noninjection pni, struct output *pop)

int	output_distortions (struct distortions psd, struct output pop)

int	output_print_data (FILE out, char titles[_MAXTITLESTRINGLENGTH_], double dataptr, int size_dataptr)

int	output_open_cl_file (struct harmonic phr, struct output pop, FILE *clfile, FileName filename, char first_line, int lmax)

int	output_one_line_of_cl (struct background pba, struct harmonic phr, struct output pop, FILE clfile, double l, double *cl, int ct_size)

int	output_open_pk_file (struct background pba, struct fourier pfo, struct output pop, FILE pkfile, FileName filename, char first_line, double z)

int	output_one_line_of_pk (FILE *pkfile, double one_k, double one_pk)

Detailed Description

Documented output module

Julien Lesgourgues, 26.08.2010

This module writes the output in files.

The following functions can be called from other modules or from the main:

output_init() (must be called after harmonic_init())
output_total_cl_at_l() (can be called even before output_init())

No memory needs to be deallocated after that, hence there is no output_free() routine like in other modules.

Function Documentation

◆ output_init()

int output_init	(	struct background *	pba,
		struct thermodynamics *	pth,
		struct perturbations *	ppt,
		struct primordial *	ppm,
		struct transfer *	ptr,
		struct harmonic *	phr,
		struct fourier *	pfo,
		struct lensing *	ple,
		struct distortions *	psd,
		struct output *	pop
	)

This routine writes the output in files.

Parameters

pba	Input: pointer to background structure (needed for calling harmonic_pk_at_z())
pth	Input: pointer to thermodynamics structure
ppt	Input: pointer perturbation structure
ppm	Input: pointer to primordial structure
ptr	Input: pointer to transfer structure
phr	Input: pointer to harmonic structure
pfo	Input: pointer to fourier structure
ple	Input: pointer to lensing structure
psd	Input: pointer to distortions structure
pop	Input: pointer to output structure

Summary:

check that we really want to output at least one file
deal with all anisotropy power spectra 's
deal with all Fourier matter power spectra P(k)'s
deal with density and matter power spectra
deal with background quantities
deal with thermodynamics quantities
deal with perturbation quantities
deal with primordial spectra
deal with heating
deal with spectral distortions

◆ output_cl()

int output_cl	(	struct background *	pba,
		struct perturbations *	ppt,
		struct harmonic *	phr,
		struct lensing *	ple,
		struct output *	pop
	)

This routines writes the output in files for anisotropy power spectra 's.

Parameters

pba	Input: pointer to background structure (needed for $T_{cmb}$ )
ppt	Input: pointer perturbation structure
phr	Input: pointer to harmonic structure
ple	Input: pointer to lensing structure
pop	Input: pointer to output structure

Summary:

define local variables
first, allocate all arrays of files and 's
second, open only the relevant files, and write a heading in each of them
third, perform loop over l. For each multipole, get all 's by calling harmonic_cl_at_l() and distribute the results to relevant files
finally, close files and free arrays of files and 's

◆ output_pk()

int output_pk	(	struct background *	pba,
		struct perturbations *	ppt,
		struct fourier *	pfo,
		struct output *	pop,
		enum pk_outputs	pk_output
	)

This routines writes the output in files for Fourier matter power spectra P(k)'s (linear or non-linear)

Parameters

pba	Input: pointer to background structure (needed for calling harmonic_pk_at_z())
ppt	Input: pointer perturbation structure
pfo	Input: pointer to fourier structure
pop	Input: pointer to output structure
pk_output	Input: pk_linear or pk_nonlinear

Summary:

define local variables
preliminary: check whether we need to output the decomposition into contributions from each initial condition
allocate arrays to store the P(k)
allocate pointer to output files
loop over pk type (_cb, _m)
loop over z
first, check that requested redshift z_pk is consistent
second, open only the relevant files and write a header in each of them
third, compute P(k) for each k
fourth, write in files
fifth, close files

◆ output_tk()

int output_tk	(	struct background *	pba,
		struct perturbations *	ppt,
		struct output *	pop
	)

This routines writes the output in files for matter transfer functions 's.

Parameters

pba	Input: pointer to background structure (needed for calling harmonic_pk_at_z())
ppt	Input: pointer perturbation structure
pop	Input: pointer to output structure

Summary:

define local variables
first, check that requested redshift z_pk is consistent
second, open only the relevant files, and write a heading in each of them
free memory and close files

◆ output_heating()

int output_heating	(	struct injection *	pin,
		struct noninjection *	pni,
		struct output *	pop
	)

Local variables

◆ output_distortions()

int output_distortions	(	struct distortions *	psd,
		struct output *	pop
	)

Local variables

◆ output_print_data()

int output_print_data	(	FILE *	out,
		char	titles[_MAXTITLESTRINGLENGTH_],
		double *	dataptr,
		int	size_dataptr
	)

Summary

First we print the titles
Then we print the data

◆ output_open_cl_file()

int output_open_cl_file	(	struct harmonic *	phr,
		struct output *	pop,
		FILE **	clfile,
		FileName	filename,
		char *	first_line,
		int	lmax
	)

This routine opens one file where some 's will be written, and writes a heading with some general information concerning its content.

Parameters

phr	Input: pointer to harmonic structure
pop	Input: pointer to output structure
clfile	Output: returned pointer to file pointer
filename	Input: name of the file
first_line	Input: text describing the content (mode, initial condition..)
lmax	Input: last multipole in the file (the first one is assumed to be 2)

Returns: the error status

Summary

First we deal with the entries that are dependent of format type
Next deal with entries that are independent of format type

◆ output_one_line_of_cl()

int output_one_line_of_cl	(	struct background *	pba,
		struct harmonic *	phr,
		struct output *	pop,
		FILE *	clfile,
		double	l,
		double *	cl,
		int	ct_size
	)

This routine write one line with l and all 's for all types (TT, TE...)

Parameters

pba	Input: pointer to background structure (needed for $T_{cmb}$ )
phr	Input: pointer to harmonic structure
pop	Input: pointer to output structure
clfile	Input: file pointer
l	Input: multipole
cl	Input: 's for all types
ct_size	Input: number of types

Returns: the error status

◆ output_open_pk_file()

int output_open_pk_file	(	struct background *	pba,
		struct fourier *	pfo,
		struct output *	pop,
		FILE **	pkfile,
		FileName	filename,
		char *	first_line,
		double	z
	)

This routine opens one file where some P(k)'s will be written, and writes a heading with some general information concerning its content.

Parameters

pba	Input: pointer to background structure (needed for h)
pfo	Input: pointer to fourier structure
pop	Input: pointer to output structure
pkfile	Output: returned pointer to file pointer
filename	Input: name of the file
first_line	Input: text describing the content (initial conditions, ...)
z	Input: redshift of the output

Returns: the error status

◆ output_one_line_of_pk()

int output_one_line_of_pk	(	FILE *	pkfile,
		double	one_k,
		double	one_pk
	)

This routine writes one line with k and P(k)

Parameters

pkfile	Input: file pointer
one_k	Input: wavenumber
one_pk	Input: matter power spectrum

Returns: the error status

Functions

Detailed Description

Function Documentation

◆ output_init()

◆ output_cl()

◆ output_pk()

◆ output_tk()

◆ output_heating()

◆ output_distortions()

◆ output_print_data()

◆ output_open_cl_file()

◆ output_one_line_of_cl()

◆ output_open_pk_file()

◆ output_one_line_of_pk()