Beckmann Sample, near normal, slope_y

Percentage Accurate: 57.2% → 98.4%

Time: 11.4s

Alternatives: 10

Speedup: 2.9×

Specification

\[\left(\left(cosTheta\_i > 0.9999 \land cosTheta\_i \leq 1\right) \land \left(2.328306437 \cdot 10^{-10} \leq u1 \land u1 \leq 1\right)\right) \land \left(2.328306437 \cdot 10^{-10} \leq u2 \land u2 \leq 1\right)\]

Math

\[\begin{array}{l} \\ \sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* (sqrt (- (log (- 1.0 u1)))) (sin (* (* 2.0 PI) u2))))

C

float code(float cosTheta_i, float u1, float u2) {
	return sqrtf(-logf((1.0f - u1))) * sinf(((2.0f * ((float) M_PI)) * u2));
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * sin(Float32(Float32(Float32(2.0) * Float32(pi)) * u2)))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = sqrt(-log((single(1.0) - u1))) * sin(((single(2.0) * single(pi)) * u2));
end

Initial Program: 57.2% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* (sqrt (- (log (- 1.0 u1)))) (sin (* (* 2.0 PI) u2))))

C

float code(float cosTheta_i, float u1, float u2) {
	return sqrtf(-logf((1.0f - u1))) * sinf(((2.0f * ((float) M_PI)) * u2));
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(sqrt(Float32(-log(Float32(Float32(1.0) - u1)))) * sin(Float32(Float32(Float32(2.0) * Float32(pi)) * u2)))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = sqrt(-log((single(1.0) - u1))) * sin(((single(2.0) * single(pi)) * u2));
end

Alternative 1: 98.4% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \sqrt[3]{{\sin \left(2 \cdot \left(\pi \cdot u2\right)\right)}^{3} \cdot {\left(-\mathsf{log1p}\left(-u1\right)\right)}^{1.5}} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (cbrt (* (pow (sin (* 2.0 (* PI u2))) 3.0) (pow (- (log1p (- u1))) 1.5))))

C

float code(float cosTheta_i, float u1, float u2) {
	return cbrtf((powf(sinf((2.0f * (((float) M_PI) * u2))), 3.0f) * powf(-log1pf(-u1), 1.5f)));
}

Julia

function code(cosTheta_i, u1, u2)
	return cbrt(Float32((sin(Float32(Float32(2.0) * Float32(Float32(pi) * u2))) ^ Float32(3.0)) * (Float32(-log1p(Float32(-u1))) ^ Float32(1.5))))
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Step-by-step derivation

   1. sub-neg 57.4%

      \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   2. log1p-define 98.5%

      \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

3. Simplified 98.5%

   \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
4. Add Preprocessing
5. Step-by-step derivation

   1. add-cbrt-cube 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)}} \cdot u2\right) \]

   2. add-cbrt-cube 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)} \cdot \color{blue}{\sqrt[3]{\left(u2 \cdot u2\right) \cdot u2}}\right) \]

   3. cbrt-unprod 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{\left(\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right)} \]

   4. pow3 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(2 \cdot \pi\right)}^{3}} \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right) \]

   5. pow3 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot \color{blue}{{u2}^{3}}}\right) \]

6. Applied egg-rr 98.5%

   \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot {u2}^{3}}\right)} \]
7. Step-by-step derivation

   1. *-commutative 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{u2}^{3} \cdot {\left(2 \cdot \pi\right)}^{3}}}\right) \]

   2. pow-prod-down 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(u2 \cdot \left(2 \cdot \pi\right)\right)}^{3}}}\right) \]

   3. associate-*l* 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\color{blue}{\left(\left(u2 \cdot 2\right) \cdot \pi\right)}}^{3}}\right) \]

   4. rem-cbrt-cube 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\left(u2 \cdot 2\right) \cdot \pi\right)} \]

   5. *-commutative 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\left(2 \cdot u2\right)} \cdot \pi\right) \]

   6. associate-*r* 98.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]

   7. sin-2 98.4%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(2 \cdot \left(\sin \left(u2 \cdot \pi\right) \cdot \cos \left(u2 \cdot \pi\right)\right)\right)} \]

   8. associate-*r* 98.4%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \sin \left(u2 \cdot \pi\right)\right) \cdot \cos \left(u2 \cdot \pi\right)\right)} \]

   9. *-commutative 98.4%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(\left(2 \cdot \sin \color{blue}{\left(\pi \cdot u2\right)}\right) \cdot \cos \left(u2 \cdot \pi\right)\right) \]

   10. *-commutative 98.4%

       \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \color{blue}{\left(\pi \cdot u2\right)}\right) \]

