--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/+util/semidiscrete2function.m	Thu Sep 17 10:12:50 2015 +0200
@@ -0,0 +1,13 @@
+function F = semidiscrete2function(D,S)
+    % TODO: Find threshhold when sparse makes sense.
+    % TODO: Why does having the spase block slow down convergence.m (or does it?)
+    if size(D,1) > 500
+        D = sparse(D);
+    end
+
+    if ~isa(S,'function_handle')
+        F = @(w,t)(D*w + S);
+    else
+        F = @(w,t)(D*w + S(t));
+    end
+end