Stuff related to the propagation matrix. More...

#include "complex.h"
#include "matpackIV.h"

Classes
class	LazyScale< base >
	Class to help with hidden temporary variables for operations of type Numeric times Class. More...

class	PropagationMatrix

class	StokesVector
	Stokes vector is as Propagation matrix but only has 4 possible values. More...

Typedefs
typedef Array< PropagationMatrix >	ArrayOfPropagationMatrix

typedef Array< ArrayOfPropagationMatrix >	ArrayOfArrayOfPropagationMatrix

typedef Array< StokesVector >	ArrayOfStokesVector

typedef Array< ArrayOfStokesVector >	ArrayOfArrayOfStokesVector

typedef Array< ArrayOfArrayOfStokesVector >	ArrayOfArrayOfArrayOfStokesVector

Functions
void	compute_transmission_matrix (Tensor3View T, const Numeric &r, const PropagationMatrix &upper_level, const PropagationMatrix &lower_level, const Index iz=0, const Index ia=0)
	Compute the matrix exponent as the transmission matrix of this propagation matrix. More...

void	compute_transmission_matrix_from_averaged_matrix_at_frequency (MatrixView T, const Numeric &r, const PropagationMatrix &averaged_propagation_matrix, const Index iv, const Index iz=0, const Index ia=0)
	Compute the matrix exponent as the transmission matrix of this propagation matrix. More...

void	compute_transmission_matrix_and_derivative (Tensor3View T, Tensor4View dT_upper_level, Tensor4View dT_lower_level, const Numeric &r, const PropagationMatrix &upper_level, const PropagationMatrix &lower_level, const Array< PropagationMatrix > &dprop_mat_upper_level, const Array< PropagationMatrix > &dprop_mat_lower_level, const Numeric &dr_dTu=0.0, const Numeric &dr_dTl=0.0, const Index it=-1, const Index iz=0, const Index ia=0)

std::ostream &	operator<< (std::ostream &os, const PropagationMatrix &pm)
	output operator More...

std::ostream &	operator<< (std::ostream &os, const ArrayOfPropagationMatrix &apm)
	output operator More...

std::ostream &	operator<< (std::ostream &os, const ArrayOfArrayOfPropagationMatrix &aapm)
	output operator More...

std::ostream &	operator<< (std::ostream &os, const StokesVector &pm)

std::ostream &	operator<< (std::ostream &os, const ArrayOfStokesVector &apm)

std::ostream &	operator<< (std::ostream &os, const ArrayOfArrayOfStokesVector &aapm)

LazyScale< PropagationMatrix >	operator* (const PropagationMatrix &pm, const Numeric &x)
	Returns a lazy multiplier. More...

LazyScale< PropagationMatrix >	operator* (const Numeric &x, const PropagationMatrix &pm)
	Returns a lazy multiplier. More...

Detailed Description

Stuff related to the propagation matrix.

Author: Richard Larsson

Date: 2017-06-23

This implementation takes advantage of symmetries to lower memory and computational costs

Some standard functions applied on the propagation matrix have been included but far from all.

Present support is for stokes dim 1-4 however please take note that circular polarization delay variable's position in the internal mdata structure moves depending on if stokes_dim is 3 or 4. It is not present elsewhere...

Definition in file propagationmatrix.h.

Typedef Documentation

◆ ArrayOfArrayOfArrayOfStokesVector

typedef Array<ArrayOfArrayOfStokesVector> ArrayOfArrayOfArrayOfStokesVector

Definition at line 1356 of file propagationmatrix.h.

◆ ArrayOfArrayOfPropagationMatrix

typedef Array<ArrayOfPropagationMatrix> ArrayOfArrayOfPropagationMatrix

Definition at line 979 of file propagationmatrix.h.

◆ ArrayOfArrayOfStokesVector

typedef Array<ArrayOfStokesVector> ArrayOfArrayOfStokesVector

Definition at line 1355 of file propagationmatrix.h.

◆ ArrayOfPropagationMatrix

typedef Array<PropagationMatrix> ArrayOfPropagationMatrix

