Add diffusion_matrix_sparse() and its tests.

[octave.git] / tests / diffusion_matrix_sparse_tests.m

diff --git a/tests/diffusion_matrix_sparse_tests.m b/tests/diffusion_matrix_sparse_tests.m
new file mode 100644 (file)
index 0000000..735c79b
--- /dev/null
+++ b/tests/diffusion_matrix_sparse_tests.m
@@ -0,0 +1,11 @@
+expected_K = [2, -1, 0, 0, -1;
+             -1, 2, -1, 0, 0;
+             0, -1, 2, -1, 0;
+             0, 0, -1, 2, -1;
+             -1, 0, 0, -1, 2];
+
+actual_K = diffusion_matrix_sparse(5);
+
+unit_test_equals("diffusion_matrix_sparse(5) looks right", ...
+                true, ...
+                isequal(actual_K, expected_K));