Reconstruction/MdcPatRec/MdcRecoUtil/MdcRecoUtil-00-01-08/MdcRecoUtil/DifNumber.h File Reference

#include <math.h>
#include <assert.h>
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "MdcRecoUtil/Code.h"
#include <iosfwd>
#include "MdcRecoUtil/DifNumber.icc"

Go to the source code of this file.

Classes
class	DifNumber

Macros
#define	MATRIX_BOUND_CHECK

Functions
DifNumber	solveQuad (const DifNumber &a, const DifNumber &b, const DifNumber &c, int pref, Code &code)

Macro Definition Documentation

MATRIX_BOUND_CHECK

#define MATRIX_BOUND_CHECK

Definition at line 27 of file Reconstruction/MdcPatRec/MdcRecoUtil/MdcRecoUtil-00-01-08/MdcRecoUtil/DifNumber.h.

Function Documentation

solveQuad()

DifNumber solveQuad	(	const DifNumber &	a,
		const DifNumber &	b,
		const DifNumber &	c,
		int	pref,
		Code &	code
	)

Definition at line 90 of file DifNumber.cxx.

{
  DifNumber descr=b*b-4.0*a*c;
  if(descr<0.0) {       // solution not real
    code.setFail(1341);
    return DifNumber(0.0);
  }
  if(a.number()==0.0){
    if(b.number()==0.0) {
      code.setFail(1342);
      return DifNumber(0.0);
    }
    code.setSuccess(40);
    return -c/b+a*c/pow(b,3);
  }
  code.setSuccess(40);
  descr=sqrt(descr);
  DifNumber s=-b;
 
  if(pref==+1) {        // positive solution
    s+=descr;
  }else if(pref==-1){       // negative solution
    s-=descr;
  }else if(pref==0) {       // smallest solution
    if(s>0.0) {s-=descr;}else {s+=descr;}
  }else {           // illegal prefrence
    code.setFail(1343);
    return DifNumber(0.0);
  }
  s/=2.0*a;
  return s;
}