8. Applied egg-rr 98.4%

   \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)} \]
9. Step-by-step derivation

   1. *-commutative 98.4%

      \[\leadsto \color{blue}{\left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right) \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)}} \]

   2. add-cbrt-cube 98.3%

      \[\leadsto \color{blue}{\sqrt[3]{\left(\left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)}} \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)} \]

   3. add-cbrt-cube 98.3%

      \[\leadsto \sqrt[3]{\left(\left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)} \cdot \color{blue}{\sqrt[3]{\left(\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)}\right) \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)}}} \]

   4. cbrt-unprod 98.3%

      \[\leadsto \color{blue}{\sqrt[3]{\left(\left(\left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)\right) \cdot \left(\left(2 \cdot \sin \left(\pi \cdot u2\right)\right) \cdot \cos \left(\pi \cdot u2\right)\right)\right) \cdot \left(\left(\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)}\right) \cdot \sqrt{-\mathsf{log1p}\left(-u1\right)}\right)}} \]

10. Applied egg-rr 98.6%

    \[\leadsto \color{blue}{\sqrt[3]{{\sin \left(2 \cdot \left(\pi \cdot u2\right)\right)}^{3} \cdot {\left(-\mathsf{log1p}\left(-u1\right)\right)}^{1.5}}} \]
11. Add Preprocessing

Alternative 2: 98.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* (sqrt (- (log1p (- u1)))) (sin (* u2 (* 2.0 PI)))))

C

float code(float cosTheta_i, float u1, float u2) {
	return sqrtf(-log1pf(-u1)) * sinf((u2 * (2.0f * ((float) M_PI))));
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(sqrt(Float32(-log1p(Float32(-u1)))) * sin(Float32(u2 * Float32(Float32(2.0) * Float32(pi)))))
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Step-by-step derivation

   1. sub-neg 57.4%

      \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   2. log1p-define 98.5%

      \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

3. Simplified 98.5%

   \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
4. Add Preprocessing

5. Final simplification 98.5%

   \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \]
6. Add Preprocessing

Alternative 3: 96.3% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := u2 \cdot \left(2 \cdot \pi\right)\\ \mathbf{if}\;t\_0 \leq 0.0031500000040978193:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot t\_0\\ \mathbf{else}:\\ \;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 - u1 \cdot -0.25\right)\right)\right)}\\ \end{array} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (let* ((t_0 (* u2 (* 2.0 PI))))
   (if (<= t_0 0.0031500000040978193)
     (* (sqrt (- (log1p (- u1)))) t_0)
     (*
      (sin t_0)
      (sqrt
       (*
        u1
        (+ 1.0 (* u1 (+ 0.5 (* u1 (- 0.3333333333333333 (* u1 -0.25))))))))))))

C

float code(float cosTheta_i, float u1, float u2) {
	float t_0 = u2 * (2.0f * ((float) M_PI));
	float tmp;
	if (t_0 <= 0.0031500000040978193f) {
		tmp = sqrtf(-log1pf(-u1)) * t_0;
	} else {
		tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * (0.5f + (u1 * (0.3333333333333333f - (u1 * -0.25f))))))));
	}
	return tmp;
}

Julia

function code(cosTheta_i, u1, u2)
	t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi)))
	tmp = Float32(0.0)
	if (t_0 <= Float32(0.0031500000040978193))
		tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * t_0);
	else
		tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) + Float32(u1 * Float32(Float32(0.3333333333333333) - Float32(u1 * Float32(-0.25))))))))));
	end
	return tmp
end

Derivation

1. Split input into 2 regimes

2. if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.00315

   1. Initial program 59.1%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Step-by-step derivation

      1. sub-neg 59.1%

         \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

      2. log1p-define 98.6%

         \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   3. Simplified 98.6%

      \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
   4. Add Preprocessing
   5. Step-by-step derivation

      1. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)}} \cdot u2\right) \]

      2. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)} \cdot \color{blue}{\sqrt[3]{\left(u2 \cdot u2\right) \cdot u2}}\right) \]

      3. cbrt-unprod 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{\left(\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right)} \]

      4. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(2 \cdot \pi\right)}^{3}} \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right) \]

      5. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot \color{blue}{{u2}^{3}}}\right) \]

   6. Applied egg-rr 98.7%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot {u2}^{3}}\right)} \]

   7. Taylor expanded in u2 around 0 98.1%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]
   8. Step-by-step derivation

      1. *-commutative 98.1%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

      2. associate-*r* 98.1%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   9. Simplified 98.1%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   if 0.00315 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)

   1. Initial program 53.8%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Add Preprocessing

   3. Taylor expanded in u1 around 0 95.6%

      \[\leadsto \sqrt{-\color{blue}{u1 \cdot \left(u1 \cdot \left(u1 \cdot \left(-0.25 \cdot u1 - 0.3333333333333333\right) - 0.5\right) - 1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 97.3%

   \[\leadsto \begin{array}{l} \mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.0031500000040978193:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 + u1 \cdot \left(0.3333333333333333 - u1 \cdot -0.25\right)\right)\right)}\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 95.7% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := u2 \cdot \left(2 \cdot \pi\right)\\ \mathbf{if}\;t\_0 \leq 0.0031500000040978193:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot t\_0\\ \mathbf{else}:\\ \;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 - u1 \cdot -0.3333333333333333\right)\right)}\\ \end{array} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (let* ((t_0 (* u2 (* 2.0 PI))))
   (if (<= t_0 0.0031500000040978193)
     (* (sqrt (- (log1p (- u1)))) t_0)
     (*
      (sin t_0)
      (sqrt (* u1 (+ 1.0 (* u1 (- 0.5 (* u1 -0.3333333333333333))))))))))

