def test_script_to(self):
for D in self.valid_args:
dmi = oc.DMI(D, crystalclass='t')
check_script(dmi._script, contains=['Oxs_DMI_T'])
dmi = oc.DMI(D, crystalclass='o')
check_script(dmi._script, contains=['Oxs_DMI_T'])
def test_skyrmion():
name = "skyrmion"
# Remove any previous simulation directories.
if os.path.exists(name):
shutil.rmtree(name)
mesh = oc.Mesh(p1=(-50e-9, -50e-9, 0),
p2=(50e-9, 50e-9, 10e-9),
cell=(5e-9, 5e-9, 5e-9))
system = oc.System(name="skyrmion")
system.hamiltonian = oc.Exchange(A=1.6e-11) + \
oc.DMI(D=4e-3, kind="interfacial") + \
oc.UniaxialAnisotropy(K1=0.51e6, K2=0.1, u=(0, 0, 1)) + \
oc.Demag() + \
oc.Zeeman(H=(0, 0, 2e5))
Ms = 1.1e6
def Ms_fun(pos):
x, y, z = pos
if (x**2 + y**2)**0.5 < 50e-9:
return Ms
else:
return 0
def m_init(pos):
x, y, z = pos
if (x**2 + y**2)**0.5 < 10e-9:
return (0, 0.1, -1)
else:
return (0, 0.1, 1)
system.m = df.Field(mesh, value=m_init, norm=Ms_fun)
md = oc.MinDriver()
md.drive(system)
# Check the magnetisation at the sample centre.
assert system.m((0, 0, 0))[2] / Ms < -0.97
# Check the magnetisation at the sample edge.
assert system.m((50e-9, 0, 0))[2] / Ms > 0
shutil.rmtree(name)
def test_valueerror(self):
with pytest.raises(ValueError):
dmi = oc.DMI(D=1, crystalclass='t')
dmi._script
def test_script_cnv_win(self):
for D in self.valid_args:
dmi = oc.DMI(D, crystalclass='cnv')
check_script(dmi._script, contains=['Oxs_DMExchange6Ngbr'])
def test_script_cnv_linux_mac(self):
for D in self.valid_args:
dmi = oc.DMI(D, crystalclass='cnv')
check_script(dmi._script, contains=['Oxs_DMI_Cnv'])
def test_script_d2d(self):
for D in self.valid_args:
dmi = oc.DMI(D, crystalclass='d2d')
check_script(dmi._script, contains=['Oxs_DMI_D2d'])