From bf219fb10e37a52fc30b780ae60909ebf3b0e58b Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Mon, 23 Jul 2012 10:28:47 +0000 Subject: [PATCH] test for finite residual git-svn-id: https://svn.dealii.org/trunk@25720 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/lac/householder.h | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/deal.II/include/deal.II/lac/householder.h b/deal.II/include/deal.II/lac/householder.h index 6dd6b64fda..948f1b443e 100644 --- a/deal.II/include/deal.II/lac/householder.h +++ b/deal.II/include/deal.II/lac/householder.h @@ -163,7 +163,7 @@ Householder::initialize(const FullMatrix& M) number2 s = (this->el(j,j) < 0) ? std::sqrt(sigma) : -std::sqrt(sigma); // number2 beta = std::sqrt(1./(sigma-s*this->el(j,j))); - + // Make column j the Householder // vector, store first entry in // diagonal @@ -233,6 +233,8 @@ Householder::least_squares (Vector& dst, number2 sum = 0.; for (unsigned int i=n ; ibackward(dst, *aux); @@ -276,6 +278,8 @@ Householder::least_squares (BlockVector& dst, number2 sum = 0.; for (unsigned int i=n ; i