#!/usr/bin/env python

import svhc
import sys
import pandas as pd
import numpy as np

if __name__=='__main__':
	
	inputfile,T,noise = sys.argv[1:]
	T,noise = int(T),float(noise)
	
	P = np.array(pd.read_csv(inputfile,sep='\t',header=None))
	
	X,comm = svhc.CREATE_DATA(P,T=T,noise=noise)
	
	pd.DataFrame(X).to_csv('dataSeries_benchmark.dat',sep='\t',index=False,header=False)
	with open('HyperNodeList_reference.dat','w') as fw:
		fw.write("\n".join([",".join(map(str,comm[i])) for i in range(len(comm))]))
	
	#pd.DataFrame(hclustval.CompressHC(comm,X.shape[0])).to_csv('CompressDendrogram_reference.dat',sep='\t',index=False,header=False)
