dpmInverse.C File Reference

#include "DataPointMapping.h"
#include "MappingInformation.h"
#include "arrayGetIndex.h"
#include <string.h>
#include "conversion.h"
#include "display.h"
#include "Inverse.h"
#include "ParallelUtility.h"

Include dependency graph for dpmInverse.C:

Macros
#define	INT_1D_ORDER_2(dr, x111, x211)   ( (1.-dr)*(x111)+dr*(x211) )
#define	INT_1D_ORDER_2_R(dr, x111, x211)   ( delta[0]*( (x211)-(x111) ) )
#define	INT_2D_ORDER_2(dr, ds, x111, x211, x121, x221)   ( (1.-ds)*((1.-dr)*(x111)+dr*(x211))+ds*((1.-dr)*(x121)+dr*(x221)) )
#define	INT_2D_ORDER_2_R(dr, ds, x111, x211, x121, x221)   ( ((1.-ds)*( (x211)-(x111) ) +ds*( (x221)-(x121) ))*delta[0] )
#define	INT_2D_ORDER_2_S(dr, ds, x111, x211, x121, x221)   ( ((1.-dr)*( (x121)-(x111) ) +dr*( (x221)-(x211) ))*delta[1] )
#define	INT_3D_ORDER_2(dr, ds, dt, x111, x211, x121, x221, x112, x212, x122, x222)   ( (1.-dt)*((1.-ds)*((1.-dr)*(x111)+dr*(x211))+ds*((1.-dr)*(x121)+dr*(x221))) + dt *((1.-ds)*((1.-dr)*(x112)+dr*(x212))+ds*((1.-dr)*(x122)+dr*(x222))) )
#define	INT_3D_ORDER_2_R(dr, ds, dt, x111, x211, x121, x221, x112, x212, x122, x222)   ( ( (1.-dt)*((1.-ds)*((x211)-(x111))+ds*((x221)-(x121))) + dt *((1.-ds)*((x212)-(x112))+ds*((x222)-(x122))) )*delta[0])
#define	INT_3D_ORDER_2_S(dr, ds, dt, x111, x211, x121, x221, x112, x212, x122, x222)   ( ( (1.-dt)*((1.-dr)*((x121)-(x111))+dr*((x221)-(x211))) + dt *((1.-dr)*((x122)-(x112))+dr*((x222)-(x212))) )*delta[1] )
#define	INT_3D_ORDER_2_T(dr, ds, dt, x111, x211, x121, x221, x112, x212, x122, x222)   ( ( (1.-dr)*((1.-ds)*((x112)-(x111))+ds*((x122)-(x121))) + dr *((1.-ds)*((x212)-(x211))+ds*((x222)-(x221))) )*delta[2] )
#define	q03(z)   (-oneSixth *((z)-1.)*((z)-2.)*((z)-3.))
#define	q13(z)   ( .5*(z) *((z)-2.)*((z)-3.))
#define	q23(z)   (-.5*(z)*((z)-1.) *((z)-3.))
#define	q33(z)   ( oneSixth*(z)*((z)-1.)*((z)-2.))
#define	q03d(z)   ( -oneSixth*(11.+(z)*(-12.+3.*(z))))
#define	q13d(z)   ( 3.+(z)*(-5.+1.5*(z)))
#define	q23d(z)   ( -1.5+(z)*(4.-1.5*(z)))
#define	q33d(z)   (oneSixth*(2.+(z)*(-6.+3.*(z))))
#define	q1x(i1, i2, i3, axis)   ( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )
#define	q2x(i1, i2, i3, axis)   ( b0*q1x(i1 ,i2 ,i3 ,axis) +b1*q1x(i1 ,i2+1,i3 ,axis) +b2*q1x(i1 ,i2+2,i3 ,axis) +b3*q1x(i1 ,i2+3,i3 ,axis) )
#define	q3x(i1, i2, i3, axis)   ( c0*q2x(i1 ,i2 ,i3 ,axis) +c1*q2x(i1 ,i2 ,i3+1,axis) +c2*q2x(i1 ,i2 ,i3+2,axis) +c3*q2x(i1 ,i2 ,i3+3,axis) )
#define	q1xr(i1, i2, i3, axis)   ( a0r*XY(i1 ,i2 ,i3 ,axis) +a1r*XY(i1+1,i2 ,i3 ,axis) +a2r*XY(i1+2,i2 ,i3 ,axis) +a3r*XY(i1+3,i2 ,i3 ,axis) )
#define	q2xr(i1, i2, i3, axis)   ( b0*q1xr(i1 ,i2 ,i3 ,axis) +b1*q1xr(i1 ,i2+1,i3 ,axis) +b2*q1xr(i1 ,i2+2,i3 ,axis) +b3*q1xr(i1 ,i2+3,i3 ,axis) )
#define	q3xr(i1, i2, i3, axis)   ( c0*q2xr(i1 ,i2 ,i3 ,axis) +c1*q2xr(i1 ,i2 ,i3+1,axis) +c2*q2xr(i1 ,i2 ,i3+2,axis) +c3*q2xr(i1 ,i2 ,i3+3,axis) )
#define	q1xs(i1, i2, i3, axis)   ( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )
#define	q2xs(i1, i2, i3, axis)   ( b0r*q1xs(i1 ,i2 ,i3 ,axis) +b1r*q1xs(i1 ,i2+1,i3 ,axis) +b2r*q1xs(i1 ,i2+2,i3 ,axis) +b3r*q1xs(i1 ,i2+3,i3 ,axis))
#define	q3xs(i1, i2, i3, axis)   ( c0*q2xs(i1 ,i2 ,i3 ,axis) +c1*q2xs(i1 ,i2 ,i3+1,axis) +c2*q2xs(i1 ,i2 ,i3+2,axis) +c3*q2xs(i1 ,i2 ,i3+3,axis) )
#define	q1xt(i1, i2, i3, axis)   ( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )
#define	q2xt(i1, i2, i3, axis)   ( b0*q1xt(i1 ,i2 ,i3 ,axis) +b1*q1xt(i1 ,i2+1,i3 ,axis) +b2*q1xt(i1 ,i2+2,i3 ,axis) +b3*q1xt(i1 ,i2+3,i3 ,axis) )
#define	q3xt(i1, i2, i3, axis)   ( c0r*q2xt(i1 ,i2 ,i3 ,axis) +c1r*q2xt(i1 ,i2 ,i3+1,axis) +c2r*q2xt(i1 ,i2 ,i3+2,axis) +c3r*q2xt(i1 ,i2 ,i3+3,axis) )
#define	XY(i0, i1, i2, i3)   xyp[i0+xyDim0*(i1+xyDim1*(i2+xyDim2*(i3)))]
#define	R(i0, i1)   rp[i0+rDim0*(i1)]
#define	R0(i0, i1)   r0p[i0+r0Dim0*(i1)]
#define	X(i0, i1)   xp[i0+xDim0*(i1)]
#define	RX(i0, i1, i2)   rxp[i0+rxDim0*(i1+rxDim1*(i2))]
#define	GRIDINDEXRANGE(i0, i1)   gridIndexRangep[i0+gridIndexRangeDim0*(i1)]
#define	dim(side, axis)   pdim[(side)+2*(axis)]
#define	singular(side, axis)   singularp[(side)+2*(axis)]

