#! /usr/bin/env python
#
def r8mat_norm_li ( m, n, a ):

#*****************************************************************************80
#
## R8MAT_NORM_LI returns the matrix Loo norm of an R8MAT.
#
#  Discussion:
#
#    The matrix L-infinity norm is defined as:
#
#      value =  max ( 1 <= I <= M ) sum ( 1 <= J <= N ) abs ( A(I,J) ).
#
#    The matrix L-infinity norm is derived from the vector L-infinity norm,
#    and satisifies:
#
#      vec_norm_li ( A * x ) <= mat_norm_li ( A ) * vec_norm_li ( x ).
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    16 February 2015
#
#  Author:
#
#    John Burkardt
#
#  Parameters:
#
#    Input, integer M, the number of rows in A.
#
#    Input, integer N, the number of columns in A.
#
#    Input, real A(M,N), the matrix whose norm is desired.
#
#    Output, real VALUE, the norm of A.
#
  value = 0.0

  for i in range ( 0, m ):
    col_sum = 0.0
    for j in range ( 0, n ):
      col_sum = col_sum + abs ( a[i,j] )
    value = max ( value, col_sum )

  return value

def r8mat_norm_li_test ( ):

#*****************************************************************************80
#
## R8MAT_NORM_LI_TEST tests R8MAT_NORM_LI.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    16 February 2015
#
#  Author:
#
#    John Burkardt
#
  import platform
  from r8mat_print import r8mat_print
  from r8mat_uniform_ab import r8mat_uniform_ab

  m = 5
  n = 4
  x1 = -5.0
  x2 = +5.0
  seed = 123456789

  a, seed = r8mat_uniform_ab ( m, n, x1, x2, seed )
  
  for j in range ( 0, n ):
    for i in range ( 0, m ):
      a[i,j] = round ( a[i,j] )

  print ( '' )
  print ( 'R8MAT_NORM_LI_TEST' )
  print ( '  Python version: %s' % ( platform.python_version ( ) ) )
  print ( '  R8MAT_NORM_LI computes the Loo norm of an R8MAT;' )

  t = r8mat_norm_li ( m, n, a )

  r8mat_print ( m, n, a, '  A:' )

  print ( '' )
  print ( '  Computed Loo norm = %g' % ( t ) )
#
#  Terminate.
#
  print ( '' )
  print ( 'R8MAT_NORM_LI_TEST' )
  print ( '  Normal end of execution.' )
  return

if ( __name__ == '__main__' ):
  from timestamp import timestamp
  timestamp ( )
  r8mat_norm_li_test ( )
  timestamp ( )