本文整理汇总了C++中VectorType::squared_magnitude方法的典型用法代码示例。如果您正苦于以下问题:C++ VectorType::squared_magnitude方法的具体用法?C++ VectorType::squared_magnitude怎么用?C++ VectorType::squared_magnitude使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类VectorType
的用法示例。
在下文中一共展示了VectorType::squared_magnitude方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: ComputeParallelCurve
AbstractModel::SufficientStatisticsVector
MultivariateModel
::GetSufficientStatistics(const Realizations &R, const Observations& Obs)
{
/// Computation of the suffisient statistics of the model
/// Basically, it is a vector (named S1 to S9) of vectors
/// S1 <- y_ij * eta_ij & S2 <- eta_ij * eta_ij
VectorType S1(m_NbTotalOfObservations), S2(m_NbTotalOfObservations);
auto itS1 = S1.begin(), itS2 = S2.begin();
for(size_t i = 0; i < m_NumberOfSubjects; ++i)
{
for(size_t j = 0; j < Obs.GetNumberOfTimePoints(i); ++j)
{
VectorType PC = ComputeParallelCurve(i, j);
*itS1 = dot_product(PC, Obs.GetSubjectCognitiveScore(i, j));
*itS2 = PC.squared_magnitude();
++itS1, ++itS2;
}
}
/// S3 <- Ksi_i & S4 <- Ksi_i * Ksi_i
VectorType S3 = R.at("Ksi");
VectorType S4 = R.at("Ksi") % R.at("Ksi");
/// S5 <- Tau_i & S6 <- Tau_i * Tau_i
VectorType S5 = R.at("Tau");
VectorType S6 = R.at("Tau") % R.at("Tau");
/// S7 <- G
VectorType S7(1, R.at("G", 0));
/// S8 <- beta_k
VectorType S8((m_ManifoldDimension-1) * m_NbIndependentSources);
ScalarType * itS8 = S8.memptr();
for(size_t i = 0; i < S8.size(); ++i)
itS8[i] = R.at("Beta#" + std::to_string(i), 0);
/// S8 <- delta_k
VectorType S9(m_ManifoldDimension - 1);
ScalarType * itS9 = S9.memptr();
for(size_t i = 1; i < S9.size(); ++i)
itS9[i] = R.at("Delta#" + std::to_string(i), 0);
return {S1, S2, S3, S4, S5, S6, S7, S8, S9};
}