#include <cstring>
#include "matpackI.h"
#include "exceptions.h"

Functions
ostream &	operator<< (ostream &os, const ConstVectorView &v)
	Output operator. More...

void	copy (ConstIterator1D origin, const ConstIterator1D &end, Iterator1D target)
	Copy data between begin and end to target. More...

void	copy (Numeric x, Iterator1D target, const Iterator1D &end)
	Copy a scalar to all elements. More...

ostream &	operator<< (ostream &os, const ConstMatrixView &v)
	Output operator. More...

void	copy (ConstIterator2D origin, const ConstIterator2D &end, Iterator2D target)
	Copy data between begin and end to target. More...

void	copy (Numeric x, Iterator2D target, const Iterator2D &end)
	Copy a scalar to all elements. More...

Numeric	operator* (const ConstVectorView &a, const ConstVectorView &b)
	Scalar product. More...

void	mult (VectorView y, const ConstMatrixView &M, const ConstVectorView &x)
	Matrix Vector multiplication. More...

void	mult (MatrixView A, const ConstMatrixView &B, const ConstMatrixView &C)
	Matrix multiplication. More...

ConstMatrixView	transpose (ConstMatrixView m)
	Const version of transpose. More...

MatrixView	transpose (MatrixView m)
	Returns the transpose. More...

MatrixView	transpose (Vector v)
	Returns the transpose. More...

void	transform (VectorView y, double(&my_func)(double), ConstVectorView x)
	A generic transform function for vectors, which can be used to implement mathematical functions operating on all elements. More...

void	transform (MatrixView y, double(&my_func)(double), ConstMatrixView x)
	A generic transform function for matrices, which can be used to implement mathematical functions operating on all elements. More...

Numeric	max (const ConstVectorView &x)
	Max function, vector version. More...

Numeric	max (const ConstMatrixView &x)
	Max function, matrix version. More...

Numeric	min (const ConstVectorView &x)
	Min function, vector version. More...

Numeric	min (const ConstMatrixView &x)
	Min function, matrix version. More...

Numeric	mean (const ConstVectorView &x)
	Mean function, vector version. More...

Numeric	mean (const ConstMatrixView &x)
	Mean function, matrix version. More...

Numeric	debug_matrixview_get_elem (MatrixView &mv, Index r, Index c)
	Helper function to access matrix elements. More...

Variables
const Joker	joker = Joker()

Detailed Description

Author: Stefan Buehler

Date: 2001-09-15

Definition in file matpackI.cc.

Function Documentation

◆ copy() [1/4]

void copy	(	ConstIterator1D	origin,
		const ConstIterator1D &	end,
		Iterator1D	target
	)

Copy data between begin and end to target.

Target must be a valid area of memory. Note that the strides in the iterators can be different, so that we can for example copy data between different kinds of subvectors.

Definition at line 597 of file matpackI.cc.

References Iterator1D::mstride, ConstIterator1D::mstride, Iterator1D::mx, and ConstIterator1D::mx.

Referenced by abs(), Sparse::insert_row(), Matrix::Matrix(), Array< base >::operator=(), MatrixView::operator=(), VectorView::operator=(), Matrix::operator=(), Vector::operator=(), and Vector::Vector().

◆ copy() [2/4]

void copy	(	ConstIterator2D	origin,
		const ConstIterator2D &	end,
		Iterator2D	target
	)

Copy data between begin and end to target.

Target must be a valid area of memory. Note that the strides in the iterators can be different, so that we can for example copy data between different kinds of subvectors.

Origin, end, and target are 2D iterators, marking rows in a matrix. For each row the 1D iterator is obtained and used to copy the elements.

Definition at line 1392 of file matpackI.cc.

References ConstVectorView::begin(), VectorView::begin(), and ConstVectorView::end().

◆ copy() [3/4]

void copy	(	Numeric	x,
		Iterator1D	target,
		const Iterator1D &	end
	)