C

float code(float cosTheta_i, float u1, float u2) {
	float t_0 = u2 * (2.0f * ((float) M_PI));
	float tmp;
	if (t_0 <= 0.0031500000040978193f) {
		tmp = sqrtf(-log1pf(-u1)) * t_0;
	} else {
		tmp = sinf(t_0) * sqrtf((u1 * (1.0f + (u1 * (0.5f - (u1 * -0.3333333333333333f))))));
	}
	return tmp;
}

Julia

function code(cosTheta_i, u1, u2)
	t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi)))
	tmp = Float32(0.0)
	if (t_0 <= Float32(0.0031500000040978193))
		tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * t_0);
	else
		tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) + Float32(u1 * Float32(Float32(0.5) - Float32(u1 * Float32(-0.3333333333333333))))))));
	end
	return tmp
end

Derivation

1. Split input into 2 regimes

2. if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.00315

   1. Initial program 59.1%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Step-by-step derivation

      1. sub-neg 59.1%

         \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

      2. log1p-define 98.6%

         \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   3. Simplified 98.6%

      \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
   4. Add Preprocessing
   5. Step-by-step derivation

      1. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)}} \cdot u2\right) \]

      2. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)} \cdot \color{blue}{\sqrt[3]{\left(u2 \cdot u2\right) \cdot u2}}\right) \]

      3. cbrt-unprod 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{\left(\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right)} \]

      4. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(2 \cdot \pi\right)}^{3}} \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right) \]

      5. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot \color{blue}{{u2}^{3}}}\right) \]

   6. Applied egg-rr 98.7%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot {u2}^{3}}\right)} \]

   7. Taylor expanded in u2 around 0 98.1%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]
   8. Step-by-step derivation

      1. *-commutative 98.1%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

      2. associate-*r* 98.1%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   9. Simplified 98.1%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   if 0.00315 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)

   1. Initial program 53.8%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Add Preprocessing

   3. Taylor expanded in u1 around 0 93.2%

      \[\leadsto \sqrt{-\color{blue}{u1 \cdot \left(u1 \cdot \left(-0.3333333333333333 \cdot u1 - 0.5\right) - 1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 96.5%

   \[\leadsto \begin{array}{l} \mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.0031500000040978193:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1 \cdot \left(1 + u1 \cdot \left(0.5 - u1 \cdot -0.3333333333333333\right)\right)}\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 94.3% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := u2 \cdot \left(2 \cdot \pi\right)\\ \mathbf{if}\;t\_0 \leq 0.0044999998062849045:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot t\_0\\ \mathbf{else}:\\ \;\;\;\;\sin t\_0 \cdot \sqrt{u1 \cdot \left(1 - u1 \cdot -0.5\right)}\\ \end{array} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (let* ((t_0 (* u2 (* 2.0 PI))))
   (if (<= t_0 0.0044999998062849045)
     (* (sqrt (- (log1p (- u1)))) t_0)
     (* (sin t_0) (sqrt (* u1 (- 1.0 (* u1 -0.5))))))))

C

float code(float cosTheta_i, float u1, float u2) {
	float t_0 = u2 * (2.0f * ((float) M_PI));
	float tmp;
	if (t_0 <= 0.0044999998062849045f) {
		tmp = sqrtf(-log1pf(-u1)) * t_0;
	} else {
		tmp = sinf(t_0) * sqrtf((u1 * (1.0f - (u1 * -0.5f))));
	}
	return tmp;
}

Julia

function code(cosTheta_i, u1, u2)
	t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi)))
	tmp = Float32(0.0)
	if (t_0 <= Float32(0.0044999998062849045))
		tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * t_0);
	else
		tmp = Float32(sin(t_0) * sqrt(Float32(u1 * Float32(Float32(1.0) - Float32(u1 * Float32(-0.5))))));
	end
	return tmp
