Add advection matrix_sparse() and a test for it.

[octave.git] / tests / advection_matrix_sparse_tests.m

diff --git a/tests/advection_matrix_sparse_tests.m b/tests/advection_matrix_sparse_tests.m
new file mode 100644 (file)
index 0000000..1880dee
--- /dev/null
+++ b/tests/advection_matrix_sparse_tests.m
@@ -0,0 +1,16 @@
+result = true;
+
+for n = [3:100]
+  ## Check the sparse matrix against the full one.
+  expected = advection_matrix(n,0,1);
+  actual = advection_matrix_sparse(n);
+
+  if (expected != actual)
+     result = false;
+  end
+end
+
+
+unit_test_equals("Sparse advection matrices agree with full ones", ...
+                true, ...
+                result);