#!/usr/bin/python -u
  
  # ./gen_dct_mat.py
  # script generates matrix with DCT transform
  #     
  # author: Karel Vesely
  #
  
  from math import *
  import sys
  
  
  from optparse import OptionParser
  
  parser = OptionParser()
  parser.add_option('--fea-dim', dest='dim', help='feature dimension')
  parser.add_option('--splice', dest='splice', help='applied splice value')
  parser.add_option('--dct-basis', dest='dct_basis', help='number of DCT basis')
  (options, args) = parser.parse_args()
  
  if(options.dim == None):
      parser.print_help()
      sys.exit(1)
  
  dim=int(options.dim)
  splice=int(options.splice)
  dct_basis=int(options.dct_basis)
  
  timeContext=2*splice+1
  
  
  #generate the DCT matrix
  M_PI = 3.1415926535897932384626433832795
  M_SQRT2 = 1.4142135623730950488016887
  
  
  #generate small DCT matrix
  print '['
  for k in range(dct_basis):
      for m in range(dim):
          for n in range(timeContext):
            if(n==0): 
                print m*'0 ',
            else: 
                print (dim-1)*'0 ',
            print str(sqrt(2.0/timeContext)*cos(M_PI/timeContext*k*(n+0.5))),
            if(n==timeContext-1):
                print (dim-m-1)*'0 ',
          print
      print 
  
  print ']'