end

Derivation

1. Split input into 2 regimes

2. if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.00449999981

   1. Initial program 59.3%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Step-by-step derivation

      1. sub-neg 59.3%

         \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

      2. log1p-define 98.6%

         \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   3. Simplified 98.6%

      \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
   4. Add Preprocessing
   5. Step-by-step derivation

      1. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)}} \cdot u2\right) \]

      2. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)} \cdot \color{blue}{\sqrt[3]{\left(u2 \cdot u2\right) \cdot u2}}\right) \]

      3. cbrt-unprod 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{\left(\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right)} \]

      4. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(2 \cdot \pi\right)}^{3}} \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right) \]

      5. pow3 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot \color{blue}{{u2}^{3}}}\right) \]

   6. Applied egg-rr 98.6%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot {u2}^{3}}\right)} \]

   7. Taylor expanded in u2 around 0 97.6%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]
   8. Step-by-step derivation

      1. *-commutative 97.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

      2. associate-*r* 97.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   9. Simplified 97.6%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   if 0.00449999981 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)

   1. Initial program 52.9%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Add Preprocessing

   3. Taylor expanded in u1 around 0 89.9%

      \[\leadsto \sqrt{-\color{blue}{u1 \cdot \left(-0.5 \cdot u1 - 1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 95.3%

   \[\leadsto \begin{array}{l} \mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.0044999998062849045:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1 \cdot \left(1 - u1 \cdot -0.5\right)}\\ \end{array} \]
5. Add Preprocessing

Alternative 6: 90.5% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := u2 \cdot \left(2 \cdot \pi\right)\\ \mathbf{if}\;t\_0 \leq 0.004600000102072954:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot t\_0\\ \mathbf{else}:\\ \;\;\;\;\sin t\_0 \cdot \sqrt{u1}\\ \end{array} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (let* ((t_0 (* u2 (* 2.0 PI))))
   (if (<= t_0 0.004600000102072954)
     (* (sqrt (- (log1p (- u1)))) t_0)
     (* (sin t_0) (sqrt u1)))))

C

float code(float cosTheta_i, float u1, float u2) {
	float t_0 = u2 * (2.0f * ((float) M_PI));
	float tmp;
	if (t_0 <= 0.004600000102072954f) {
		tmp = sqrtf(-log1pf(-u1)) * t_0;
	} else {
		tmp = sinf(t_0) * sqrtf(u1);
	}
	return tmp;
}

Julia

function code(cosTheta_i, u1, u2)
	t_0 = Float32(u2 * Float32(Float32(2.0) * Float32(pi)))
	tmp = Float32(0.0)
	if (t_0 <= Float32(0.004600000102072954))
		tmp = Float32(sqrt(Float32(-log1p(Float32(-u1)))) * t_0);
	else
		tmp = Float32(sin(t_0) * sqrt(u1));
	end
	return tmp
end

Derivation

1. Split input into 2 regimes

2. if (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2) < 0.0046000001

   1. Initial program 59.5%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Step-by-step derivation

      1. sub-neg 59.5%

         \[\leadsto \sqrt{-\log \color{blue}{\left(1 + \left(-u1\right)\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

      2. log1p-define 98.6%

         \[\leadsto \sqrt{-\color{blue}{\mathsf{log1p}\left(-u1\right)}} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]

   3. Simplified 98.6%

      \[\leadsto \color{blue}{\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \]
   4. Add Preprocessing
   5. Step-by-step derivation

      1. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\color{blue}{\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)}} \cdot u2\right) \]

      2. add-cbrt-cube 98.6%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)} \cdot \color{blue}{\sqrt[3]{\left(u2 \cdot u2\right) \cdot u2}}\right) \]

      3. cbrt-unprod 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{\left(\left(\left(2 \cdot \pi\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(2 \cdot \pi\right)\right) \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right)} \]

      4. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{\color{blue}{{\left(2 \cdot \pi\right)}^{3}} \cdot \left(\left(u2 \cdot u2\right) \cdot u2\right)}\right) \]

      5. pow3 98.7%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot \color{blue}{{u2}^{3}}}\right) \]

   6. Applied egg-rr 98.7%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \sin \color{blue}{\left(\sqrt[3]{{\left(2 \cdot \pi\right)}^{3} \cdot {u2}^{3}}\right)} \]

   7. Taylor expanded in u2 around 0 97.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]
   8. Step-by-step derivation

      1. *-commutative 97.5%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

      2. associate-*r* 97.5%

         \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   9. Simplified 97.5%

      \[\leadsto \sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)} \]

   if 0.0046000001 < (*.f32 (*.f32 #s(literal 2 binary32) (PI.f32)) u2)

   1. Initial program 52.3%

      \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
   2. Add Preprocessing

   3. Taylor expanded in u1 around 0 -0.0%

      \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}\right)} \]
   4. Step-by-step derivation

      1. associate-*r* -0.0%

         \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}} \]

      2. unpow2 -0.0%

         \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{\left(\sqrt{-1} \cdot \sqrt{-1}\right)} \]

      3. rem-square-sqrt 3.6%

         \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{-1} \]

      4. *-commutative 3.6%

         \[\leadsto \color{blue}{\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1}\right)} \cdot -1 \]

      5. associate-*l* 3.6%

         \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(\sqrt{u1} \cdot -1\right)} \]

      6. *-commutative 3.6%

         \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-1 \cdot \sqrt{u1}\right)} \]

      7. mul-1-neg 3.6%

         \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-\sqrt{u1}\right)} \]

   5. Simplified 3.6%

      \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)} \]
   6. Step-by-step derivation

      1. expm1-log1p-u 3.6%

         \[\leadsto \color{blue}{\mathsf{expm1}\left(\mathsf{log1p}\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)\right)\right)} \]

      2. expm1-undefine 3.6%

         \[\leadsto \color{blue}{e^{\mathsf{log1p}\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)\right)} - 1} \]

      3. *-commutative 3.6%

         \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\left(-\sqrt{u1}\right) \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)}\right)} - 1 \]

      4. add-sqr-sqrt -0.0%

         \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\left(\sqrt{-\sqrt{u1}} \cdot \sqrt{-\sqrt{u1}}\right)} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

      5. sqrt-unprod 57.8%

         \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\sqrt{\left(-\sqrt{u1}\right) \cdot \left(-\sqrt{u1}\right)}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

      6. sqr-neg 57.8%

         \[\leadsto e^{\mathsf{log1p}\left(\sqrt{\color{blue}{\sqrt{u1} \cdot \sqrt{u1}}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

      7. add-sqr-sqrt 57.8%

         \[\leadsto e^{\mathsf{log1p}\left(\sqrt{\color{blue}{u1}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

      8. *-commutative 57.8%

         \[\leadsto e^{\mathsf{log1p}\left(\sqrt{u1} \cdot \sin \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right)\right)} - 1 \]

   7. Applied egg-rr 57.8%

      \[\leadsto \color{blue}{e^{\mathsf{log1p}\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)} - 1} \]
   8. Step-by-step derivation

      1. log1p-undefine 57.8%

         \[\leadsto e^{\color{blue}{\log \left(1 + \sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)}} - 1 \]

      2. rem-exp-log 57.8%

         \[\leadsto \color{blue}{\left(1 + \sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)} - 1 \]

      3. +-commutative 57.8%

         \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right) + 1\right)} - 1 \]

      4. associate--l+ 79.6%

         \[\leadsto \color{blue}{\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right) + \left(1 - 1\right)} \]

      5. *-commutative 79.6%

         \[\leadsto \color{blue}{\sin \left(2 \cdot \left(\pi \cdot u2\right)\right) \cdot \sqrt{u1}} + \left(1 - 1\right) \]

      6. *-commutative 79.6%

         \[\leadsto \sin \left(2 \cdot \color{blue}{\left(u2 \cdot \pi\right)}\right) \cdot \sqrt{u1} + \left(1 - 1\right) \]

      7. metadata-eval 79.6%

         \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1} + \color{blue}{0} \]

      8. mul0-lft 79.6%

         \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1} + \color{blue}{0 \cdot \sqrt{u1}} \]

      9. distribute-rgt-in 79.6%

         \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) + 0\right)} \]

      10. +-rgt-identity 79.6%

          \[\leadsto \sqrt{u1} \cdot \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]

      11. *-commutative 79.6%

          \[\leadsto \sqrt{u1} \cdot \sin \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

      12. count-2 79.6%

          \[\leadsto \sqrt{u1} \cdot \sin \color{blue}{\left(\pi \cdot u2 + \pi \cdot u2\right)} \]

      13. distribute-rgt-out 79.6%

          \[\leadsto \sqrt{u1} \cdot \sin \color{blue}{\left(u2 \cdot \left(\pi + \pi\right)\right)} \]

      14. count-2 79.6%

          \[\leadsto \sqrt{u1} \cdot \sin \left(u2 \cdot \color{blue}{\left(2 \cdot \pi\right)}\right) \]

   9. Simplified 79.6%

      \[\leadsto \color{blue}{\sqrt{u1} \cdot \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 92.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;u2 \cdot \left(2 \cdot \pi\right) \leq 0.004600000102072954:\\ \;\;\;\;\sqrt{-\mathsf{log1p}\left(-u1\right)} \cdot \left(u2 \cdot \left(2 \cdot \pi\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1}\\ \end{array} \]
5. Add Preprocessing

Alternative 7: 76.9% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1} \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* (sin (* u2 (* 2.0 PI))) (sqrt u1)))

