A class that implements a search for prime numbers. More...

#include <src/CPUCruncher.h>

Inheritance diagram for CPUCruncher:

Collaboration diagram for CPUCruncher:

Public Types
typedef tbb::concurrent_hash_map< std::string, unsigned int >	CHM

Public Types inherited from Algorithm
typedef Gaudi::PluginService::Factory< IAlgorithm , const std::string &, ISvcLocator >	Factory

Public Types inherited from CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > >
using	base_class = CommonMessaging

Public Member Functions
virtual StatusCode	execute ()
	the execution of the algorithm More...

virtual StatusCode	initialize ()
	Its initialization. More...

virtual StatusCode	finalize ()
	the finalization of the algorithm More...

double	get_runtime () const

	CPUCruncher (const std::string &name, ISvcLocator *pSvc)

virtual	~CPUCruncher ()
	virtual & protected desctrustor More...

Public Member Functions inherited from GaudiAlgorithm
StatusCode	initialize () override
	standard initialization method More...

StatusCode	execute () override
	standard execution method More...

StatusCode	finalize () override
	standard finalization method More...

StatusCode	sysExecute () override
	the generic actions for the execution. More...

void	put (IDataProviderSvc svc, DataObject object, const std::string &address, const bool useRootInTES=true) const
	Register a data object or container into Gaudi Event Transient Store. More...

