Exemple #1
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
Exemple #2
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
Exemple #3
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
Exemple #4
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
Exemple #5
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
Exemple #6
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
Exemple #7
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
Exemple #8
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