#Solver
from scipyAndersonOZsolver import ScipyAndersonOZsolver
#plot
import matplotlib.pyplot as plt
import pylab as pl
if __name__ == '__main__':
solver = ScipyAndersonOZsolver(port = 0)
solver.setPotentialByName('HardSphere')
#solver.setNumberOfIterations(10*solver.getNumberOfIterations())
#Relaxing ConvergenceCriterion by a factor 1000
solver.setConvergenceCriterion(1.0e3*solver.getConvergenceCriterion())
#solver.setPotentialByName('LennardJones', 0.5)
startDensity = 0.58
targetDensity = 0.63
totalDeltaDensity = targetDensity - startDensity
density = startDensity
deltaDensity = totalDeltaDensity
decayFactor = 2.5
deltaDensity = 0.01
print "start approaching target density", targetDensity
while density < targetDensity:
#deltaDensity /= decayFactor
