]> https://gitweb.dealii.org/ - dealii.git/commitdiff
The matrix_norm_square function forgot to conjugate its argument. Adjust the document...
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 9 Dec 2007 22:36:37 +0000 (22:36 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 9 Dec 2007 22:36:37 +0000 (22:36 +0000)
git-svn-id: https://svn.dealii.org/trunk@15587 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/full_matrix.h
deal.II/lac/include/lac/full_matrix.templates.h

index 9903ea99284fd67ef3579669b1f5b09587b30cac..976398c612f0b318bfdf2f884dcdc4d358a4cede 100644 (file)
@@ -440,8 +440,7 @@ class FullMatrix : public Table<2,number>
 
                                     /**
                                      * Return the square of the norm
-                                     * of the vector <tt>v</tt> with
-                                     * respect to the norm induced by
+                                     * of the vector <tt>v</tt> induced by
                                      * this matrix,
                                      * i.e. <i>(v,Mv)</i>. This is
                                      * useful, e.g. in the finite
index 32629554145ebd4686f270e505bec9bd5945e8cb..e9a48242ea3755686c0be8a2254a8e89a3a7c526 100644 (file)
@@ -620,7 +620,7 @@ FullMatrix<number>::matrix_norm_square (const Vector<number2> &v) const
       while (val_ptr != val_end_of_row)
        s += number(*val_ptr++) * number(*v_ptr++);
 
-      sum += s* number(v(row));
+      sum += s * number(numbers::NumberTraits<number2>::conjugate(v(row)));
     }
 
   return sum;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.