\[1.000000000000000006295358232172963997211 \cdot 10^{-150} \lt \left|x\right| \lt 9.999999999999999808355961724373745905731 \cdot 10^{149}\]

\[\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{\left(4 \cdot p\right) \cdot p + x \cdot x}}\right)}\]

\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{\left(4 \cdot p\right) \cdot p + x \cdot x}}\right)}

double f(double p, double x) {
        double r9077 = 0.5;
        double r9078 = 1.0;
        double r9079 = x;
        double r9080 = 4.0;
        double r9081 = p;
        double r9082 = r9080 * r9081;
        double r9083 = r9082 * r9081;
        double r9084 = r9079 * r9079;
        double r9085 = r9083 + r9084;
        double r9086 = sqrt(r9085);
        double r9087 = r9079 / r9086;
        double r9088 = r9078 + r9087;
        double r9089 = r9077 * r9088;
        double r9090 = sqrt(r9089);
        return r9090;
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2019315 +o rules:numerics
(FPCore (p x)
  :name "Given's Rotation SVD example"
  :precision binary64
  :pre (< 1.00000000000000001e-150 (fabs x) 9.99999999999999981e149)

  :herbie-target
  (sqrt (+ 0.5 (/ (copysign 0.5 x) (hypot 1 (/ (* 2 p) x)))))

  (sqrt (* 0.5 (+ 1 (/ x (sqrt (+ (* (* 4 p) p) (* x x))))))))

Backtrace

match: no matching clause for 'binary64		L	C
match:error	/opt/racket-7.2/collects/racket/match/runtime.rkt	24	0
for-loop	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt	246	3
make-multiplication-subnode	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt	244	0
make-multiplication-node	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt	234	0
parse-loop54	/opt/racket-7.2/collects/racket/match/compiler.rkt	418	15
map	/opt/racket-7.2/collects/racket/private/map.rkt	35	13
f45	/opt/racket-7.2/collects/racket/match/compiler.rkt	507	40
simplify	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/reduce.rkt	16	0
approximate9	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/core/taylor.rkt	12	0
(unnamed)	/opt/racket-7.2/collects/racket/contract/private/arrow-higher-order.rkt	360	33
location-do	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt	64	0
location-do	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt	64	0
location-do	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt	64	0
location-do	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt	64	0
location-do	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/programs.rkt	64	0
(unnamed)	/opt/racket-7.2/collects/racket/contract/private/arrow-higher-order.rkt	360	33
for-loop	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	151	6
taylor-alt	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	145	0
for-loop	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	166	7
gen-series!	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	159	0
run-iter!	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	326	0
for-loop	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	364	4
run-improve47	/data/pavpan/nightlies/herbie/fix-interface-bugs/src/mainloop.rkt	345	0
(unnamed)	/opt/racket-7.2/collects/racket/private/more-scheme.rkt	261	28
run	/opt/racket-7.2/share/pkgs/profile-lib/main.rkt	39	2
profile-thunk16	/opt/racket-7.2/share/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/opt/racket-7.2/collects/racket/private/more-scheme.rkt	261	28