Expose even more terms to Python from C Struct

Author: rpep

fidimag/common/dipolar/dipolar.pyx

@@ -18,6 +18,9 @@ cdef extern from "dipolar.h":
         double *tensor_yy
         double *tensor_yz
         double *tensor_zz
+        double *hx
+        double *hy
+        double *hz
         complex *Nxx
         complex *Nxy
         complex *Nxz
@@ -47,12 +50,12 @@ cdef class FFTDemag(object):
     cdef fft_demag_plan *_c_plan
     cdef int total_length
     cdef np.float64_t[:] tensor_xx_p, tensor_xy_p, tensor_xz_p, tensor_yy_p, \
-                         tensor_yz_p, tensor_zz_p
+                         tensor_yz_p, tensor_zz_p, hx_p, hy_p, hz_p
     cdef np.complex128_t[:] Nxx_p, Nxy_p, Nxz_p, Nyy_p, Nyz_p, Nzz_p, Hx_p, \
                             Hy_p, Hz_p, Mx_p, My_p, Mz_p
     cdef public np.ndarray tensor_xx, tensor_xy, tensor_xz, tensor_yy, \
                          tensor_yz, tensor_zz, Nxx, Nxy, Nxz, Nyy, Nyz, Nzz, \
-                         Mx, My, Mz, Hx, Hy, Hz
+                         Mx, My, Mz, Hx, Hy, Hz, hx, hy, hz
     #tensor_type could be 'dipolar', 'demag' or '2d_pbc'
     def __cinit__(self, dx, dy, dz, nx, ny, nz, tensor_type='dipolar'):
         self._c_plan = create_plan()
@@ -103,6 +106,13 @@ cdef class FFTDemag(object):
         self.Hz_p = <np.complex128_t[:self.total_length]> self._c_plan.Hz
         self.Hz = np.asarray(self.Hz_p)
 
+        self.hx_p = <np.float64_t[:self.total_length]> self._c_plan.hx
+        self.hx = np.asarray(self.hx_p)
+        self.hy_p = <np.float64_t[:self.total_length]> self._c_plan.hy
+        self.hy = np.asarray(self.hy_p)
+        self.hz_p = <np.float64_t[:self.total_length]> self._c_plan.hz
+        self.hz = np.asarray(self.hz_p)
+
 
         self.tensor_xy_p = <np.float64_t[:self.total_length]> self._c_plan.tensor_xy
         self.tensor_xy = np.asarray(self.tensor_xy_p)