#! /usr/bin/python3
#
def ppm_example ( width, height ):

#*****************************************************************************80
#
## ppm_example() sets up sample RGB data suitable for a PPM file.
#
#  Licensing:
#
#    This code is distributed under the MIT license.
#
#  Modified:
#
#    14 September 2018
#
#  Author:
#
#    John Burkardt
#
#  Input:
#
#    integer WIDTH, HEIGHT, the number of columns and rows of data.
#    Reasonable values would be WIDTH = 200, HEIGHT = 600.
#
#  Output:
#
#    integer RGB[WIDTH*HEIGHT,3], the RGB data.
#
  import numpy as np

  rgb = np.zeros ( [width * height, 3 ], dtype = np.int32 )

  k = 0

  for i in range ( 0, height ):

    y = float ( height - 1 - i ) / float ( height - 1 )

    for j in range ( 0, width ):

      x = float ( j ) / float ( width - 1 )

      f1 = 4.0 * ( x - 0.5 ) ** 2
      f2 = np.sin ( np.pi * x )
      f3 = x

      if ( y <= f1 ):
        rgb[k,0] = int ( 255.0 * f1 )
      else:
        rgb[k,0] = 50

      if ( y <= f2 ):
        rgb[k,1] = int ( 255.0 * f2 )
      else:
        rgb[k,1] = 150

      if ( y <= f3 ):
        rgb[k,2] = int ( 255.0 * f3 )
      else:
        rgb[k,2] = 250

      k = k + 1

  return rgb

def ppmb_write ( file_name, width, height, rgb ):

#*****************************************************************************80
#
## ppmb_write() writes a binary PPM graphics file.
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    14 September 2018
#
#  Author:
#
#    John Burkardt
#
#  Input:
#
#    string FILE_NAME, the name of the file.
#
#    integer WIDTH, HEIGHT, the width and height of the graphics image.
#
#    integer RGB[3,WIDTH*HEIGHT], a list of WIDTH*HEIGHT triples
#    of R, G, B values between 0 and 255.
#
  import numpy as np
  import struct

  maxval = 255

  file_handle = open ( file_name, 'wb' )
#
#  Set up the header.
#
  ppm_header = f'P6 {width} {height} {maxval}\n'
  file_handle.write ( bytearray ( ppm_header, 'ascii' ) ) 
#
#  Convert 2D array to 1D vector.
#  Depending on how RESHAPE flattens this array, I may need the transpose.
#
  rgbV = np.reshape ( rgb, 3 * width * height )
#
#  Pack entries of vector into a string of bytes, replacing each integer
#  as an unsigned 1 byte character.
#
  rgbB = struct.pack ( '%sB' % len(rgbV), *rgbV )
  file_handle.write ( rgbB )

  file_handle.close ( )

  return

def ppmb_write_test ( ):

#*****************************************************************************80
#
## ppmb_write_test() tests ppmb_write().
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    14 September 2018
#
#  Author:
#
#    John Burkardt
#
  print ( '' )
  print ( 'ppmb_write_test():' )
  print ( '  ppmb_write() writes a binary PPM graphics file.' )

  file_name = 'ppmb_io.ppm'
#
#  Tiny example.
#
  width = 3
  height = 2

  rgb = \
  [ \
    [ 255,   0,   0 ], \
    [   0, 255,   0 ], \
    [   0,   0, 255 ], \
    [ 255, 255,   0 ], \
    [ 255, 255, 255 ], \
    [   0,   0,   0 ] \
  ]
#
#  Large example.
#
  width = 200
  height = 600
  rgb = ppm_example ( width, height )

  ppmb_write ( file_name, width, height, rgb )

  print ( '' )
  print ( '  Graphics data stored in file "%s".' % ( file_name ) )
#
#  Terminate.
#
  print ( '' )
  print ( 'ppmb_write_test():' )
  print ( '  Normal end of execution.' )
  return

def timestamp ( ):

#*****************************************************************************80
#
## timestamp() prints the date as a timestamp.
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    21 August 2019
#
#  Author:
#
#    John Burkardt
#
  import time

  t = time.time ( )
  print ( time.ctime ( t ) )

  return

if ( __name__ == '__main__' ):
  timestamp ( )
  ppmb_write_test ( )
  timestamp ( )