####################################################################################

#

# A file that takes in a .dat file (a time series that has been discretiezed)

# and outputs

import cssr_interface

import sys

import pylab

import numpy

import ipdb

# Use CSSR to generate the CSM files

fname = 'timeseries_synthetic/Bernoulli/sample1'

historyLength = 5

is_multiline = False

cssr_interface.run_CSSR(filename = fname, L = historyLength, savefiles = True, showdot = True, is_multiline = is_multiline, showCSSRoutput = True)

hist_length, Cmu, hmu, num_states = cssr_interface.parseResultFile(fname)

print 'C_mu: {0}\nh_mu: {1}\nNumber of States: {2}'.format(Cmu, hmu, num_states)