bug fix

Author: Weiwei Wang

examples/mc/skx.py

@@ -12,12 +12,15 @@ def init_m(pos):
 
     return (0,0,1)
 
+def random_m(pos):
+    return np.random.random(3) - 0.5
+
 
 def run(mesh):
 
     mc = MonteCarlo(mesh, name='test1')
     mc.set_m(init_m)
-    mc.set_options(H=[0,0,0.5], J=50.0, D=0.27*50, T=5.0)
+    mc.set_options(H=[0,0,0.0], J=50.0, D=0.27*50, T=5.0)
     mc.run(steps=200000, save_m_steps=None, save_vtk_steps=1000, save_data_steps=10)
 
 if __name__=='__main__':

fidimag/atomistic/lib/mc.c

@@ -2,7 +2,7 @@
 
 void uniform_random_sphere(double *phi, double *ct, int n){
 	for (int i=0;i<n;i++){
-		phi[i] = single_random()*2*M_PI;
+		phi[i] = single_random()*2*WIDE_PI;
 		ct[i] = 2*single_random()-1;
 	}
 }
@@ -13,7 +13,7 @@ void uniform_random_sphere(double *phi, double *ct, int n){
 void random_spin_uniform(double *spin, int n){
     for (int i=0;i<n;i++){
         int j=3*i;
-        double phi= single_random()*2*M_PI;
+        double phi= single_random()*2*WIDE_PI;
         double ct = 2*single_random()-1;
         double st = sqrt(1-ct*ct);
         spin[j] = st*cos(phi); //mx = sin(theta)*cos(phi)