#include <G4DNASmoluchowskiDiffusion.hh>

Public Types
enum	PreviousAction { IncreaseProba, DecreaseProba, Undefined }

Public Member Functions
	BoundingBox (double xmin, double xmax, double toleranceY)

void	Print ()

bool	Propose (double proposedXValue, double proposedProba, double nextProba, double &returnedValue)

Public Attributes
double	fXmax

double	fXmin

double	fXmaxDef

double	fXminDef

double	fToleranceY

double	fSum

double	fIncreasingCumulativeFunction

PreviousAction	fPreviousAction

Detailed Description

Definition at line 138 of file G4DNASmoluchowskiDiffusion.hh.

Member Enumeration Documentation

enum G4DNASmoluchowskiDiffusion::BoundingBox::PreviousAction

Enumerator
IncreaseProba
DecreaseProba
Undefined

Definition at line 148 of file G4DNASmoluchowskiDiffusion.hh.

     {
       IncreaseProba,
       DecreaseProba,
       Undefined
     };

Constructor & Destructor Documentation

G4DNASmoluchowskiDiffusion::BoundingBox::BoundingBox	(	double	xmin,
		double	xmax,
		double	toleranceY
	)

inline

Definition at line 157 of file G4DNASmoluchowskiDiffusion.hh.

                                    :
      fXmax(xmax), fXmin(xmin),
      fToleranceY(toleranceY),
      fSum(0)
     {
       if(fXmax < fXmin)
       {
         double tmp = fXmin;
         fXmin = fXmax;
         fXmax = tmp;
       }
       
       fXminDef = fXmin;
       fXmaxDef = fXmax;
       fPreviousAction = BoundingBox::Undefined;
       fIncreasingCumulativeFunction = (GetCumulativeProbability(fXmax) - GetCumulativeProbability(fXmin))/(fXmax-fXmin);
     }

Here is the call graph for this function:

Member Function Documentation

void G4DNASmoluchowskiDiffusion::BoundingBox::Print ( )

inline

Definition at line 177 of file G4DNASmoluchowskiDiffusion.hh.

     {
       G4cout << "fXmin: " << fXmin << " | fXmax: " << fXmax << G4endl;
     }

bool G4DNASmoluchowskiDiffusion::BoundingBox::Propose	(	double	proposedXValue,
		double	proposedProba,
		double	nextProba,
		double &	returnedValue
	)

inline

Definition at line 182 of file G4DNASmoluchowskiDiffusion.hh.

     {
 //      G4cout << "---------------------------" << G4endl;
 //      G4cout << "Proposed x value: " << proposedXValue
 //          << "| proposedProba: " << proposedProba
 //          << "| nextProba: " << nextProba
 //          << " | fXmin: " << fXmin << " (" << G4DNASmoluchowskiDiffusion::GetCumulativeProbability(fXmin) <<")"
 //          << " | fXmax: " << fXmax << " (" << G4DNASmoluchowskiDiffusion::GetCumulativeProbability(fXmax) <<")"
 //          << G4endl;
 
       bool returnFlag = false;
       
       if(proposedProba < nextProba-fToleranceY) // proba trop petite ==> augmente
       {
         // G4cout << "proposedProba < nextProba-fToleranceY" << G4endl;
 
         if(fIncreasingCumulativeFunction > 0) // croissant
         {
           if(proposedXValue > fXmin)
             fXmin = proposedXValue;
         }
         else if(fIncreasingCumulativeFunction < 0) // decroissant
         {
           if(proposedXValue < fXmax)
             fXmax = proposedXValue;
         }
         
         returnedValue = (fXmax + fXmin)/2;
         returnFlag = false;
         fPreviousAction = BoundingBox::IncreaseProba;
       }
       else if(proposedProba > nextProba+fToleranceY) // proba trop grande
       {
         // G4cout << "proposedProba > nextProba+fToleranceY" << G4endl;
 
         if(fIncreasingCumulativeFunction>0)
         {
           if(proposedXValue < fXmax)
             fXmax = proposedXValue;
         }
         else if(fIncreasingCumulativeFunction<0)
         {
           if(proposedXValue > fXmin)
           {
             fXmin = proposedXValue;
           }
         }
         
         returnedValue = (fXmax + fXmin)/2;
         returnFlag = false;
         fPreviousAction = BoundingBox::DecreaseProba;
       }
       else
       {
         // G4cout << "IN THE INTERVAL !! : " << nextProba << G4endl;
         fSum = proposedProba;
         
         // Assuming search for next proba is increasing
         if(fIncreasingCumulativeFunction<0)
         {
          fXmin = fXminDef;
          fXmax = proposedXValue;
         }
         else if(fIncreasingCumulativeFunction>0)
         {
           fXmin = proposedXValue;
           fXmax = fXmaxDef;
         }
         returnFlag = true;
         fPreviousAction = BoundingBox::Undefined;
       }
       
       return returnFlag;
     }