Opened 13 years ago
Closed 13 years ago
#266 closed defect (fixed)
Corrections to ropp_abel_tool etc
Reported by: | Ian Culverwell | Owned by: | Ian Culverwell |
---|---|---|---|
Priority: | normal | Milestone: | 6.0 |
Component: | ropp_pp | Version: | 5.0 |
Keywords: | abel, invert | Cc: |
Description
ROPP6.0 beta testing uncovered some errors in the ropp_pp_{abel,invert}_{exp,lin}.f90 subroutines. Most of these appeared because the beta tester's compiler failed on SQRT(<0), whereas all ours produced NaNs which were processed OK.
These have been fixed, and the beta tester has run the code successfully on his machine.
As well as removing the SQRT(<0) bugs by generating ropp_MDFV instead of NaNs, other improvements have been made:
- Cleaning out of the test2.nc file, which contains the results of the _LIN routines, instead of continually appending to it;
- Consistent use of linear extrapolation to estimate the starting bangle and dlogn/dx in the _LIN routines;
- Correction to the upper bounds of some integrals (this led to many of the SQRT(<0)s);
- Correction to the expressions for the bangle and refrac (use of variables at lower bound: failure to do this can lead to discontinuities in bangle curve);
- General tidying up where needed.
Change history (4)
comment:1 by , 13 years ago
comment:4 by , 13 years ago
Resolution: | → fixed |
---|---|
Status: | new → closed |
Note:
See TracTickets
for help on using tickets.
Passes make test tests.