Copy a scalar to all elements.

Definition at line 611 of file matpackI.cc.

◆ copy() [4/4]

void copy	(	Numeric	x,
		Iterator2D	target,
		const Iterator2D &	end
	)

Copy a scalar to all elements.

Definition at line 1407 of file matpackI.cc.

References VectorView::begin(), and VectorView::end().

◆ debug_matrixview_get_elem()

Numeric debug_matrixview_get_elem	(	MatrixView &	mv,
		Index	r,
		Index	c
	)

Helper function to access matrix elements.

Because of function inlining the operator() is not accessible from the debuggger. This function helps to access Matrix elements from within the debugger.

Parameters

mv	MatrixView
r	Row index
c	Column index

Author: Oliver Lemke

Date: 2004-05-10

Definition at line 1944 of file matpackI.cc.

◆ max() [1/2]

Numeric max ( const ConstMatrixView & x )

Max function, matrix version.

Definition at line 1800 of file matpackI.cc.

References ConstVectorView::begin(), ConstMatrixView::begin(), ConstVectorView::end(), ConstMatrixView::end(), and max().

◆ max() [2/2]

Numeric max ( const ConstVectorView & x )

Max function, vector version.

Definition at line 1782 of file matpackI.cc.

References ConstVectorView::begin(), ConstVectorView::end(), and max().

Referenced by max().

◆ mean() [1/2]

Numeric mean ( const ConstMatrixView & x )

Mean function, matrix version.

Definition at line 1878 of file matpackI.cc.

References ConstVectorView::begin(), ConstMatrixView::begin(), ConstVectorView::end(), ConstMatrixView::end(), mean(), ConstMatrixView::ncols(), and ConstMatrixView::nrows().

◆ mean() [2/2]

Numeric mean ( const ConstVectorView & x )

Mean function, vector version.

Definition at line 1862 of file matpackI.cc.

References ConstVectorView::begin(), ConstVectorView::end(), mean(), and ConstVectorView::nelem().

Referenced by abs_lookupSetup(), abs_lookupTestAccMC(), mean(), and polynomial_basis_func().

◆ min() [1/2]

Numeric min ( const ConstMatrixView & x )

Min function, matrix version.

Definition at line 1840 of file matpackI.cc.

References ConstVectorView::begin(), ConstMatrixView::begin(), ConstVectorView::end(), ConstMatrixView::end(), and min().

◆ min() [2/2]

Numeric min ( const ConstVectorView & x )

Min function, vector version.

Definition at line 1822 of file matpackI.cc.

References ConstVectorView::begin(), ConstVectorView::end(), and min().

Referenced by min().

◆ mult() [1/2]

void mult	(	MatrixView	A,
		const ConstMatrixView &	B,
		const ConstMatrixView &	C
	)

Matrix multiplication.

A = B*C. Note that the order is different from MTL, output comes first! Dimensions of A, B, and C must match. No memory reallocation takes place, only the data is copied. Using this function on overlapping MatrixViews belonging to the same Matrix will lead to unpredictable results. In particular, this means that A and B must not be the same matrix!

Definition at line 1652 of file matpackI.cc.

References VectorView::begin(), ConstMatrixView::begin(), MatrixView::begin(), VectorView::end(), MatrixView::end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and transpose().

◆ mult() [2/2]

void mult	(	VectorView	y,
		const ConstMatrixView &	M,
		const ConstVectorView &	x
	)

Matrix Vector multiplication.

y = M*x. Note that the order is different from MTL, output comes first! Dimensions of y, M, and x must match. No memory reallocation takes place, only the data is copied. Using this function on overlapping MatrixViews belonging to the same Matrix will lead to unpredictable results. In particular, this means that A and B must not be the same matrix! The implementation here is different from the other multiplication routines. It does not use iterators but a more drastic approach to gain maximum performance.

Definition at line 1607 of file matpackI.cc.

