/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/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
Defines
#define	MATRIX_BOUND_CHECK
Functions
DifNumber	solveQuad (const DifNumber &a, const DifNumber &b, const DifNumber &c, int pref, Code &code)

---

Define Documentation

#define MATRIX_BOUND_CHECK

Definition at line 27 of file DifNumber.h.

---

Function Documentation

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

Definition at line 91 of file DifNumber.cxx.

References DifNumber::number(), s, Code::setFail(), and Code::setSuccess().

{
  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;
}

---