#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 {0}

double	fIncreasingCumulativeFunction

PreviousAction	fPreviousAction

Detailed Description

Definition at line 138 of file G4DNASmoluchowskiDiffusion.hh.

Member Enumeration Documentation

◆ PreviousAction

enum G4DNASmoluchowskiDiffusion::BoundingBox::PreviousAction

Enumerator
IncreaseProba
DecreaseProba
Undefined

Definition at line 148 of file G4DNASmoluchowskiDiffusion.hh.

    {
      IncreaseProba,
      DecreaseProba,
      Undefined
    };

Constructor & Destructor Documentation

◆ BoundingBox()

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

inline

Definition at line 157 of file G4DNASmoluchowskiDiffusion.hh.

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

Member Function Documentation

◆ Print()

void G4DNASmoluchowskiDiffusion::BoundingBox::Print ( )

inline

Definition at line 176 of file G4DNASmoluchowskiDiffusion.hh.

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

◆ Propose()

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

inline

Definition at line 181 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;
    }