options = {
'MaxIters': 5000,
'MaxFunEvals': 5000,
'NormGradTol': 1e-6,
'FunValDiff': 1e-0,
'StepLength': 2,
'StepTol': 1e-14,
'Output': True
}
moptions = scftmodel2d_options(nspecies=2,
nblend=1,
nblock=2,
ndeg=100,
fA=0.2,
chiAB=0.25,
dim=2,
T0=20,
T1=80,
nupdate=1,
order=1,
rdir=sys.argv[3])
if sys.argv[1] == 'quadtree':
mesh = init_mesh(n=5, h=12)
mesh = mesh.to_pmesh()
elif sys.argv[1] == 'halfedge':
mesh = halfedgemesh(n=6, h=12)
elif sys.argv[1] == 'complex':
mesh = complex_mesh(r=20, filename=sys.argv[2])
Halftest = HalfEdgeAVEMTest(mesh,
fieldstype=3,
options = {
'MaxIters': 5000,
'MaxFunEvals': 5000,
'NormGradTol': 1e-6,
'FunValDiff': 1e-6,
'StepLength': 2,
'StepTol': 1e-14,
'etarefTol': 0.1,
'Output': True
}
moptions = scftmodel2d_options(nspecies=2,
nblend=1,
nblock=2,
ndeg=100,
fA=0.5,
chiAB=0.20,
dim=2,
T0=40,
T1=160,
nupdate=1,
order=2,
rdir=sys.argv[1])
Halftest = HalfEdgeAVEMTest(fieldstype=4,
moptions=moptions,
optoptions=options)
Halftest.run()
options = {
'MaxIters': 5000,
'MaxFunEvals': 5000,
'NormGradTol': 1e-7,
'FunValDiff': 1e-6,
'StepLength': 2,
'StepTol': 1e-14,
'Output': True
}
moptions = scftmodel2d_options(nspecies=2,
nblend=1,
nblock=2,
ndeg=100,
fA=0.2,
chiAB=0.25,
dim=2,
T0=10,
T1=40,
nupdate=1,
order=2)
mesh = halfedgemesh(n=1)
Halftest = HalfEdgeAVEMTest(mesh,
fieldstype=3,
moptions=moptions,
optoptions=options)
Halftest.run()