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
