math_funcs.h File Reference

#include "matpackI.h"
#include "matpackIII.h"

Go to the source code of this file.

Functions
Numeric	fac (const Index n)
	fac More...
Index	integer_div (const Index &x, const Index &y)
	integer_div More...
Numeric	last (ConstVectorView x)
	last More...
Index	last (const ArrayOfIndex &x)
	last More...
void	linspace (Vector &x, const Numeric start, const Numeric stop, const Numeric step)
	linspace More...
void	nlinspace (Vector &x, const Numeric start, const Numeric stop, const Index n)
	nlinspace More...
void	nlinspace (VectorView x, const Numeric start, const Numeric stop, const Index n)
void	nlogspace (Vector &x, const Numeric start, const Numeric stop, const Index n)
	nlogspace More...
Numeric	AngIntegrate_trapezoid (ConstMatrixView Integrand, ConstVectorView za_grid, ConstVectorView aa_grid)
	AngIntegrate_trapezoid. More...
Numeric	AngIntegrate_trapezoid (ConstVectorView Integrand, ConstVectorView za_grid)
	AngIntegrate_trapezoid. More...
Numeric	AngIntegrate_trapezoid_opti (ConstMatrixView Integrand, ConstVectorView za_grid, ConstVectorView aa_grid, ConstVectorView grid_stepsize)
	AngIntegrate_trapezoid_opti. More...
Numeric	sign (const Numeric &x)
	sign More...
void	mgd (VectorView psd, const Vector &x, const Numeric &n0, const Numeric &mu, const Numeric &la, const Numeric &ga)
void	mgd_with_derivatives (VectorView psd, MatrixView jac_data, const Vector &x, const Numeric &n0, const Numeric &mu, const Numeric &la, const Numeric &ga, const bool &do_n0_jac, const bool &do_mu_jac, const bool &do_la_jac, const bool &do_ga_jac)
void	delanoe_shape_with_derivative (VectorView psd, MatrixView jac_data, const Vector &x, const Numeric &alpha, const Numeric &beta)
	! Shape functions for normalized PSD. More...
Numeric	mod_gamma_dist (Numeric x, Numeric N0, Numeric Lambda, Numeric mu, Numeric gamma)
	Generalized Modified Gamma Distribution. More...
void	unitl (Vector &x)
	unitl More...
void	flat (VectorView x, ConstMatrixView X)
	flat More...
void	flat (VectorView x, ConstTensor3View X)
	flat More...
void	reshape (MatrixView X, ConstVectorView x)
	reshape More...
void	reshape (Tensor3View X, ConstVectorView x)
	reshape More...
void	calculate_weights_linear (Vector &x, Vector &w, const Index nph)
	calculate_weights_linear More...

Detailed Description

Author

Patrick Eriksson Patri.nosp@m.ck.E.nosp@m.rikss.nosp@m.on@c.nosp@m.halme.nosp@m.rs.s.nosp@m.e

Date

2000-09-18

Contains the decleration of the functions in math_funcs.cc.

Definition in file math_funcs.h.

Function Documentation

AngIntegrate_trapezoid() [1/2]

Numeric AngIntegrate_trapezoid	(	ConstMatrixView	Integrand,
		ConstVectorView	za_grid,
		ConstVectorView	aa_grid
	)

AngIntegrate_trapezoid.

Performs an integration of a matrix over all directions defined in angular grids using the trapezoidal integration method.

Parameters

Integrand	The Matrix to be integrated
za_grid	The zenith angle grid
aa_grid	The azimuth angle grid

Returns

The resulting integral

Definition at line 296 of file math_funcs.cc.

References ARTS::Var::aa_grid(), DEG2RAD, is_size(), and ARTS::Var::za_grid().

Referenced by AngIntegrate_trapezoid_opti(), doit_scat_fieldCalc(), doit_scat_fieldCalcLimb(), doit_scat_fieldNormalize(), get_stepwise_scattersky_source(), scat_dataCheck(), and scat_dataReduceT().

AngIntegrate_trapezoid() [2/2]

Numeric AngIntegrate_trapezoid	(	ConstVectorView	Integrand,
		ConstVectorView	za_grid
	)

AngIntegrate_trapezoid.

Performs an integration of a matrix over all directions defined in angular grids using the trapezoidal integration method. The integrand is independant of the azimuth angle. The integration over the azimuth angle gives a 2*PI

Parameters

Integrand	Input : The vector to be integrated
za_grid	Input : The zenith angle grid

Author

Claas Teichmann

Date

2003-05-13

Returns

The resulting integral

Definition at line 393 of file math_funcs.cc.

References DEG2RAD, is_size(), PI, and ARTS::Var::za_grid().

