G4SimplexDownhill< T > Class Template Reference

#include <G4SimplexDownhill.hh>

Public Member Functions
	G4SimplexDownhill (T *tp, G4int n)
	~G4SimplexDownhill ()
G4double	GetMinimum ()
std::vector< G4double >	GetMinimumPoint ()

---

Detailed Description

template<class T>
class G4SimplexDownhill< T >

Definition at line 49 of file G4SimplexDownhill.hh.

---

Constructor & Destructor Documentation

template<class T>

G4SimplexDownhill< T >::G4SimplexDownhill	(	T *	tp,
		G4int	n
	)			[inline]

Definition at line 54 of file G4SimplexDownhill.hh.

        : currentValue(0.), target(tp), numberOfVariable(n)
      { init(); }

template<class T>

G4SimplexDownhill< T >::~G4SimplexDownhill

(

)

---

Member Function Documentation

template<class T>

G4double G4SimplexDownhill< T >::GetMinimum

(

)

Definition at line 64 of file G4SimplexDownhill.icc.

Referenced by G4SimplexDownhill< T >::GetMinimumPoint().

{

   initialize();

// First Tryal;
  
   //G4cout << "Begin First Trials" << G4endl;
   doDownhill();
   //G4cout << "End First Trials" << G4endl;

   std::vector< G4double >::iterator it_minh =
      std::min_element( currentHeights.begin() , currentHeights.end() );
   G4int imin = -1;
   G4int i = 0;
   for ( std::vector< G4double >::iterator it = currentHeights.begin();
         it != currentHeights.end(); it++ )
   {
      if ( it == it_minh )
      {
         imin = i;
      }
      i++;
   }
   minimumPoint = currentSimplex[ imin ];

// Second Trial

   //std::vector< G4double > minimumPoint = currentSimplex[ 0 ];
   initialize();

   currentSimplex[ numberOfVariable ] = minimumPoint;

   //G4cout << "Begin Second Trials" << G4endl;
   doDownhill();
   //G4cout << "End Second Trials" << G4endl;
   
   G4double sum = std::accumulate( currentHeights.begin() ,
                                   currentHeights.end() , 0.0 );
   G4double average = sum/(numberOfVariable+1); 
   G4double minimum = average;

   minimized = true;

   return minimum;

}

template<class T>

std::vector< G4double > G4SimplexDownhill< T >::GetMinimumPoint

(

)

Definition at line 419 of file G4SimplexDownhill.icc.

References G4SimplexDownhill< T >::GetMinimum().

{
   if ( minimized != true )
   {
      GetMinimum();
   }

   std::vector< G4double >::iterator it_minh =
     std::min_element( currentHeights.begin() , currentHeights.end() );;
   G4int imin = -1;
   G4int i = 0;
   for ( std::vector< G4double >::iterator 
         it = currentHeights.begin(); it != currentHeights.end(); it++ )
   {
      if ( it == it_minh )
      {
         imin = i;
      }
      i++;
   }
   minimumPoint = currentSimplex[ imin ];

   return minimumPoint; 
}

---

The documentation for this class was generated from the following files:

• G4SimplexDownhill.hh
• G4SimplexDownhill.icc

---