Inheritance diagram for torus:

Collaboration diagram for torus:

Public Member Functions
	torus (double x=0, double y=0, double z=0, double aa=0.3, double cc=1.0)

virtual	~torus ()

void	project_points2geom (int dim, double oldcoords, double newcoords, double *derivs) const

Public Member Functions inherited from geomObject
virtual	~geomObject ()

virtual double	compute_projecterror (int dim, double *oldcoords) const

virtual void	compute_projecterror (int dim, int nverts, double *oldcoords, double &l1err, double &l2err, double &linferr) const

double	Twonorm (int dim, double *vec) const

Private Attributes
double	centerx

double	centery

double	centerz

double	a

double	c

Detailed Description

Definition at line 125 of file geomObject.cpp.

Constructor & Destructor Documentation

◆ torus()

torus::torus	(	double	x = `0`,
		double	y = `0`,
		double	z = `0`,
		double	aa = `0.3`,
		double	cc = `1.0`
	)

inline

Definition at line 128 of file geomObject.cpp.

  : centerx( x ), centery( y ), centerz( z ), a( aa ), c( cc )
  {
  assert( aa > 0 && cc > 0 );
  }

◆ ~torus()

virtual torus::~torus ( )

inlinevirtual

Definition at line 133 of file geomObject.cpp.

133 {}

Member Function Documentation

◆ project_points2geom()

void torus::project_points2geom	(	int	dim,
		double *	oldcoords,
		double *	newcoords,
		double *	derivs
	)		const

inlinevirtual

Implements geomObject.

Definition at line 134 of file geomObject.cpp.

  {
  assert( dim == 3 && oldcoords && newcoords );
  double transfer[3] = { oldcoords[0] - centerx, oldcoords[1] - centery, oldcoords[2] - centerz };
  double twodnrm = geomObject::Twonorm( 2, transfer );
  double tubecenter[3] = { c * transfer[0] / twodnrm + centerx, c * transfer[1] / twodnrm + centery, centerz };
  double direction[3] = { oldcoords[0] - tubecenter[0], oldcoords[1] - tubecenter[1],
  oldcoords[2] - tubecenter[2] };
  double len = geomObject::Twonorm( 3, direction );
  assert( len > 0 );
  direction[0] /= len;
  direction[1] /= len;
  direction[2] /= len;
  
  if( derivs )
  {
  derivs[0] = direction[0];
  derivs[1] = direction[1];
  derivs[2] = direction[2];
  }
  
  newcoords[0] = tubecenter[0] + a * direction[0];
  newcoords[1] = tubecenter[1] + a * direction[1];
  newcoords[2] = tubecenter[2] + a * direction[2];
  }