VERB_code_2.3
rroots.h File Reference
#include <iostream>
#include <iomanip>
#include <math.h>
#include <stdlib.h>

Go to the source code of this file.

Macros

#define maxiter   5000
 maximum number of iterations
 
#define DBL_EPSILON   1e-15
 some other epsilon and stuff
 

Functions

int roots (double *a, int n, double *wr, double *wi)
 Extract individual real or complex roots from list of quadratic factors.
 
void get_quads (double *a, int n, double *quad, double *x)
 

Detailed Description

Finds all roots of polynomial by first finding quadratic factors using Bairstow's method, then extracting roots from quadratics. Implements new algorithm for managing multiple roots.

Date
(C) 2002, 2003,
Author
C. Bond. All rights reserved.

Definition in file rroots.h.

Function Documentation

void get_quads ( double *  a,
int  n,
double *  quad,
double *  x 
)

Top level routine to manage the determination of all roots of the given polynomial 'a', returning the quadratic factors (and possibly one linear factor) in 'x'.

Todo:
Should gives error in case roots not founded. But that is commented.

Definition at line 201 of file rroots.cpp.