#!/usr/bin/python
from socal_caller import *
from cvxopt import matrix
from optparse import OptionParser
from random import shuffle
from socal_trainer import *
import sys
import math
import time
c_aa = matrix([0,0],(2,1))
E_aa = matrix([1,0,0,1],(2,2))
c_ab = matrix([1,1],(2,1))
E_ab = matrix([1,0,0,1],(2,2))
c_bb = matrix([2,2],(2,1))
E_bb = matrix([1,0,0,1],(2,2))
sc = socal_caller(c_aa,E_aa,c_ab,E_ab,c_bb,E_ab,32)
xs = [matrix([0,0],(2,1)),matrix([1,1],(2,1)),matrix([2,2],(2,1))]
for x in xs:
print sc.mindist(x)
c_aa = None
E_aa = None
c_ab = None
E_ab = None
c_bb = None
E_bb = None
if(e_aa != None):
c_aa = e_aa['c']
E_aa = e_aa['E']
if(e_ab != None):
c_ab = e_ab['c']
E_ab = e_ab['E']
if(e_bb != None):
c_bb = e_bb['c']
E_bb = e_bb['E']
sc = socal_caller(c_aa,E_aa,c_ab,E_ab,c_bb,E_bb,c4,e_aa,e_ab,e_bb)
# run the caller
for i in xrange(len(indv_list)):
info = cols[i+3].split(':')
geno = info[1]
sensea_int = float(info[0].split('/')[0])
senseb_int = float(info[0].split('/')[1])
# do calling
x = matrix([sensea_int, senseb_int],(2,1))
result = sc.mindist(x)
call = result[0]
score = result[1]
# write out to file