References M, ConstVectorView::mdata, Range::mextent, ConstVectorView::mrange, Range::mstart, and Range::mstride.

Referenced by cloud_RT_surface(), get_ppath_cloudrtvars(), interpTArray(), iy_transmission_mult(), iyBeerLambertStandardClearsky(), iyBeerLambertStandardCloudbox(), iyEmissionStandardClearsky(), iyFOS(), jacobianCalcAbsSpeciesPerturbations(), jacobianCalcFreqShiftAndStretchInterp(), jacobianCalcPointingZaInterp(), jacobianCalcPointingZaRecalc(), jacobianCalcTemperaturePerturbations(), map_daa(), matrix_exp(), MatrixMatrixMultiply(), MCGeneral(), MCIPA(), mcPathTraceGeneral(), mcPathTraceIPA(), rte_step_std(), sensor_responseAntenna(), sensor_responseBackend(), sensor_responseBeamSwitching(), sensor_responseFillFgrid(), sensor_responseFrequencySwitching(), sensor_responseMixer(), sensor_responsePolarisation(), sensor_responseWMRF(), SparseSparseMultiply(), surface_calc(), test32(), test4(), test40(), test42(), test45(), test5(), test_lusolve4D(), trans_step_std(), VectorMatrixMultiply(), yCalc(), and yCalc2().

◆ operator*()

Numeric operator*	(	const ConstVectorView &	a,
		const ConstVectorView &	b
	)

Scalar product.

The two vectors may be identical.

Definition at line 1582 of file matpackI.cc.

References ConstVectorView::begin(), ConstVectorView::end(), and ConstVectorView::nelem().

◆ operator<<() [1/2]

ostream& operator<<	(	ostream &	os,
		const ConstMatrixView &	v
	)

Output operator.

This demonstrates how iterators can be used to traverse the matrix. The iterators know which part of the matrix is ‘active’, and also the strides in both directions. This function is a bit more complicated than necessary to illustrate the concept, because the formating should look nice. This means that the first row, and the first element in each row, have to be treated individually.

Definition at line 914 of file matpackI.cc.

References ConstVectorView::begin(), ConstMatrixView::begin(), ConstVectorView::end(), and ConstMatrixView::end().

◆ operator<<() [2/2]

ostream& operator<<	(	ostream &	os,
		const ConstVectorView &	v
	)

Output operator.

This demonstrates how iterators can be used to traverse the vector. The iterators know which part of the vector is ‘active’, and also the stride.

Definition at line 275 of file matpackI.cc.

References ConstVectorView::begin(), and ConstVectorView::end().

◆ transform() [1/2]

void transform	(	MatrixView	y,
		double(&)(double)	my_func,
		ConstMatrixView	x
	)

A generic transform function for matrices, which can be used to implement mathematical functions operating on all elements.

Because we have this, we don't need explicit functions like sqrt for matrices! The type of the mathematical function is double (&my_func)(double). Numeric would not work here, since mathematical functions for float do not exist!

transform(y,sin,x) computes y = sin(x)

This function can also be used for Vectors, because there is a conversion to MatrixView.

The two Matrix views may be the same one, in which case the conversion happens in place.

Parameters

y	Output: The results of the function acting on each element of x.
my_func	A function (e.g., sqrt).
x	A matrix.

Definition at line 1760 of file matpackI.cc.

References ConstVectorView::begin(), VectorView::begin(), ConstMatrixView::begin(), MatrixView::begin(), ConstVectorView::end(), ConstMatrixView::end(), ConstMatrixView::ncols(), and ConstMatrixView::nrows().

◆ transform() [2/2]

void transform	(	VectorView	y,
		double(&)(double)	my_func,
		ConstVectorView	x
	)

A generic transform function for vectors, which can be used to implement mathematical functions operating on all elements.

Because we have this, we don't need explicit functions like sqrt for matrices! The type of the mathematical function is double (&my_func)(double). Numeric would not work here, since mathematical functions for float do not exist!

