#!/usr/bin/python

import pylab as pl
from indep import indep
from move import move

def t1calc(S, Srows, T1, T1rows, sel, pr, Frows, Fbcols, abar, gap=None):
    b = pl.array([1])
    sel = sel.squeeze()
    while (T1rows < Frows - Fbcols) and pl.any(sel==1) and b.size!=0:
        b = indep(S[pl.logical_and(sel, sel),:], gap)
        if b.size != 0:
            b = move(b, pl.find(sel), Srows)
            b = b.squeeze()
            j = b.shape[-1]-1
            while not (b[j] and pl.any(abar==j)):
                j -= 1
                if j == -1:
                   print 'Message from t1calc, called from left_prm'
                   print 'Denominator is noninvertible'
                   exit()
            aa = pl.remainder(pr, Frows) == pl.remainder(j, Frows)
            bb = j<pr
            if not pl.any(bb * aa):
                if T1.size is 0:
                    T1 = b.copy()
                else:
                    T1 = pl.row_stack((T1, b))
                T1rows += 1
            if pr.size is 0:
                pr = pl.array([j])
            else:
                pr = pl.vstack((pr, j))
            while j < Srows:
                sel[j] = 0
                j += Frows
    return T1, T1rows, sel, pr

if __name__== "__main__":
    pass