Negate the poisson matrix during its generation.
authorMichael Orlitzky <michael@orlitzky.com>
Tue, 16 Oct 2012 02:16:13 +0000 (22:16 -0400)
committerMichael Orlitzky <michael@orlitzky.com>
Tue, 16 Oct 2012 02:16:13 +0000 (22:16 -0400)
poisson_matrix.m

index 696185ff0e6418649aecb60d47e05d274e9b33f3..4171374f0cb03799d3312bc29758218101717173 100644 (file)
@@ -60,6 +60,7 @@ function A = poisson_matrix(integerN, x0, xN)
     A = cat(1, A, u_row);
   end
 
-  ## Finally, append the last row for xN.
-  A = cat(1, A, u_xN_coeffs);
+  ## Finally, append the last row for xN and negate the whole thing (see
+  ## the definition of the poisson problem).
+  A = - cat(1, A, u_xN_coeffs);
 end