acom v2 convorder experiments 120-->80

tests/acom_v2_convorder_tps1_theta1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+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 \
+  , T_start=0.0, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+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_V2_ConvOrder_TPS1_theta1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v2_convorder_tps1_theta5e-1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+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 \
+  , T_start=0.0, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+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_V2_ConvOrder_TPS1_theta5e-1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v2_convorder_tps1pf_theta1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+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 \
+  , T_start=0.0, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+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_V2_ConvOrder_TPS1PF_theta1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v2_convorder_tps1pf_theta5e-1.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+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 \
+  , T_start=0.0, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+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_V2_ConvOrder_TPS1PF_theta5e-1_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()

tests/acom_v2_convorder_tps2.py

+from commix import *
+from commix.addons import ConvergenceOrderWrapper
+
+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 \
+  , T_start=0.0, T_end=0.91e-09, timeStepSize=1e-13/64)
+
+from math import sqrt
+e = 1.0/100.0
+par.m0 = (e, -e, sqrt(1-2*e*e))
+
+out = OutputOptions()
+out.video = False
+
+wrapper = ConvergenceOrderWrapper(referenceScheme=TPS2AB, otherScheme=TPS2AB, nrSimulations=5, geometry=geo, parameters=par, output=out, numthreads=4)
+
+wrapper.SetResultsFolder("ACOM_V2_ConvOrder_TPS2_Structured_Result")
+
+wrapper.GetReady()
+wrapper.Do()