Merge branch 'sparsecholesky_fix' into 'master'

Another fix for size 0 matrices See merge request jschoeberl/ngsolve!340
NGSolve · Mar 27, 2018 · 878f236d23286543bc7dfd2361adf5e0e25af928 · 878f236
2 parents e7c5de7 + d6828a5
commit 878f236d23286543bc7dfd2361adf5e0e25af928
Unified Split

Showing with 6 additions and 1 deletion.

+6 −1 linalg/sparsecholesky.cpp
diff --git a/linalg/sparsecholesky.cpp b/linalg/sparsecholesky.cpp
@@ -396,7 +396,12 @@ namespace ngla
         first_microtask.Append (microtasks.Size());
 
         // int nb = extdofs.Size() / 256 + 1;
-        int nb = (extdofs.Size()+255) / 256;
+        // int nb = (extdofs.Size()+255) / 256;
+	int nb = 0;
+	if(BlockDofs(i).Size()) {
+	  auto extdofs = BlockExtDofs (i);
+	  nb = (extdofs.Size()+255) / 256;
+	}
 
         if (nb == 1)
           // if (false)