def compute_field_uniform(self, t=0, spin=None): m = spin if spin is not None else self.spin clib.compute_exchange_field(m, self.field, self.mu_s_inv, self.energy, self.Jx, self.Jy, self.Jz, self.neighbours, self.n, self.n_ngbs) return self.field
def compute_field(self, t=0, spin=None): if spin is not None: m = spin else: m = self.spin clib.compute_exchange_field(m, self.field, self.energy, self.Jx, self.Jy, self.Jz, self.neighbours, self.n) return self.field * self.mu_s_inv