transform(y,sin,x) computes y = sin(x)

Although the matrix version of this can also be used for vectors, thanks to the automatic interpretation of a vector as a one column matrix, this one is slightly more efficient. However, the difference is very small (only a few percent).

The two views may be the same one, in which case the conversion happens in place.

Parameters

y	Output: The results of the function acting on each element of x.
my_func	A function (e.g., sqrt).
x	A vector.

Definition at line 1728 of file matpackI.cc.

References ConstVectorView::begin(), VectorView::begin(), ConstVectorView::end(), and ConstVectorView::nelem().

Referenced by abs_lookupSetup(), abs_lookupSetupBatch(), abs_lookupSetupWide(), GasAbsLookup::Adapt(), atm_fields_compactFromMatrix(), atm_fields_compactFromMatrixChevalAll(), AtmFieldsRefinePgrid(), itw2p(), p2gridpos(), p2gridpos_poly(), test2(), test31(), test6(), test7(), and VectorLogSpace().

◆ transpose() [1/3]

ConstMatrixView transpose ( ConstMatrixView m )

Const version of transpose.

Definition at line 1689 of file matpackI.cc.

References ConstMatrixView::mcr, ConstMatrixView::mdata, and ConstMatrixView::mrr.

Referenced by abs_fieldCalc(), atm_fields_compactFromMatrix(), atm_fields_compactFromMatrixChevalAll(), iyBeerLambertStandardClearsky(), iyEmissionStandardClearsky(), mult(), test1(), test33(), test41(), test42(), test46(), test5(), transpose(), WMRFSelectChannels(), and y_unitApply().

◆ transpose() [2/3]

MatrixView transpose ( MatrixView m )

Returns the transpose.

This creates a special MatrixView for the transpose. The original is not changed!

Definition at line 1696 of file matpackI.cc.

References ConstMatrixView::mcr, ConstMatrixView::mdata, and ConstMatrixView::mrr.

◆ transpose() [3/3]

MatrixView transpose ( Vector v )

Returns the transpose.

This creates a special MatrixView for the transpose. The original is not changed!

Definition at line 1703 of file matpackI.cc.

References transpose().

Variable Documentation

◆ joker

const Joker joker = Joker()

extern

