# HG changeset patch
# User Britton Smith <brittonsmith@gmail.com>
# Date 1247594311 21600
# Branch raytrace
# Node ID 8d36a87a55153e4f449c50c3800d4058f5973965
# Parent  f0df96c8e8874fe1dd5526adb6f4254a6e35ec6c
Adding testing print statements.

diff -r f0df96c8e8874fe1dd5526adb6f4254a6e35ec6c -r 8d36a87a55153e4f449c50c3800d4058f5973965 yt/lagos/RTIntegrator.pyx
--- a/yt/lagos/RTIntegrator.pyx	Tue Jul 14 10:25:58 2009 -0600
+++ b/yt/lagos/RTIntegrator.pyx	Tue Jul 14 11:58:31 2009 -0600
@@ -117,7 +117,8 @@
     for i in range(3):
         cur_ind[i] = np.floor((intersect[i] - left_edge[i])/dx[i])
         tmax[i] = (((cur_ind[i]+step[i])*dx[i])+left_edge[i]-u[i])/v[i]
-        if step[i] < 0: cur_ind[i] -= 1
+        if step[i] < 0:
+            cur_ind[i] -= 1
         tdelta[i] = abs(dx[i]/v[i])
     # The variable intersect contains the point we first pierce the grid
     enter_t = intersect_t
@@ -127,13 +128,13 @@
            not (0 <= cur_ind[1] < grid_mask.shape[1]) or \
            not (0 <= cur_ind[2] < grid_mask.shape[2]):
             break
-#        else:
-#            grid_mask[cur_ind[0], cur_ind[1], cur_ind[2]] = 1
         # Note that we are calculating t on the fly, but we get *negative* t
         # values from what they should be.
         # If we've reached t = 1, we are done.
         if tmax[0] >= 1 and tmax[1] >= 1 and tmax[2] >= 1:
             grid_t[cur_ind[0], cur_ind[1], cur_ind[2]] = 1.0 - enter_t
+            print "Ending index: ", cur_ind[0],cur_ind[1],cur_ind[2]
+            print "tmax: ",tmax[0],tmax[1],tmax[2]
             if grid_t[cur_ind[0], cur_ind[1], cur_ind[2]] > 0:
                 grid_mask[cur_ind[0], cur_ind[1], cur_ind[2]] = 1
             break