Definition at line 978 of file propagationmatrix.h.

◆ ArrayOfStokesVector

typedef Array<StokesVector> ArrayOfStokesVector

Definition at line 1354 of file propagationmatrix.h.

Function Documentation

◆ compute_transmission_matrix()

void compute_transmission_matrix	(	Tensor3View	T,
		const Numeric &	r,
		const PropagationMatrix &	upper_level,
		const PropagationMatrix &	lower_level,
		const Index	iz = `0`,
		const Index	ia = `0`
	)

Compute the matrix exponent as the transmission matrix of this propagation matrix.

The propagation matrix is multiplied by -r and level-averaged before exponent is applied.

upper_level and lower_level propagation matrices should thus be the level matrices and r the distance between these levels. The same is true for the derivative matrices.

Stokes dim 1 and 4 have been tested more. Stokes dim 2 and 3 have been found to work but could still have hidden errors for uncommon cases

Parameters

[in,out]	T	transmission tensor with outmost dimension being frequency
[in]	r	the distance over which the propagation matrix causes the transmission
[in]	upper_level	The upper level propagation matrix
[in]	lower_level	The lower level propagation matrix
[in]	iz	Zenith index
[in]	ia	Azimuth index

Definition at line 33 of file propagationmatrix.cc.

References F, i, joker, PropagationMatrix::K12(), PropagationMatrix::Kjj(), PropagationMatrix::NumberOfFrequencies(), and PropagationMatrix::StokesDimensions().

Referenced by get_stepwise_transmission_matrix().

◆ compute_transmission_matrix_and_derivative()

void compute_transmission_matrix_and_derivative	(	Tensor3View	T,
		Tensor4View	dT_upper_level,
		Tensor4View	dT_lower_level,
		const Numeric &	r,
		const PropagationMatrix &	upper_level,
		const PropagationMatrix &	lower_level,
		const Array< PropagationMatrix > &	dprop_mat_upper_level,
		const Array< PropagationMatrix > &	dprop_mat_lower_level,
		const Numeric &	dr_dTu = `0.0`,
		const Numeric &	dr_dTl = `0.0`,
		const Index	it = `-1`,
		const Index	iz = `0`,
		const Index	ia = `0`
	)

Compute the matrix exponent as the transmission matrix of this propagation matrix

The propagation matrix is multiplied by -r and level-averaged before exponent is applied.

upper_level and lower_level propagation matrices should thus be the level matrices and r the distance between these levels. The same is true for the derivative matrices.

Stokes dim 1 and 4 have been tested more. Stokes dim 2 and 3 have been found to work but could still have hidden errors for uncommon cases

Parameters

[in,out]	T	transmission tensor with outmost dimension being frequency
[in,out]	dT_upp	transmission tensors derivative with respect to derivatives of the propagation matrix for upper level
[in,out]	dT_low	transmission tensors derivative with respect to derivatives of the propagation matrix for lower level
[in]	r	the distance over which the propagation matrix causes the transmission
[in]	upper_level	The upper level propagation matrix
[in]	lower_level	The lower level propagation matrix
[in]	dprop_mat_upp	derivatives of the upper propagation matrix with respect to some parameter (is multiplied by -0.5 r)
[in]	dprop_mat_low	derivatives of the lower propagation matrix with respect to some parameter (is multiplied by -0.5 r)
[in]	dr_dTu	Distance temperature derivative for upper level
[in]	dr_dTl	Distance temperature derivative for lower level
[in]	it	Position of temperature derivative (ignored at -1)
[in]	iz	Zenith index
[in]	ia	Azimuth index

Definition at line 478 of file propagationmatrix.cc.

References F, i, joker, PropagationMatrix::K12(), PropagationMatrix::Kjj(), Array< base >::nelem(), PropagationMatrix::NumberOfFrequencies(), and PropagationMatrix::StokesDimensions().

Referenced by get_stepwise_transmission_matrix().

◆ compute_transmission_matrix_from_averaged_matrix_at_frequency()

