def test_sim_spins(do_plot=False):
mesh = CuboidMesh(nx=10, ny=5, nz=1)
sim = Sim(mesh, name='10spin')
alpha = 0.1
gamma = 2.21e5
sim.alpha = alpha
sim.gamma = gamma
sim.mu_s = 1.0
sim.set_m((1, 0, 0))
print sim.spin
H0 = 1e5
sim.add(Zeeman((0, 0, H0)))
ts = np.linspace(0, 1e-9, 101)
mx = []
my = []
mz = []
real_ts = []
for t in ts:
sim.run_until(t)
real_ts.append(sim.t)
#print sim.t, abs(sim.spin_length()[0] - 1)
av = sim.compute_average()
mx.append(av[0])
my.append(av[1])
mz.append(av[2])
#sim.save_vtk()
mz = np.array(mz)
# print mz
a_mx, a_my, a_mz = single_spin(alpha, gamma, H0, ts)
print sim.stat()
if do_plot:
plot(real_ts,
mx,
my,
mz,
a_mx,
a_my,
a_mz,
name='spins.pdf',
title='integrating spins')
print("Max Deviation = {0}".format(np.max(np.abs(mz - a_mz))))
assert np.max(np.abs(mz - a_mz)) < 5e-7
def test_sim_spins(do_plot=False):
mesh = CuboidMesh(nx=10, ny=5, nz=1)
sim = Sim(mesh, name='10spin')
alpha = 0.1
gamma = 2.21e5
sim.alpha = alpha
sim.gamma = gamma
sim.mu_s = 1.0
sim.set_m((1, 0, 0))
print sim.spin
H0 = 1e5
sim.add(Zeeman((0, 0, H0)))
ts = np.linspace(0, 1e-9, 101)
mx = []
my = []
mz = []
real_ts = []
for t in ts:
sim.run_until(t)
real_ts.append(sim.t)
#print sim.t, abs(sim.spin_length()[0] - 1)
av = sim.compute_average()
mx.append(av[0])
my.append(av[1])
mz.append(av[2])
#sim.save_vtk()
mz = np.array(mz)
# print mz
a_mx, a_my, a_mz = single_spin(alpha, gamma, H0, ts)
print sim.stat()
if do_plot:
plot(real_ts, mx, my, mz, a_mx, a_my, a_mz, name='spins.pdf', title='integrating spins')
print("Max Deviation = {0}".format(
np.max(np.abs(mz - a_mz))))
assert np.max(np.abs(mz - a_mz)) < 5e-7
def test_sim_single_spin_vode(do_plot=False):
mesh = CuboidMesh(nx=1, ny=1, nz=1)
sim = Sim(mesh, name='spin')
alpha = 0.1
gamma = 2.21e5
sim.alpha = alpha
sim.gamma = gamma
sim.mu_s = 1.0
sim.set_m((1, 0, 0))
H0 = 1e5
sim.add(Zeeman((0, 0, H0)))
ts = np.linspace(0, 1e-9, 101)
mx = []
my = []
mz = []
real_ts = []
for t in ts:
sim.run_until(t)
real_ts.append(sim.t)
#print sim.t, abs(sim.spin_length()[0] - 1)
mx.append(sim.spin[0])
my.append(sim.spin[1])
mz.append(sim.spin[2])
mz = np.array(mz)
a_mx, a_my, a_mz = single_spin(alpha, gamma, H0, ts)
print sim.stat()
if do_plot:
plot(real_ts, mx, my, mz, a_mx, a_my, a_mz)
print("Max Deviation = {0}".format(
np.max(np.abs(mz - a_mz))))
assert np.max(np.abs(mz - a_mz)) < 5e-7
def test_sim_single_spin_vode(do_plot=False):
mesh = CuboidMesh(nx=1, ny=1, nz=1)
sim = Sim(mesh, name='spin')
alpha = 0.1
gamma = 2.21e5
sim.alpha = alpha
sim.gamma = gamma
sim.mu_s = 1.0
sim.set_m((1, 0, 0))
H0 = 1e5
sim.add(Zeeman((0, 0, H0)))
ts = np.linspace(0, 1e-9, 101)
mx = []
my = []
mz = []
real_ts = []
for t in ts:
sim.run_until(t)
real_ts.append(sim.t)
#print sim.t, abs(sim.spin_length()[0] - 1)
mx.append(sim.spin[0])
my.append(sim.spin[1])
mz.append(sim.spin[2])
mz = np.array(mz)
a_mx, a_my, a_mz = single_spin(alpha, gamma, H0, ts)
print sim.stat()
if do_plot:
plot(real_ts, mx, my, mz, a_mx, a_my, a_mz)
print("Max Deviation = {0}".format(np.max(np.abs(mz - a_mz))))
assert np.max(np.abs(mz - a_mz)) < 5e-7