#!/usr/bin/env python
# Copyright (c) 2013, Ryan Bressler and the Inspectra Contributors 
# All rights reserved.

# Redistribution and use in source and binary forms, with or without
# modification, are permitted provided that the following conditions are met:
#     * Redistributions of source code must retain the above copyright
#       notice, this list of conditions and the following disclaimer.
#     * Redistributions in binary form must reproduce the above copyright
#       notice, this list of conditions and the following disclaimer in the
#       documentation and/or other materials provided with the distribution.
#     * Neither the name of The Institute for Systems Biology, GraphSpectrometer nor the
#       names of its contributors may be used to endorse or promote products
#       derived from this software without specific prior written permission.

# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
# ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
# DISCLAIMED. IN NO EVENT SHALL <COPYRIGHT HOLDER> BE LIABLE FOR ANY
# DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
# (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
# LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
# ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
# (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
# SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
'''
Usage splitfm.py featurematrix.fm targetfeature
'''
import sys
import csv

def main():
	fm = sys.argv[1]
	target = sys.argv[2]

	fo = open(fm)
	featurereader=csv.reader(fo, delimiter='\t')
	clasis={}
	for row in featurereader:
		if row[0]==target:
			for i,v in enumerate(row):
				if i > 0:
					if not v in clasis:
						clasis[v]=[0]
					clasis[v].append(i)

	fo.close()

	outfos = {}

	for key in clasis:
		outfos[key]=open(fm+"."+key,"w")
	
	fo = open(fm)
	featurereader=csv.reader(fo, delimiter='\t')
	for row in featurereader:
		for key in clasis:
			outfos[key].write("\t".join((row[i] for i in clasis[key]))+"\n")

	fo.close()

	for key in outfos:
		outfos[key].close()




	
if __name__ == '__main__':
	main()