x + \left(y \cdot z\right) \cdot \left(\tanh \left(\frac{t}{y}\right) - \tanh \left(\frac{x}{y}\right)\right)\left(y \cdot \tanh \left(\frac{t}{y}\right) + y \cdot \left(-\tanh \left(\frac{x}{y}\right)\right)\right) \cdot z + xdouble f(double x, double y, double z, double t) {
double r243151 = x;
double r243152 = y;
double r243153 = z;
double r243154 = r243152 * r243153;
double r243155 = t;
double r243156 = r243155 / r243152;
double r243157 = tanh(r243156);
double r243158 = r243151 / r243152;
double r243159 = tanh(r243158);
double r243160 = r243157 - r243159;
double r243161 = r243154 * r243160;
double r243162 = r243151 + r243161;
return r243162;
}
double f(double x, double y, double z, double t) {
double r243163 = y;
double r243164 = t;
double r243165 = r243164 / r243163;
double r243166 = tanh(r243165);
double r243167 = r243163 * r243166;
double r243168 = x;
double r243169 = r243168 / r243163;
double r243170 = tanh(r243169);
double r243171 = -r243170;
double r243172 = r243163 * r243171;
double r243173 = r243167 + r243172;
double r243174 = z;
double r243175 = r243173 * r243174;
double r243176 = r243175 + r243168;
return r243176;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 4.6 |
|---|---|
| Target | 2.1 |
| Herbie | 1.5 |
Initial program 4.6
Simplified2.1
rmApplied add-cube-cbrt2.5
Applied associate-*l*2.5
rmApplied fma-udef2.5
Simplified1.5
rmApplied sub-neg1.5
Applied distribute-lft-in1.5
Final simplification1.5
herbie shell --seed 2020089 +o rules:numerics
(FPCore (x y z t)
:name "SynthBasics:moogVCF from YampaSynth-0.2"
:precision binary64
:herbie-target
(+ x (* y (* z (- (tanh (/ t y)) (tanh (/ x y))))))
(+ x (* (* y z) (- (tanh (/ t y)) (tanh (/ x y))))))