double code(double x, double y, double z, double t, double a) {
return ((double) (((double) (((double) (x * y)) * z)) / ((double) sqrt(((double) (((double) (z * z)) - ((double) (t * a))))))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((z <= -3.38525079024628e+116)) {
VAR = ((double) (y * ((double) (x / ((double) (((double) (((double) (0.5 * ((double) (((double) (a / z)) * t)))) - z)) / z))))));
} else {
double VAR_1;
if ((z <= 1.2268572723557027e-234)) {
VAR_1 = ((double) (y * ((double) (z * ((double) (x / ((double) sqrt(((double) (((double) (z * z)) - ((double) (a * t))))))))))));
} else {
double VAR_2;
if (((z <= 1.3233303272720813e-182) || !(z <= 1.584046164603484e+132))) {
VAR_2 = ((double) (y * ((double) (x / ((double) (((double) (z + ((double) (((double) (((double) (a / z)) * t)) * -0.5)))) / z))))));
} else {
VAR_2 = ((double) (y * ((double) (x / ((double) (((double) sqrt(((double) (((double) (z * z)) - ((double) (a * t)))))) / z))))));
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.9 |
|---|---|
| Target | 7.8 |
| Herbie | 7.2 |
if z < -3.3852507902462798e116Initial program 46.3
Simplified44.0
rmApplied add-cube-cbrt44.2
Applied associate-*l*44.2
Simplified44.2
rmApplied pow144.2
Applied pow144.2
Applied pow-prod-down44.2
Applied pow144.2
Applied pow-prod-down44.2
Applied pow144.2
Applied pow144.2
Applied pow-prod-down44.2
Applied pow-prod-down44.2
Simplified44.0
Taylor expanded around -inf 6.0
Simplified1.6
if -3.3852507902462798e116 < z < 1.22685727235570267e-234Initial program 12.2
Simplified10.7
rmApplied add-cube-cbrt11.3
Applied associate-*l*11.3
Simplified10.8
rmApplied pow110.8
Applied pow110.8
Applied pow-prod-down10.8
Applied pow110.8
Applied pow-prod-down10.8
Applied pow110.8
Applied pow110.8
Applied pow-prod-down10.8
Applied pow-prod-down10.8
Simplified10.7
rmApplied *-un-lft-identity10.7
Applied *-un-lft-identity10.7
Applied sqrt-prod10.7
Applied times-frac10.7
Applied *-un-lft-identity10.7
Applied times-frac10.7
Simplified10.7
Simplified12.0
if 1.22685727235570267e-234 < z < 1.32333032727208128e-182 or 1.58404616460348404e132 < z Initial program 44.7
Simplified43.7
rmApplied add-cube-cbrt43.9
Applied associate-*l*43.9
Simplified43.7
rmApplied pow143.7
Applied pow143.7
Applied pow-prod-down43.7
Applied pow143.7
Applied pow-prod-down43.7
Applied pow143.7
Applied pow143.7
Applied pow-prod-down43.7
Applied pow-prod-down43.7
Simplified43.7
Taylor expanded around inf 10.3
Simplified6.6
if 1.32333032727208128e-182 < z < 1.58404616460348404e132Initial program 9.0
Simplified5.7
rmApplied add-cube-cbrt6.6
Applied associate-*l*6.6
Simplified5.6
rmApplied pow15.6
Applied pow15.6
Applied pow-prod-down5.6
Applied pow15.6
Applied pow-prod-down5.6
Applied pow15.6
Applied pow15.6
Applied pow-prod-down5.6
Applied pow-prod-down5.6
Simplified5.6
Final simplification7.2
herbie shell --seed 2020191
(FPCore (x y z t a)
:name "Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< z -3.1921305903852764e+46) (neg (* y x)) (if (< z 5.976268120920894e+90) (/ (* x z) (/ (sqrt (- (* z z) (* a t))) y)) (* y x)))
(/ (* (* x y) z) (sqrt (- (* z z) (* t a)))))