#!/usr/bin/python

import pylab as pl

def makezero(B, gap=None):

    B = pl.array(B)   # B is a vector!

    B = B.squeeze()

    if len(B.shape)!=1:

        print "B must be a vector!"

        exit()

    if gap == None:

        gap = 1e8

    temp = B[pl.find(B)]

    temp = -pl.sort(-abs(temp))

    ratio = temp[:-1]/pl.float32(temp[1:])

    if pl.find(ratio>gap).size != 0:

        min_ind = min(pl.find(ratio>gap))

        our_eps = temp[min_ind+1]

        zeroind = pl.find(abs(B)<=our_eps)

        B[zeroind] = 0

    return pl.atleast_2d(B) 

if __name__== "__main__":

    B = pl.array([2,3,.000000000004,5,.00006])

    print B

    B = makezero(B)

    print B