Sample trimmed logistic on [-pi, pi]

Average Error: 0.3 → 0.3

Time: 18.3s

Precision: binary32

Cost: 16800

\[\left(2.328306437 \cdot 10^{-10} \leq u \land u \leq 1\right) \land \left(0 \leq s \land s \leq 1.0651631\right)\]

Math

\[\left(-s\right) \cdot \log \left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\pi}{s}}} - \frac{1}{1 + e^{\frac{\pi}{s}}}\right) + \frac{1}{1 + e^{\frac{\pi}{s}}}} - 1\right) \]

↓

\[s \cdot \left(-\log \left(-1 + \frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} + \frac{1 - u}{1 + e^{\frac{1}{\frac{s}{\pi}}}}}\right)\right) \]

FPCore

(FPCore (u s)
 :precision binary32
 (*
  (- s)
  (log
   (-
    (/
     1.0
     (+
      (*
       u
       (- (/ 1.0 (+ 1.0 (exp (/ (- PI) s)))) (/ 1.0 (+ 1.0 (exp (/ PI s))))))
      (/ 1.0 (+ 1.0 (exp (/ PI s))))))
    1.0))))

↓

(FPCore (u s)
 :precision binary32
 (*
  s
  (-
   (log
    (+
     -1.0
     (/
      1.0
      (+
       (/ u (+ 1.0 (exp (/ (- PI) s))))
       (/ (- 1.0 u) (+ 1.0 (exp (/ 1.0 (/ s PI))))))))))))

C

float code(float u, float s) {
	return -s * logf(((1.0f / ((u * ((1.0f / (1.0f + expf((-((float) M_PI) / s)))) - (1.0f / (1.0f + expf((((float) M_PI) / s)))))) + (1.0f / (1.0f + expf((((float) M_PI) / s)))))) - 1.0f));
}

↓

float code(float u, float s) {
	return s * -logf((-1.0f + (1.0f / ((u / (1.0f + expf((-((float) M_PI) / s)))) + ((1.0f - u) / (1.0f + expf((1.0f / (s / ((float) M_PI))))))))));
}

Julia

function code(u, s)
	return Float32(Float32(-s) * log(Float32(Float32(Float32(1.0) / Float32(Float32(u * Float32(Float32(Float32(1.0) / Float32(Float32(1.0) + exp(Float32(Float32(-Float32(pi)) / s)))) - Float32(Float32(1.0) / Float32(Float32(1.0) + exp(Float32(Float32(pi) / s)))))) + Float32(Float32(1.0) / Float32(Float32(1.0) + exp(Float32(Float32(pi) / s)))))) - Float32(1.0))))
end

↓

function code(u, s)
	return Float32(s * Float32(-log(Float32(Float32(-1.0) + Float32(Float32(1.0) / Float32(Float32(u / Float32(Float32(1.0) + exp(Float32(Float32(-Float32(pi)) / s)))) + Float32(Float32(Float32(1.0) - u) / Float32(Float32(1.0) + exp(Float32(Float32(1.0) / Float32(s / Float32(pi))))))))))))
end

MATLAB

function tmp = code(u, s)
	tmp = -s * log(((single(1.0) / ((u * ((single(1.0) / (single(1.0) + exp((-single(pi) / s)))) - (single(1.0) / (single(1.0) + exp((single(pi) / s)))))) + (single(1.0) / (single(1.0) + exp((single(pi) / s)))))) - single(1.0)));
end

↓

function tmp = code(u, s)
	tmp = s * -log((single(-1.0) + (single(1.0) / ((u / (single(1.0) + exp((-single(pi) / s)))) + ((single(1.0) - u) / (single(1.0) + exp((single(1.0) / (s / single(pi))))))))));
end

Derivation

1. Initial program 0.3

   \[\left(-s\right) \cdot \log \left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\pi}{s}}} - \frac{1}{1 + e^{\frac{\pi}{s}}}\right) + \frac{1}{1 + e^{\frac{\pi}{s}}}} - 1\right) \]

2. Simplified 0.3

   \[\leadsto \color{blue}{\left(-s\right) \cdot \log \left(\frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} - \frac{u + -1}{1 + e^{\frac{\pi}{s}}}} + -1\right)} \]
   Proof

   [Start] 0.3	\[ \left(-s\right) \cdot \log \left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\pi}{s}}} - \frac{1}{1 + e^{\frac{\pi}{s}}}\right) + \frac{1}{1 + e^{\frac{\pi}{s}}}} - 1\right) \]
   sub-neg [=>] 0.3	\[ \left(-s\right) \cdot \log \color{blue}{\left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\pi}{s}}} - \frac{1}{1 + e^{\frac{\pi}{s}}}\right) + \frac{1}{1 + e^{\frac{\pi}{s}}}} + \left(-1\right)\right)} \]

3. Applied egg-rr 0.3

   \[\leadsto \left(-s\right) \cdot \log \left(\frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} - \frac{u + -1}{1 + e^{\color{blue}{{\left(\frac{s}{\pi}\right)}^{-1}}}}} + -1\right) \]

4. Applied egg-rr 0.3

   \[\leadsto \left(-s\right) \cdot \log \left(\frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} - \frac{u + -1}{1 + e^{\color{blue}{\frac{1}{\frac{s}{\pi}}}}}} + -1\right) \]

5. Final simplification 0.3

   \[\leadsto s \cdot \left(-\log \left(-1 + \frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} + \frac{1 - u}{1 + e^{\frac{1}{\frac{s}{\pi}}}}}\right)\right) \]

Alternatives

Alternative 1
Error	0.3
Cost	16736

\[\left(-s\right) \cdot \log \left(-1 + \frac{1}{\frac{u}{1 + e^{\frac{-\pi}{s}}} + \frac{1 - u}{1 + e^{\frac{\pi}{s}}}}\right) \]

Alternative 2
Error	0.8
Cost	9984

\[s \cdot \left(-\log \left(-1 + \frac{1 + e^{\frac{-\pi}{s}}}{u}\right)\right) \]

Alternative 3
Error	7.7
Cost	9920

\[\left(-s\right) \cdot \log \left(\frac{1 + e^{\frac{-\pi}{s}}}{u}\right) \]

Alternative 4
Error	20.2
Cost	9888

\[s \cdot \left(-{\left(\sqrt[3]{\log \left(-1 + \frac{2}{u}\right)}\right)}^{3}\right) \]

Alternative 5
Error	20.2
Cost	9856

\[\left(-s\right) \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\log \left(-1 + \frac{2}{u}\right)\right)\right) \]

Alternative 6
Error	20.2
Cost	6816

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

Alternative 7
Error	20.2
Cost	3776

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

Alternative 8
Error	20.2
Cost	3712

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

Alternative 9
Error	20.2
Cost	3584

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

Alternative 10
Error	20.2
Cost	3456

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

Alternative 11
Error	20.2
Cost	3456

\[\left(-s\right) \cdot \mathsf{log1p}\left(\frac{2}{u} + -2\right) \]

Alternative 12
Error	20.2
Cost	3392

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

Alternative 13
Error	28.3
Cost	3232

\[-\pi \]

Reproduce

herbie shell --seed 2023073 
(FPCore (u s)
  :name "Sample trimmed logistic on [-pi, pi]"
  :precision binary32
  :pre (and (and (<= 2.328306437e-10 u) (<= u 1.0)) (and (<= 0.0 s) (<= s 1.0651631)))
  (* (- s) (log (- (/ 1.0 (+ (* u (- (/ 1.0 (+ 1.0 (exp (/ (- PI) s)))) (/ 1.0 (+ 1.0 (exp (/ PI s)))))) (/ 1.0 (+ 1.0 (exp (/ PI s)))))) 1.0))))