void	put (DataObject *object, const std::string &address, const bool useRootInTES=true) const
	Register a data object or container into Gaudi Event Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	get (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Templated access to the data in Gaudi Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	getIfExists (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Templated access to the data in Gaudi Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	get (const std::string &location, const bool useRootInTES=true) const
	Templated access to the data from Gaudi Event Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	getIfExists (const std::string &location, const bool useRootInTES=true) const
	Templated access to the data in Gaudi Transient Store. More...

template<class TYPE >
TYPE *	getDet (IDataProviderSvc *svc, const std::string &location) const
	Templated access to the detector data from the Gaudi Detector Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	getDetIfExists (IDataProviderSvc *svc, const std::string &location) const
	Templated access to the detector data from the Gaudi Detector Transient Store. More...

template<class TYPE >
TYPE *	getDet (const std::string &location) const
	Templated access to the detector data from the Gaudi Detector Transient Store. More...

template<class TYPE >
Gaudi::Utils::GetData< TYPE >::return_type	getDetIfExists (const std::string &location) const
	Templated access to the detector data from the Gaudi Detector Transient Store. More...

template<class TYPE >
bool	exist (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Check the existence of a data object or container in the Gaudi Transient Event Store. More...

template<class TYPE >
bool	exist (const std::string &location, const bool useRootInTES=true) const
	Check the existence of a data object or container in the Gaudi Transient Event Store. More...

template<class TYPE >
bool	existDet (IDataProviderSvc *svc, const std::string &location) const
	Check the existence of detector objects in the Gaudi Transient Detector Store. More...

template<class TYPE >
bool	existDet (const std::string &location) const
	Check the existence of detector objects in the Gaudi Transient Detector Store. More...

template<class TYPE , class TYPE2 >
Gaudi::Utils::GetData< TYPE >::return_type	getOrCreate (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Get the existing data object from Gaudi Event Transient store. More...

template<class TYPE , class TYPE2 >
Gaudi::Utils::GetData< TYPE >::return_type	getOrCreate (const std::string &location, const bool useRootInTES=true) const
	Get the existing data object from Gaudi Event Transient store. More...

	GaudiAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
	Standard constructor (protected) More...

	~GaudiAlgorithm () override=default
	destructor, virtual and protected More...

SmartIF< INTupleSvc > &	evtColSvc () const
	Access the standard event collection service. More...

	GaudiAlgorithm (const GaudiAlgorithm &)=delete

GaudiAlgorithm &	operator= (const GaudiAlgorithm &)=delete

template<class T >
Property *	declareProperty (const std::string &name, DataObjectHandle< T > &hndl, const std::string &doc="none") const

Public Member Functions inherited from GaudiCommon< Algorithm >
Gaudi::Utils::GetData< TYPE >::return_type	get (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Templated access to the data in Gaudi Transient Store. More...

Gaudi::Utils::GetData< TYPE >::return_type	getIfExists (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Quicker version of the get function which bypasses the check on the retrieved data. More...

bool	exist (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Check the existence of a data object or container in the Gaudi Transient Event Store. More...

Gaudi::Utils::GetData< TYPE >::return_type	getOrCreate (IDataProviderSvc *svc, const std::string &location, const bool useRootInTES=true) const
	Get the existing data object from Gaudi Event Transient store. More...

DataObject *	put (IDataProviderSvc svc, DataObject object, const std::string &location, const bool useRootInTES=true) const
	Register a data object or container into Gaudi Event Transient Store. More...

TOOL *	tool (const std::string &type, const std::string &name, const IInterface *parent=0, bool create=true) const
	Useful method for the easy location of tools. More...

TOOL *	tool (const std::string &type, const IInterface *parent=0, bool create=true) const
	A useful method for the easy location of tools. More...

SmartIF< SERVICE >	svc (const std::string &name, const bool create=true) const
	A useful method for the easy location of services. More...

IUpdateManagerSvc *	updMgrSvc () const
	Short-cut to locate the Update Manager Service. More...

StatusCode	Error (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
	Print the error message and return with the given StatusCode. More...

StatusCode	Warning (const std::string &msg, const StatusCode st=StatusCode::FAILURE, const size_t mx=10) const
	Print the warning message and return with the given StatusCode. More...

StatusCode	Info (const std::string &msg, const StatusCode st=StatusCode::SUCCESS, const size_t mx=10) const
	Print the info message and return with the given StatusCode. More...

StatusCode	Print (const std::string &msg, const StatusCode st=StatusCode::SUCCESS, const MSG::Level lev=MSG::INFO) const
	Print the message and return with the given StatusCode. More...

void	Assert (const bool ok, const std::string &message="", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
	Assertion - throw exception if the given condition is not fulfilled. More...

void	Assert (const bool ok, const char *message, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
	Assertion - throw exception if the given condition is not fulfilled. More...

void	Exception (const std::string &msg, const GaudiException &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
	Create and (re)-throw a given GaudiException. More...

void	Exception (const std::string &msg, const std::exception &exc, const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
	Create and (re)-throw a given exception. More...

void	Exception (const std::string &msg="no message", const StatusCode sc=StatusCode(StatusCode::FAILURE, true)) const
	Create and throw an exception with the given message. More...

const Statistics &	counters () const
	accessor to all counters More...

StatEntity &	counter (const std::string &tag) const
	accessor to certain counter by name More...

bool	typePrint () const
	Insert the actual C++ type of the algorithm/tool in the messages ? More...

bool	propsPrint () const
	Print properties at initialization ? More...

bool	statPrint () const
	Print statistical counters at finalization ? More...

bool	errorsPrint () const
	Print error counters at finalization ? More...

long	printStat (const MSG::Level level=MSG::ALWAYS) const
	perform the actual printout of statistical counters More...

long	printErrors (const MSG::Level level=MSG::ALWAYS) const
	perform the actual printout of error counters More...

long	printProps (const MSG::Level level=MSG::ALWAYS) const
	perform the actual printout of properties More...

void	registerCondition (const std::string &condition, StatusCode(CallerClass::*mf)()=nullptr)
	register the current instance to the UpdateManagerSvc as a consumer for a condition. More...

void	registerCondition (const std::string &condition, CondType &condPtrDest, StatusCode(CallerClass::mf)()=NULL)
	register the current instance to the UpdateManagerSvc as a consumer for a condition. More...

void	registerCondition (char condition, StatusCode(CallerClass::mf)()=NULL)
	just to avoid conflicts with the version using a pointer to a template class. More...

void	registerCondition (TargetClass condition, StatusCode(CallerClass::mf)()=NULL)
	register the current instance to the UpdateManagerSvc as a consumer for a condition. More...

StatusCode	runUpdate ()
	asks the UpdateManagerSvc to perform an update of the instance (if needed) without waiting the next BeginEvent incident. More...

	GaudiCommon (const std::string &name, ISvcLocator *pSvcLocator)
	Algorithm constructor. More...

	GaudiCommon (const std::string &type, const std::string &name, const IInterface *parent)
	Tool constructor. More...

StatusCode	initialize () override
	standard initialization method More...

StatusCode	finalize () override
	standard finalization method More...

StatusCode	release (const IInterface *interface) const
	Manual forced (and 'safe') release of the active tool or service. More...

const Services &	services () const
	get the list of aquired services More...

const std::string &	context () const
	Returns the "context" string. Used to identify different processing states. More...

const std::string &	rootInTES () const
	Returns the "rootInTES" string. More...

const std::string	fullTESLocation (const std::string &location, const bool useRootInTES) const
	Returns the full correct event location given the rootInTes settings. More...

Public Member Functions inherited from Algorithm
	Algorithm (const std::string &name, ISvcLocator *svcloc, const std::string &version=PACKAGE_VERSION)
	Constructor. More...

	~Algorithm () override=default
	Destructor. More...

StatusCode	sysStart () override
	Reinitialization method invoked by the framework. More...

StatusCode	sysInitialize () override
	Initialization method invoked by the framework. More...

StatusCode	sysReinitialize () override
	Reinitialization method invoked by the framework. More...

StatusCode	sysRestart () override
	Restart method invoked by the framework. More...

StatusCode	sysExecute () override
	The actions to be performed by the algorithm on an event. More...

StatusCode	sysStop () override
	System stop. More...

StatusCode	sysFinalize () override
	System finalization. More...

StatusCode	sysBeginRun () override
	beginRun method invoked by the framework. More...

StatusCode	sysEndRun () override
	endRun method invoked by the framework. More...

const std::string &	name () const override
	The identifying name of the algorithm object. More...

const Gaudi::StringKey &	nameKey () const override

const std::string &	type () const override
	The type of the algorithm object. More...

void	setType (const std::string &type) override

const std::string &	version () const override

unsigned int	index () const override

StatusCode	configure () override
	Dummy implementation of IStateful::configure() method. More...

StatusCode	terminate () override
	Dummy implementation of IStateful::terminate() method. More...

StatusCode	initialize () override
	the default (empty) implementation of IStateful::initialize() method More...

StatusCode	start () override
	the default (empty) implementation of IStateful::start() method More...

StatusCode	stop () override
	the default (empty) implementation of IStateful::stop() method More...

StatusCode	finalize () override
	the default (empty) implementation of IStateful::finalize() method More...

StatusCode	reinitialize () override
	the default (empty) implementation of IStateful::reinitialize() method More...

StatusCode	restart () override
	the default (empty) implementation of IStateful::restart() method More...

Gaudi::StateMachine::State	FSMState () const override
	returns the current state of the algorithm More...

Gaudi::StateMachine::State	targetFSMState () const override
	returns the state the algorithm will be in after the ongoing transition More...

bool	isExecuted () const override
	Has this algorithm been executed since the last reset? More...

void	setExecuted (bool state) override
	Set the executed flag to the specified state. More...

void	resetExecuted () override
	Reset the executed state of the Algorithm for the duration of the current event. More...

StatusCode	beginRun () override
	Algorithm begin run. More...

StatusCode	endRun () override
	Algorithm end run. This method is called at the end of the event loop. More...

bool	isEnabled () const override
	Is this algorithm enabled or disabled? More...

bool	filterPassed () const override
	Did this algorithm pass or fail its filter criterion for the last event? More...

void	setFilterPassed (bool state) override
	Set the filter passed flag to the specified state. More...

int	errorCount () const
	Get the number of failures of the algorithm. More...

template<class T >
StatusCode	service (const std::string &name, T *&psvc, bool createIf=true) const
	Access a service by name, creating it if it doesn't already exist. More...

template<class T >
StatusCode	service (const std::string &svcType, const std::string &svcName, T *&psvc) const
	Access a service by name and type, creating it if it doesn't already exist. More...

SmartIF< IService >	service (const std::string &name, const bool createIf=true, const bool quiet=false) const
	Return a pointer to the service identified by name (or "type/name") More...

template<class T >
SmartIF< T >	service (const std::string &name, bool createIf=true, bool quiet=false) const

SmartIF< IAuditorSvc > &	auditorSvc () const
	The standard auditor service.May not be invoked before sysInitialize() has been invoked. More...

SmartIF< IChronoStatSvc > &	chronoSvc () const
	The standard Chrono & Stat service, Return a pointer to the service if present. More...

SmartIF< IChronoStatSvc > &	chronoStatService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IDataProviderSvc > &	detSvc () const
	The standard detector data service. More...

SmartIF< IDataProviderSvc > &	detDataService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IConversionSvc > &	detCnvSvc () const
	The standard detector data persistency conversion service. More...

SmartIF< IConversionSvc > &	detDataCnvService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IDataProviderSvc > &	eventSvc () const
	The standard event data service. More...

SmartIF< IDataProviderSvc > &	evtSvc () const
	shortcut for method eventSvc More...

SmartIF< IDataProviderSvc > &	eventDataService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IConversionSvc > &	eventCnvSvc () const
	The standard event data persistency conversion service. More...

SmartIF< IConversionSvc > &	eventDataCnvService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IHistogramSvc > &	histoSvc () const
	The standard histogram service. More...

SmartIF< IHistogramSvc > &	histogramDataService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< INTupleSvc > &	ntupleSvc () const
	The standard N tuple service. More...

SmartIF< INTupleSvc > &	ntupleService () const
	Obsoleted name, kept due to the backwards compatibility. More...

SmartIF< IRndmGenSvc > &	randSvc () const
	AIDA-based NTuple service Returns a pointer to the AIDATuple service if present. More...

SmartIF< IToolSvc > &	toolSvc () const
	The standard ToolSvc service, Return a pointer to the service if present. More...

SmartIF< IExceptionSvc > &	exceptionSvc () const
	Get the exception Service. More...

SmartIF< IAlgContextSvc > &	contextSvc () const
	get Algorithm Context Service More...

SmartIF< ITimelineSvc > &	timelineSvc () const

SmartIF< ISvcLocator > &	serviceLocator () const override
	The standard service locator. More...

SmartIF< ISvcLocator > &	svcLoc () const
	shortcut for method serviceLocator More...

SmartIF< IHiveWhiteBoard > &	whiteboard () const

bool	registerContext () const
	register for Algorithm Context Service? More...

StatusCode	createSubAlgorithm (const std::string &type, const std::string &name, Algorithm *&pSubAlg)
	Create a sub algorithm. More...

const std::vector< Algorithm * > *	subAlgorithms () const
	List of sub-algorithms. Returns a pointer to a vector of (sub) Algorithms. More...

std::vector< Algorithm * > *	subAlgorithms ()
	List of sub-algorithms. Returns a pointer to a vector of (sub) Algorithms. More...

StatusCode	setProperty (const Property &p) override
	Implementation of IProperty::setProperty. More...

StatusCode	setProperty (const std::string &s) override
	Implementation of IProperty::setProperty. More...

StatusCode	setProperty (const std::string &n, const std::string &v) override
	Implementation of IProperty::setProperty. More...

StatusCode	getProperty (Property *p) const override
	Implementation of IProperty::getProperty. More...

const Property &	getProperty (const std::string &name) const override
	Implementation of IProperty::getProperty. More...

StatusCode	getProperty (const std::string &n, std::string &v) const override
	Implementation of IProperty::getProperty. More...

const std::vector< Property * > &	getProperties () const override
	Implementation of IProperty::getProperties. More...

bool	hasProperty (const std::string &name) const override
	Implementation of IProperty::hasProperty. More...

PropertyMgr *	getPropertyMgr ()

StatusCode	setProperties ()
	Set the algorithm's properties. More...

template<class T >
Property *	declareProperty (const std::string &name, T &property, const std::string &doc="none") const
	Declare the named property. More...

Property *	declareRemoteProperty (const std::string &name, IProperty *rsvc, const std::string &rname="") const
	Declare remote named properties. More...

template<class T >
Property *	declareProperty (const std::string &name, ToolHandle< T > &hndl, const std::string &doc="none") const

template<class T >
Property *	declareProperty (const std::string &name, ToolHandleArray< T > &hndlArr, const std::string &doc="none") const

SmartIF< IMonitorSvc > &	monitorSvc () const
	Access the monitor service. More...

template<class T >
void	declareInfo (const std::string &name, const T &var, const std::string &desc) const
	Declare monitoring information. More...

void	declareInfo (const std::string &name, const std::string &format, const void *var, int size, const std::string &desc) const
	Declare monitoring information (special case) More...

template<class TYPE >
StatusCode	setProperty (const std::string &name, const TYPE &value)
	set the property form the value More...

EventContext *	getContext () const
	get the context More...

void	setContext (EventContext *context)
	set the context More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override

virtual const DataObjIDColl &	extraInputDeps () const override

virtual const DataObjIDColl &	extraOutputDeps () const override

virtual void	acceptDHVisitor (IDataHandleVisitor *) const override

const DataObjIDColl &	inputDataObjs () const

const DataObjIDColl &	outputDataObjs () const

void	commitHandles () override

void	registerTool (IAlgTool *tool) const

void	deregisterTool (IAlgTool *tool) const

template<class T >
StatusCode	declareTool (ToolHandle< T > &handle, std::string toolTypeAndName="", bool createIf=true)

const std::vector< IAlgTool * > &	tools () const

bool	isClonable () const override
	Specifies the clonability of the algorithm. More...

unsigned int	cardinality () const override
	Return the cardinality. More...

const std::vector< std::string > &	neededResources () const override

Public Member Functions inherited from CommonMessagingBase
virtual	~CommonMessagingBase ()=default
	Virtual destructor. More...

SmartIF< IMessageSvc > &	msgSvc () const
	The standard message service. More...

MsgStream &	msgStream () const
	Return an uninitialized MsgStream. More...

MsgStream &	msgStream (const MSG::Level level) const
	Predefined configurable message stream for the efficient printouts. More...

MsgStream &	always () const
	shortcut for the method msgStream(MSG::ALWAYS) More...

MsgStream &	fatal () const
	shortcut for the method msgStream(MSG::FATAL) More...

MsgStream &	err () const
	shortcut for the method msgStream(MSG::ERROR) More...

MsgStream &	error () const
	shortcut for the method msgStream(MSG::ERROR) More...

MsgStream &	warning () const
	shortcut for the method msgStream(MSG::WARNING) More...

MsgStream &	info () const
	shortcut for the method msgStream(MSG::INFO) More...

MsgStream &	debug () const
	shortcut for the method msgStream(MSG::DEBUG) More...

MsgStream &	verbose () const
	shortcut for the method msgStream(MSG::VERBOSE) More...

MsgStream &	msg () const
	shortcut for the method msgStream(MSG::INFO) More...

MSG::Level	msgLevel () const
	get the output level from the embedded MsgStream More...

MSG::Level	outputLevel () const __attribute__((deprecated))
	Backward compatibility function for getting the output level. More...

bool	msgLevel (MSG::Level lvl) const
	get the output level from the embedded MsgStream More...

Public Attributes
int	evtslots = 10

tuple	whiteboard

tuple	eventloopmgr

tuple	a1

tuple	a2 = CPUCruncher("A2")

tuple	a3 = CPUCruncher("A3")

tuple	a4 = CPUCruncher("A4")

string	EvtSel = 'NONE'

list	ExtSvc = [whiteboard]

	EventLoop = eventloopmgr,

list	TopAlg = [a1,a2,a3,a4]

Private Member Functions
	CPUCruncher ()
	the default constructor is disabled More...

	CPUCruncher (const CPUCruncher &)
	the copy constructor is disabled More...

CPUCruncher &	operator= (const CPUCruncher &)
	the assignement operator is disabled More...

void	findPrimes (const unsigned long int)
	The CPU intensive function. More...

void	calibrate ()
	Calibrate. More...

long unsigned int	getNCaliIters (double)

Private Attributes
double	m_avg_runtime

double	m_var_runtime

bool	m_local_rndm_gen

bool	m_shortCalib

const uint	MAX_INPUTS = 40

const uint	MAX_OUTPUTS = 10

std::vector< DataObjectHandle< DataObject > * >	m_inputHandles

std::vector< DataObjectHandle< DataObject > * >	m_outputHandles

std::vector< std::string >	m_inpKeys

std::vector< std::string >	m_outKeys

unsigned int	m_rwRepetitions

bool	m_sleepyExecution

Static Private Attributes
static std::vector< unsigned int >	m_niters_vect

static std::vector< double >	m_times_vect

static CHM	m_name_ncopies_map

Additional Inherited Members
Protected Types inherited from GaudiCommon< Algorithm >
using	base_class = Algorithm

typedef std::map< std::string, StatEntity >	Statistics
	the actual type of general counters More...

typedef std::map< std::string, unsigned int >	Counter
	the actual type error/warning counter More...

typedef std::vector< IAlgTool * >	AlgTools
	storage for active tools More...

typedef std::vector< SmartIF< IService > >	Services
	storage for active services More...

Protected Member Functions inherited from GaudiCommon< Algorithm >
	~GaudiCommon () override=default
	Destructor. More...

StatusCode	releaseTool (const IAlgTool *tool) const
	manual forced (and 'safe') release of the tool More...

StatusCode	releaseSvc (const IInterface *svc) const
	manual forced (and 'safe') release of the service More...

Protected Member Functions inherited from Algorithm
virtual void	declareInput (Gaudi::DataHandle *im) override

virtual void	declareOutput (Gaudi::DataHandle *im) override

std::vector< IAlgTool * > &	tools ()

bool	isInitialized () const override
	Has the Algorithm already been initialized? More...

bool	isFinalized () const override
	Has the Algorithm already been finalized? More...

void	setIndex (const unsigned int &idx) override
	set instantiation index of Alg More...

Protected Member Functions inherited from CommonMessaging< implements< IAlgorithm, IDataHandleHolder, IProperty, IStateful > >
void	updateMsgStreamOutputLevel (int level)
	Update the output level of the cached MsgStream. More...

Protected Attributes inherited from GaudiCommon< Algorithm >
SmartIF< ICounterSummarySvc >	m_counterSummarySvc
	a pointer to the CounterSummarySvc More...

std::vector< std::string >	m_counterList
	list of counters to declare. Set by property CounterList. This can be a regular expression. More...

std::vector< std::string >	m_statEntityList

Protected Attributes inherited from Algorithm
EventContext *	m_event_context
	Event specific data for multiple event processing. More...

SmartIF< PropertyMgr >	m_propertyMgr
	For management of properties. More...

std::unique_ptr< IDataHandleVisitor >	m_updateDataHandles
	Hook for for derived classes to provide a custom visitor for data handles. More...

Static Protected Attributes inherited from GaudiCommon< Algorithm >
static const bool	IgnoreRootInTES
	Simple definition to be used with the new useRootInTES argument get<TYPE> and put methods. More...

static const bool	UseRootInTES
	Simple definition to be used with the new useRootInTES argument get<TYPE> and put methods. More...

Detailed Description

A class that implements a search for prime numbers.

A lot of allocations are performed on purpose not to stress only the FP unit of the CPU.

Definition at line 19 of file CPUCruncher.h.

Member Typedef Documentation

typedef tbb::concurrent_hash_map<std::string,unsigned int> CPUCruncher::CHM

Definition at line 24 of file CPUCruncher.h.

Constructor & Destructor Documentation

CPUCruncher::CPUCruncher	(	const std::string &	name,
		ISvcLocator *	pSvc
	)

Definition at line 19 of file CPUCruncher.cpp.

                           :
         GaudiAlgorithm(name, pSvc), m_avg_runtime(1.), m_var_runtime(.01), m_shortCalib(
                 false) {
 
   declareProperty("inpKeys", m_inpKeys);
   declareProperty("outKeys", m_outKeys);
 
 
     declareProperty("avgRuntime", m_avg_runtime,
             "Average runtime of the module.");
     declareProperty("varRuntime", m_var_runtime,
             "Variance of the runtime of the module.");
     declareProperty("localRndm", m_local_rndm_gen = true,
             "Decide if the local random generator is to be used");
     declareProperty("NIterationsVect", m_niters_vect,
             "Number of iterations for the calibration.");
     declareProperty("NTimesVect", m_times_vect,
             "Number of seconds for the calibration.");
     declareProperty("shortCalib", m_shortCalib = false,
             "Enable coarse grained calibration");
     declareProperty("RwRepetitions", m_rwRepetitions = 1,
             "Increase access to the WB");
     declareProperty("SleepyExecution", m_sleepyExecution = false,
                     "Sleep during execution instead of crunching");
 
     // Register the algo in the static concurrent hash map in order to
     // monitor the # of copies
     CHM::accessor name_ninstances;
     m_name_ncopies_map.insert(name_ninstances, name);
     name_ninstances->second += 1;
 }

CPUCruncher::~CPUCruncher ( )

virtual

virtual & protected desctrustor

Definition at line 52 of file CPUCruncher.cpp.

                           {
   for (uint i = 0; i < m_inputHandles.size(); ++i) {
         delete m_inputHandles[i];
     }
 
     m_outputHandles.resize(MAX_OUTPUTS);
   for (uint i = 0; i < m_outputHandles.size(); ++i) {
         delete m_outputHandles[i];
     }
 }

CPUCruncher::CPUCruncher ( )

private

the default constructor is disabled

CPUCruncher::CPUCruncher ( const CPUCruncher & )

private

the copy constructor is disabled

Member Function Documentation

void CPUCruncher::calibrate ( )

private

Calibrate.

Definition at line 97 of file CPUCruncher.cpp.

                            {
 
   MsgStream log(msgSvc(), name());
   m_niters_vect.push_back(0);
   m_niters_vect.push_back(500);
   m_niters_vect.push_back(600);
   m_niters_vect.push_back(700);
   m_niters_vect.push_back(800);
   m_niters_vect.push_back(1000);
   m_niters_vect.push_back(1300);
   m_niters_vect.push_back(1600);
   m_niters_vect.push_back(2000);
   m_niters_vect.push_back(2300);
   m_niters_vect.push_back(2600);
   m_niters_vect.push_back(3000);
   m_niters_vect.push_back(3300);
   m_niters_vect.push_back(3500);
   m_niters_vect.push_back(3900);
   m_niters_vect.push_back(4200);
   m_niters_vect.push_back(5000);
   m_niters_vect.push_back(6000);
   m_niters_vect.push_back(8000);
   m_niters_vect.push_back(10000);
   m_niters_vect.push_back(12000);
   m_niters_vect.push_back(15000);
   m_niters_vect.push_back(17000);
   m_niters_vect.push_back(20000);
   m_niters_vect.push_back(25000);
   m_niters_vect.push_back(30000);
   m_niters_vect.push_back(35000);
   m_niters_vect.push_back(40000);
   m_niters_vect.push_back(60000);
   if (!m_shortCalib){
     m_niters_vect.push_back(100000);
     m_niters_vect.push_back(200000);
     }
 
 
   m_times_vect.resize(m_niters_vect.size());
   m_times_vect[0]=0.;
 
 
   log << MSG::INFO << "Starting calibration..." << endmsg;
   for (unsigned int i=1;i<m_niters_vect.size();++i){
    unsigned long niters=m_niters_vect[i];
    unsigned int trials = 30;
    do{
     auto start_cali=tbb::tick_count::now();
     findPrimes(niters);
     auto stop_cali=tbb::tick_count::now();
     double deltat = (stop_cali-start_cali).seconds();
     m_times_vect[i]=deltat;
     log << MSG::DEBUG << "Calibration: # iters = " << niters << " => " << deltat << endmsg;
     trials--;
     } while(trials > 0 and m_times_vect[i]<m_times_vect[i-1]); // make sure that they are monotonic
    }
    log << MSG::INFO << "Calibration finished!" << endmsg;
 }

StatusCode CPUCruncher::execute ( )

virtual

the execution of the algorithm

Definition at line 248 of file CPUCruncher.cpp.

 {
 
   MsgStream logstream(msgSvc(), name());
 
   if (m_sleepyExecution) {
     logstream  << MSG::DEBUG << "Going to dream for: "<< m_avg_runtime << endmsg;
     std::chrono::duration<double> dreamtime( m_avg_runtime );
 
     tbb::tick_count starttbb=tbb::tick_count::now();
     std::this_thread::sleep_for(dreamtime);
     tbb::tick_count endtbb=tbb::tick_count::now();
     // actual sleeping time can be longer due to scheduling or resource contention delays
     const double actualDreamTime=(endtbb-starttbb).seconds();
 
     logstream  << MSG::DEBUG << "Actual dreaming time was: "<< actualDreamTime << "s" << endmsg;
 
     return StatusCode::SUCCESS;
   }
 
   float runtime;
 
   if (m_local_rndm_gen){
   /* This will disappear with a thread safe random number generator svc
    * Use box mueller to generate gaussian randoms
    * The quality is not good for in depth study given that the generator is a
    * linear congruent.
    * Throw away basically a free number: we are in a cpu cruncher after all.
    * The seed is taken from the clock, but we could assign a seed per module to
    * ensure reproducibility.
    *
    * This is not an overkill but rather an exercise towards a thread safe
    * random number generation.
    */
 
   auto getGausRandom = [] (double mean, double sigma) -> double {
 
     unsigned int seed = std::clock();
 
     auto getUnifRandom = [] (unsigned int & seed) ->double {
       // from numerical recipies
       constexpr unsigned int m = 232;
       constexpr unsigned int a = 1664525;
       constexpr unsigned int c = 1013904223;
       seed = (a * seed + c) % m;
       const double unif = double(seed) / m;
       return unif;
       };
 
     const double unif1 = getUnifRandom(seed);
     const double unif2 = getUnifRandom(seed);
     const double normal = sqrt(-2.*log(unif1))*cos(2*M_PI*unif2);
     return normal*sigma + mean;
     };
     runtime = fabs(getGausRandom( m_avg_runtime , m_var_runtime ));
     //End Of temp block
   } else {
     // Should be a member.
     HiveRndm::HiveNumbers rndmgaus(randSvc(), Rndm::Gauss( m_avg_runtime , m_var_runtime ));
     runtime = std::fabs(rndmgaus());
   }
 
   tbb::tick_count starttbb=tbb::tick_count::now();
   logstream  << MSG::DEBUG << "Runtime will be: "<< runtime << endmsg;
   if (getContext())
     logstream  << MSG::DEBUG << "Start event " <<  getContext()->evt()
            << " in slot " << getContext()->slot()
            << " on pthreadID " << std::hex << pthread_self() << std::dec
            << endmsg;
 
   for (auto & inputHandle: m_inputHandles){
       if(!inputHandle->isValid())
         continue;
 
     DataObject* obj = nullptr;
     for (unsigned int i=0; i<m_rwRepetitions;++i){
       obj = inputHandle->get();
     }
     if (obj == nullptr)
       logstream << MSG::ERROR << "A read object was a null pointer." << endmsg;
   }
 
   const unsigned long n_iters= getNCaliIters(runtime);
   findPrimes( n_iters );
 
   for (auto & outputHandle: m_outputHandles){
       if(!outputHandle->isValid())
         continue;
 
       outputHandle->put(new DataObject());
   }
 
   for (auto & inputHandle: m_inputHandles){
       if(!inputHandle->isValid())
         continue;
 
     for (unsigned int i=1; i<m_rwRepetitions;++i){
       inputHandle->get();
     }
   }
 
   tbb::tick_count endtbb=tbb::tick_count::now();
 
   const double actualRuntime=(endtbb-starttbb).seconds();
 
   if (getContext())
     logstream << MSG::DEBUG << "Finish event " <<  getContext()->evt()
       //         << " on pthreadID " << getContext()->m_thread_id
          << " in " << actualRuntime  << " seconds" << endmsg;
 
   logstream << MSG::DEBUG << "Timing: ExpectedRuntime= " << runtime
                          << " ActualRuntime= " << actualRuntime
                          << " Ratio= " << runtime/actualRuntime
                          << " Niters= " << n_iters << endmsg;
 
 
   return StatusCode::SUCCESS ;
 }

StatusCode CPUCruncher::finalize ( )

virtual

the finalization of the algorithm

Definition at line 369 of file CPUCruncher.cpp.

 {
   MsgStream log(msgSvc(), name());
 
   unsigned int ninstances;
 
   {
     CHM::const_accessor const_name_ninstances;
     m_name_ncopies_map.find(const_name_ninstances,name());
     ninstances=const_name_ninstances->second;
   }
 
   constexpr double s2ms=1000.;
     // do not show repetitions
     if (ninstances!=0){
       log << MSG::INFO << "Summary: name= "<< name()
                        <<"\t avg_runtime= " << m_avg_runtime*s2ms
                        << "\t n_clones= " << ninstances << endmsg;
 
       CHM::accessor name_ninstances;
       m_name_ncopies_map.find(name_ninstances,name());
       name_ninstances->second=0;
       }
 
   return GaudiAlgorithm::finalize () ;
 }

void CPUCruncher::findPrimes ( const unsigned long int n_iterations )

private

The CPU intensive function.

Definition at line 191 of file CPUCruncher.cpp.

                                                                    {
 
 
   MsgStream log(msgSvc(), name());
 
   // Flag to trigger the allocation
   bool is_prime;
 
   // Let's prepare the material for the allocations
   unsigned int primes_size=1;
   unsigned long* primes = new unsigned long[primes_size];
   primes[0]=2;
 
   unsigned long i = 2;
 
   // Loop on numbers
   for (unsigned long int iiter=0;iiter<n_iterations;iiter++ ){
     // Once at max, it returns to 0
     i+=1;
 
     // Check if it can be divided by the smaller ones
     is_prime = true;
     for (unsigned long j=2;j<i && is_prime;++j){
       if (i%j == 0){
         is_prime = false;
       }
     }// end loop on numbers < than tested one
     if (is_prime){
       // copy the array of primes (INEFFICIENT ON PURPOSE!)
       unsigned int new_primes_size = 1 + primes_size;
       unsigned long* new_primes = new unsigned long[new_primes_size];
 
       for (unsigned int prime_index=0; prime_index<primes_size;prime_index++){
         new_primes[prime_index]=primes[prime_index];
       }
       // attach the last prime
       new_primes[primes_size]=i;
 
       // Update primes array
       delete[] primes;
       primes = new_primes;
       primes_size=new_primes_size;
     } // end is prime
 
   } // end of while loop
 
   // Fool Compiler optimisations:
   for (unsigned int prime_index=0; prime_index<primes_size;prime_index++)
     if (primes[prime_index] == 4)
       log << "This does never happen, but it's necessary too fool aggressive compiler optimisations!"<< endmsg ;
 
   delete[] primes;
 
 }

double CPUCruncher::get_runtime ( ) const

inline

Definition at line 33 of file CPUCruncher.h.

33 { return m_avg_runtime; };

CPUCruncher::m_avg_runtime

double m_avg_runtime

Definition: CPUCruncher.h:57

unsigned long CPUCruncher::getNCaliIters ( double runtime )

private

Definition at line 156 of file CPUCruncher.cpp.

                                                       {
 
   unsigned int smaller_i=0;
   double time=0.;
   bool found=false;
   // We know that the first entry is 0, so we start to iterate from 1
   for (unsigned int i=1;i<m_times_vect.size();i++){
     time = m_times_vect[i];
     if (time>runtime){
       smaller_i=i-1;
       found=true;
       break;
     }
   }
 
   // Case 1: we are outside the interpolation range, we take the last 2 points
   if (not found)
     smaller_i=m_times_vect.size()-2;
 
   // Case 2: we maeke a linear interpolation
   // y=mx+q
   const double x0=m_times_vect[smaller_i];
   const double x1=m_times_vect[smaller_i+1];
   const double y0=m_niters_vect[smaller_i];
   const double y1=m_niters_vect[smaller_i+1];
   const double m=(y1-y0)/(x1-x0);
   const double q=y0-m*x0;
 
   const unsigned long nCaliIters =  m * runtime + q ;
   //always() << x0 << "<" << runtime << "<" << x1 << " Corresponding to " << nCaliIters << " iterations" << endmsg;
 
   return nCaliIters ;
 }

StatusCode CPUCruncher::initialize ( )

virtual

Its initialization.

Definition at line 63 of file CPUCruncher.cpp.

                                   {
   if (m_times_vect.size()==0){
     calibrate();
   }
 
   // This is a bit ugly. There is no way to declare a vector of DataObjectHandles, so
   // we need to wait until initialize when we've read in the input and output key 
   // properties, and know their size, and then turn them
   // into Handles and register them with the framework by calling declareProperty. We
   // could call declareInput/declareOutput on them too.
 
   int i=0;
   for (auto k: m_inpKeys) {
     debug() << "adding input key " << k << endmsg;
     m_inputHandles.push_back( new DataObjectHandle<DataObject>( k, Gaudi::DataHandle::Reader, this ));
     declareProperty("dummy_in_" + std::to_string(i), *(m_inputHandles.back()) );
     i++;
   }
 
   i = 0;
   for (auto k: m_outKeys) {
     debug() << "adding output key " << k << endmsg;
     m_outputHandles.push_back( new DataObjectHandle<DataObject>( k, Gaudi::DataHandle::Writer, this ));
     declareProperty("dummy_out_" + std::to_string(i), *(m_outputHandles.back()) );
     i++;
   }
 
   return StatusCode::SUCCESS ;
 }

CPUCruncher& CPUCruncher::operator= ( const CPUCruncher & )

private

the assignement operator is disabled

Member Data Documentation

tuple CPUCruncher.a1

Initial value:

 = CPUCruncher("A1",
                  varRuntime=.1,
                  avgRuntime=.5,
                  shortCalib = True )

Definition at line 16 of file CPUCruncher.py.

tuple CPUCruncher.a2 = CPUCruncher("A2")

Definition at line 22 of file CPUCruncher.py.

tuple CPUCruncher.a3 = CPUCruncher("A3")

Definition at line 26 of file CPUCruncher.py.

tuple CPUCruncher.a4 = CPUCruncher("A4")

Definition at line 30 of file CPUCruncher.py.

CPUCruncher.EventLoop = eventloopmgr,

Definition at line 37 of file CPUCruncher.py.

tuple CPUCruncher.eventloopmgr

Initial value:

 = HiveEventLoopMgr(MaxEventsParallel = evtslots,
                                 MaxAlgosParallel  = 20,
                                 NumThreads = 8,
                                 AlgosDependencies = [[],['a1'],['a1'],['a2','a3']])

Definition at line 11 of file CPUCruncher.py.

string CPUCruncher.EvtSel = 'NONE'

Definition at line 35 of file CPUCruncher.py.

int CPUCruncher.evtslots = 10

Definition at line 6 of file CPUCruncher.py.

list CPUCruncher.ExtSvc = [whiteboard]

Definition at line 36 of file CPUCruncher.py.

double CPUCruncher::m_avg_runtime

private

Definition at line 57 of file CPUCruncher.h.

std::vector<std::string> CPUCruncher::m_inpKeys

private

Definition at line 74 of file CPUCruncher.h.

std::vector<DataObjectHandle<DataObject> *> CPUCruncher::m_inputHandles

private

Definition at line 71 of file CPUCruncher.h.

bool CPUCruncher::m_local_rndm_gen

private

Definition at line 59 of file CPUCruncher.h.

CPUCruncher::CHM CPUCruncher::m_name_ncopies_map

staticprivate

Definition at line 78 of file CPUCruncher.h.

std::vector< unsigned int > CPUCruncher::m_niters_vect

staticprivate

Definition at line 63 of file CPUCruncher.h.

std::vector<std::string> CPUCruncher::m_outKeys

private

Definition at line 74 of file CPUCruncher.h.

std::vector<DataObjectHandle<DataObject> *> CPUCruncher::m_outputHandles

private

Definition at line 72 of file CPUCruncher.h.

unsigned int CPUCruncher::m_rwRepetitions

private

Definition at line 76 of file CPUCruncher.h.

bool CPUCruncher::m_shortCalib

private

Definition at line 60 of file CPUCruncher.h.

bool CPUCruncher::m_sleepyExecution

private

Definition at line 81 of file CPUCruncher.h.

std::vector< double > CPUCruncher::m_times_vect

staticprivate

Definition at line 64 of file CPUCruncher.h.

double CPUCruncher::m_var_runtime

private

Definition at line 58 of file CPUCruncher.h.

const uint CPUCruncher::MAX_INPUTS = 40

private

Definition at line 68 of file CPUCruncher.h.

const uint CPUCruncher::MAX_OUTPUTS = 10

private

Definition at line 69 of file CPUCruncher.h.

list CPUCruncher.TopAlg = [a1,a2,a3,a4]

Definition at line 38 of file CPUCruncher.py.

tuple CPUCruncher.whiteboard

Initial value:

1 = HiveWhiteBoard("EventDataSvc",

2 EventSlots = evtslots)

Definition at line 8 of file CPUCruncher.py.

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

GaudiHive/src/CPUCruncher.h
GaudiHive/options/CPUCruncher.py
GaudiHive/src/CPUCruncher.cpp

Public Types

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation