acom v3 convorder experiments --> start at 5ps

tests/acom_v3_convorder_tps1_theta1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+geo = Geometry(geometry="data_ACOM/acom_meshStructured2.5nmCubes.vol", scaling=1.0e-09)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , T_start=5e-12, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+par.m0 = "data_ACOM/acom_tilted_5ps.vtk"
+
+out = OutputOptions()
+out.video = False
+
+otherPar = par.DeepCopy()
+otherPar.theta = 1.0
+otherPar.dmImplicitTreatment = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS1, nrSimulations=5, geometry=geo, parameters=par, otherParameters=otherPar, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V3_ConvOrder_TPS1_theta1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v3_convorder_tps1_theta5e-1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+geo = Geometry(geometry="data_ACOM/acom_meshStructured2.5nmCubes.vol", scaling=1.0e-09)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , T_start=5e-12, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+par.m0 = "data_ACOM/acom_tilted_5ps.vtk"
+
+out = OutputOptions()
+out.video = False
+
+otherPar = par.DeepCopy()
+otherPar.theta = 0.5
+otherPar.dmImplicitTreatment = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS1, nrSimulations=5, geometry=geo, parameters=par, otherParameters=otherPar, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V3_ConvOrder_TPS1_theta5e-1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v3_convorder_tps1pf_theta1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+geo = Geometry(geometry="data_ACOM/acom_meshStructured2.5nmCubes.vol", scaling=1.0e-09)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , T_start=5e-12, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+par.m0 = "data_ACOM/acom_tilted_5ps.vtk"
+
+out = OutputOptions()
+out.video = False
+
+otherPar = par.DeepCopy()
+otherPar.theta = 1.0
+otherPar.dmImplicitTreatment = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS1PF, nrSimulations=5, geometry=geo, parameters=par, otherParameters=otherPar, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V3_ConvOrder_TPS1PF_theta1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v3_convorder_tps1pf_theta5e-1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+geo = Geometry(geometry="data_ACOM/acom_meshStructured2.5nmCubes.vol", scaling=1.0e-09)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , T_start=5e-12, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+par.m0 = "data_ACOM/acom_tilted_5ps.vtk"
+
+out = OutputOptions()
+out.video = False
+
+otherPar = par.DeepCopy()
+otherPar.theta = 0.5
+otherPar.dmImplicitTreatment = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS1PF, nrSimulations=5, geometry=geo, parameters=par, otherParameters=otherPar, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V3_ConvOrder_TPS1PF_theta5e-1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v3_convorder_tps2.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+geo = Geometry(geometry="data_ACOM/acom_meshStructured2.5nmCubes.vol", scaling=1.0e-09)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , T_start=5e-12, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+par.m0 = "data_ACOM/acom_tilted_5ps.vtk"
+
+out = OutputOptions()
+out.video = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS2AB, nrSimulations=5, geometry=geo, parameters=par, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V3_ConvOrder_TPS2_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v3_initial_condition.py

+from commix import *
+
+h = 10 / 4
+
+geo = Geometry(geometry=Cuboid(80, 80, 10), meshSize=(h, h, h), scaling=1.0e-09, structuredMesh=True)
+
+from math import pi
+par = Parameters(A=pi*2.0e-11, Ms=1.0e6, alpha=0.08, K=0.0, D=pi*2.0e-3 \
+  , maxwellCoupling = "none", dmCoupling="bulk", TOL=1e-12, dmImplicitTreatment=True \
+  , timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+exp = Integrator(scheme=TPS2AB, geometry=geo, parameters=par, numthreads=4)
+exp.SetResultsFolder("ACOM_V3_Convorder_InitialCondition_Result")
+exp.RecordMagnetization(record=True, maxwell=False, every=0.1e-12)
+
+exp.Integrate(duration=5e-12)