diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2021-08-04 15:18:08 -0400 |
| commit | 70fd14e98bc476a84aef352eaf063cfcba44ab83 (patch) | |
| tree | 2737f5c1c21c76824f49c6a1b9c7061952b646f5 /Invert.f | |
| parent | c396a4bc7a3282de9d2191de0e1076fec9df0783 (diff) | |
| download | moog-70fd14e98bc476a84aef352eaf063cfcba44ab83.tar.gz | |
Use modern equality operators
Diffstat (limited to 'Invert.f')
| -rwxr-xr-x | Invert.f | 10 |
1 files changed, 5 insertions, 5 deletions
@@ -29,18 +29,18 @@ c****************************************************************************** c here begins the inversion do n=1,order - if (n .eq. order) go to 244 + if (n == order) go to 244 g = a(n,n) n1 = n+1 ng = n do m=n1,order f = a(m,n) - if (f .gt. g) then + if (f > g) then g = f ng = m endif enddo - if (ng .eq. n) go to 244 + if (ng == n) go to 244 do k=1,order d = i(n,k) e = i(ng,k) @@ -52,7 +52,7 @@ c here begins the inversion a(n,k) = g enddo 244 g = a(n,n) - if (g .eq. 0.0) then + if (g == 0.0) then write (nf1out,1001) return endif @@ -60,7 +60,7 @@ c here begins the inversion a(n,k) = a(n,k)/g i(n,k) = i(n,k)/g enddo - if (n .eq. order) go to 27 + if (n == order) go to 27 do k=n1,order f = -a(k,n) do j=1,order |