AngIntegrate_trapezoid_opti()

Numeric AngIntegrate_trapezoid_opti	(	ConstMatrixView	Integrand,
		ConstVectorView	za_grid,
		ConstVectorView	aa_grid,
		ConstVectorView	grid_stepsize
	)

AngIntegrate_trapezoid_opti.

Performs an integration of a matrix over all directions defined in angular grids using the trapezoidal integration method.

In addition to the "old fashined" integration method, it checks whether the stepsize is constant. If it is, it uses a faster method, if not, it uses the old one.

Parameters

Integrand	Input : The Matrix to be integrated
za_grid	Input : The zenith angle grid
aa_grid	Input : The azimuth angle grid
grid_stepsize	Input : stepsize of the grid

Returns

The resulting integral

Author

Claas Teichmann claas.nosp@m.@sat.nosp@m..phys.nosp@m.ik.u.nosp@m.ni-br.nosp@m.emen.nosp@m..de

Date

2003/05/28

Definition at line 340 of file math_funcs.cc.

References ARTS::Var::aa_grid(), AngIntegrate_trapezoid(), DEG2RAD, is_size(), temp, and ARTS::Var::za_grid().

Referenced by doit_scat_fieldCalc(), doit_scat_fieldCalcLimb(), and test_AngIntegrate_trapezoid_opti().

calculate_weights_linear()

void calculate_weights_linear	(	Vector &	x,
		Vector &	w,
		const Index	nph
	)

calculate_weights_linear

Function to set the evaluation points and the corresponding weights for numerical integration on the domain from [-1,1]

Parameters:

Parameters

[out]	x	evaluation points \paran[out] w integration weights
[in]	nph	number of evaluation points per hemisphere

Definition at line 813 of file math_funcs.cc.

References N, nlinspace(), w(), and ARTS::Var::x().

Referenced by AngularGridsSetFluxCalc().

delanoe_shape_with_derivative()

void delanoe_shape_with_derivative	(	VectorView	psd,
		MatrixView	jac_data,
		const Vector &	x,
		const Numeric &	alpha,
		const Numeric &	beta
	)

! Shape functions for normalized PSD.

This function implements the shape function F(X, alpha, beta) from as proposed by Delanoe et al. in "Normalized particle size distribution for remote sensing application".

Parameters

[OUT]	psd On return contains the values of F corresponding to the values in x.
[OUT]	jac_data On return contains the first derivative of F w.r.t x evaluated at the values in x.
[IN]	x The values at which to evaluate the shape functions and and derivatives.
[IN]	alpha The alpha parameter of the shape function.
[IN]	beta The beta parameter of the shape function.

Definition at line 619 of file math_funcs.cc.

fac()

Numeric fac

(

const Index

n

)

fac

Calculates the factorial.

The function asserts that n must be >= 0

Returns

The factorial

Parameters

n

Nominator

Author

Oliver Lemke

Date

2003-08-15

Definition at line 63 of file math_funcs.cc.

Referenced by Linefunctions::apply_VVW_scaling(), g_legendre_poly_norm_schmidt(), g_legendre_poly_norm_schmidt_deriv(), g_legendre_poly_norm_schmidt_deriv1(), g_legendre_poly_norm_schmidt_deriv2(), g_legendre_poly_norm_schmidt_deriv3(), g_legendre_poly_norm_schmidt_deriv4(), legendre_poly_norm_schmidt(), legendre_poly_norm_schmidt_deriv(), Absorption::PredefinedModel::makarov2020_o2_lines_mpm(), MCRadar(), MPM93_N2_continuum(), Linefunctions::set_voigt(), specular_losCalc(), and ze_cfac().

flat() [1/2]

void flat	(	VectorView	x,
		ConstMatrixView	X
	)

flat

Flattens a matrix to a vector

The matrix is read from front, i.e. rows are looped first. In Matlab this equals x=X(:).

Parameters

[out]	x	The vector. Should already be sized
[in]	X	The matrix.

Author

Patrick Eriksson

Date

2015-09-09

Definition at line 707 of file math_funcs.cc.

References ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().

Referenced by xaStandard().

flat() [2/2]

void flat	(	VectorView	x,
		ConstTensor3View	X
	)

flat

Converts Tensor3 to a vector

The matrix is read from front, i.e. pages are looped first, followed by rows. In Matlab this equals x=X(:).

Parameters

[out]	x	The vector. Should already be sized
[in]	X	The tensor.

Author

Patrick Eriksson

Date

2015-09-09

Definition at line 733 of file math_funcs.cc.

References ConstTensor3View::ncols(), ConstTensor3View::npages(), ConstTensor3View::nrows(), and ARTS::Var::x().