C

float code(float cosTheta_i, float u1, float u2) {
	return sinf((u2 * (2.0f * ((float) M_PI)))) * sqrtf(u1);
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(sin(Float32(u2 * Float32(Float32(2.0) * Float32(pi)))) * sqrt(u1))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = sin((u2 * (single(2.0) * single(pi)))) * sqrt(u1);
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Add Preprocessing

3. Taylor expanded in u1 around 0 -0.0%

   \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}\right)} \]
4. Step-by-step derivation

   1. associate-*r* -0.0%

      \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}} \]

   2. unpow2 -0.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{\left(\sqrt{-1} \cdot \sqrt{-1}\right)} \]

   3. rem-square-sqrt 4.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{-1} \]

   4. *-commutative 4.0%

      \[\leadsto \color{blue}{\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1}\right)} \cdot -1 \]

   5. associate-*l* 4.0%

      \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(\sqrt{u1} \cdot -1\right)} \]

   6. *-commutative 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-1 \cdot \sqrt{u1}\right)} \]

   7. mul-1-neg 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-\sqrt{u1}\right)} \]

5. Simplified 4.0%

   \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)} \]
6. Step-by-step derivation

   1. expm1-log1p-u 4.0%

      \[\leadsto \color{blue}{\mathsf{expm1}\left(\mathsf{log1p}\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)\right)\right)} \]

   2. expm1-undefine 5.3%

      \[\leadsto \color{blue}{e^{\mathsf{log1p}\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)\right)} - 1} \]

   3. *-commutative 5.3%

      \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\left(-\sqrt{u1}\right) \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)}\right)} - 1 \]

   4. add-sqr-sqrt -0.0%

      \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\left(\sqrt{-\sqrt{u1}} \cdot \sqrt{-\sqrt{u1}}\right)} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

   5. sqrt-unprod 30.5%

      \[\leadsto e^{\mathsf{log1p}\left(\color{blue}{\sqrt{\left(-\sqrt{u1}\right) \cdot \left(-\sqrt{u1}\right)}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

   6. sqr-neg 30.5%

      \[\leadsto e^{\mathsf{log1p}\left(\sqrt{\color{blue}{\sqrt{u1} \cdot \sqrt{u1}}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

   7. add-sqr-sqrt 30.5%

      \[\leadsto e^{\mathsf{log1p}\left(\sqrt{\color{blue}{u1}} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right)} - 1 \]

   8. *-commutative 30.5%

      \[\leadsto e^{\mathsf{log1p}\left(\sqrt{u1} \cdot \sin \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right)\right)} - 1 \]

7. Applied egg-rr 30.5%

   \[\leadsto \color{blue}{e^{\mathsf{log1p}\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)} - 1} \]
8. Step-by-step derivation

   1. log1p-undefine 30.5%

      \[\leadsto e^{\color{blue}{\log \left(1 + \sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)}} - 1 \]

   2. rem-exp-log 30.5%

      \[\leadsto \color{blue}{\left(1 + \sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right)\right)} - 1 \]

   3. +-commutative 30.5%

      \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right) + 1\right)} - 1 \]

   4. associate--l+ 76.8%

      \[\leadsto \color{blue}{\sqrt{u1} \cdot \sin \left(2 \cdot \left(\pi \cdot u2\right)\right) + \left(1 - 1\right)} \]

   5. *-commutative 76.8%

      \[\leadsto \color{blue}{\sin \left(2 \cdot \left(\pi \cdot u2\right)\right) \cdot \sqrt{u1}} + \left(1 - 1\right) \]

   6. *-commutative 76.8%

      \[\leadsto \sin \left(2 \cdot \color{blue}{\left(u2 \cdot \pi\right)}\right) \cdot \sqrt{u1} + \left(1 - 1\right) \]

   7. metadata-eval 76.8%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1} + \color{blue}{0} \]

   8. mul0-lft 76.8%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1} + \color{blue}{0 \cdot \sqrt{u1}} \]

   9. distribute-rgt-in 76.8%

      \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) + 0\right)} \]

   10. +-rgt-identity 76.8%

       \[\leadsto \sqrt{u1} \cdot \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)} \]

   11. *-commutative 76.8%

       \[\leadsto \sqrt{u1} \cdot \sin \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right) \]

   12. count-2 76.8%

       \[\leadsto \sqrt{u1} \cdot \sin \color{blue}{\left(\pi \cdot u2 + \pi \cdot u2\right)} \]

   13. distribute-rgt-out 76.8%

       \[\leadsto \sqrt{u1} \cdot \sin \color{blue}{\left(u2 \cdot \left(\pi + \pi\right)\right)} \]

   14. count-2 76.8%

       \[\leadsto \sqrt{u1} \cdot \sin \left(u2 \cdot \color{blue}{\left(2 \cdot \pi\right)}\right) \]

