weights_1d fails
Currently, the routine get_upper_source_index_desc_eq
fails, when the following arguments are passed:
- targetBound -12.0000000000000
- vals(cachedIndex+1_i8) -16.0000000000000
- cachedIndex 2000
- step -1.999999999995339E-003
Note that step is calculated and should really be -2.0E-3.
it calculates
- advance 2000.00000000466
- floored 2000
and thus
- abs(real(floored, kind=dp) - advance) 4.664798325393349E-009
- maxTolerance 1.000000000000000E-009
and returns 0
, but should 1
.