function matrix_vector_multiplication_test ( )

%*****************************************************************************80
%
%% matrix_vector_multiplication_test() computes a matrix vector product using einsum().
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    09 August 2026
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '  matrix_vector_multiplication_test():\n' );

  A = [ ...
    11, 12, 13; ...
    21, 22, 23 ];
  fprintf ( 1, '\n' );
  fprintf ( 1, '  A:\n' );
  disp ( A );

  v = [ 1, 2 ];
  fprintf ( 1, '  v:\n' );
  disp ( v );

  Av = einsum ( 'ij,i -> j', A, v );
  fprintf ( 1, '  Av = einsum ( "ij -> j", A, v ):\n');
  disp ( Av );

  return
end

