#include <Trapezoid2RectangleMappingX.h>

Public Member Functions
void	rectangle (double xtrap, double ytrap, double &xrec, double &yrec) const

void	trapezoid (double xrec, double yrec, double &xtrap, double &ytrap) const

	Trapezoid2RectangleMappingX ()

	Trapezoid2RectangleMappingX (double x0, double y0, double bovera, double h)
	normal trapezoid case, a != b More...

	Trapezoid2RectangleMappingX (double x0, double y0, double k)
	special parallelogram case, a == b. The meaning of k changes. More...

Private Attributes
double	k_

bool	parallel_

double	x0_

double	y0_

Detailed Description

Maps a trapezoidal coordinate system into a cartesian one. It is assumed that x is the coordinate along the trapezoid bases while y is along the trapezoid height.

Date:: 2010/10/18 22:34:16

Revision:: 1.3

Author: T. Todorov

Definition at line 22 of file Trapezoid2RectangleMappingX.h.

Constructor & Destructor Documentation

Trapezoid2RectangleMappingX::Trapezoid2RectangleMappingX ( )

inline

Definition at line 25 of file Trapezoid2RectangleMappingX.h.

25 {}

Trapezoid2RectangleMappingX::Trapezoid2RectangleMappingX	(	double	x0,
		double	y0,
		double	bovera,
		double	h
	)

inline

normal trapezoid case, a != b

Definition at line 28 of file Trapezoid2RectangleMappingX.h.

References gather_cfg::cout, and k_.

                                                                               :
     x0_(x0), y0_(y0), parallel_(false)
   {
     k_ = 2/h * (bovera-1.) / (bovera+1.);
 
 #ifdef DEBUG_GRID_TRM
     std::cout << "Trapezoid2RectangleMappingX constructed with x0,y0 " << x0 << " " << y0 
      << " b/a= " << bovera << " h= " << h << std::endl;
 #endif
   }

Trapezoid2RectangleMappingX::Trapezoid2RectangleMappingX	(	double	x0,
		double	y0,
		double	k
	)

inline

special parallelogram case, a == b. The meaning of k changes.

Definition at line 40 of file Trapezoid2RectangleMappingX.h.

References gather_cfg::cout.

                                                                :
     x0_(x0), y0_(y0), k_(k), parallel_(true)
   {
 #ifdef DEBUG_GRID_TRM
     std::cout << "Trapezoid2RectangleMappingX constructed with x0,y0 " << x0 << " " << y0 
      << " k= " << k << std::endl;
 #endif
   }

Member Function Documentation

void Trapezoid2RectangleMappingX::rectangle	(	double	xtrap,
		double	ytrap,
		double &	xrec,
		double &	yrec
	)		const

inline

Definition at line 49 of file Trapezoid2RectangleMappingX.h.

References gather_cfg::cout, k_, parallel_, x0_, and y0_.

Referenced by TrapezoidalCylindricalMFGrid::toGridFrame(), TrapezoidalCartesianMFGrid::toGridFrame(), TrapezoidalCartesianMFGrid::TrapezoidalCartesianMFGrid(), TrapezoidalCylindricalMFGrid::TrapezoidalCylindricalMFGrid(), and TrapezoidalCartesianMFGrid::uncheckedValueInTesla().

                                             {
 
     yrec = ytrap - y0_;
     if (!parallel_) xrec = (xtrap - x0_) / (1.+ yrec*k_);
     else            xrec = xtrap - x0_ + k_*yrec;
 
 #ifdef DEBUG_GRID
     std::cout << xtrap << " " << ytrap << " transformed to rectangle " << xrec << " " << yrec << std::endl;
 #endif
   }

void Trapezoid2RectangleMappingX::trapezoid	(	double	xrec,
		double	yrec,
		double &	xtrap,
		double &	ytrap
	)		const

inline

Definition at line 61 of file Trapezoid2RectangleMappingX.h.

References gather_cfg::cout, k_, parallel_, x0_, and y0_.

Referenced by TrapezoidalCylindricalMFGrid::fromGridFrame(), and TrapezoidalCartesianMFGrid::fromGridFrame().

                                               {
     if (!parallel_) xtrap = x0_ + xrec * (1.+ yrec*k_);
     else            xtrap = x0_ + xrec - k_*yrec;
     ytrap = y0_ + yrec;
 
 #ifdef DEBUG_GRID
     std::cout << xrec << " " << yrec << " transformed to trapezoid " << xtrap << " " << ytrap << std::endl;
 #endif
   }