Changeset 243 for trunk/src/pivlab.m

Timestamp:

Apr 19, 2011, 5:41:35 PM (13 years ago)

Author:

sommeria

Message:

bug repair in uvmat

File:

• trunk/src/pivlab.m (modified) (2 diffs)

trunk/src/pivlab.m

@@ -136 +136 @@
     if y <= npy-1 && y >= 1
         f0 = log(result_conv(y,x));
-        f1 = log(result_conv(y-1,x));
-        f2 = log(result_conv(y+1,x));
+        f1 = real(log(result_conv(y-1,x)));
+        f2 = real(log(result_conv(y+1,x)));
         peaky = peaky+ (f1-f2)/(2*f1-4*f0+2*f2);
     else
@@ -145 +145 @@
     if x <= npx-1 && x >= 1
         f0 = log(result_conv(y,x));
-        f1 = log(result_conv(y,x-1));
-        f2 = log(result_conv(y,x+1));
+        f1 = real(log(result_conv(y,x-1)));
+        f2 = real(log(result_conv(y,x+1)));
         peakx = peakx+ (f1-f2)/(2*f1-4*f0+2*f2);
     else