MATH_QuadraticPolynom.hpp

Go to the documentation of this file.

#ifndef MATH_QuadraticPolynom_HPP
#define MATH_QuadraticPolynom_HPP


template<class T>
void MATH_QuadraticPolynom::solve(const tReal& a,const tReal& b,const tReal& c,
                                  int& n, T solutions[3])  {
  
    n=1;
    if (a==0) {
        if (b==0) {
            if (c==0) {
                n=1;
                solutions[0]=0;
            } else {
                n=0;
            }
        } else {
            n=1;
            solutions[0]=-c/b;
        }
        return ;
    }
    double delta=b*b-4*a*c;
    if (delta<0) {
        n=0;
    } else if (delta>0) {
        delta=sqrt(delta);
        n=2;
        if (a>0) {
            solutions[0]=(-b-delta)/(2*a);
            solutions[1]=(-b+delta)/(2*a);
        } else {
            solutions[1]=(-b-delta)/(2*a);
            solutions[0]=(-b+delta)/(2*a);
        }
    } else {
        n=1;
        solutions[0]=-b/(2*a);
    }
}

#endif