def _log_current_status(contigs,cluster_prob,p,z,run):
clustering = _get_current_clustering(contigs,p,z)
Output.write_clustering_result(clustering,cluster_prob ,p,arguments=None,tmpfile=True,tmpfile_suffix=run)
#Prep output settings
#=============================
#TODO: Needs to be fixed to allow only coverage etc.
Output.set_output_path(outfile,args)
#=============================
#Calling clustering
#=============================
(clusters,clust_prob,centroids) = main(cluster_func, contigs, p, args.cluster_count,args.epsilon,args.iterations, \
args.runs,args.verbose,args.serial, **params)
#=============================
#Printing Results
#=============================
Output.write_clustering_result(clusters,clust_prob,centroids,idx,args)
#=============================
#Preprocess timeseries data for coverage
#=============================
elif args.script == 'preprocess':
if args.output:
args.output = open(args.output,'w+')
else:
args.output = sys.stdout
if not args.contigs:
print "Contigs file not correctly supplied, will now exit"
sys.exit(-1)
main_preprocess(args)
else:
pass
#Prep output settings
#=============================
#TODO: Needs to be fixed to allow only coverage etc.
Output.set_output_path(outfile,args)
#=============================
#Calling clustering
#=============================
(clusters,clust_prob,responsibilities, probability_parameters) = main(cluster_func, contigs, p, args.cluster_count,args.epsilon,args.iterations, \
args.runs,args.verbose,args.serial, **params)
#=============================
#Printing Results
#=============================
Output.write_clustering_result(clusters,clust_prob,responsibilities, probability_parameters,idx,args)
#=============================
#Preprocess timeseries data for coverage
#=============================
elif args.script == 'preprocess':
if args.output:
args.output = open(args.output,'w+')
else:
args.output = sys.stdout
if not args.contigs:
print "Contigs file not correctly supplied, will now exit"
sys.exit(-1)
main_preprocess(args)
else:
pass