Macro Definition Documentation

#define dim	(		side,
			axis
	)		pdim[(side)+2*(axis)]

#define GRIDINDEXRANGE	(		i0,
			i1
	)		gridIndexRangep[i0+gridIndexRangeDim0*(i1)]

Referenced by DataPointMapping::basicInverseS().

#define INT_1D_ORDER_2	(		dr,
			x111,
			x211
	)		( (1.-dr)*(x111)+dr*(x211) )

#define INT_1D_ORDER_2_R	(		dr,
			x111,
			x211
	)		( delta[0]*( (x211)-(x111) ) )

#define INT_2D_ORDER_2	(		dr,
			ds,
			x111,
			x211,
			x121,
			x221
	)		( (1.-ds)*((1.-dr)*(x111)+dr*(x211))+ds*((1.-dr)*(x121)+dr*(x221)) )

Referenced by DataPointMapping::basicInverseS().

#define INT_2D_ORDER_2_R	(		dr,
			ds,
			x111,
			x211,
			x121,
			x221
	)		( ((1.-ds)*( (x211)-(x111) ) +ds*( (x221)-(x121) ))*delta[0] )

Referenced by DataPointMapping::basicInverseS().

#define INT_2D_ORDER_2_S	(		dr,
			ds,
			x111,
			x211,
			x121,
			x221
	)		( ((1.-dr)*( (x121)-(x111) ) +dr*( (x221)-(x211) ))*delta[1] )

