#! /usr/bin/env python
#
def c8mat_identity ( n ):

#*****************************************************************************80
#
## C8MAT_IDENTITY returns the identity matrix.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license. 
#
#  Modified:
#
#    14 February 2015
#
#  Author:
#
#    John Burkardt
#
#  Parameters:
#
#    Input, integer N, the order of the matrix.
#
#    Output, complex C(N,N), the identity matrix.
#
  import numpy as np

  c = np.zeros ( ( n, n ), 'complex' )

  for i in range ( 0, n ): 
    c[i][i] = complex ( 1.0, 0.0 )

  return c

def c8mat_identity_test ( ):

#*****************************************************************************80
#
## C8MAT_IDENTITY_TEST tests C8MAT_IDENTITY.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license. 
#
#  Modified:
#
#    14 February 2015
#
#  Author:
#
#    John Burkardt
#
  import numpy as np
  import platform
  from c8mat_print import c8mat_print

  m = 4
  n = 4

  print ( '' )
  print ( 'C8MAT_IDENTITY_TEST' )
  print ( '  Python version: %s' % ( platform.python_version ( ) ) )
  print ( '  C8MAT_IDENTITY returns the complex identity matrix.' )

  c = c8mat_identity ( n )

  c8mat_print ( m, n, c, '  Identity matrix:' )
#
#  Terminate.
#
  print ( '' )
  print ( 'C8MAT_IDENTITY_TEST:' )
  print ( '  Normal end of execution.' )
  return

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