#include <G4JTPolynomialSolver.hh>

Public Member Functions
	G4JTPolynomialSolver ()

	~G4JTPolynomialSolver ()

G4int	FindRoots (G4double op, G4int degree, G4double zeror, G4double *zeroi)

Private Member Functions
void	Quadratic (G4double a, G4double b1, G4double c, G4double sr, G4double si, G4double lr, G4double li)

void	ComputeFixedShiftPolynomial (G4int l2, G4int *nz)

void	QuadraticPolynomialIteration (G4double uu, G4double vv, G4int *nz)

void	RealPolynomialIteration (G4double sss, G4int nz, G4int *iflag)

void	ComputeScalarFactors (G4int *type)

void	ComputeNextPolynomial (G4int *type)

void	ComputeNewEstimate (G4int type, G4double uu, G4double vv)

void	QuadraticSyntheticDivision (G4int n, G4double u, G4double v, std::vector< G4double > &p, std::vector< G4double > &q, G4double a, G4double b)

Private Attributes
std::vector< G4double >	p

std::vector< G4double >	qp

std::vector< G4double >	k

std::vector< G4double >	qk

std::vector< G4double >	svk

G4double	sr

G4double	si

G4double	u

G4double	v

G4double	a

G4double	b

G4double	c

G4double	d

G4double	a1

G4double	a3

G4double	a7

G4double	e

G4double	f

G4double	g

G4double	h

G4double	szr

G4double	szi

G4double	lzr

G4double	lzi

G4int	n

Static Private Attributes
static const G4double	base = 2

static const G4double	eta = DBL_EPSILON

static const G4double	infin = DBL_MAX

static const G4double	smalno = DBL_MIN

static const G4double	are = DBL_EPSILON

static const G4double	mre = DBL_EPSILON

static const G4double	lo = DBL_MIN/DBL_EPSILON

Detailed Description

Definition at line 70 of file G4JTPolynomialSolver.hh.

Constructor & Destructor Documentation

G4JTPolynomialSolver::G4JTPolynomialSolver ( )

Definition at line 48 of file G4JTPolynomialSolver.cc.

G4JTPolynomialSolver::~G4JTPolynomialSolver ( )

Definition at line 58 of file G4JTPolynomialSolver.cc.

Member Function Documentation

void G4JTPolynomialSolver::ComputeFixedShiftPolynomial	(	G4int	l2,
		G4int *	nz
	)

private

Definition at line 311 of file G4JTPolynomialSolver.cc.

References a, b, ComputeNewEstimate(), ComputeNextPolynomial(), ComputeScalarFactors(), k, n, p, qp, QuadraticPolynomialIteration(), QuadraticSyntheticDivision(), RealPolynomialIteration(), sr, svk, u, and v.

Referenced by FindRoots().

Here is the call graph for this function:

Here is the caller graph for this function:

void G4JTPolynomialSolver::ComputeNewEstimate	(	G4int	type,
		G4double *	uu,
		G4double *	vv
	)

private

Definition at line 779 of file G4JTPolynomialSolver.cc.

References a, a1, a3, a4, a5, a7, b, b1, b2, c, c1, c2, c3, c4, d, f, g, h, k, n, p, u, and v.

Referenced by ComputeFixedShiftPolynomial(), and QuadraticPolynomialIteration().

Here is the caller graph for this function:

void G4JTPolynomialSolver::ComputeNextPolynomial ( G4int * type )

private

Definition at line 734 of file G4JTPolynomialSolver.cc.

References a, a1, a3, a7, b, eta, k, n, qk, and qp.

Referenced by ComputeFixedShiftPolynomial(), and QuadraticPolynomialIteration().

Here is the caller graph for this function:

void G4JTPolynomialSolver::ComputeScalarFactors ( G4int * type )

private

Definition at line 692 of file G4JTPolynomialSolver.cc.

References a, a1, a3, a7, b, c, d, e, eta, f, g, h, k, n, qk, QuadraticSyntheticDivision(), u, and v.

Referenced by ComputeFixedShiftPolynomial(), and QuadraticPolynomialIteration().

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4JTPolynomialSolver::FindRoots	(	G4double *	op,
		G4int	degree,
		G4double *	zeror,
		G4double *	zeroi
	)

Definition at line 62 of file G4JTPolynomialSolver.cc.

References base, bb, ComputeFixedShiftPolynomial(), deg, eta, factor, G4Exp(), G4Log(), G4Pow::GetInstance(), infin, k, lo, lzi, lzr, G4INCL::Math::max(), G4INCL::Math::min(), n, p, G4Pow::powN(), qk, qp, Quadratic(), si, smalno, sr, svk, szi, szr, u, and v.

Referenced by G4TwistBoxSide::DistanceToSurface(), G4TwistTrapParallelSide::DistanceToSurface(), G4TwistTrapAlphaSide::DistanceToSurface(), and G4Torus::TorusRootsJT().

Here is the call graph for this function:

Here is the caller graph for this function:

void G4JTPolynomialSolver::Quadratic	(	G4double	a,
		G4double	b1,
		G4double	c,
		G4double *	sr,
		G4double *	si,
		G4double *	lr,
		G4double *	li
	)

private

Definition at line 848 of file G4JTPolynomialSolver.cc.

References b1, and bb.

Referenced by FindRoots(), and QuadraticPolynomialIteration().

Here is the caller graph for this function:

void G4JTPolynomialSolver::QuadraticPolynomialIteration	(	G4double *	uu,
		G4double *	vv,
		G4int *	nz
	)

private

Definition at line 480 of file G4JTPolynomialSolver.cc.

References a, are, b, ComputeNewEstimate(), ComputeNextPolynomial(), ComputeScalarFactors(), eta, lzi, lzr, mre, n, p, qp, Quadratic(), QuadraticSyntheticDivision(), szi, szr, u, and v.

Referenced by ComputeFixedShiftPolynomial().

Here is the call graph for this function:

Here is the caller graph for this function:

void G4JTPolynomialSolver::QuadraticSyntheticDivision	(	G4int	n,
		G4double *	u,
		G4double *	v,
		std::vector< G4double > &	p,
		std::vector< G4double > &	q,
		G4double *	a,
		G4double *	b
	)