matrix_matrix_multiplication_hadamard_test <- function ( )

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

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

B = matrix ( c (
    1, 2, 3,
    4, 5, 6 ), nrow = 2, ncol = 3, byrow = TRUE )
cat ( '\n' )
cat ( '  B:\n' )
print ( B )

AB = einsum ( 'ij,ij -> ij', A, B )
cat ( '\n' )
cat ( '  AB = einsum ( \'ij,ij -> ij\', A, B ):\n' )
print ( AB )