void GenericMatrixCovarianceFunction<P_V,P_M,Q_V,Q_M>::covMatrix(const P_V& positionVector1, const P_V& positionVector2, Q_M& imageMatrix) const { if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) { *m_env.subDisplayFile() << "Entering GenericMatrixCovarianceFunction<P_V,P_M,Q_V,Q_M>::covMatrix()" << std::endl; } unsigned int matrixOrder = m_imageSet.vectorSpace().dimLocal(); queso_require_equal_to_msg(imageMatrix.numRowsLocal(), matrixOrder, "imageMatrix has invalid number of rows"); queso_require_equal_to_msg(imageMatrix.numCols(), matrixOrder, "imageMatrix has invalid number of columns"); queso_require_msg(m_covRoutinePtr, "m_covRoutinePtr = NULL"); m_covRoutinePtr(positionVector1, positionVector2, m_routineDataPtr, imageMatrix); if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) { *m_env.subDisplayFile() << "Leaving GenericMatrixCovarianceFunction<P_V,P_M,Q_V,Q_M>::covMatrix()" << std::endl; } return; }
double GenericScalarCovarianceFunction<V,M>::value(const V& positionVector1, const V& positionVector2) const { if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) { *m_env.subDisplayFile() << "Entering GenericScalarCovarianceFunction<V,M>::value()" << std::endl; } UQ_FATAL_TEST_MACRO(m_covRoutinePtr == NULL, m_env.worldRank(), "GenericScalarCovarianceFunction<V,M>::value()", "m_covRoutinePtr = NULL"); double result = m_covRoutinePtr(positionVector1, positionVector2, m_routineDataPtr); if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) { *m_env.subDisplayFile() << "Leaving GenericScalarCovarianceFunction<V,M>::value()" << std::endl; } return result; }