From a523040ed8c9b39256174b727100ac9296bfaa4a Mon Sep 17 00:00:00 2001
From: Matthias Frey <matthias.frey@psi.ch>
Date: Wed, 15 Jul 2020 14:51:58 +0200
Subject: [PATCH] SigmaGenerator: move initialization
---
src/Distribution/SigmaGenerator.cpp | 4 +---
1 file changed, 1 insertion(+), 3 deletions(-)
diff --git a/src/Distribution/SigmaGenerator.cpp b/src/Distribution/SigmaGenerator.cpp
index 63e861cbf..57ee37f2f 100644
--- a/src/Distribution/SigmaGenerator.cpp
+++ b/src/Distribution/SigmaGenerator.cpp
@@ -772,8 +772,6 @@ SigmaGenerator::updateInitialSigma(const matrix_t& M,
void SigmaGenerator::updateSigma(const std::vector<matrix_t>& Mscs,
const std::vector<matrix_t>& Mcycs)
{
- matrix_t M = boost::numeric::ublas::matrix<double>(6,6);
-
std::ofstream writeSigma;
if (write_m) {
@@ -795,7 +793,7 @@ void SigmaGenerator::updateSigma(const std::vector<matrix_t>& Mscs,
for (unsigned int i = 1; i < nSteps_m; ++i) {
// transfer matrix for one angle
- M = boost::numeric::ublas::prod(Mscs[i - 1],Mcycs[i - 1]);
+ matrix_t M = boost::numeric::ublas::prod(Mscs[i - 1],Mcycs[i - 1]);
// transfer the matrix sigma
sigmas_m[i] = matt_boost::gemmm<matrix_t>(M,sigmas_m[i - 1],
boost::numeric::ublas::trans(M));
--
GitLab