def sundials_jtn(self, mp, Jmp, t, m, fy): # we can not copy mp to self.spin since m and self.spin is one object. # self.spin[:] = mp[:] print "NO jac..........." self.compute_effective_field_jac(t, mp) clib.compute_llg_jtimes( Jmp, m, fy, mp, self.field, self.alpha, self._pins, self.gamma, self.n, self.do_procession, self.default_c ) return 0
def sundials_jtn(self, mp, Jmp, t, m, fy): # we can not copy mp to self.spin since m and self.spin is one object. #self.spin[:] = mp[:] print('NO jac...........') self.compute_effective_field_jac(t, mp) clib.compute_llg_jtimes(Jmp, m, fy, mp, self.field, self.alpha, self._pins, self.gamma, self.n, self.do_precession, self.default_c) return 0
def sundials_jtimes(self, mp, Jmp, t, m, fy): self.compute_effective_field_jac(t, mp) clib.compute_llg_jtimes(Jmp, m, fy, mp, self.field, self.alpha, self._pins, self.gamma, self.n, self.do_precession, self.default_c) return 0