ExperimentStorage<S_V,S_M,D_V,D_M>::ExperimentStorage(
const VectorSpace<S_V,S_M>& scenarioSpace,
unsigned int numExperiments)
:
m_env (scenarioSpace.env()),
m_scenarioSpace (scenarioSpace),
m_paper_n (numExperiments),
m_paper_n_ys_transformed (m_paper_n,0),
m_paper_n_y (0),
m_addId (0),
m_scenarioVecs_standard (m_paper_n, (S_V*) NULL),
m_dataVecs_transformed (m_paper_n, (D_V*) NULL),
m_covMats_transformed_inv(m_paper_n, (D_M*) NULL),
m_y_space (NULL),
m_yVec_transformed (NULL),
m_Wy (NULL)
{
if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
*m_env.subDisplayFile() << "Entering ExperimentStorage<S_V,S_M,D_V,D_M>::constructor()"
<< "\n m_paper_n = " << m_paper_n
<< std::endl;
}
if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 2)) {
*m_env.subDisplayFile() << "Leaving ExperimentStorage<S_V,S_M,D_V,D_M>::constructor()"
<< std::endl;
}
}
VectorSpace<V,M>::VectorSpace(const VectorSpace<V,M>& aux)
: VectorSet<V,M>(aux.env(),((std::string)(aux.m_prefix)).c_str(),INFINITY),
m_dimGlobal(aux.m_dimGlobal),
m_map(newMap()),
m_dimLocal(m_map->NumMyElements()),
m_componentsNamesArray(NULL),
m_componentsNamesVec(NULL),
m_emptyComponentName(""),
m_zeroVector(new V(m_env,*m_map))
{
if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 5)) {
*m_env.subDisplayFile() << "Entering VectorSpace<V,M>::constructor(2)"
<< ": aux.m_componentsNamesArray = " << aux.m_componentsNamesArray
<< ", aux.m_componentsNamesVec = " << aux.m_componentsNamesVec
<< std::endl;
}
if (aux.m_componentsNamesArray != NULL) {
m_componentsNamesArray = new DistArray<std::string>(*(aux.m_componentsNamesArray));
}
if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 5)) {
*m_env.subDisplayFile() << "Leaving VectorSpace<V,M>::constructor(2)"
<< std::endl;
}
}
ArrayOfOneDTables<V,M>::ArrayOfOneDTables(
const char* prefix,
const VectorSpace<V,M>& rowSpace)
:
m_env (rowSpace.env()),
m_prefix ((std::string)(prefix)+""),
m_rowSpace (rowSpace ),
m_oneDTables(m_rowSpace.map(),1)
{
for (unsigned int i = 0; i < (unsigned int) m_oneDTables.MyLength(); ++i) {
m_oneDTables(i,0) = NULL;
}
}
ArrayOfOneDGrids<V,M>::ArrayOfOneDGrids(
const char* prefix,
const VectorSpace<V,M>& rowSpace)
:
m_env (rowSpace.env() ),
m_prefix ((std::string)(prefix)+""),
m_rowSpace (rowSpace ),
m_oneDGrids (m_rowSpace.map(),1),
m_sizes (NULL),
m_minPositions(NULL),
m_maxPositions(NULL)
{
for (unsigned int i = 0; i < (unsigned int) m_oneDGrids.MyLength(); ++i) {
m_oneDGrids(i,0) = NULL;
}
}
BaseTKGroup<V,M>::BaseTKGroup(
const char* prefix,
const VectorSpace<V,M>& vectorSpace,
const std::vector<double>& scales)
:
m_emptyEnv (NULL),
m_env (vectorSpace.env()),
m_prefix (prefix),
m_vectorSpace (&vectorSpace),
m_scales (scales.size(),1.),
m_preComputingPositions(scales.size()+1,NULL), // Yes, +1
m_rvs (scales.size(),NULL) // IMPORTANT: it stays like this for scaledTK, but it will be overwritten to '+1' by hessianTK constructor
{
for (unsigned int i = 0; i < m_scales.size(); ++i) {
m_scales[i] = scales[i];
}
}