9. Simplified 76.8%

   \[\leadsto \color{blue}{\sqrt{u1} \cdot \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right)} \]

10. Final simplification 76.8%

    \[\leadsto \sin \left(u2 \cdot \left(2 \cdot \pi\right)\right) \cdot \sqrt{u1} \]
11. Add Preprocessing

Alternative 8: 66.6% accurate, 2.9× speedup

Math

\[\begin{array}{l} \\ 2 \cdot \left(\sqrt{u1} \cdot \left(\pi \cdot u2\right)\right) \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2) :precision binary32 (* 2.0 (* (sqrt u1) (* PI u2))))

C

float code(float cosTheta_i, float u1, float u2) {
	return 2.0f * (sqrtf(u1) * (((float) M_PI) * u2));
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(Float32(2.0) * Float32(sqrt(u1) * Float32(Float32(pi) * u2)))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = single(2.0) * (sqrt(u1) * (single(pi) * u2));
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Add Preprocessing
3. Step-by-step derivation

   1. add-sqr-sqrt 56.9%

      \[\leadsto \color{blue}{\sqrt{\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)} \cdot \sqrt{\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)}} \]

   2. sqrt-unprod 56.9%

      \[\leadsto \color{blue}{\sqrt{\left(\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)\right) \cdot \left(\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)\right)}} \]

   3. swap-sqr 56.9%

      \[\leadsto \sqrt{\color{blue}{\left(\sqrt{-\log \left(1 - u1\right)} \cdot \sqrt{-\log \left(1 - u1\right)}\right) \cdot \left(\sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right)\right)}} \]

