Bugfix: Initialise sumMagToDiag result to zero. Henrik Rusche

2016-06-10 10:25:12 +01:00 · 2016-06-10 10:25:12 +01:00 · 378b24ba22
commit 378b24ba22
parent f87887df69
1 changed files with 2 additions and 0 deletions
--- a/src/foam/primitives/VectorN/expandContract/ExpandTensorNI.H
+++ b/src/foam/primitives/VectorN/expandContract/ExpandTensorNI.H
@ -358,6 +358,8 @@ inline void sumMagToDiag
    const TensorN<Cmpt, length>& t
 )
 {
+    result = VectorN<Cmpt, length>::zero;
+
    for (register int i = 0; i < TensorN<Cmpt, length>::rowLength; i++)
    {
        for (register int j = 0; j < TensorN<Cmpt, length>::rowLength; j++)