62 #ifndef G4JTPOLYNOMIALSOLVER_HH
63 #define G4JTPOLYNOMIALSOLVER_HH
83 std::vector<G4double>
p;
84 std::vector<G4double>
qp;
85 std::vector<G4double>
k;
86 std::vector<G4double>
qk;
87 std::vector<G4double>
svk;
118 std::vector<G4double> &p,
119 std::vector<G4double> &q,
static const G4double are
void ComputeFixedShiftPolynomial(G4int l2, G4int *nz)
std::vector< G4double > svk
std::vector< G4double > qk
static const G4double mre
void ComputeScalarFactors(G4int *type)
std::vector< G4double > p
static const G4double smalno
std::vector< G4double > qp
void ComputeNewEstimate(G4int type, G4double *uu, G4double *vv)
static const G4double eta
std::vector< G4double > k
static const G4double infin
void ComputeNextPolynomial(G4int *type)
void RealPolynomialIteration(G4double *sss, G4int *nz, G4int *iflag)
static const G4double base
static const char sss[MAX_N_PAR+2]
void QuadraticSyntheticDivision(G4int n, G4double *u, G4double *v, std::vector< G4double > &p, std::vector< G4double > &q, G4double *a, G4double *b)
static const double degree
void QuadraticPolynomialIteration(G4double *uu, G4double *vv, G4int *nz)
G4int FindRoots(G4double *op, G4int degree, G4double *zeror, G4double *zeroi)
void Quadratic(G4double a, G4double b1, G4double c, G4double *sr, G4double *si, G4double *lr, G4double *li)