From 777fd25870365dc984b416dc0b9f5cf544d6a3bb Mon Sep 17 00:00:00 2001
From: Matthew Emmett <memmett@gmail.com>
Date: Tue, 8 Jul 2014 16:32:10 -0700
Subject: [PATCH] quadrature: Put zeros in interp matrix where appropriate.

This fixes the MLSDC examples.  Guess how long that took to track down?

Signed-off-by: Matthew Emmett <memmett@gmail.com>
---
 include/pfasst/quadrature.hpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/include/pfasst/quadrature.hpp b/include/pfasst/quadrature.hpp
index ff2b7d63..301ca72b 100644
--- a/include/pfasst/quadrature.hpp
+++ b/include/pfasst/quadrature.hpp
@@ -281,7 +281,9 @@ namespace pfasst
 
         if (abs(num) > 1e-32) {
           mat(i, j) = num / den;
-        }
+        } else {
+	  mat(i, j) = 0.0;
+	}
       }
     }