Referenced by DataPointMapping::basicInverseS().

#define INT_3D_ORDER_2	(		dr,
			ds,
			dt,
			x111,
			x211,
			x121,
			x221,
			x112,
			x212,
			x122,
			x222
	)		( (1.-dt)*((1.-ds)*((1.-dr)*(x111)+dr*(x211))+ds*((1.-dr)*(x121)+dr*(x221))) + dt *((1.-ds)*((1.-dr)*(x112)+dr*(x212))+ds*((1.-dr)*(x122)+dr*(x222))) )

Referenced by DataPointMapping::basicInverseS().

#define INT_3D_ORDER_2_R	(		dr,
			ds,
			dt,
			x111,
			x211,
			x121,
			x221,
			x112,
			x212,
			x122,
			x222
	)		( ( (1.-dt)*((1.-ds)*((x211)-(x111))+ds*((x221)-(x121))) + dt *((1.-ds)*((x212)-(x112))+ds*((x222)-(x122))) )*delta[0])

Referenced by DataPointMapping::basicInverseS().

#define INT_3D_ORDER_2_S	(		dr,
			ds,
			dt,
			x111,
			x211,
			x121,
			x221,
			x112,
			x212,
			x122,
			x222
	)		( ( (1.-dt)*((1.-dr)*((x121)-(x111))+dr*((x221)-(x211))) + dt *((1.-dr)*((x122)-(x112))+dr*((x222)-(x212))) )*delta[1] )

Referenced by DataPointMapping::basicInverseS().

#define INT_3D_ORDER_2_T	(		dr,
			ds,
			dt,
			x111,
			x211,
			x121,
			x221,
			x112,
			x212,
			x122,
			x222
	)		( ( (1.-dr)*((1.-ds)*((x112)-(x111))+ds*((x122)-(x121))) + dr *((1.-ds)*((x212)-(x211))+ds*((x222)-(x221))) )*delta[2] )

Referenced by DataPointMapping::basicInverseS().

#define q03

(

z

)

(-oneSixth *((z)-1.)*((z)-2.)*((z)-3.))

Referenced by DataPointMapping::basicInverseS().

#define q03d

(

z

)

( -oneSixth*(11.+(z)*(-12.+3.*(z))))

Referenced by DataPointMapping::basicInverseS().

#define q13

(

z

)

( .5*(z) *((z)-2.)*((z)-3.))

Referenced by DataPointMapping::basicInverseS().

#define q13d

(

z

)

( 3.+(z)*(-5.+1.5*(z)))

Referenced by DataPointMapping::basicInverseS().

#define q1x	(		i1,
			i2,
			i3,
			axis
	)		( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )

#define q1xr	(		i1,
			i2,
			i3,
			axis
	)		( a0r*XY(i1 ,i2 ,i3 ,axis) +a1r*XY(i1+1,i2 ,i3 ,axis) +a2r*XY(i1+2,i2 ,i3 ,axis) +a3r*XY(i1+3,i2 ,i3 ,axis) )

#define q1xs	(		i1,
			i2,
			i3,
			axis
	)		( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )

#define q1xt	(		i1,
			i2,
			i3,
			axis
	)		( a0* XY(i1 ,i2 ,i3 ,axis) +a1* XY(i1+1,i2 ,i3 ,axis) +a2* XY(i1+2,i2 ,i3 ,axis) +a3* XY(i1+3,i2 ,i3 ,axis) )

#define q23

(

z

)

(-.5*(z)*((z)-1.) *((z)-3.))

Referenced by DataPointMapping::basicInverseS().

#define q23d

(

z

)

( -1.5+(z)*(4.-1.5*(z)))

Referenced by DataPointMapping::basicInverseS().

#define q2x	(		i1,
			i2,
			i3,
			axis
	)		( b0*q1x(i1 ,i2 ,i3 ,axis) +b1*q1x(i1 ,i2+1,i3 ,axis) +b2*q1x(i1 ,i2+2,i3 ,axis) +b3*q1x(i1 ,i2+3,i3 ,axis) )

Referenced by DataPointMapping::basicInverseS().

#define q2xr	(		i1,
			i2,
			i3,
			axis
	)		( b0*q1xr(i1 ,i2 ,i3 ,axis) +b1*q1xr(i1 ,i2+1,i3 ,axis) +b2*q1xr(i1 ,i2+2,i3 ,axis) +b3*q1xr(i1 ,i2+3,i3 ,axis) )

