Changeset 8578 for trunk/MagicSoft/Mars/mbase
- Timestamp:
- 06/18/07 15:29:32 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mbase/MMath.cc
r8563 r8578 1 1 /* ======================================================================== *\ 2 ! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.3 3 2007-06-16 21:50:26tbretz Exp $2 ! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.34 2007-06-18 14:29:32 tbretz Exp $ 3 3 ! -------------------------------------------------------------------------- 4 4 ! … … 691 691 const Double_t sqrtd = TMath::Sqrt(D); 692 692 693 const Double_t S = MMath::Sqrt3(R + sqrtd); 694 const Double_t T = MMath::Sqrt3(R - sqrtd); 695 696 x1 = (S+T) - a/3.; // real root 693 const Double_t A = TMah::Sign(1., R)*MMath::Sqrt3(TMath::Abs(R)+sqrtd); 694 695 // The case A==0 cannot happen. This would imply D==0 696 // if (A==0) 697 // { 698 // x1 = -a/3; 699 // return 1; 700 // } 701 702 x1 = (A+Q/A)-a/3; 703 704 //const Double_t S = MMath::Sqrt3(R + sqrtd); 705 //const Double_t T = MMath::Sqrt3(R - sqrtd); 706 //x1 = (S+T) - a/3.; // real root 697 707 698 708 return 1;
Note:
See TracChangeset
for help on using the changeset viewer.