function wathen_matrix_test ( )

%*****************************************************************************80
%
%% wathen_matrix_test() tests wathen_matrix().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 February 2019
%
%  Author:
%
%    John Burkardt
%
  addpath ( '../wathen_matrix' );

  timestamp ( );
  fprintf ( 1, '\n' );
  fprintf ( 1, 'wathen_matrix_test():\n' );
  fprintf ( 1, '  MATLAB/Octave version %s.\n', version ( ) );
  fprintf ( 1, '  Test wathen_matrix().\n' );
%
%  Set up and solve a single system with FULL, BANDED, SPARSE and DAVIS.
%
  wathen_test01 ( );
  wathen_test02 ( );
  wathen_test03 ( );
  wathen_test04 ( );
%
%  Demonstrate storage for a sequence of problems.
%
  wathen_test05 ( );
%
%  For various problem sizes:
%  * time FULL
%  * time BANDED
%  * time FULL, BANDED, SPARSE, DAVIS.
%  * time SPARSE, DAVIS (larger systems).
%
  wathen_test06 ( );
  wathen_test07 ( );
  wathen_test08 ( );
  wathen_test09 ( );
%
%  Use the conjugate gradient solver
%  * with FULL storage
%  * with ST storage
%  * with BANDED storage.
%  * with sparse storage.
%
  wathen_test10 ( );
  wathen_test11 ( );
  wathen_test115 ( );
  wathen_test12 ( );
%
%  Use spy() to display the sparsity of the matrix.
%
  wathen_spy_test ( );
%
%  Get XY coordinates of nodes and plot the grid.
%
  wathen_xy_test ( );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'wathen_matrix_test():\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  fprintf ( 1, '\n' );
  timestamp ( );

  rmpath ( '../wathen_matrix' );

  return
end
function timestamp ( )

%*****************************************************************************80
%
%% timestamp() prints the current YMDHMS date as a timestamp.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 February 2003
%
%  Author:
%
%    John Burkardt
%
  t = now;
  c = datevec ( t );
  s = datestr ( c, 0 );
  fprintf ( 1, '%s\n', s );

  return
end