4. Applied egg-rr 73.4%

   \[\leadsto \color{blue}{\sqrt{\mathsf{log1p}\left(u1\right) \cdot {\sin \left(\left(u2 \cdot 2\right) \cdot \pi\right)}^{2}}} \]

5. Taylor expanded in u1 around 0 75.5%

   \[\leadsto \sqrt{\color{blue}{u1 \cdot {\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)}^{2}}} \]
6. Step-by-step derivation

   1. *-commutative 75.5%

      \[\leadsto \sqrt{u1 \cdot {\sin \left(2 \cdot \color{blue}{\left(\pi \cdot u2\right)}\right)}^{2}} \]

   2. associate-*r* 75.5%

      \[\leadsto \sqrt{u1 \cdot {\sin \color{blue}{\left(\left(2 \cdot \pi\right) \cdot u2\right)}}^{2}} \]

7. Simplified 75.5%

   \[\leadsto \sqrt{\color{blue}{u1 \cdot {\sin \left(\left(2 \cdot \pi\right) \cdot u2\right)}^{2}}} \]

8. Taylor expanded in u2 around 0 67.6%

   \[\leadsto \color{blue}{2 \cdot \left(\sqrt{u1} \cdot \left(u2 \cdot \pi\right)\right)} \]

9. Final simplification 67.6%

   \[\leadsto 2 \cdot \left(\sqrt{u1} \cdot \left(\pi \cdot u2\right)\right) \]
10. Add Preprocessing

Alternative 9: 4.6% accurate, 2.9× speedup

Math

\[\begin{array}{l} \\ -2 \cdot \left(\pi \cdot \left(u2 \cdot \sqrt{u1}\right)\right) \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* -2.0 (* PI (* u2 (sqrt u1)))))

C

float code(float cosTheta_i, float u1, float u2) {
	return -2.0f * (((float) M_PI) * (u2 * sqrtf(u1)));
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(Float32(-2.0) * Float32(Float32(pi) * Float32(u2 * sqrt(u1))))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = single(-2.0) * (single(pi) * (u2 * sqrt(u1)));
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Add Preprocessing

3. Taylor expanded in u1 around 0 -0.0%

   \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}\right)} \]
