On 10/01/2017 16:52, Black Michael wrote:
try your loopback test using McLeod's estimator.
Just replace peakup.f90 with this:


subroutine peakup(ym,y0,yp,dx)
! McLeod's estimator
  delta = (ym-yp)/(2*y0+ym+yp)
dx=(sqrt(1+8*delta*delta)-1)/4*delta;
  return
end subroutine peakup

Hi Mike,

here you go:

16:58:29 660 1 1500 1500.000 0.000 -39.4 78.9 16:58:43 660 1 1501 1501.088 0.088 -39.5 78.116:58:57 660 1 1502 1501.953 -0.047 -30.2 68.216:59:16 880 1 1503 1503.039 0.039 -27.5 63.016:59:39 1210 1 1504 1503.902 -0.098 -39.5 77.916:59:50 1210 1 1505 1504.991 -0.009 -39.5 79.017:00:11 2500 1 1506 1506.078 0.078 -39.6 78.317:00:23 2500 1 1507 1506.944 -0.056 -35.8 74.817:00:41 3330 1 1508 1508.030 0.030 -39.4 78.817:00:53 3330 1 1509 1508.893 -0.107 -39.6 77.817:01:11 5000 1 1510 1509.983 -0.017 -39.5 79.0

Not an improvement, actually worse than the orginal non-windowed DFT plus and peakup function.

73
Bill
G4WJS.

------------------------------------------------------------------------------
Developer Access Program for Intel Xeon Phi Processors
Access to Intel Xeon Phi processor-based developer platforms.
With one year of Intel Parallel Studio XE.
Training and support from Colfax.
Order your platform today. http://sdm.link/xeonphi
_______________________________________________
wsjt-devel mailing list
wsjt-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/wsjt-devel

Reply via email to