tests/advection_matrix_sparse_tests.m

   1 result = true;
   2
   3 for n = [3:100]
   4   ## Check the sparse matrix against the full one.
   5   expected = advection_matrix(n,0,1);
   6   actual = advection_matrix_sparse(n);
   7
   8   if (expected != actual)
   9      result = false;
  10   end
  11 end
  12
  13
  14 unit_test_equals("Sparse advection matrices agree with full ones", ...
  15                  true, ...
  16                  result);