DeterministicAnnealing Class Reference

#include <DeterministicAnnealing.h>

Inheritance diagram for DeterministicAnnealing:

Classes
class	A

Public Member Functions
void	anneal ()
DeterministicAnnealing *	clone () const
double	currentTemp () const
double	cutoff () const
void	debug () const
	DeterministicAnnealing (float cutoff=3.0)
	DeterministicAnnealing (const std::vector< float > &sched, float cutoff=3.0)
double	initialTemp () const
bool	isAnnealed () const
double	phi (double chi2) const
void	resetAnnealing ()
double	weight (double chi2) const
Public Member Functions inherited from AnnealingSchedule
virtual	~AnnealingSchedule ()

Private Attributes
double	theChi2cut
unsigned int	theIndex
bool	theIsAnnealed
std::vector< float >	theTemperatures

Detailed Description

Definition at line 7 of file DeterministicAnnealing.h.

Constructor & Destructor Documentation

DeterministicAnnealing::DeterministicAnnealing

(

float

cutoff = 3.0

)

Definition at line 9 of file DeterministicAnnealing.cc.

Referenced by clone().

                                                              :
  theTemperatures({256,64,16,4,2,1}),

DeterministicAnnealing::DeterministicAnnealing	(	const std::vector< float > &	sched,
		float	cutoff = 3.0
	)

Definition at line 15 of file DeterministicAnnealing.cc.

                                                                                             :
  theTemperatures(sched),theIndex(0), theChi2cut ( cutoff*cutoff ), theIsAnnealed ( false )
{
}

Member Function Documentation

void DeterministicAnnealing::anneal ( )

virtual

Implements AnnealingSchedule.

Definition at line 20 of file DeterministicAnnealing.cc.

References theIndex, theIsAnnealed, and theTemperatures.

{
  if ( theIndex < ( theTemperatures.size() - 1 ) )
  {
    theIndex++; 
  } else {
    theIsAnnealed = true;
  };
}

DeterministicAnnealing* DeterministicAnnealing::clone ( void  ) const

inlinevirtual

Implements AnnealingSchedule.

Definition at line 51 of file DeterministicAnnealing.h.

References DeterministicAnnealing().

  {
    return new DeterministicAnnealing ( * this );
  };

double DeterministicAnnealing::currentTemp ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 63 of file DeterministicAnnealing.cc.

References theIndex, and theTemperatures.

{
  return theTemperatures[theIndex];
}

double DeterministicAnnealing::cutoff ( ) const

virtual

Returns phi(chi2) / ( phi(cutoff^2) + sum_i { phi(chi2s[i]) } )

Implements AnnealingSchedule.

Definition at line 58 of file DeterministicAnnealing.cc.

References mathSSE::sqrt(), and theChi2cut.

{
  return sqrt(theChi2cut);
}

void DeterministicAnnealing::debug ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 78 of file DeterministicAnnealing.cc.

References gather_cfg::cout, i, and theTemperatures.

{
  cout << "[DeterministicAnnealing] schedule=";
  for ( vector< float >::const_iterator i=theTemperatures.begin(); 
        i!=theTemperatures.end() ; ++i )
  {
    cout << *i << " ";
  };
  cout << endl;
}

double DeterministicAnnealing::initialTemp ( ) const

virtual

Implements AnnealingSchedule.

Definition at line 68 of file DeterministicAnnealing.cc.

References theTemperatures.

{
  return theTemperatures[0];
}

bool DeterministicAnnealing::isAnnealed ( ) const

virtual

is it annealed yet?

Implements AnnealingSchedule.

Definition at line 73 of file DeterministicAnnealing.cc.

References theIsAnnealed.

{
  return theIsAnnealed;
}

double DeterministicAnnealing::phi ( double  chi2 ) const

inlinevirtual

phi ( chi2 ) = e^( -.5*chi2 / T )

Implements AnnealingSchedule.

Definition at line 53 of file DeterministicAnnealing.cc.

References create_public_lumi_plots::exp, theIndex, and theTemperatures.

Referenced by weight().

{
  return exp ( -.5 * chi2 / theTemperatures[theIndex] );
}

void DeterministicAnnealing::resetAnnealing ( )

virtual

Implements AnnealingSchedule.

Definition at line 47 of file DeterministicAnnealing.cc.

References theIndex, and theIsAnnealed.

{
  theIndex=0;
  theIsAnnealed = false;
}

double DeterministicAnnealing::weight ( double  chi2 ) const

virtual

Returns phi(chi2) / ( phi(cutoff^2) + phi(chi2) ),

Implements AnnealingSchedule.

Definition at line 30 of file DeterministicAnnealing.cc.

References CommonUtils::isinf(), phi(), and theChi2cut.

Referenced by cuy.ValElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), and cuy.graphElement::__init__().

{
  long double mphi = phi ( chi2 );
  /*
  if ( mphi < std::numeric_limits<double>::epsilon() ) return 0.;
  return 1. / ( 1. + phi ( theChi2cut * theChi2cut ) / mphi );
  */
  // return mphi / ( mphi + phi ( theChi2cut ) );
  long double newtmp = mphi / ( mphi + phi ( theChi2cut ) );
  if ( std::isinf(newtmp ) )
  {
    if ( chi2 < theChi2cut ) newtmp=1.;
    else newtmp=0.;
  }
  return newtmp;
}

Member Data Documentation

double DeterministicAnnealing::theChi2cut

private

Definition at line 59 of file DeterministicAnnealing.h.

Referenced by cutoff(), and weight().

unsigned int DeterministicAnnealing::theIndex

private

Definition at line 58 of file DeterministicAnnealing.h.

Referenced by anneal(), currentTemp(), phi(), and resetAnnealing().

bool DeterministicAnnealing::theIsAnnealed

private

Definition at line 60 of file DeterministicAnnealing.h.

Referenced by anneal(), isAnnealed(), and resetAnnealing().

std::vector<float> DeterministicAnnealing::theTemperatures

private

Definition at line 54 of file DeterministicAnnealing.h.

Referenced by anneal(), currentTemp(), debug(), initialTemp(), and phi().