x \cdot \cos y - z \cdot \sin y
\mathsf{expm1}\left(\mathsf{log1p}\left(\sqrt[3]{\cos y}\right)\right) \cdot \left({\left(\cos y \cdot \cos y\right)}^{\frac{1}{3}} \cdot x\right) - z \cdot \sin ydouble f(double x, double y, double z) {
double r12224363 = x;
double r12224364 = y;
double r12224365 = cos(r12224364);
double r12224366 = r12224363 * r12224365;
double r12224367 = z;
double r12224368 = sin(r12224364);
double r12224369 = r12224367 * r12224368;
double r12224370 = r12224366 - r12224369;
return r12224370;
}
double f(double x, double y, double z) {
double r12224371 = y;
double r12224372 = cos(r12224371);
double r12224373 = cbrt(r12224372);
double r12224374 = log1p(r12224373);
double r12224375 = expm1(r12224374);
double r12224376 = r12224372 * r12224372;
double r12224377 = 0.3333333333333333;
double r12224378 = pow(r12224376, r12224377);
double r12224379 = x;
double r12224380 = r12224378 * r12224379;
double r12224381 = r12224375 * r12224380;
double r12224382 = z;
double r12224383 = sin(r12224371);
double r12224384 = r12224382 * r12224383;
double r12224385 = r12224381 - r12224384;
return r12224385;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.1
rmApplied add-cube-cbrt0.4
Applied associate-*r*0.4
rmApplied pow1/316.4
Applied pow1/316.3
Applied pow-prod-down0.2
rmApplied expm1-log1p-u0.3
Final simplification0.3
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x y z)
:name "Diagrams.ThreeD.Transform:aboutX from diagrams-lib-1.3.0.3, A"
(- (* x (cos y)) (* z (sin y))))