Disney BSSRDF, PDF of scattering profile

Specification

?
\[\left(0 \leq s \land s \leq 256\right) \land \left(10^{-6} < r \land r < 1000000\right)\]
\[\begin{array}{l} \\ \frac{0.25 \cdot e^{\frac{-r}{s}}}{\left(\left(2 \cdot \pi\right) \cdot s\right) \cdot r} + \frac{0.75 \cdot e^{\frac{-r}{3 \cdot s}}}{\left(\left(6 \cdot \pi\right) \cdot s\right) \cdot r} \end{array} \]
(FPCore (s r)
 :precision binary32
 (+
  (/ (* 0.25 (exp (/ (- r) s))) (* (* (* 2.0 PI) s) r))
  (/ (* 0.75 (exp (/ (- r) (* 3.0 s)))) (* (* (* 6.0 PI) s) r))))
float code(float s, float r) {
	return ((0.25f * expf((-r / s))) / (((2.0f * ((float) M_PI)) * s) * r)) + ((0.75f * expf((-r / (3.0f * s)))) / (((6.0f * ((float) M_PI)) * s) * r));
}
function code(s, r)
	return Float32(Float32(Float32(Float32(0.25) * exp(Float32(Float32(-r) / s))) / Float32(Float32(Float32(Float32(2.0) * Float32(pi)) * s) * r)) + Float32(Float32(Float32(0.75) * exp(Float32(Float32(-r) / Float32(Float32(3.0) * s)))) / Float32(Float32(Float32(Float32(6.0) * Float32(pi)) * s) * r)))
end
function tmp = code(s, r)
	tmp = ((single(0.25) * exp((-r / s))) / (((single(2.0) * single(pi)) * s) * r)) + ((single(0.75) * exp((-r / (single(3.0) * s)))) / (((single(6.0) * single(pi)) * s) * r));
end
\begin{array}{l}

\\
\frac{0.25 \cdot e^{\frac{-r}{s}}}{\left(\left(2 \cdot \pi\right) \cdot s\right) \cdot r} + \frac{0.75 \cdot e^{\frac{-r}{3 \cdot s}}}{\left(\left(6 \cdot \pi\right) \cdot s\right) \cdot r}
\end{array}

Reproduce

?
herbie shell --seed 2024055 
(FPCore (s r)
  :name "Disney BSSRDF, PDF of scattering profile"
  :precision binary32
  :pre (and (and (<= 0.0 s) (<= s 256.0)) (and (< 1e-6 r) (< r 1000000.0)))
  (+ (/ (* 0.25 (exp (/ (- r) s))) (* (* (* 2.0 PI) s) r)) (/ (* 0.75 (exp (/ (- r) (* 3.0 s)))) (* (* (* 6.0 PI) s) r))))

Please file a bug report with this information.

Backtrace

vector-ref: index is out of range index: 34 valid range: [0, 33] vector: '#(133 113 133 113 93 133 113 133 113 93 73 133 173 153 133 113 173 173 153 133 113 93 133 193 173 153 133 113 173 153 133 113 93 73)LC
(unnamed)/home/nightlies/herbie/artem-minus-and-trig-tuning/src/compiler.rkt456
(unnamed)/home/nightlies/herbie/artem-minus-and-trig-tuning/src/ground-truth.rkt472
ival-eval/home/nightlies/herbie/artem-minus-and-trig-tuning/src/ground-truth.rkt690
batch-prepare-points/home/nightlies/herbie/artem-minus-and-trig-tuning/src/sampling.rkt1130
sample-points/home/nightlies/herbie/artem-minus-and-trig-tuning/src/sampling.rkt1610
setup-context!/home/nightlies/herbie/artem-minus-and-trig-tuning/src/mainloop.rkt3000
get-alternatives/report/home/nightlies/herbie/artem-minus-and-trig-tuning/src/sandbox.rkt1690
(unnamed)/home/nightlies/herbie/artem-minus-and-trig-tuning/src/sandbox.rkt2626
(unnamed)/usr/local/racket/share/pkgs/profile-lib/main.rkt4010
profile-thunk/usr/local/racket/share/pkgs/profile-lib/main.rkt90
(unnamed)/usr/local/racket/collects/racket/engine.rkt4224