diff --git a/ChangeLog b/ChangeLog
index f6983934..c1fc2702 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -7,7 +7,7 @@
 - updated README
 - don't use tables for bilinear on float data for a small speedup (thanks
   Nicolas Robidoux)
-- no tables for uchar either, about a 15% speedup (thank Nicolas)
+- no tables for uchar either, about a 15% speedup (thanks Nicolas)
 
 30/11/10 started 7.24.0
 - bump for new stable
diff --git a/libvips/resample/interpolate.c b/libvips/resample/interpolate.c
index 5b2b6834..cf461403 100644
--- a/libvips/resample/interpolate.c
+++ b/libvips/resample/interpolate.c
@@ -10,6 +10,8 @@
  * 7/1/11
  * 	- don't use tables for bilinear on float data for a small speedup
  * 	  (thanks Nicolas)
+ * 12/1/11
+ * 	- faster, more accuarate uchar bilinear (thanks Nicolas)
  */
 
 /*
@@ -335,9 +337,9 @@ G_DEFINE_TYPE( VipsInterpolateBilinear, vips_interpolate_bilinear,
  */
 #define BILINEAR_INT( TYPE ) { \
 	TYPE *tq = (TYPE *) out; \
- 	\
+	\
 	const int X = (x - ix) * VIPS_INTERPOLATE_SCALE; \
-	const int Y = (y - iy) * VIPS_INTERPOLATE_SCALE; \
+	const int Y = (iy - y) * VIPS_INTERPOLATE_SCALE; \
 	\
 	const TYPE *tp1 = (TYPE *) p1; \
 	const TYPE *tp2 = (TYPE *) p2; \
@@ -350,7 +352,7 @@ G_DEFINE_TYPE( VipsInterpolateBilinear, vips_interpolate_bilinear,
 		const int bot = tp3[z] + \
 			((X * (tp4[z] - tp3[z])) >> VIPS_INTERPOLATE_SHIFT); \
 		\
-		tq[z] = top + ((Y * (bot - top)) >> VIPS_INTERPOLATE_SHIFT); \
+		tq[z] = top - ((Y * (bot - top)) >> VIPS_INTERPOLATE_SHIFT); \
 	} \
 }
 
@@ -381,6 +383,27 @@ G_DEFINE_TYPE( VipsInterpolateBilinear, vips_interpolate_bilinear,
 			c3 * tp3[z] + c4 * tp4[z]; \
 }
 
+/* Float arithmetic, used for int32 and float types, no tables.
+ */
+#define BILINEAR_FLOAT2( TYPE ) { \
+	TYPE *tq = (TYPE *) out; \
+	\
+	const double X = x - ix; \
+	const double Y = iy - y; \
+	\
+	const TYPE *tp1 = (TYPE *) p1; \
+	const TYPE *tp2 = (TYPE *) p2; \
+	const TYPE *tp3 = (TYPE *) p3; \
+	const TYPE *tp4 = (TYPE *) p4; \
+	\
+	for( z = 0; z < b; z++ ) { \
+		const double top = tp1[z] + X * (tp2[z] - tp1[z]); \
+		const double bot = tp3[z] + X * (tp4[z] - tp3[z]); \
+		\
+		tq[z] = top - Y * (bot - top); \
+	} \
+}
+
 /* Expand for band types. with a fixed-point interpolator and a float
  * interpolator.
  */