4. Step-by-step derivation

   1. associate-*r* -0.0%

      \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}} \]

   2. unpow2 -0.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{\left(\sqrt{-1} \cdot \sqrt{-1}\right)} \]

   3. rem-square-sqrt 4.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{-1} \]

   4. *-commutative 4.0%

      \[\leadsto \color{blue}{\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1}\right)} \cdot -1 \]

   5. associate-*l* 4.0%

      \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(\sqrt{u1} \cdot -1\right)} \]

   6. *-commutative 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-1 \cdot \sqrt{u1}\right)} \]

   7. mul-1-neg 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-\sqrt{u1}\right)} \]

5. Simplified 4.0%

   \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)} \]

6. Taylor expanded in u2 around 0 4.3%

   \[\leadsto \color{blue}{-2 \cdot \left(\sqrt{u1} \cdot \left(u2 \cdot \pi\right)\right)} \]
7. Step-by-step derivation

   1. associate-*r* 4.3%

      \[\leadsto -2 \cdot \color{blue}{\left(\left(\sqrt{u1} \cdot u2\right) \cdot \pi\right)} \]

8. Simplified 4.3%

   \[\leadsto \color{blue}{-2 \cdot \left(\left(\sqrt{u1} \cdot u2\right) \cdot \pi\right)} \]

9. Final simplification 4.3%

   \[\leadsto -2 \cdot \left(\pi \cdot \left(u2 \cdot \sqrt{u1}\right)\right) \]
10. Add Preprocessing

Alternative 10: 4.6% accurate, 2.9× speedup

Math

\[\begin{array}{l} \\ \left(\sqrt{u1} \cdot \left(\pi \cdot u2\right)\right) \cdot -2 \end{array} \]

FPCore

(FPCore (cosTheta_i u1 u2)
 :precision binary32
 (* (* (sqrt u1) (* PI u2)) -2.0))

C

float code(float cosTheta_i, float u1, float u2) {
	return (sqrtf(u1) * (((float) M_PI) * u2)) * -2.0f;
}

Julia

function code(cosTheta_i, u1, u2)
	return Float32(Float32(sqrt(u1) * Float32(Float32(pi) * u2)) * Float32(-2.0))
end

MATLAB

function tmp = code(cosTheta_i, u1, u2)
	tmp = (sqrt(u1) * (single(pi) * u2)) * single(-2.0);
end

Derivation

1. Initial program 57.4%

   \[\sqrt{-\log \left(1 - u1\right)} \cdot \sin \left(\left(2 \cdot \pi\right) \cdot u2\right) \]
2. Add Preprocessing

3. Taylor expanded in u1 around 0 -0.0%

   \[\leadsto \color{blue}{\sqrt{u1} \cdot \left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}\right)} \]
4. Step-by-step derivation

   1. associate-*r* -0.0%

      \[\leadsto \color{blue}{\left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot {\left(\sqrt{-1}\right)}^{2}} \]

   2. unpow2 -0.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{\left(\sqrt{-1} \cdot \sqrt{-1}\right)} \]

   3. rem-square-sqrt 4.0%

      \[\leadsto \left(\sqrt{u1} \cdot \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right)\right) \cdot \color{blue}{-1} \]

   4. *-commutative 4.0%

      \[\leadsto \color{blue}{\left(\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \sqrt{u1}\right)} \cdot -1 \]

   5. associate-*l* 4.0%

      \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(\sqrt{u1} \cdot -1\right)} \]

   6. *-commutative 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-1 \cdot \sqrt{u1}\right)} \]

   7. mul-1-neg 4.0%

      \[\leadsto \sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \color{blue}{\left(-\sqrt{u1}\right)} \]

5. Simplified 4.0%

   \[\leadsto \color{blue}{\sin \left(2 \cdot \left(u2 \cdot \pi\right)\right) \cdot \left(-\sqrt{u1}\right)} \]

6. Taylor expanded in u2 around 0 4.3%

   \[\leadsto \color{blue}{-2 \cdot \left(\sqrt{u1} \cdot \left(u2 \cdot \pi\right)\right)} \]

7. Final simplification 4.3%

   \[\leadsto \left(\sqrt{u1} \cdot \left(\pi \cdot u2\right)\right) \cdot -2 \]
8. Add Preprocessing

Reproduce

herbie shell --seed 2024097 
(FPCore (cosTheta_i u1 u2)
  :name "Beckmann Sample, near normal, slope_y"
  :precision binary32
  :pre (and (and (and (> cosTheta_i 0.9999) (<= cosTheta_i 1.0)) (and (<= 2.328306437e-10 u1) (<= u1 1.0))) (and (<= 2.328306437e-10 u2) (<= u2 1.0)))
  (* (sqrt (- (log (- 1.0 u1)))) (sin (* (* 2.0 PI) u2))))