Merge pull request #36 from memmett/feature/fix-interp

Fix interpolation matrix.

Author: memmett

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;
+	}
       }
     }