Changeset 490 in Sophya for trunk/SophyaLib/NTools/generalfit.cc


Ignore:
Timestamp:
Oct 21, 1999, 5:25:53 PM (26 years ago)
Author:
ansari
Message:

Merge avec PEIDA++ (~V 3.8) et nettoyage pour nouveau PPersist Reza+cmv 21/10/99

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaLib/NTools/generalfit.cc

    r307 r490  
    77#include <string>
    88
    9 #if defined(__KCC__)
    10 using std::string ;
    11 #endif
    12 
    139#include "pexceptions.h"
    1410#include "generalfit.h"
    … …  
    4137  : mNVar(nVar), mNPar(nPar)
    4238{
    43  DBASSERT( nVar > 0 && nPar > 0 );
     39 ASSERT( nVar > 0 && nPar > 0 );
    4440 deltaParm = new double[nPar];
    4541 tmpParm   = new double[nPar];
    … …  
    8682//--
    8783{
    88  DBASSERT(numPar >= 0 && numPar < mNPar);
     84 ASSERT(numPar >= 0 && numPar < mNPar);
    8985 deltaParm[numPar] = d;
    9086}
    … …  
    226222  : mNPar(nPar)
    227223{
    228  DBASSERT( nPar>0 );
     224 ASSERT( nPar>0 );
    229225 deltaParm = new double[nPar];
    230226 END_CONSTRUCTOR
    … …  
    307303//--
    308304{
    309  DBASSERT(numPar >= 0 && numPar < mNPar);
     305 ASSERT(numPar >= 0 && numPar < mNPar);
    310306 
    311307 deltaParm[numPar] = d;
    … …  
    375371    D             (f->NPar())
    376372{
    377  DBASSERT(mNVar>0 && mNPar>0);
    378  DBASSERT(mNPar<1000000);
     373 ASSERT(mNVar>0 && mNPar>0);
     374 ASSERT(mNPar<1000000);
    379375
    380376 TRY {
    … …  
    418414    D             (f->NPar())
    419415{
    420  DBASSERT( mNPar>0 );
    421  DBASSERT( mNPar < 1000000 );
     416 ASSERT( mNPar>0 );
     417 ASSERT( mNPar < 1000000 );
    422418
    423419 TRY {
    … …  
    573569//--
    574570{
    575  DBASSERT(n>=0 && n<mNPar);
     571 ASSERT(n>=0 && n<mNPar);
    576572 Eps(n) = (ep<=0.) ? EPS_FIT_MIN: ep;
    577573 if(debugLevel>0) cout<<"SetEps("<<n<<") = "<<Eps(n)<<endl;
    … …  
    616612//--
    617613{
    618  DBASSERT( mFuncXi2  == NULL );
    619  DBASSERT( f != NULL );
    620  DBASSERT( f->NVar() == mNVar );
    621  DBASSERT( f->NPar() == mNPar );
     614 ASSERT( mFuncXi2  == NULL );
     615 ASSERT( f != NULL );
     616 ASSERT( f->NVar() == mNVar );
     617 ASSERT( f->NPar() == mNPar );
    622618 mFunction = f;
    623619 if(debugLevel>0) cout<<"SetFunction "<<mFunction<<endl;
    … …  
    632628//--
    633629{
    634  DBASSERT( mFunction == NULL );
    635  DBASSERT( f != NULL );
    636  DBASSERT( f->NPar() == mNPar );
     630 ASSERT( mFunction == NULL );
     631 ASSERT( f != NULL );
     632 ASSERT( f->NPar() == mNPar );
    637633 mFuncXi2  = f;
    638634 if(debugLevel>0) cout<<"SetFuncXi2 "<<mFuncXi2<<endl;
    … …  
    646642//--
    647643{
    648  DBASSERT( data->NVar()==mNVar || mFunction==NULL );
     644 ASSERT( data->NVar()==mNVar );
    649645 mData = data;
    650646 mNddl = mData->NDataGood() - mNParFree;
    … …  
    661657//--
    662658{
    663  DBASSERT(n>=0 && n<mNPar);
     659 ASSERT(n>=0 && n<mNPar);
    664660
    665661 Param(n)     = value;
    … …  
    688684//--
    689685{
    690  DBASSERT(n>=0 && n<mNPar);
     686 ASSERT(n>=0 && n<mNPar);
    691687 SetParam(n,value,step,min,max);
    692688 nameParam[n] = name;
    … …  
    700696//--
    701697{
    702  DBASSERT(n>=0 && n<mNPar);
     698 ASSERT(n>=0 && n<mNPar);
    703699 Param(n) = value;
    704700 if(debugLevel) {cout<<"Set_Param "; PrintParm(n);}
    … …  
    713709//--
    714710{
    715  DBASSERT(n>=0 && n<mNPar);
     711 ASSERT(n>=0 && n<mNPar);
    716712 if(step>0.) {
    717713   if( fixParam[n] ) { fixParam[n]=0; mNParFree++;}
    … …  
    734730//--
    735731{
    736  DBASSERT(i>=0 && i<mNPar);
     732 ASSERT(i>=0 && i<mNPar);
    737733 if(val<0.) minStepDeriv(i) = Eps(i);
    738734   else     minStepDeriv(i) = val;
    … …  
    759755//--
    760756{
    761  DBASSERT(n>=0 && n<mNPar && max>min);
     757 ASSERT(n>=0 && n<mNPar && max>min);
    762758
    763759 minParam(n)  = min;
    … …  
    778774//--
    779775{
    780  DBASSERT(n>=0 && n<mNPar && maxParam(n)>minParam(n));
     776 ASSERT(n>=0 && n<mNPar && maxParam(n)>minParam(n));
    781777 SetBound(n,minParam(n),maxParam(n));
    782778}
    … …  
    788784//--
    789785{
    790  DBASSERT(n>=0 && n<mNPar);
     786 ASSERT(n>=0 && n<mNPar);
    791787
    792788 if( boundParam[n] ) {
    … …  
    814810//--
    815811{
    816  DBASSERT(n>=0 && n<mNPar);
     812 ASSERT(n>=0 && n<mNPar);
    817813
    818814 Param(n) = v;
    … …  
    833829//--
    834830{
    835  DBASSERT(n>=0 && n<mNPar);
     831 ASSERT(n>=0 && n<mNPar);
    836832 SetFix(n,Param(n));
    837833}
    … …  
    843839//--
    844840{
    845  DBASSERT(n>=0 && n<mNPar);
     841 ASSERT(n>=0 && n<mNPar);
    846842
    847843 if( fixParam[n] ) {
    … …  
    873869//--
    874870{
    875  DBASSERT(n>=0 && n<mNPar);
     871 ASSERT(n>=0 && n<mNPar);
    876872 return Param(n);
    877873}
    … …  
    892888//--
    893889{
    894  DBASSERT(n>=0 && n<mNPar);
     890 ASSERT(n>=0 && n<mNPar);
    895891 return errParam(n);
    896892}
    … …  
    902898//--
    903899{
    904  DBASSERT(i>=0 && i<mNPar && j>=0 && j<mNPar);
     900 ASSERT(i>=0 && i<mNPar && j>=0 && j<mNPar);
    905901 return ATGA(i,j);
    906902}
    … …  
    913909//--
    914910{
    915  DBASSERT(n>=0 && n<mNPar);
     911 ASSERT(n>=0 && n<mNPar);
    916912 return stepParam(n);
    917913}
    … …  
    923919//--
    924920{
    925  DBASSERT(n>=0 && n<mNPar);
     921 ASSERT(n>=0 && n<mNPar);
    926922 return maxParam(n);
    927923}
    … …  
    933929//--
    934930{
    935  DBASSERT(n>=0 && n<mNPar);
     931 ASSERT(n>=0 && n<mNPar);
    936932 return minParam(n);
    937933}
    … …  
    986982//--
    987983{
    988  DBASSERT(n>=0 && n<mNPar);
     984 ASSERT(n>=0 && n<mNPar);
    989985
    990986 cout<<"Par["<<n<<"] "<<nameParam[n]
    … …  
    13961392//--
    13971393{
    1398 double c2 = -.1;
     1394double c2 = -1.;
    13991395if(par==NULL) par = Param.Data();
    14001396if( mData->NData() <= 0 ) {nddl = -100; return 0.;}
    … …  
    16221618void GeneralFit::CheckSanity()
    16231619{
    1624   DBASSERT( mData != NULL );
    1625   DBASSERT( mFunction != NULL || mFuncXi2 != NULL );
     1620  ASSERT( mData != NULL );
     1621  ASSERT( mFunction != NULL || mFuncXi2 != NULL );
    16261622  if( mFunction != NULL ) {
    1627     DBASSERT( mFunction->NVar() == mNVar );
    1628     DBASSERT( mData->NVar() == mNVar );
     1623    ASSERT( mFunction->NVar() == mNVar );
     1624    ASSERT( mData->NVar() == mNVar );
    16291625  }
    1630   DBASSERT( mNParFree > 0 && mNParFree <= mNPar );
    1631   DBASSERT( mNParBound >= 0 && mNParBound <= mNPar );
    1632   DBASSERT( mNParFree <= mData->NDataGood() );
     1626  ASSERT( mNParFree > 0 && mNParFree <= mNPar );
     1627  ASSERT( mNParBound >= 0 && mNParBound <= mNPar );
     1628  ASSERT( mNParFree <= mData->NDataGood() );
    16331629}
    16341630
    … …  
    16381634// D = (max-min)/Pi
    16391635{
    1640   // DBASSERT(i>=0 && i<mNPar);
     1636  // ASSERT(i>=0 && i<mNPar);
    16411637 C(i) = D(i) = 0.;
    16421638 if( !boundParam[i] || fixParam[i] ) return;
    … …  
    16611657// tr = tan( (p-C)/D )
    16621658{
    1663  // DBASSERT(i>=0 && i<mNPar);
     1659 // ASSERT(i>=0 && i<mNPar);
    16641660 double tr = p;
    16651661 if(boundParam[i]) tr = tan((p-C(i))/D(i));
    … …  
    16921688// p = C+D*atan(tr)
    16931689{
    1694  // DBASSERT(i>=0 && i<mNPar);
     1690 // ASSERT(i>=0 && i<mNPar);
    16951691 double p = tr;
    16961692 if(boundParam[i]) p = C(i)+D(i)*atan(tr);
    … …  
    17241720// attention: df/dp = (1+tr**2)/D * dF/dtr = coeff * dF/dtr
    17251721{
    1726  // DBASSERT(i>=0 && i<mNPar);
     1722 // ASSERT(i>=0 && i<mNPar);
    17271723 double coeff = 1.;
    17281724 if(boundParam[i]) coeff = (1.+tr*tr)/D(i);
    … …  
    17561752// attention: df/dtr = D/(1+tr**2) * dF/dp = coeff * dF/dp
    17571753{
    1758   // DBASSERT(i>=0 && i<mNPar);
     1754  // ASSERT(i>=0 && i<mNPar);
    17591755 double coeff = 1.;
    17601756 if(boundParam[i]) coeff = D(i)/(1.+tr*tr);
Note: See TracChangeset for help on using the changeset viewer.