diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-10-16 08:04:38 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-10-16 08:04:38 +0000 |
commit | 362c486e4731e0d6bf5185d2b2f55b73aa289462 (patch) | |
tree | c45096d6ca96c8c63d8464b8272fd521ea43dd6b /src | |
parent | 1743a158fb0fb7883e9cf6d92d9bf572650203e3 (diff) |
Small optimization to flow algorithm (to the adaptive stepsize). Should
be about 10% faster (at least on the testsuite). It breaks the testsuite
so I will commit a new one in a minute.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@128 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src')
-rw-r--r-- | src/AHFinder_flow.F | 7 |
1 files changed, 1 insertions, 6 deletions
diff --git a/src/AHFinder_flow.F b/src/AHFinder_flow.F index 0a862ed..84d853a 100644 --- a/src/AHFinder_flow.F +++ b/src/AHFinder_flow.F @@ -310,11 +310,10 @@ reldiff = dabs(ahfdiff - ahfdiff_old) . / dabs(ahfdiff + ahfdiff_old) - 100 new = .false. ! If the relative difference is too large, we reduce the -! stepsize and ignore this iteration. +! stepsize for the next iteration. if ((reldiff.gt.maxchange).or. . (ahfdiff.gt.maxchange*ahfsum)) then @@ -322,10 +321,6 @@ new = .true. A = 0.5D0*A - c0 = c0_old - cc = cc_old - cs = cs_old - ! If the diference is too small, we can safely increase ! the stepsize for the next iteration. |