Referenced by abs_coefCalcSaveMemory(), abs_fieldCalc(), abs_h2oSet(), abs_lookupSetup(), abs_lookupSetupBatch(), abs_lookupSetupWide(), abs_lookupTestAccMC(), abs_lookupTestAccuracy(), abs_n2Set(), abs_scalar_gasExtractFromLookup(), abs_scalar_gasFromAbsCoef(), abs_vecAddGas(), abs_vecAddPart(), abs_vecTransform(), abs_xsec_per_speciesAddConts(), abs_xsec_per_speciesAddLines(), AbsInputFromAtmFields(), GasAbsLookup::Adapt(), antenna1d_matrix(), antenna_responseGaussian(), AntennaMultiBeamsToPencilBeams(), atm_fields_compactAddConstant(), atm_fields_compactAddSpecies(), atm_fields_compactExpand(), atm_fields_compactFromMatrix(), atm_fields_compactFromMatrixChevalAll(), AtmFieldsCalc(), AtmFieldsFromCompact(), AtmFieldsFromCompactChevalAll(), AtmFieldsRefinePgrid(), basics_checkedCalc(), calc_lookup_error(), clear_rt_vars_at_gp(), cloud_atm_vars_by_gp(), cloud_fieldsCalc(), cloud_ppath_update1D_planeparallel(), cloud_ppath_update3D(), cloud_RT_no_background(), cloud_RT_surface(), CloudboxGetIncoming(), CloudboxGetIncoming1DAtm(), cloudboxSetAutomatically(), cloudy_rt_vars_at_gp(), diy_from_path_to_rgrids(), doit_i_fieldSetClearsky(), doit_i_fieldUpdateSeq3D(), doit_scat_fieldCalc(), doit_scat_fieldCalcLimb(), DoitCloudboxFieldPut(), DoitScatteringDataPrepare(), dtauc_ssalbCalc(), ext_matAddGas(), ext_matAddPart(), ext_matTransform(), Extract(), GasAbsLookup::Extract(), findZ11max(), fos_yStandard(), get_iy_of_background(), get_ppath_atmvars(), get_ppath_cloudrtvars(), get_ppath_pnd(), get_ppath_rtvars(), get_refr_index_1d(), get_refr_index_2d(), get_refr_index_3d(), interp(), interp_atmfield_by_itw(), interp_atmsurface_by_itw(), interp_cloud_coeff1D(), interp_gfield3(), interp_scat_angle_temperature(), interpolate_scat_angle(), interpolate_scat_angleDOIT(), interpTArray(), iwp_cloud_opt_pathCalc(), iy_interp_cloudbox_field(), iy_transmission_mult(), iy_transmission_mult_scalar_tau(), iyb_calc(), iyBeerLambertStandardClearsky(), iyBeerLambertStandardCloudbox(), iyEmissionStandardClearsky(), iyFOS(), iyMC(), jacobianCalcAbsSpeciesPerturbations(), jacobianCalcPointingZaRecalc(), jacobianCalcTemperaturePerturbations(), Matrix1ColFromVector(), Matrix1RowFromVector(), Matrix2ColFromVectors(), Matrix2RowFromVectors(), Matrix3ColFromVectors(), Matrix3RowFromVectors(), matrix_exp(), mc_IWP_cloud_opt_pathCalc(), MCGeneral(), MCIPA(), mcPathTraceGeneral(), mcPathTraceIPA(), opt_prop_sptFromData(), opt_prop_sptFromMonoData(), opt_propExtract(), pha_mat_singleExtract(), pha_mat_sptFromData(), pha_mat_sptFromDataDOITOpt(), pha_mat_sptFromMonoData(), pha_matTransform(), pmomCalc(), pnd_fieldCalc(), pnd_fieldSetup(), ppath_calc(), ppath_copy(), ppath_start_stepping(), ppath_stepGeometric(), ppath_stepRefractionEuler(), rte_posAddRgeoid(), scat_data_monoCalc(), scat_data_rawCheck(), Select(), sensor_posAddRgeoid(), sensor_responseAntenna(), sensor_responseSimpleAMSU(), surface_calc(), surfaceFlatReflectivity(), surfaceFlatRefractiveIndex(), surfaceFlatVaryingEmissivity(), surfacetilt(), test01(), test04(), test06(), test1(), test31(), test36(), test4(), test40(), test6(), test7(), test9(), VectorExtractFromMatrix(), VectorZtanToZaRefr1D(), xsec_continuum_tag(), xsec_species(), y_unitApply(), ybatchCalc(), ybatchMetProfiles(), ybatchMetProfilesClear(), yCalc(), yCalc2(), z_at_lat_2d(), z_at_latlon(), z_fieldFromHSE(), and za_gridOpt().

Functions

Variables

Detailed Description

Function Documentation

◆ copy() [1/4]

◆ copy() [2/4]

◆ copy() [3/4]

◆ copy() [4/4]

◆ debug_matrixview_get_elem()

◆ max() [1/2]

◆ max() [2/2]

◆ mean() [1/2]

◆ mean() [2/2]

◆ min() [1/2]

◆ min() [2/2]

◆ mult() [1/2]

◆ mult() [2/2]

◆ operator*()

◆ operator<<() [1/2]

◆ operator<<() [2/2]

◆ transform() [1/2]

◆ transform() [2/2]

◆ transpose() [1/3]

◆ transpose() [2/3]

◆ transpose() [3/3]

Variable Documentation

◆ joker