matrix_vector_multiplication_test <- function ( )

#*****************************************************************************80
#
## matrix_vector_multiplication_test() computes a matrix vector product using einsum().
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    07 August 2026
#
#  Author:
#
#    John Burkardt
#
cat ( '\n' )
cat ( '  matrix_vector_multiplication_test():\n' )

A = matrix ( c (
    11, 12, 13,
    21, 22, 23 ), nrow = 2, ncol = 3, byrow = TRUE )
cat ( '  A:\n' )
print ( A )

v = c ( 1, 2 )
cat ( '  v:\n' )
print ( v )

Av = einsum ( 'ij,i -> j', A, v )
cat ( '  Av = einsum ( \'ij -> j\', A, v ):\n')
print ( Av )