"""
Author: A.Alexeevski
Computes Fourier coefficients for  function
"""

from fourier  import *
import optparse
import sys

if len(sys.argv) == 1:
    print( "Computes Fourier coefficients for  function")    
    print( "Use option -h for parameters")
    exit()
####################################################################################################

parser = optparse.OptionParser()
parser.add_option("-i", "--infile", help="File (default func.txt) with function:  ", default="func.txt")
parser.add_option("-o", "--outfile", help="Outfile (default: func_ft.txt) for function", default="func_ft.txt")
parser.add_option("-F", "--noice_F", help="Random F noice (default: 0) in percents: \"-n 5\" ", default=0, type=float)
parser.add_option("-P", "--noice_phi", help="Random phi noice (default: 0) in percents: \"-n 5\" ", default=0, type=float)
opt, args = parser.parse_args()

ff = Fourier()

f = open(opt.infile, 'r')
X,Y = ff.read_func( f )
f.close()   

## noice

ft = ff.fourier(Y)

if abs(opt.noice_F)>= 0.001:
    f=ff.randomize_F(ft,opt.noice_F)
if abs(opt.noice_phi)>= 0.001:
    f=ff.randomize_phi(ft,opt.noice_phi)

g = open(opt.outfile, 'w')
ff.write_f(ft, g)
g.close()

print"..Done"