polynomialRootSolution.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           polynomialRootSolution.h
// author(s):      Jan Novak
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 This program is distributed in the hope that it will be useful,
//                 but WITHOUT ANY WARRANTY; without even the implied warranty of
//                 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//                 GNU Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************
#ifndef MUMECH_POL_ROOTS_SOLU_H
#define MUMECH_POL_ROOTS_SOLU_H

#include <complex>


namespace mumech {

//macros and variables definitions
#define _MACHINE_EPS_   1E-100
#define _sgrt3div2_     0.8660254037844386

//definition of the complex type
typedef struct{
  double  a, b;
} _COMPLEX;

class polynomialRootSolution
{
 public:
  polynomialRootSolution () {}
  virtual ~polynomialRootSolution () {}
  
  
  std::complex<double> * GetCubicPolyRoots( double a, double b, double c, std::complex<double> * roots );
  
  
  std::complex<double>* GetCubicPolyRoots( std::complex<double> a, std::complex<double> b, std::complex<double> c,
                       std::complex<double>* roots );
  
  
  static void GauLegF(double *x, double *w, const int nPoints, const double minLim, const double maxLim);
  static void GauHerQ(double *x, double *w, const int nPoints);
  
};

} // end of namespace mumech

#endif // end of MUMECH_POL_ROOTS_SOLU_H

/*end of file*/