interpolation_poly.h File Reference

Header file for interpolation_poly.cc. More...

#include "matpackI.h"
#include "interpolation.h"

Include dependency graph for interpolation_poly.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes

struct GridPosPoly

Structure to store a grid position for higher order interpolation. More...

Typedefs

typedef Array< GridPosPoly > ArrayOfGridPosPoly

An Array of grid positions.

Functions

ostream & operator<< (ostream &os, const GridPosPoly &gp)

Output operator for GridPosPoly.

void gridpos_poly (ArrayOfGridPosPoly &gp, ConstVectorView old_grid, ConstVectorView new_grid, const Index order, const Numeric &extpolfac=0.5)

Set up grid positions for higher order interpolation.

void gridpos_poly (GridPosPoly &gp, ConstVectorView old_grid, const Numeric &new_grid, const Index order, const Numeric &extpolfac=0.5)

gridpos_poly

void interpweights (VectorView itw, const GridPosPoly &tc)

Red 1D interpolation weights.

void interpweights (VectorView itw, const GridPosPoly &tr, const GridPosPoly &tc)

Red 2D interpolation weights.

void interpweights (VectorView itw, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 3D interpolation weights.

void interpweights (VectorView itw, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 4D interpolation weights.

void interpweights (VectorView itw, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 5D interpolation weights.

void interpweights (VectorView itw, const GridPosPoly &tv, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 6D interpolation weights.

Numeric interp (ConstVectorView itw, ConstVectorView a, const GridPosPoly &tc)

Red 1D Interpolate.

Numeric interp (ConstVectorView itw, ConstMatrixView a, const GridPosPoly &tr, const GridPosPoly &tc)

Red 2D Interpolate.

Numeric interp (ConstVectorView itw, ConstTensor3View a, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 3D Interpolate.

Numeric interp (ConstVectorView itw, ConstTensor4View a, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 4D Interpolate.

Numeric interp (ConstVectorView itw, ConstTensor5View a, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 5D Interpolate.

Numeric interp (ConstVectorView itw, ConstTensor6View a, const GridPosPoly &tv, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)

Red 6D Interpolate.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &cgp)

Compute 1D interpolation weights.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 2D interpolation weights for a sequence of positions.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 3D interpolation weights for a sequence of positions.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 4D interpolation weights for a sequence of positions.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 5D interpolation weights for a sequence of positions.

void interpweights (MatrixView itw, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 6D interpolation weights for a sequence of positions.

void interp (VectorView ia, ConstMatrixView itw, ConstVectorView a, const ArrayOfGridPosPoly &cgp)

Interpolate 1D field.

void interp (VectorView ia, ConstMatrixView itw, ConstMatrixView a, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 2D field to a sequence of positions.

void interp (VectorView ia, ConstMatrixView itw, ConstTensor3View a, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 3D field to a sequence of positions.

void interp (VectorView ia, ConstMatrixView itw, ConstTensor4View a, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 4D field to a sequence of positions.

void interp (VectorView ia, ConstMatrixView itw, ConstTensor5View a, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 5D field to a sequence of positions.

void interp (VectorView ia, ConstMatrixView itw, ConstTensor6View a, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 6D field to a sequence of positions.

void interpweights (Tensor3View itw, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 2D interpolation weights for an entire field.

void interpweights (Tensor4View itw, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 3D interpolation weights for an entire field.

void interpweights (Tensor5View itw, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 4D interpolation weights for an entire field.

void interpweights (Tensor6View itw, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 5D interpolation weights for an entire field.

void interpweights (Tensor7View itw, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Compute 6D interpolation weights for an entire field.

void interp (MatrixView ia, ConstTensor3View itw, ConstMatrixView a, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 2D field to another 2D field.

void interp (Tensor3View ia, ConstTensor4View itw, ConstTensor3View a, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 3D field to another 3D field.

void interp (Tensor4View ia, ConstTensor5View itw, ConstTensor4View a, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 4D field to another 4D field.

void interp (Tensor5View ia, ConstTensor6View itw, ConstTensor5View a, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 5D field to another 5D field.

void interp (Tensor6View ia, ConstTensor7View itw, ConstTensor6View a, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)

Interpolate 6D field to another 6D field.

Detailed Description

Header file for interpolation_poly.cc.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-02-04

Definition in file interpolation_poly.h.

Typedef Documentation

typedef Array<GridPosPoly> ArrayOfGridPosPoly

An Array of grid positions.

See GridPosPoly for details.

Definition at line 78 of file interpolation_poly.h.

Function Documentation

void gridpos_poly	(	GridPosPoly &	gp,
		ConstVectorView	old_grid,
		const Numeric &	new_grid,
		const Index	order,
		const Numeric &	extpolfac
	)

gridpos_poly

Creates a grid position structure for higher order interpolation.

This is a gateway to the function for arrays of GridPosPoly, to be used for e.g. "red interpolation".

Return values:

gp

The GridPos structure.

Parameters:

	old_grid	The original grid.
	new_grid	The position where we want to have the interpolated value.
	order	Interpolation order. 1 = linear, 2 = quadratic, etc.. The number of points used in the interpolation will be order+1.
	extpolfac	Extrapolation factor. Default value is 0.5, which means that extrapolation of half of the last grid distance is allowed. You don't have to specify this.

Author:: Stefan Buehler

Date:: 2008-03-06

Definition at line 196 of file interpolation_poly.cc.

References gridpos_poly().

void gridpos_poly	(	ArrayOfGridPosPoly &	gp,
		ConstVectorView	old_grid,
		ConstVectorView	new_grid,
		const Index	order,
		const Numeric &	extpolfac
	)

Set up grid positions for higher order interpolation.

This function performs the same task as gridpos, but for arbitrary orders of interpolation. (Linear is also possible, as a special case.)

The formula for calculating the weights w is taken from Numerical Recipes, 2nd edition, section 3.1, eq. 3.1.1.

Parameters:

	gp	Output: An array of grid positions.
	old_grid	Original grid.
	new_grid	New grid.
	order	Interpolation order. 1 = linear, 2 = quadratic, etc.. The number of points used in the interpolation will be order+1.
	extpolfac	Extrapolation fraction. Should normally not be specified, then the default of 0.5 is used.

Definition at line 89 of file interpolation_poly.cc.

References gridpos(), IMAX, IMIN, is_size(), and ConstVectorView::nelem().

Referenced by abs_lookupSetupBatch(), choose_abs_nls_pert(), choose_abs_t_pert(), GasAbsLookup::Extract(), gridpos_poly(), test07(), and test08().

void interp	(	Tensor6View	ia,
		ConstTensor7View	itw,
		ConstTensor6View	a,
		const ArrayOfGridPosPoly &	vgp,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 6D field to another 6D field.

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2727 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), ns, and sum_check_epsilon.

void interp	(	Tensor5View	ia,
		ConstTensor6View	itw,
		ConstTensor5View	a,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 5D field to another 5D field.

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2629 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), ns, and sum_check_epsilon.

void interp	(	Tensor4View	ia,
		ConstTensor5View	itw,
		ConstTensor4View	a,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 4D field to another 4D field.

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2540 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	Tensor3View	ia,
		ConstTensor4View	itw,
		ConstTensor3View	a,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 3D field to another 3D field.

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2460 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	MatrixView	ia,
		ConstTensor3View	itw,
		ConstMatrixView	a,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 2D field to another 2D field.

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

This function can be easily distinguished from the other interpolation function (that creates a sequence of interpolated values), because of the dimension of ia and itw.

The size of ia and itw in all dimensions must be consistent with the grid position Arrays.

Return values:

ia

Interpolated field.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2390 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstTensor6View	a,
		const ArrayOfGridPosPoly &	vgp,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 6D field to a sequence of positions.

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1953 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstTensor5View	a,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 5D field to a sequence of positions.

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1870 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstTensor4View	a,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 4D field to a sequence of positions.

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1793 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstTensor3View	a,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 3D field to a sequence of positions.

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1722 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstMatrixView	a,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 2D field to a sequence of positions.

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

This performs a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The output vector ia also must have the same length. And the dimension of itw must be consistent with this.

This function can be easily distinguished from the other interpolation function (that creates an entire field of interpolated values), because of the dimension of ia and itw.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1657 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

void interp	(	VectorView	ia,
		ConstMatrixView	itw,
		ConstVectorView	a,
		const ArrayOfGridPosPoly &	cgp
	)

Interpolate 1D field.

For this 1D case there is no distinction between "blue" and "green" type interpolation.

The output vector ia must have the same length as the grid position vector cgp. And the dimension of itw must be consistent with this.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

For this 1D case there is no distinction between "blue" and "green" type interpolation.

The output vector ia must have the same length as the grid position vector cgp. And the dimension of itw must be consistent with this.

Return values:

ia

Vector containing the interpolated field values.

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1598 of file interpolation.cc.

References GridPos::idx, is_same_within_epsilon(), is_size(), joker, LOOPIDX, Array< base >::nelem(), and sum_check_epsilon.

Numeric interp	(	ConstVectorView	itw,
		ConstTensor6View	a,
		const GridPosPoly &	tv,
		const GridPosPoly &	ts,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 6D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	tv	The grid position for the vitrine dimension.
	ts	The grid position for the shelf dimension.
	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 791 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

Numeric interp	(	ConstVectorView	itw,
		ConstTensor5View	a,
		const GridPosPoly &	ts,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 5D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	ts	The grid position for the shelf dimension.
	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 728 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

Numeric interp	(	ConstVectorView	itw,
		ConstTensor4View	a,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 4D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 670 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

Numeric interp	(	ConstVectorView	itw,
		ConstTensor3View	a,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 3D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 617 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

Numeric interp	(	ConstVectorView	itw,
		ConstMatrixView	a,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 2D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 569 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

Numeric interp	(	ConstVectorView	itw,
		ConstVectorView	a,
		const GridPosPoly &	tc
	)

Red 1D Interpolate.

"Red" interpolation returns just a scalar.

The dimension of itw must be consistent with the dimension of the interpolation (2^n).

Parameters:

	itw	Interpolation weights.
	a	The field to interpolate.
	tc	The grid position for the column dimension.

Returns:: Interpolated value.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 527 of file interpolation_poly.cc.

References is_same_within_epsilon(), is_size(), LOOPIDX, ConstVectorView::nelem(), ConstVectorView::sum(), sum_check_epsilon, and GridPosPoly::w.

void interpweights	(	Tensor7View	itw,
		const ArrayOfGridPosPoly &	vgp,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 6D interpolation weights for an entire field.

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2310 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, Array< base >::nelem(), and ns.

void interpweights	(	Tensor6View	itw,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 5D interpolation weights for an entire field.

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2233 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, Array< base >::nelem(), and ns.

void interpweights	(	Tensor5View	itw,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 4D interpolation weights for an entire field.

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2164 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	Tensor4View	itw,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 3D interpolation weights for an entire field.

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2103 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	Tensor3View	itw,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 2D interpolation weights for an entire field.

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "green" type interpolation of the field, that means that the grid position Arrays are interpreted as defining the grids for the interpolated field.

The dimensions of itw must be consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Tensor with one more dimension than the number of grid position Arrays.

Return values:

itw

Interpolation weights

Parameters:

	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 2041 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	vgp,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 6D interpolation weights for a sequence of positions.

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	vgp	The grid position Array for the vitrine dimension.
	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1539 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	sgp,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 5D interpolation weights for a sequence of positions.

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	sgp	The grid position Array for the shelf dimension.
	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1475 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	bgp,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 4D interpolation weights for a sequence of positions.

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	bgp	The grid position Array for the book dimension.
	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1416 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	pgp,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 3D interpolation weights for a sequence of positions.

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	pgp	The grid position Array for the page dimension.
	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1362 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	rgp,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 2D interpolation weights for a sequence of positions.

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

Compute the weights for a "blue" type interpolation of the field, that means that the grid position Arrays are interpreted as defining a sequence of positions. ALL GRID POSITION ARRAYS MUST HAVE THE SAME LENGTH!

The dimensions of itw must be also consistent with this.

Note that we still do not need the actual field for this step.

This function can be easily distinguished from the other interpweights function (for "green" interpolation), because the output is a Matrix, whereas in the other case it is a Tensor with one more dimension than there are input grid position Arrays.

Return values:

itw

Interpolation weights.

Parameters:

	rgp	The grid position Array for the row dimension.
	cgp	The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1305 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	MatrixView	itw,
		const ArrayOfGridPosPoly &	cgp
	)

Compute 1D interpolation weights.

For this 1D case there is no distinction between "blue" and "green" type interpolation.

The dimensions of itw must be consistent with cgp.

Note that we still do not need the actual field for this step.

Return values:

itw

Interpolation weights.

Parameters:

cgp

The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler@ltu.se>

Date:: Fri May 3 08:55:51 2002

For this 1D case there is no distinction between "blue" and "green" type interpolation.

The dimensions of itw must be consistent with cgp.

Note that we still do not need the actual field for this step.

Return values:

itw

Interpolation weights.

Parameters:

cgp

The grid position Array for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 1232 of file interpolation.cc.

References is_size(), LOOPIT, LOOPW, and Array< base >::nelem().

void interpweights	(	VectorView	itw,
		const GridPosPoly &	tv,
		const GridPosPoly &	ts,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 6D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

	tv	The grid position for the vitrine dimension.
	ts	The grid position for the shelf dimension.
	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 481 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

void interpweights	(	VectorView	itw,
		const GridPosPoly &	ts,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 5D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

	ts	The grid position for the shelf dimension.
	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 435 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

void interpweights	(	VectorView	itw,
		const GridPosPoly &	tb,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 4D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

	tb	The grid position for the book dimension.
	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 393 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

void interpweights	(	VectorView	itw,
		const GridPosPoly &	tp,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 3D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

	tp	The grid position for the page dimension.
	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 355 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

void interpweights	(	VectorView	itw,
		const GridPosPoly &	tr,
		const GridPosPoly &	tc
	)

Red 2D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

	tr	The grid position for the row dimension.
	tc	The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 322 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

void interpweights	(	VectorView	itw,
		const GridPosPoly &	tc
	)

Red 1D interpolation weights.

"Red" interpolation returns just a scalar, so the weights are stored in a Vector.

The length of itw must be consistent with the dimension of the field to be interpolated (2^n).

Return values:

itw

Interpolation weights.

Parameters:

tc

The grid position for the column dimension.

Author:: Stefan Buehler <sbuehler(at)ltu.se>

Date:: 2008-03-06

Definition at line 288 of file interpolation_poly.cc.

References is_size(), LOOPW, ConstVectorView::nelem(), and GridPosPoly::w.

ostream& operator<<	(	ostream &	os,
		const GridPosPoly &	gp
	)

Output operator for GridPosPoly.

This is just intended for testing and debugging.

Parameters:

	os	Output stream.
	gp	Grid position.

Returns:: The output stream.

Definition at line 252 of file interpolation_poly.cc.

References GridPosPoly::idx, and GridPosPoly::w.


Classes
struct	GridPosPoly
	Structure to store a grid position for higher order interpolation. More...
Typedefs
typedef Array< GridPosPoly >	ArrayOfGridPosPoly
	An Array of grid positions.
Functions
ostream &	operator<< (ostream &os, const GridPosPoly &gp)
	Output operator for GridPosPoly.
void	gridpos_poly (ArrayOfGridPosPoly &gp, ConstVectorView old_grid, ConstVectorView new_grid, const Index order, const Numeric &extpolfac=0.5)
	Set up grid positions for higher order interpolation.
void	gridpos_poly (GridPosPoly &gp, ConstVectorView old_grid, const Numeric &new_grid, const Index order, const Numeric &extpolfac=0.5)
	gridpos_poly
void	interpweights (VectorView itw, const GridPosPoly &tc)
	Red 1D interpolation weights.
void	interpweights (VectorView itw, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 2D interpolation weights.
void	interpweights (VectorView itw, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 3D interpolation weights.
void	interpweights (VectorView itw, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 4D interpolation weights.
void	interpweights (VectorView itw, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 5D interpolation weights.
void	interpweights (VectorView itw, const GridPosPoly &tv, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 6D interpolation weights.
Numeric	interp (ConstVectorView itw, ConstVectorView a, const GridPosPoly &tc)
	Red 1D Interpolate.
Numeric	interp (ConstVectorView itw, ConstMatrixView a, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 2D Interpolate.
Numeric	interp (ConstVectorView itw, ConstTensor3View a, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 3D Interpolate.
Numeric	interp (ConstVectorView itw, ConstTensor4View a, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 4D Interpolate.
Numeric	interp (ConstVectorView itw, ConstTensor5View a, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 5D Interpolate.
Numeric	interp (ConstVectorView itw, ConstTensor6View a, const GridPosPoly &tv, const GridPosPoly &ts, const GridPosPoly &tb, const GridPosPoly &tp, const GridPosPoly &tr, const GridPosPoly &tc)
	Red 6D Interpolate.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &cgp)
	Compute 1D interpolation weights.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 2D interpolation weights for a sequence of positions.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 3D interpolation weights for a sequence of positions.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 4D interpolation weights for a sequence of positions.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 5D interpolation weights for a sequence of positions.
void	interpweights (MatrixView itw, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 6D interpolation weights for a sequence of positions.
void	interp (VectorView ia, ConstMatrixView itw, ConstVectorView a, const ArrayOfGridPosPoly &cgp)
	Interpolate 1D field.
void	interp (VectorView ia, ConstMatrixView itw, ConstMatrixView a, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 2D field to a sequence of positions.
void	interp (VectorView ia, ConstMatrixView itw, ConstTensor3View a, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 3D field to a sequence of positions.
void	interp (VectorView ia, ConstMatrixView itw, ConstTensor4View a, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 4D field to a sequence of positions.
void	interp (VectorView ia, ConstMatrixView itw, ConstTensor5View a, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 5D field to a sequence of positions.
void	interp (VectorView ia, ConstMatrixView itw, ConstTensor6View a, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 6D field to a sequence of positions.
void	interpweights (Tensor3View itw, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 2D interpolation weights for an entire field.
void	interpweights (Tensor4View itw, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 3D interpolation weights for an entire field.
void	interpweights (Tensor5View itw, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 4D interpolation weights for an entire field.
void	interpweights (Tensor6View itw, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 5D interpolation weights for an entire field.
void	interpweights (Tensor7View itw, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Compute 6D interpolation weights for an entire field.
void	interp (MatrixView ia, ConstTensor3View itw, ConstMatrixView a, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 2D field to another 2D field.
void	interp (Tensor3View ia, ConstTensor4View itw, ConstTensor3View a, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 3D field to another 3D field.
void	interp (Tensor4View ia, ConstTensor5View itw, ConstTensor4View a, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 4D field to another 4D field.
void	interp (Tensor5View ia, ConstTensor6View itw, ConstTensor5View a, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 5D field to another 5D field.
void	interp (Tensor6View ia, ConstTensor7View itw, ConstTensor6View a, const ArrayOfGridPosPoly &vgp, const ArrayOfGridPosPoly &sgp, const ArrayOfGridPosPoly &bgp, const ArrayOfGridPosPoly &pgp, const ArrayOfGridPosPoly &rgp, const ArrayOfGridPosPoly &cgp)
	Interpolate 6D field to another 6D field.