def main(argv=None):
(opts,args)= parser.parse_args()
mydata=load_file(args[0])
claims=mydata.claims
if opts.method==NF3:
strat=gen_strat(mydata.ideal)
for c in claims:
proof(to_if_then(c),strat)
del strat
try:
del c
except NameError:
pass
else:
if opts.method==LINEAR_LEAD_NOREDSB:
reductors=ll_encode(mydata.ideal)
for c in claims:
proofll(to_if_then(c),reductors,redsb=False)
del reductors
try:
del c
except NameError:
pass
else:
reductors=ll_encode(mydata.ideal, reduce=True)
for c in claims:
proofll(to_if_then(c),reductors)
del reductors
try:
del c
except NameError:
pass
return 0
def main(argv=None):
(opts, args) = parser.parse_args()
mydata = load_file(args[0])
claims = mydata.claims
if opts.method == NF3:
strat = gen_strat(mydata.ideal)
for c in claims:
proof(to_if_then(c), strat)
del strat
try:
del c
except NameError:
pass
else:
if opts.method == LINEAR_LEAD_NOREDSB:
reductors = ll_encode(mydata.ideal)
for c in claims:
proofll(to_if_then(c), reductors, redsb=False)
del reductors
try:
del c
except NameError:
pass
else:
reductors = ll_encode(mydata.ideal, reduce=True)
for c in claims:
proofll(to_if_then(c), reductors)
del reductors
try:
del c
except NameError:
pass
return 0
msg = '\nNumber of cycles to repeat: '
RC, k = getInteger(msg)
if ( RC != 0 ):
fpIn.close()
fpOut.close()
sys.exit(3)
random.seed() # Initialize random number generator
###################################################
# Main code #
###################################################
PK, maxNV = readPK(fpIn) # Convert file to PolyBoRi format
maxNV += 1 # Not an index but a quantity
writePolyForPolyBoRi(PK, maxNV, "forPB") # Generate file for Groebner Basis
data = load_file("forPB") # Import file to use in GB
acumm = 0.0
####### Computes Groebner Basis for a set of polynomials equal to a random "y"
for i in range( k-1 ):
polys = copy.copy(data.ideal) # Copy original data ( = 0 )
y = randVect( len(data.ideal) ) # Generates random "y"
polys = equalZero(polys, y) # Set of "polys" - "y" = 0
startTime = time.time()
gb = groebner_basis( polys ) # Computes Groebner Basis
endTime = time.time()
acumm += (endTime - startTime)
print ".",
sys.stdout.flush()
