def relax_system(mesh):
sim = Sim(mesh, name='relax')
sim.set_tols(rtol=1e-6, atol=1e-6)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.6e5
sim.do_procession = False
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
dmi = DMI(D=1.3e-3)
sim.add(dmi)
anis = UniaxialAnisotropy(-3.25e4, axis=(0, 0, 1))
sim.add(anis)
zeeman = Zeeman((0, 0, 6.014576e4))
sim.add(zeeman, save_field=True)
sim.relax(dt=1e-13,
stopping_dmdt=1e-2,
max_steps=5000,
save_m_steps=None,
save_vtk_steps=50)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
sim = Sim(mesh, name='relax')
sim.set_tols(rtol=1e-10, atol=1e-10)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.0e5
sim.do_procession = False
sim.set_m((1, 0.25, 0.1))
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
sim.relax(dt=1e-13,
stopping_dmdt=0.01,
max_steps=5000,
save_m_steps=100,
save_vtk_steps=50)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
sim = Sim(mesh, name='relax')
sim.set_tols(rtol=1e-10, atol=1e-14)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = spatial_Ms
sim.do_procession = False
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
mT = 795.7747154594767
ONE_DEGREE_PER_NS = 17453292.52
sim.relax(dt=1e-13, stopping_dmdt=0.01, max_steps=5000,
save_m_steps=100, save_vtk_steps=50)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
sim = Sim(mesh, name='relax')
sim.set_tols(rtol=1e-10, atol=1e-14)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = spatial_Ms
sim.do_procession = False
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
mT = 795.7747154594767
ONE_DEGREE_PER_NS = 17453292.52
sim.relax(dt=1e-13,
stopping_dmdt=0.01,
max_steps=5000,
save_m_steps=100,
save_vtk_steps=50)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
# Only relaxation
sim = Sim(mesh, name='relax')
# Simulation parameters
sim.set_tols(rtol=1e-8, atol=1e-10)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.6e5
sim.do_procession = False
# The initial state passed as a function
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
# Energies
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
anis = UniaxialAnisotropy(5e4)
sim.add(anis)
# dmi = DMI(D=8e-4)
# sim.add(dmi)
# mT = 795.7747154594767
# ONE_DEGREE_PER_NS = 17453292.52
# Start relaxation and save the state in m0.npy
sim.relax(dt=1e-14, stopping_dmdt=0.00001, max_steps=5000,
save_m_steps=None, save_vtk_steps=None)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
sim = Sim(mesh, name="relax")
sim.set_tols(rtol=1e-6, atol=1e-6)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.6e5
sim.do_procession = False
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
dmi = DMI(D=1.3e-3)
sim.add(dmi)
anis = UniaxialAnisotropy(-3.25e4, axis=(0, 0, 1))
sim.add(anis)
zeeman = Zeeman((0, 0, 6.014576e4))
sim.add(zeeman, save_field=True)
sim.relax(dt=1e-13, stopping_dmdt=1e-2, max_steps=5000, save_m_steps=None, save_vtk_steps=50)
np.save("m0.npy", sim.spin)
def relax_system(mesh):
# Only relaxation
sim = Sim(mesh, name='relax')
# Simulation parameters
sim.set_tols(rtol=1e-8, atol=1e-10)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.6e5
sim.do_procession = False
# The initial state passed as a function
sim.set_m(init_m)
# sim.set_m(np.load('m0.npy'))
# Energies
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
anis = UniaxialAnisotropy(5e4)
sim.add(anis)
# dmi = DMI(D=8e-4)
# sim.add(dmi)
# mT = 795.7747154594767
# ONE_DEGREE_PER_NS = 17453292.52
# Start relaxation and save the state in m0.npy
sim.relax(dt=1e-14,
stopping_dmdt=0.00001,
max_steps=5000,
save_m_steps=None,
save_vtk_steps=None)
np.save('m0.npy', sim.spin)
def relax_system(mesh):
sim = Sim(mesh, name="relax")
sim.set_tols(rtol=1e-10, atol=1e-10)
sim.alpha = 0.5
sim.gamma = 2.211e5
sim.Ms = 8.0e5
sim.do_procession = False
sim.set_m((1, 0.25, 0.1))
# sim.set_m(np.load('m0.npy'))
A = 1.3e-11
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
sim.relax(dt=1e-13, stopping_dmdt=0.01, max_steps=5000, save_m_steps=100, save_vtk_steps=50)
np.save("m0.npy", sim.spin)
