Explorar o código

findpeaks: pcurv can be zero, avoid that...

Fernando Pablo Lopez-Lezcano %!s(int64=8) %!d(string=hai) anos
pai
achega
a786ab3018
Modificáronse 1 ficheiros con 5 adicións e 1 borrados
  1. 5 1
      calibration/src/sphear_findpeaks.m

+ 5 - 1
calibration/src/sphear_findpeaks.m

@@ -99,7 +99,11 @@ function [peakamps,peaklocs,peakwidths,resid] = ...
       ipeak = ipeak + 1;
       peaklocs(ipeak) = ploc;
       peakamps(ipeak) = pamp;
-      peakwidths(ipeak) = - 1/pcurv;
+      if (any(pcurv == 0))
+	peakwidths(ipeak) = 0;
+      else
+	peakwidths(ipeak) = - 1/pcurv;
+      end
       nrej = 0;
       if pdebug>1
 	peaksarr(plocq) = pamp;