Commit 568983fc authored by Fuchsi*'s avatar Fuchsi*

Merge branch 'development' of hemioGit:xerus/xerus into development

parents 97294718 96ceef03
Pipeline #1225 failed with stages
in 11 minutes and 17 seconds
......@@ -43,6 +43,7 @@ namespace xerus { namespace uq { namespace impl_uqRaAdf {
const size_t N;
const size_t d;
double initialResiduum = 0.0;
const double targetResidual;
const size_t maxRank = 50;
const double minRankEps = 1e-10;
......@@ -456,7 +457,7 @@ namespace xerus { namespace uq { namespace impl_uqRaAdf {
}
}
LOG(uqADF, "Residual decrease from " << std::scientific << 0.0 /* TODO */ << " to " << std::scientific << residuals.back() << " in " << _iteration << " iterations.");
LOG(uqADF, "Residual decrease from " << std::scientific << initialResiduum << " to " << std::scientific << residuals.back() << " in " << _iteration << " iterations.");
}
......@@ -469,6 +470,11 @@ namespace xerus { namespace uq { namespace impl_uqRaAdf {
calc_right_stack(corePosition);
}
double optResidual, testResidual;
std::vector<double> setResiduals;
std::tie(optResidual, testResidual, setResiduals) = calc_residuals(0);
initialResiduum = optResidual; /* TODO: inefficient */
for(size_t iteration = 0; maxIterations == 0 || iteration < maxIterations; ++iteration) {
double optResidual, testResidual;
std::vector<double> setResiduals;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment