Average Error: 12.4 → 0.4
Time: 9.2s
Precision: binary32
\[0 \leq s \land s \leq 256 \land 2.328306437 \cdot 10^{-10} \leq u \land u \leq 0.25\]
\[s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\]
↓
\[\begin{array}{l} \mathbf{if}\;4 \cdot u \leq 0.031028371304273605:\\ \;\;\;\;s \cdot \left(4 \cdot u + \left(u \cdot u\right) \cdot \left(u \cdot 21.333333333333332 + \left(8 + \left(u \cdot u\right) \cdot 64\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\\ \end{array}\]
s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)
↓
\begin{array}{l}
\mathbf{if}\;4 \cdot u \leq 0.031028371304273605:\\
\;\;\;\;s \cdot \left(4 \cdot u + \left(u \cdot u\right) \cdot \left(u \cdot 21.333333333333332 + \left(8 + \left(u \cdot u\right) \cdot 64\right)\right)\right)\\

\mathbf{else}:\\
\;\;\;\;s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\\

\end{array}
(FPCore (s u) :precision binary32 (* s (log (/ 1.0 (- 1.0 (* 4.0 u))))))
↓
(FPCore (s u)
 :precision binary32
 (if (<= (* 4.0 u) 0.031028371304273605)
   (*
    s
    (+
     (* 4.0 u)
     (* (* u u) (+ (* u 21.333333333333332) (+ 8.0 (* (* u u) 64.0))))))
   (* s (log (/ 1.0 (- 1.0 (* 4.0 u)))))))
float code(float s, float u) {
	return s * logf(1.0f / (1.0f - (4.0f * u)));
}
↓
float code(float s, float u) {
	float tmp;
	if ((4.0f * u) <= 0.031028371304273605f) {
		tmp = s * ((4.0f * u) + ((u * u) * ((u * 21.333333333333332f) + (8.0f + ((u * u) * 64.0f)))));
	} else {
		tmp = s * logf(1.0f / (1.0f - (4.0f * u)));
	}
	return tmp;
}

Error

Bits error versus s

Bits error versus u

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes

  2. if (*.f32 4 u) < 0.0310283713

    1. Initial program 14.5

      \[s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\]

    2. Simplified13.6

      \[\leadsto \color{blue}{s \cdot \left(-\log \left(1 - 4 \cdot u\right)\right)}\]

    3. Taylor expanded around 0 0.4

      \[\leadsto \color{blue}{21.333333333333332 \cdot \left(s \cdot {u}^{3}\right) + \left(4 \cdot \left(s \cdot u\right) + \left(8 \cdot \left(s \cdot {u}^{2}\right) + 64 \cdot \left(s \cdot {u}^{4}\right)\right)\right)}\]

    4. Simplified0.3

      \[\leadsto \color{blue}{s \cdot \left(u \cdot \left(4 + u \cdot 8\right) + \left(21.333333333333332 \cdot {u}^{3} + 64 \cdot {u}^{4}\right)\right)}\]
    5. Using strategy rm

    6. Applied distribute-rgt-in_binary320.2

      \[\leadsto s \cdot \left(\color{blue}{\left(4 \cdot u + \left(u \cdot 8\right) \cdot u\right)} + \left(21.333333333333332 \cdot {u}^{3} + 64 \cdot {u}^{4}\right)\right)\]

    7. Simplified0.2

      \[\leadsto s \cdot \left(\left(\color{blue}{u \cdot 4} + \left(u \cdot 8\right) \cdot u\right) + \left(21.333333333333332 \cdot {u}^{3} + 64 \cdot {u}^{4}\right)\right)\]

    8. Simplified0.2

      \[\leadsto s \cdot \left(\left(u \cdot 4 + \color{blue}{8 \cdot \left(u \cdot u\right)}\right) + \left(21.333333333333332 \cdot {u}^{3} + 64 \cdot {u}^{4}\right)\right)\]
    9. Using strategy rm

    10. Applied associate-+l+_binary320.2

      \[\leadsto s \cdot \color{blue}{\left(u \cdot 4 + \left(8 \cdot \left(u \cdot u\right) + \left(21.333333333333332 \cdot {u}^{3} + 64 \cdot {u}^{4}\right)\right)\right)}\]

    11. Simplified0.2

      \[\leadsto s \cdot \left(u \cdot 4 + \color{blue}{\left(u \cdot u\right) \cdot \left(21.333333333333332 \cdot u + \left(8 + 64 \cdot \left(u \cdot u\right)\right)\right)}\right)\]

    if 0.0310283713 < (*.f32 4 u)

    1. Initial program 1.5

      \[s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\]
  3. Recombined 2 regimes into one program.

  4. Final simplification0.4

    \[\leadsto \begin{array}{l} \mathbf{if}\;4 \cdot u \leq 0.031028371304273605:\\ \;\;\;\;s \cdot \left(4 \cdot u + \left(u \cdot u\right) \cdot \left(u \cdot 21.333333333333332 + \left(8 + \left(u \cdot u\right) \cdot 64\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2021175 
(FPCore (s u)
  :name "Disney BSSRDF, sample scattering profile, lower"
  :precision binary32
  :pre (and (<= 0.0 s 256.0) (<= 2.328306437e-10 u 0.25))
  (* s (log (/ 1.0 (- 1.0 (* 4.0 u))))))