w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2024-11-22 04:11:16 -05:00
Code Issues Packages Projects Releases Wiki Activity

Fix the loop exit criterion in polfit.f90.

Browse Source
This commit is contained in:
Joe Taylor 2021-07-01 11:14:46 -04:00
parent e0c06168cf
commit 525c791a16
2 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
map65/libm65/polfit.f90
View File

@ -56,15 +56,16 @@ subroutine polfit(y,npts,a)
deltaa(j)=deltaa(j)*fn/3.
! write(*,4000) iter,j,a,deltaa,chisq2
!4000 format(2i2,2(2x,3f8.2),f12.5)
enddo
enddo ! j=1,nterms
chisqr=fchisq_pol(y,npts,a)
! write(*,4000) 0,0,a,chisqr
if(chisqr.lt.1.0) exit
if(deltaa(1).lt.0.01*(a(2)-a(1)) .and. deltaa(2).lt.0.01*(a(2)-a(1)) &
.and. deltaa(3).lt.1.0) exit
if(chisqr/chisqr0.gt.0.99) exit
a(3)=mod(a(3)+360.0,180.0)
chisqr0=chisqr
enddo
enddo ! iter
a(3)=mod(a(3)+360.0,180.0)
return
end subroutine polfit