#Solver
from picardOZsolver import PicardOZsolver
from andersonOZsolver import AndersonOZsolver
from scipyAndersonOZsolver import ScipyAndersonOZsolver
from scipyNewtonKrylovOZsolver import ScipyNewtonKrylovOZsolver
if __name__ == '__main__':
#Instantiate solver class (Picard iteration, etc.)
#For star potential, ScipyAndersonOZsolver gives best result
#s = PicardOZsolver(port = 0)
#s = AndersonOZsolver(port = 0)
s = ScipyAndersonOZsolver(port = 0)
#s = ScipyNewtonKrylovOZsolver(port = 0)
s.printPotentialSetterArguments()
#s.setPotentialByName('HardSphere')
#s.setPotentialByName('LennardJones', 0.2)
s.setPotentialByName('Star', 20)
#Define the density range to scan
densityRange = np.arange(0.3, 0.4, 0.1)
#densityRange = np.asarray([0.3])
#densityRange = np.asarray([0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5])
densityRDFdictionary = {}
#Start loop over densities and measure wall clock time
t_start = time.time()
