solnihv:: set result[3] to zero, as otherwise it may become infinite

2019-10-26 18:24:45 +02:00 · 2019-10-26 18:24:45 +02:00 · a9d3efe02c
parent fde37d56a9
commit a9d3efe02c
1 changed files with 2 additions and 0 deletions
--- a/nonisotropic.cpp
+++ b/nonisotropic.cpp
@ -110,6 +110,8 @@ EX namespace solnihv {
  
      for(int t=0; t<3; t++)
        res[t] = S2(ax) * (bx-ix) + S2(bx) * (ix-ax);
+      
+      res[3] = 0;
      }
    
    return res;