Set the indexes of lattice sites with zero magnetic moment to -1 in the neighbours matrix (both in micro and atomistic simulations)

Author: davidcorteso

fidimag/atomistic/llg.py

@@ -153,6 +153,9 @@ def set_mu_s(self, value):
             if self._mu_s[i] == 0.0:
                 self._pins[i] = 1
 
+                # Set the neighbour index to -1 for sites with mu_s = 0
+                self.mesh.neighbours[self.mesh.neighbours == i] = -1
+
     mu_s = property(get_mu_s, set_mu_s)
 
     def add(self, interaction, save_field=False):

fidimag/micro/llg.py

@@ -161,6 +161,9 @@ def set_Ms(self, value):
             if self._Ms[i] == 0.0:
                 self._pins[i] = 1
 
+                # Set the neighbour index to -1 for sites with Ms = 0
+                self.mesh.neighbours[self.mesh.neighbours == i] = -1
+
         self.Ms_const = np.max(self._Ms)
 
     Ms = property(get_Ms, set_Ms)