Polynomial.h

Go to the documentation of this file.

//            QMcBeaver
//
//         Constructed by 
//
//     Michael Todd Feldmann 
//              and 
//   David Randall "Chip" Kent IV
//
// Copyright 2000-2.  All rights reserved.
//
// drkent@users.sourceforge.net mtfeldmann@users.sourceforge.net

#ifndef Polynomial_H
#define Polynomial_H

#include "Array1D.h"
#include "FunctionR1toR1.h"
#include "Complex.h"
#include "Exception.h"

class Polynomial : public FunctionR1toR1
{
private:
  Array1D<double> coefficients;

  Array1D<double> firstDerivativeCoefficients;

  Array1D<double> secondDerivativeCoefficients;

  Array1D<double> thirdDerivativeCoefficients;

  double f;

  double df;

  double d2f;

  double d3f;

  double x;

  bool evaluatedF;

  bool evaluatedDF;

  bool evaluatedD2F;

public:
  Polynomial();

  Polynomial(Array1D<double> & coeffs);

  void initialize(Array1D<double> & coeffs);

  void evaluate(double x);
  double getFunctionValue();

  double getFunctionValue(double x);

  double getFirstDerivativeValue();
  double getSecondDerivativeValue();
  double getThirdDerivativeValue();

  double get_p_a(int ai);
  double get_p2_xa(int ai);
  double get_p3_xxa(int ai);

  Array1D<double> getCoefficients();

  Array1D<Complex> getRoots();

  void operator=(Polynomial rhs);

  int getNumberCoefficients();

  double getCoefficient(int i);

  void print(ostream & strm);

private:
  void initialize();

  double evaluate(double x, Array1D<double> & coeffs);

  void evaluateAll(double x, Array1D<double> & coeffs);

  void laguer(Array1D<Complex> &a, int m, Complex &x, int *its, bool *calcOK);

  Array1D<Complex> zroots(Array1D<Complex> & a, bool polish, bool *calcOK);
};

#endif

---