def __init__(self, mesh, chi=1e-3, name='unnamed'):
self.chi = chi
super(LLBarFull, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.add(Relaxation(chi))
self.beta = 0
def __init__(self, mesh, name='unnamed'):
super(LLBar, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.field_perp = np.zeros(3 * self.n, dtype=np.float)
self.beta = 0
class LLBarFull(LLG):
def __init__(self, mesh, chi=1e-3, name='unnamed'):
self.chi = chi
super(LLBarFull, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.add(Relaxation(chi))
self.beta = 0
def sundials_rhs(self, t, y, ydot):
self.t = t
# already synchronized when call this funciton
# self.spin[:]=y[:]
self.compute_effective_field(t)
delta_h = self.lap.compute_laplace_field(self.field, self._Ms)
clib.compute_llg_rhs_baryakhtar(ydot, self.spin, self.field, delta_h,
self.alpha, self.beta, self._pins,
self.gamma, self.n, self.do_procession)
#ydot[:] = self.dm_dt[:]
return 0
class LLBar(LLG):
def __init__(self, mesh, name='unnamed'):
super(LLBar, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.field_perp = np.zeros(3 * self.n, dtype=np.float)
self.beta = 0
def sundials_rhs(self, t, y, ydot):
self.t = t
# already synchronized when call this funciton
# self.spin[:]=y[:]
self.compute_effective_field(t)
clib.compute_perp_field(self.spin, self.field, self.field_perp, self.n)
delta_h = self.lap.compute_laplace_field(self.field_perp, self._Ms)
clib.compute_llg_rhs_baryakhtar_reduced(ydot, self.spin, self.field,
delta_h, self.alpha, self.beta,
self._pins, self.gamma, self.n,
self.do_procession,
self.default_c)
#ydot[:] = self.dm_dt[:]
return 0
class LLBarFull(LLG):
def __init__(self, mesh, chi=1e-3, name='unnamed'):
self.chi = chi
super(LLBarFull, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.add(Relaxation(chi))
self.beta = 0
def sundials_rhs(self, t, y, ydot):
self.t = t
# already synchronized when call this funciton
# self.spin[:]=y[:]
self.compute_effective_field(t)
delta_h = self.lap.compute_laplace_field(self.field, self._Ms)
clib.compute_llg_rhs_baryakhtar(ydot,
self.spin,
self.field,
delta_h,
self.alpha,
self.beta,
self._pins,
self.gamma,
self.n,
self.do_precession)
#ydot[:] = self.dm_dt[:]
return 0
def __init__(self, mesh, chi=1e-3, name='unnamed'):
self.chi = chi
super(LLBarFull, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.add(Relaxation(chi))
self.beta = 0
def __init__(self, mesh, name='unnamed'):
super(LLBar, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.field_perp = np.zeros(3 * self.n, dtype=np.float)
self.beta = 0
class LLBar(LLG):
def __init__(self, mesh, name='unnamed'):
super(LLBar, self).__init__(mesh, name=name)
self.lap = Laplace(mesh)
self.field_perp = np.zeros(3 * self.n, dtype=np.float)
self.beta = 0
def sundials_rhs(self, t, y, ydot):
self.t = t
# already synchronized when call this funciton
# self.spin[:]=y[:]
self.compute_effective_field(t)
clib.compute_perp_field(
self.spin, self.field, self.field_perp, self.n)
delta_h = self.lap.compute_laplace_field(self.field_perp, self._Ms)
clib.compute_llg_rhs_baryakhtar_reduced(ydot,
self.spin,
self.field,
delta_h,
self.alpha,
self.beta,
self._pins,
self.gamma,
self.n,
self.do_precession,
self.default_c)
#ydot[:] = self.dm_dt[:]
return 0
