checkpoint --- there are NANs in result

src/fluid_solver.cpp

@@ -171,7 +171,7 @@ fluid_solver<R>::fluid_solver( \
  \
     /* ``physical'' parameters etc */ \
  \
-    this->nu = 0.01; \
+    this->nu = 0.1; \
     this->dkx = DKX; \
     this->dky = DKY; \
     this->dkz = DKZ; \

test.py

@@ -61,10 +61,11 @@ Kdata0[..., 1] = Kdata01
 Kdata0[..., 2] = Kdata02
 Kdata0.tofile("Kdata0")
 run_test('test_step')
-#Kdata1 = np.fromfile('Kdata1', dtype = np.complex64).reshape(Kdata0.shape)
-#
+Kdata1 = np.fromfile('Kdata1', dtype = np.complex64).reshape(Kdata0.shape)
+
 #print np.max(np.abs(Kdata0 - Kdata1))
-#print np.max(np.abs(Kdata0))
+print np.max(np.abs(Kdata0))
+print np.max(np.abs(Kdata1))
 #
 #fig = plt.figure(figsize=(12, 6))
 #a = fig.add_subplot(121)