def test_exch_field_oommf(A=1e-11, Ms=2.6e5):
"""
Compare the exchange field from Fidimag with an equivalent
OOMMF simulation. OOMMF field data is taken from an OVF file.
"""
mesh = CuboidMesh(nx=10, ny=3, nz=2, dx=0.5, unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = Ms
exch = UniformExchange(A=A)
sim.add(exch)
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
field = exch.compute_field()
# An equivalent initial magnetisation for OOMMF
# The spatial variables are rescale since they are in nm
init_m0 = (r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] '
+ r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] '
+ r'0 '
+ r'] ')
field_oommf = compute_exch_field(mesh, Ms=Ms, init_m0=init_m0, A=A)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
# Test if the maximum relative errors between both simulations
# is small enough, for every field component
assert max([mx0, mx1, mx2]) < 1e-12
def test_exch_field_oommf(A=1e-11, Ms=2.6e5):
"""
Compare the exchange field from Fidimag with an equivalent
OOMMF simulation. OOMMF field data is taken from an OVF file.
"""
mesh = CuboidMesh(nx=10, ny=3, nz=2, dx=0.5, unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = Ms
exch = UniformExchange(A=A)
sim.add(exch)
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
field = exch.compute_field()
# An equivalent initial magnetisation for OOMMF
# The spatial variables are rescale since they are in nm
init_m0 = (
r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] ' +
r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] ' + r'0 ' + r'] ')
field_oommf = compute_exch_field(mesh, Ms=Ms, init_m0=init_m0, A=A)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
# Test if the maximum relative errors between both simulations
# is small enough, for every field component
assert max([mx0, mx1, mx2]) < 1e-12
def test_with_oommf_spatial_Ms(A=1e-11):
def spatial_Ms(pos):
x, y = pos[0], pos[1]
if x**2 + y**2 < 5**2:
return 2e4
else:
return 0
init_m0 = (
r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] ' +
r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] ' + r'0 ' + r'] ')
init_Ms = """
if { ($x * $x + $y * $y) < 5e-9 * 5e-9 } {
return 2e4
} else {
return 0
}
"""
mesh = CuboidMesh(nx=12, ny=10, nz=2, dx=0.5, unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = spatial_Ms
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
field = exch.compute_field()
field_oommf = compute_exch_field(mesh,
init_m0=init_m0,
A=A,
spatial_Ms=init_Ms)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
assert max([mx0, mx1, mx2]) < 1e-12
field = demag.compute_field()
field_oommf = compute_demag_field(mesh,
spatial_Ms=init_Ms,
init_m0=init_m0)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
assert max([mx0, mx1, mx2]) < 1e-11
def test_with_oommf_spatial_Ms(A=1e-11):
def spatial_Ms(pos):
x, y = pos[0], pos[1]
if x ** 2 + y ** 2 < 5 ** 2:
return 2e4
else:
return 0
init_m0 = (r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] '
+ r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] '
+ r'0 '
+ r'] ')
init_Ms = """
if { ($x * $x + $y * $y) < 5e-9 * 5e-9 } {
return 2e4
} else {
return 0
}
"""
mesh = CuboidMesh(nx=12, ny=10, nz=2, dx=0.5, unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = spatial_Ms
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
field = exch.compute_field()
field_oommf = compute_exch_field(
mesh, init_m0=init_m0, A=A, spatial_Ms=init_Ms)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
assert max([mx0, mx1, mx2]) < 1e-12
field = demag.compute_field()
field_oommf = compute_demag_field(
mesh, spatial_Ms=init_Ms, init_m0=init_m0)
mx0, mx1, mx2 = compare_fields(field_oommf, field)
assert max([mx0, mx1, mx2]) < 1e-11
def test_demag_field_oommf_large(Ms=8e5, A=1.3e-11):
mesh = CuboidMesh(nx=150,
ny=50,
nz=1,
dx=2.5,
dy=2.5,
dz=3,
unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = Ms
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
demag_field = demag.compute_field()
exch_field = exch.compute_field()
#exact = demag.compute_exact()
init_m0 = (
r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] ' +
r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] ' + r'0 ' + r'] ')
demag_oommf = compute_demag_field(mesh, Ms=Ms, init_m0=init_m0)
exch_oommf = compute_exch_field(mesh, Ms=Ms, init_m0=init_m0, A=A)
mx0, mx1, mx2 = compare_fields(demag_oommf, demag_field)
#print mx0, mx1, mx2
assert max([mx0, mx1, mx2]) < 5e-10
mx0, mx1, mx2 = compare_fields(exch_oommf, exch_field)
#print mx0, mx1, mx2
assert max([mx0, mx1, mx2]) < 1e-11
def test_demag_field_oommf_large(Ms=8e5, A=1.3e-11):
mesh = CuboidMesh(nx=150, ny=50, nz=1, dx=2.5, dy=2.5, dz=3, unit_length=1e-9)
sim = Sim(mesh)
sim.Ms = Ms
exch = UniformExchange(A=A)
sim.add(exch)
demag = Demag()
sim.add(demag)
def init_m(pos):
x, y, z = pos
return (np.sin(x) + y + 2.3 * z, np.cos(x) + y + 1.3 * z, 0)
sim.set_m(init_m)
demag_field = demag.compute_field()
exch_field = exch.compute_field()
#exact = demag.compute_exact()
init_m0 = (r'return [list [expr {sin($x * 1e9) + $y * 1e9 + $z * 2.3e9}] '
+ r' [expr {cos($x * 1e9) + $y * 1e9 + $z * 1.3e9}] '
+ r'0 '
+ r'] ')
demag_oommf = compute_demag_field(mesh, Ms=Ms, init_m0=init_m0)
exch_oommf = compute_exch_field(mesh, Ms=Ms, init_m0=init_m0, A=A)
mx0, mx1, mx2 = compare_fields(demag_oommf, demag_field)
#print mx0, mx1, mx2
assert max([mx0, mx1, mx2]) < 5e-10
mx0, mx1, mx2 = compare_fields(exch_oommf, exch_field)
#print mx0, mx1, mx2
assert max([mx0, mx1, mx2]) < 1e-11