Changeset 508 in Sophya for trunk/SophyaLib/NTools/generalfit.cc
- Timestamp:
- Oct 25, 1999, 12:36:22 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/NTools/generalfit.cc
r490 r508 874 874 875 875 //++ 876 Vector GeneralFit::GetParm()876 OVector GeneralFit::GetParm() 877 877 // 878 878 // Retourne les valeurs des parametres dans un vecteur. … … 1096 1096 { 1097 1097 volatile double oldChi2; 1098 Matrix COVAR(mNPar,mNPar);1099 Vector DA(mNPar);1100 Vector dparam(mNPar);1101 Vector paramTry(mNPar);1102 Vector param_tr(mNPar);1103 Vector paramTry_tr(mNPar);1104 Vector step_tr(mNPar);1098 OMatrix COVAR(mNPar,mNPar); 1099 OVector DA(mNPar); 1100 OVector dparam(mNPar); 1101 OVector paramTry(mNPar); 1102 OVector param_tr(mNPar); 1103 OVector paramTry_tr(mNPar); 1104 OVector step_tr(mNPar); 1105 1105 nStop = nStopL = nStep = 0; 1106 1106 Chi2 = oldChi2 = 0.; … … 1515 1515 1516 1516 ////////////////////////////////////////////////////////////////////// 1517 void GeneralFit::write_in_step(double ci2, Vector& par)1517 void GeneralFit::write_in_step(double ci2,OVector& par) 1518 1518 { 1519 1519 if(FileStep==NULL) return; … … 1524 1524 1525 1525 ////////////////////////////////////////////////////////////////////// 1526 void GeneralFit::TryFunc( Vector& par,Vector& par_tr)1526 void GeneralFit::TryFunc(OVector& par,OVector& par_tr) 1527 1527 { 1528 1528 BETA_Try = 0; 1529 1529 ATGA_Try = 0; 1530 1530 Chi2 = 0; 1531 Vector deriv(mNPar);1532 Vector derivtr(mNPar);1531 OVector deriv(mNPar); 1532 OVector derivtr(mNPar); 1533 1533 double result; 1534 1534 … … 1564 1564 1565 1565 ////////////////////////////////////////////////////////////////////// 1566 void GeneralFit::TryXi2( Vector& par,Vector& par_tr)1566 void GeneralFit::TryXi2(OVector& par,OVector& par_tr) 1567 1567 { 1568 1568 double c, *parloc; … … 1664 1664 1665 1665 ////////////////////////////////////////////////////////////////////// 1666 Vector GeneralFit::p_vers_tr(Vector const& p)1667 { 1668 Vector tr(p);1666 OVector GeneralFit::p_vers_tr(OVector const& p) 1667 { 1668 OVector tr(p); 1669 1669 for(int i=0;i<mNPar;i++) { 1670 1670 if( fixParam[i] || ! boundParam[i] ) continue; … … 1675 1675 1676 1676 ////////////////////////////////////////////////////////////////////// 1677 void GeneralFit::p_vers_tr( Vector const& p,Vector& tr)1677 void GeneralFit::p_vers_tr(OVector const& p,OVector& tr) 1678 1678 { 1679 1679 for(int i=0;i<mNPar;i++) { … … 1695 1695 1696 1696 ////////////////////////////////////////////////////////////////////// 1697 Vector GeneralFit::tr_vers_p(Vector const& tr)1698 { 1699 Vector p(tr);1697 OVector GeneralFit::tr_vers_p(OVector const& tr) 1698 { 1699 OVector p(tr); 1700 1700 for(int i=0;i<mNPar;i++) { 1701 1701 if( fixParam[i] || ! boundParam[i] ) continue; … … 1706 1706 1707 1707 ////////////////////////////////////////////////////////////////////// 1708 void GeneralFit::tr_vers_p( Vector const& tr,Vector& p)1708 void GeneralFit::tr_vers_p(OVector const& tr,OVector& p) 1709 1709 { 1710 1710 for(int i=0;i<mNPar;i++) { … … 1727 1727 1728 1728 ////////////////////////////////////////////////////////////////////// 1729 Vector GeneralFit::dp_vers_dtr(Vector const& dp,Vector const& tr)1730 { 1731 Vector dtr(dp);1729 OVector GeneralFit::dp_vers_dtr(OVector const& dp,OVector const& tr) 1730 { 1731 OVector dtr(dp); 1732 1732 for(int i=0;i<mNPar;i++) { 1733 1733 if( fixParam[i] || ! boundParam[i] ) continue; … … 1738 1738 1739 1739 ////////////////////////////////////////////////////////////////////// 1740 void GeneralFit::dp_vers_dtr( Vector const& dp,Vector const& tr,Vector& dtr)1740 void GeneralFit::dp_vers_dtr(OVector const& dp,OVector const& tr,OVector& dtr) 1741 1741 { 1742 1742 for(int i=0;i<mNPar;i++) { … … 1759 1759 1760 1760 ////////////////////////////////////////////////////////////////////// 1761 Vector GeneralFit::dtr_vers_dp(Vector const& dtr,Vector const& tr)1762 { 1763 Vector dp(dtr);1761 OVector GeneralFit::dtr_vers_dp(OVector const& dtr,OVector const& tr) 1762 { 1763 OVector dp(dtr); 1764 1764 for(int i=0;i<mNPar;i++) { 1765 1765 if( fixParam[i] || ! boundParam[i] ) continue; … … 1771 1771 ////////////////////////////////////////////////////////////////////// 1772 1772 // inline fonction pour aller + vite dans le try() 1773 //void GeneralFit::dtr_vers_dp( Vector const& dtr,Vector const& tr,Vector& dp)1774 1775 ////////////////////////////////////////////////////////////////////// 1776 int GeneralFit::put_in_limits_for_deriv( Vector const& p,Vector& dp,double dist)1773 //void GeneralFit::dtr_vers_dp(OVector const& dtr,OVector const& tr,OVector& dp) 1774 1775 ////////////////////////////////////////////////////////////////////// 1776 int GeneralFit::put_in_limits_for_deriv(OVector const& p,OVector& dp,double dist) 1777 1777 // 1-/ Redefinit dp pour qu'il soit superieur a minStepDeriv 1778 1778 // 2-/ Redefinit dp pour que p+/-dp reste dans les limites (parametre borne)
Note:
See TracChangeset
for help on using the changeset viewer.