Referenced by DataPointMapping::basicInverseS().

#define q2xs	(		i1,
			i2,
			i3,
			axis
	)		( b0r*q1xs(i1 ,i2 ,i3 ,axis) +b1r*q1xs(i1 ,i2+1,i3 ,axis) +b2r*q1xs(i1 ,i2+2,i3 ,axis) +b3r*q1xs(i1 ,i2+3,i3 ,axis))

Referenced by DataPointMapping::basicInverseS().

#define q2xt	(		i1,
			i2,
			i3,
			axis
	)		( b0*q1xt(i1 ,i2 ,i3 ,axis) +b1*q1xt(i1 ,i2+1,i3 ,axis) +b2*q1xt(i1 ,i2+2,i3 ,axis) +b3*q1xt(i1 ,i2+3,i3 ,axis) )

#define q33

(

z

)

( oneSixth*(z)*((z)-1.)*((z)-2.))

Referenced by DataPointMapping::basicInverseS().

#define q33d

(

z

)

(oneSixth*(2.+(z)*(-6.+3.*(z))))

Referenced by DataPointMapping::basicInverseS().

#define q3x	(		i1,
			i2,
			i3,
			axis
	)		( c0*q2x(i1 ,i2 ,i3 ,axis) +c1*q2x(i1 ,i2 ,i3+1,axis) +c2*q2x(i1 ,i2 ,i3+2,axis) +c3*q2x(i1 ,i2 ,i3+3,axis) )

Referenced by DataPointMapping::basicInverseS().

#define q3xr	(		i1,
			i2,
			i3,
			axis
	)		( c0*q2xr(i1 ,i2 ,i3 ,axis) +c1*q2xr(i1 ,i2 ,i3+1,axis) +c2*q2xr(i1 ,i2 ,i3+2,axis) +c3*q2xr(i1 ,i2 ,i3+3,axis) )

Referenced by DataPointMapping::basicInverseS().

#define q3xs	(		i1,
			i2,
			i3,
			axis
	)		( c0*q2xs(i1 ,i2 ,i3 ,axis) +c1*q2xs(i1 ,i2 ,i3+1,axis) +c2*q2xs(i1 ,i2 ,i3+2,axis) +c3*q2xs(i1 ,i2 ,i3+3,axis) )

Referenced by DataPointMapping::basicInverseS().

#define q3xt	(		i1,
			i2,
			i3,
			axis
	)		( c0r*q2xt(i1 ,i2 ,i3 ,axis) +c1r*q2xt(i1 ,i2 ,i3+1,axis) +c2r*q2xt(i1 ,i2 ,i3+2,axis) +c3r*q2xt(i1 ,i2 ,i3+3,axis) )

Referenced by DataPointMapping::basicInverseS().

#define R	(		i0,
			i1
	)		rp[i0+rDim0*(i1)]

#define R0	(		i0,
			i1
	)		r0p[i0+r0Dim0*(i1)]

Referenced by RevolutionMapping::basicInverse(), RevolutionMapping::basicInverseS(), DataPointMapping::basicInverseS(), ParallelOverlappingGridInterpolator::computeError(), HyperbolicMapping::createCurveOnSurface(), display(), Oges::formRhsAndSolutionVectors(), Ogen::interpolateAll(), Interpolate::interpolateCoarseFromFine(), interpolatePoints(), RevolutionMapping::map(), NormalMapping::map(), Mapping::mapGrid(), NormalMapping::mapS(), RevolutionMapping::mapS(), Regrid::merge(), Ogen::preInterpolateGrids(), Regrid::regridAligned(), Regrid::regridRotated(), Oges::storeSolutionIntoGridFunction(), and PlotIt::streamLines3d().

#define RX	(		i0,
			i1,
			i2
	)		rxp[i0+rxDim0*(i1+rxDim1*(i2))]

#define singular	(		side,
			axis
	)		singularp[(side)+2*(axis)]

Referenced by DataPointMapping::basicInverseS(), Mapping::checkMapping(), NormalMapping::initialize(), and CrossSectionMapping::mapGeneralWithPolarSingularity().

#define X	(		i0,
			i1
	)		xp[i0+xDim0*(i1)]

#define XY	(		i0,
			i1,
			i2,
			i3
	)		xyp[i0+xyDim0*(i1+xyDim1*(i2+xyDim2*(i3)))]