integer_div()

Index integer_div	(	const Index &	x,
		const Index &	y
	)

integer_div

Performs an integer division.

The function asserts that the reminder of the division x/y is 0.

Returns

The quotient

Parameters

x	Nominator
y	Denominator

Author

Patrick Eriksson

Date

2002-08-11

Definition at line 87 of file math_funcs.cc.

References is_multiple(), ARTS::Var::x(), and ARTS::Var::y().

last() [1/2]

Index last

(

const ArrayOfIndex &

x

)

last

Returns the last value of an index array.

Returns

The last value of x.

Parameters

x	An index array.

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 180 of file math_funcs.cc.

References ARTS::Var::x().

last() [2/2]

Numeric last

(

ConstVectorView

x

)

last

Returns the last value of a vector.

Returns

The last value of x.

Parameters

x

A vector.

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 165 of file math_funcs.cc.

References ARTS::Var::x().

Referenced by chk_scat_data(), find_effective_channel_boundaries(), mixer_matrix(), pos2refell_r(), scat_data_singleTmatrix(), ScatSpeciesExtendTemperature(), sensor_responseAntenna(), sensor_responseBackend(), sensor_responseMixer(), and ySimpleSpectrometer().

linspace()

void linspace	(	Vector &	x,
		const Numeric	start,
		const Numeric	stop,
		const Numeric	step
	)

linspace

Linearly spaced vector with specified spacing.

The first element of x is always start. The next value is start+step etc. Note that the last value can deviate from stop. The step can be both positive and negative. (in Matlab notation: start:step:stop)

Size of result is adjusted within this function!

Parameters

x	Output: linearly spaced vector
start	first value in x
stop	last value of x <= stop
step	distance between values in x

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 204 of file math_funcs.cc.

References Zeeman::start(), and ARTS::Var::x().

mgd()

void mgd	(	VectorView	psd,
		const Vector &	x,
		const Numeric &	n0,
		const Numeric &	mu,
		const Numeric &	la,
		const Numeric &	ga
	)

Modified gamma distribution

Uses all four free parameters (n0, mu, la, ga) to calculate psd(D) = n0 * D^mu * exp( -la * x^ga )

Reference: Eq 1 of Petty & Huang, JAS, (2011).

Parameters

psd	Particle number density per x-interval. Sizing of vector should be done before calling the function.
x	Mass or size.
n0	See above.
mu	See above.
la	See above.
ga	See above.

Author

Jana Mendrok, Patrick Eriksson

Date

2017-06-07

Definition at line 451 of file math_funcs.cc.

References ConstVectorView::nelem(), pow(), and ARTS::Var::x().

mgd_with_derivatives()

void mgd_with_derivatives	(	VectorView	psd,
		MatrixView	jac_data,
		const Vector &	x,
		const Numeric &	n0,
		const Numeric &	mu,
		const Numeric &	la,
		const Numeric &	ga,
		const bool &	do_n0_jac,
		const bool &	do_mu_jac,
		const bool &	do_la_jac,
		const bool &	do_ga_jac
	)

Modified gamma distribution, and derivatives

As mgd, but this version can also return the derivate of psd with respect to the four parameters.

Parameters

psd	Particle number density per x-interval. Sizing of vector should be done before calling the function.
jac_data	Container for returning jacobian data. Shall be a matrix with four rows, where the rows match n0, mu, la and ga. Number of columns same as length of psd.
x	Mass or size.
n0	See above.
mu	See above.
la	See above.
ga	See above.
do_n0_jac	Flag to actually calculate d_psd/d_n0
do_mu_jac	Flag to actually calculate d_psd/d_mu
do_la_jac	Flag to actually calculate d_psd/d_la
do_ga_jac	Flag to actually calculate d_psd/d_ga

Author

Patrick Eriksson

Date

2017-06-07

Definition at line 530 of file math_funcs.cc.

References ConstMatrixView::ncols(), ConstVectorView::nelem(), ConstMatrixView::nrows(), pow(), and ARTS::Var::x().

Referenced by psd_mgd_mass_and_something(), psdModifiedGamma(), and psdModifiedGammaMass().

mod_gamma_dist()

Numeric mod_gamma_dist	(	Numeric	x,
		Numeric	N0,
		Numeric	Lambda,
		Numeric	mu,
		Numeric	gamma
	)

Generalized Modified Gamma Distribution.

Returns number density per unit of 'x' as function of 'x'.

Returns

dN Number density as function of x.

Parameters

x	Numeric
N0	Numeric, Scaling parameter
Lambda	Numeric, Shape parameter
mu	Numeric, Shape parameter
gamma	Numeric, Shape parameter

Author

Manfred Brath

