#!/usr/bin/python
"""This python script uses the function thinfilm_simulation within the
file simulation.py in order to perform a parameter analysis of
hlib. Before the simulations can be run one has to create the
corresponding meshes by using the shell script setup.sh."""
import simulation
quadrature=[i for i in range(3,6)]
leafsize = [i*10 for i in range(1,5)]
polynome = [i for i in range(3,6)]
geo_criterion = [i*1.0 for i in range(1,5)]
accuracy = [10**i for i in range(-6,-2)]
#quadrature=[i for i in range(4,6)]
#leafsize = [i*5 for i in range(4,5)]
#polynome = [i for i in range(4,5)]
#geo_criterion = [i*0.5 for i in range(4,5)]
#accuracy = [10**i for i in range(-4,-3)]
for nfdeg in quadrature:
for nmin in leafsize:
for p in polynome:
for eta in geo_criterion:
for eps in accuracy:
simulation.iterate_hlib_simulations(nfdeg=3,nmin=20,p=4,eta=2.0,eps=0.0001,dirlist=['thinfilm50_50_5'])
print "Simulation has finished.\n"
for algorithm in algorithmlist:
for nfdeg in nfdeglist:
for nmin in nminlist:
for eta in etalist:
for eps_aca in eps_acalist:
for eps in epslist:
for p in plist:
for kmax in kmaxlist:
simulation.iterate_hlib_simulations(
algorithm=algorithm,
nfdeg=nfdeg,
nmin=nmin,
eta=eta,
eps_aca=eps_aca,
eps=eps,
p=p,
kmax=kmax,
dirlist=directories,
)
if algorithm > 1 and algorithm != 5:
break
if algorithm <= 1 or algorithm == 5:
break
if algorithm == 5 or algorithm == 6:
break
if algorithm == 2 or algorithm == 7:
break
