In this paper, we propose a structure preserving method using a Crank-Nicolson's type method with an implicit Gauss-Seidel fractional iteration. Such a method is of first-order accuracy in time and second-order accuracy in space, stable and length preserving. Such a proposed method brings great benefits for the theoretical analysis. The numerical accuracy, norm preserving and stability are verified for 1D and 3D tests.