Date

2015-01-19

Definition at line 653 of file math_funcs.cc.

References pow(), and ARTS::Var::x().

nlinspace() [1/2]

void nlinspace	(	Vector &	x,
		const Numeric	start,
		const Numeric	stop,
		const Index	n
	)

nlinspace

Linearly spaced vector with specified length.

Returns a vector equally and linearly spaced between start and stop of length n. (equals the Matlab function linspace)

The length must be > 1.

Parameters

x	Output: linearly spaced vector
start	first value in x
stop	last value of x <= stop
n	length of x

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 231 of file math_funcs.cc.

References Zeeman::start(), and ARTS::Var::x().

Referenced by AngularGridsSetFluxCalc(), AntennaConstantGaussian1D(), calc_ssp_fixed_test(), calc_ssp_random_test(), calculate_weights_linear(), DOAngularGridsSet(), doit_scat_fieldCalcLimb(), DoitScatteringDataPrepare(), f_gridFromSensorAMSU(), f_gridFromSensorHIRS(), fillSpeciesAuxDataWithPartitionFunctionsFromSpeciesData(), FouComp_1ScatElem(), gaussian_response_autogrid(), get_angs(), iyHybrid(), iyHybrid2(), pha_mat_sptFromDataDOITOpt(), pha_mat_sptFromMonoData(), pnd_fieldZero(), sensor_responseFillFgrid(), test46(), test_ecs20(), test_mpm20(), VectorNLinSpace(), and xml_read_from_stream().

nlinspace() [2/2]

void nlinspace	(	VectorView	x,
		const Numeric	start,
		const Numeric	stop,
		const Index	n
	)

Definition at line 241 of file math_funcs.cc.

References Zeeman::start(), and ARTS::Var::x().

nlogspace()

void nlogspace	(	Vector &	x,
		const Numeric	start,
		const Numeric	stop,
		const Index	n
	)

nlogspace

Logarithmically spaced vector with specified length.

Returns a vector logarithmically spaced vector between start and stop of length n (equals the Matlab function logspace)

The length must be > 1.

Parameters

x	Output: logarithmically spaced vector
start	first value in x
stop	last value of x <= stop
n	length of x

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 267 of file math_funcs.cc.

References Zeeman::start(), and ARTS::Var::x().

Referenced by VectorNLogSpace().

reshape() [1/2]

void reshape	(	MatrixView	X,
		ConstVectorView	x
	)

reshape

Converts vector to Matrix

The matrix is filled from front, i.e. rows are looped first, followed by cols. In Matlab this equals X = reshape( x, [ X.nrows(), X.ncols() ]

Parameters

[out]	X	The matrix. Should already be sized
[in]	x	The vector.

Author

Patrick Eriksson

Date

2015-09-10

Definition at line 789 of file math_funcs.cc.

References ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().

reshape() [2/2]

void reshape	(	Tensor3View	X,
		ConstVectorView	x
	)

reshape

Converts vector to Tensor3

The tensor is filled from front, i.e. pages are looped first, followed by rows. In Matlab this equals X = reshape( x, [ X.npages(), X.nrows(), X.ncols() ]

Parameters

[out]	X	The tensor. Should already be sized
[in]	x	The vector.

Author

Patrick Eriksson

Date

2015-09-10

Definition at line 761 of file math_funcs.cc.

References ConstTensor3View::ncols(), ConstTensor3View::npages(), ConstTensor3View::nrows(), and ARTS::Var::x().

Referenced by AtmFieldPerturb(), and x2artsAtmAndSurf().

sign()

Numeric sign

(

const Numeric &

x

)

sign

Returns the sign of a numeric value.

The function returns 1 if the value is greater than zero, 0 if it equals zero and -1 if it is less than zero.

Returns

The sign of x (see above).

Parameters

x	A Numeric.

Author

Patrick Eriksson

Date

2000-06-27

Definition at line 423 of file math_funcs.cc.

References ARTS::Var::x().

Referenced by error_if_limb_ppath(), binostream::float2ieee_double(), binostream::float2ieee_single(), geompath_from_r1_to_r2(), binistream::ieee_double2float(), binistream::ieee_single2float(), poslos2cart(), r_crossing_2d(), raytrace_2d_linear_basic(), specular_losCalc(), specular_losCalcNoTopography(), and surfaceSemiSpecularBy3beams().

unitl()

void unitl

(

Vector &

x

)

unitl

Normalises a vector to have unit length.

The standard Euclidean norm is used (2-norm).

param x In/Out: A vector.

Author

Patrick Eriksson

Date

2012-02-12

Definition at line 687 of file math_funcs.cc.

References sqrt(), and ARTS::Var::x().