Update the new code to fit the eigen transition, fix more tests.

src/kbmod/fake_data_creator.py

@@ -110,8 +110,8 @@ def make_fake_ImageStack(self):
         for i in range(self.num_times):
             img = LayeredImage(
                 ("%06i" % i),
-                self.width,
                 self.height,
+                self.width,
                 self.noise_level,
                 self.noise_level**2,
                 self.times[i],
@@ -142,7 +142,7 @@ def insert_object(self, trj):
             # re-set the image. This last step needs to be done
             # explicitly because of how pybind handles references.
             current = self.stack.get_single_image(i)
-            add_fake_object(current, px, py, trj.flux, current.get_psf())
+            add_fake_object(current, py, px, trj.flux, current.get_psf())
 
         # Save the trajectory into the internal list.
         self.trajectories.append(trj)

src/kbmod/search/bindings.cpp

@@ -32,6 +32,7 @@ PYBIND11_MODULE(search, m) {
   search::layered_image_bindings(m);
   search::image_stack_bindings(m);
   search::stack_search_bindings(m);
+  search::stamp_creator_bindings(m);
   search::trajectory_bindings(m);
   search::pixel_pos_bindings(m);
   search::image_moments_bindings(m);

src/kbmod/search/raw_image.cpp

@@ -297,12 +297,12 @@ RawImage& RawImage::operator=(RawImage&& source) {
       for(int j = 0; j < height; ++j) {
         for(int i = 0; i < width; ++i){
           if (bitmask(j, i) == -1){
-            if (((j-1 > 0) && (bitmask(j-1, i) == itr-1)) ||
-                ((i-1 > 0) && (bitmask(j, i-1) == itr-1)) ||
+            if (((j-1 >= 0) && (bitmask(j-1, i) == itr-1)) ||
+                ((i-1 >= 0) && (bitmask(j, i-1) == itr-1)) ||
                 ((j+1 < height) && (bitmask(j+1, i) == itr-1)) ||
                 ((i+1 < width) && (bitmask(j, i+1) == itr-1))){
               bitmask(j, i) = itr;
-            }
+            } 
           }
         } // for i
       } // for j
@@ -351,10 +351,10 @@ RawImage& RawImage::operator=(RawImage&& source) {
             dist2 = new_dist2;
           }
         }
-    }
-
+      } // for x
+    } // for y
     return result;
-}
+  }
 
 
   // Find the basic image moments in order to test if stamps have a gaussian shape.
@@ -581,8 +581,8 @@ RawImage& RawImage::operator=(RawImage&& source) {
           // and value based filtering.
           result(y, x) = 0.0;
         }
-    }
-
+      } // for x
+    } // for y
     return RawImage(result);
   }