void compute_transmission_matrix_from_averaged_matrix_at_frequency	(	MatrixView	T,
		const Numeric &	r,
		const PropagationMatrix &	averaged_propagation_matrix,
		const Index	iv,
		const Index	iz = `0`,
		const Index	ia = `0`
	)

Compute the matrix exponent as the transmission matrix of this propagation matrix.

The propagation matrix is multiplied by -r and level-averaged before exponent is applied.

Stokes dim 1 and 4 have been tested more. Stokes dim 2 and 3 have been found to work but could still have hidden errors for uncommon cases

Parameters

[in,out]	T	transmission matrix at the position
[in]	r	the distance over which the propagation matrix causes the transmission
[in]	averaged_propagation_matrix	The propagation matrix
[in]	iv	Frequency index
[in]	iz	Zenith index
[in]	ia	Azimuth index

Definition at line 272 of file propagationmatrix.cc.

References PropagationMatrix::K12(), PropagationMatrix::Kjj(), sqrt(), and PropagationMatrix::StokesDimensions().

Referenced by rte_step_doit_replacement().

◆ operator*() [1/2]

LazyScale<PropagationMatrix> operator*	(	const PropagationMatrix &	pm,
		const Numeric &	x
	)

inline

Returns a lazy multiplier.

Parameters

[in]	pm	Propagation matrix
[in]	x	Scale

Returns: LazyScale<PropagationMatrix> A lazy multiplier

Definition at line 1369 of file propagationmatrix.h.

References pm.

◆ operator*() [2/2]

LazyScale<PropagationMatrix> operator*	(	const Numeric &	x,
		const PropagationMatrix &	pm
	)

inline

Returns a lazy multiplier.

Parameters

[in]	x	Scale
[in]	pm	Propagation matrix

Returns: LazyScale<PropagationMatrix> A lazy multiplier

Definition at line 1380 of file propagationmatrix.h.

References pm.

◆ operator<<() [1/6]

std::ostream& operator<<	(	std::ostream &	os,
		const PropagationMatrix &	pm
	)

output operator

Definition at line 2036 of file propagationmatrix.cc.

References PropagationMatrix::Data().

◆ operator<<() [2/6]

std::ostream& operator<<	(	std::ostream &	os,
		const ArrayOfPropagationMatrix &	apm
	)

output operator

Definition at line 2041 of file propagationmatrix.cc.

References pm.

◆ operator<<() [3/6]

std::ostream& operator<<	(	std::ostream &	os,
		const ArrayOfArrayOfPropagationMatrix &	aapm
	)

output operator

Definition at line 2047 of file propagationmatrix.cc.

◆ operator<<() [4/6]

std::ostream& operator<<	(	std::ostream &	os,
		const StokesVector &	pm
	)

Definition at line 2054 of file propagationmatrix.cc.

References PropagationMatrix::Data().

◆ operator<<() [5/6]

std::ostream& operator<<	(	std::ostream &	os,
		const ArrayOfStokesVector &	apm
	)

Definition at line 2059 of file propagationmatrix.cc.

◆ operator<<() [6/6]

std::ostream& operator<<	(	std::ostream &	os,
		const ArrayOfArrayOfStokesVector &	aapm
	)

Definition at line 2064 of file propagationmatrix.cc.

Classes

Typedefs

Functions

Detailed Description

Typedef Documentation

◆ ArrayOfArrayOfArrayOfStokesVector

◆ ArrayOfArrayOfPropagationMatrix

◆ ArrayOfArrayOfStokesVector

◆ ArrayOfPropagationMatrix

◆ ArrayOfStokesVector

Function Documentation

◆ compute_transmission_matrix()

◆ compute_transmission_matrix_and_derivative()

◆ compute_transmission_matrix_from_averaged_matrix_at_frequency()

◆ operator*() [1/2]

◆ operator*() [2/2]

◆ operator<<() [1/6]

◆ operator<<() [2/6]

◆ operator<<() [3/6]

◆ operator<<() [4/6]

◆ operator<<() [5/6]

◆ operator<<() [6/6]