\[\frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 \cdot \left(\left(1 + \frac{\sqrt{5} - 1}{2} \cdot \cos x\right) + \frac{3 - \sqrt{5}}{2} \cdot \cos y\right)}
\]
↓
\[\frac{\mathsf{fma}\left(\sqrt{2} \cdot \left(\sin x + \sin y \cdot -0.0625\right), \left(\cos x - \cos y\right) \cdot \left(\sin y + \sin x \cdot -0.0625\right), 2\right)}{3 \cdot \left(\mathsf{fma}\left(\cos x, -0.5 + \sqrt{5} \cdot 0.5, 1\right) + \frac{\cos y}{1.5 + \sqrt{1.25}}\right)}
\]
double code(double x, double y) {
return (2.0 + (((sqrt(2.0) * (sin(x) - (sin(y) / 16.0))) * (sin(y) - (sin(x) / 16.0))) * (cos(x) - cos(y)))) / (3.0 * ((1.0 + (((sqrt(5.0) - 1.0) / 2.0) * cos(x))) + (((3.0 - sqrt(5.0)) / 2.0) * cos(y))));
}
↓
double code(double x, double y) {
return fma((sqrt(2.0) * (sin(x) + (sin(y) * -0.0625))), ((cos(x) - cos(y)) * (sin(y) + (sin(x) * -0.0625))), 2.0) / (3.0 * (fma(cos(x), (-0.5 + (sqrt(5.0) * 0.5)), 1.0) + (cos(y) / (1.5 + sqrt(1.25)))));
}
function code(x, y)
return Float64(Float64(2.0 + Float64(Float64(Float64(sqrt(2.0) * Float64(sin(x) - Float64(sin(y) / 16.0))) * Float64(sin(y) - Float64(sin(x) / 16.0))) * Float64(cos(x) - cos(y)))) / Float64(3.0 * Float64(Float64(1.0 + Float64(Float64(Float64(sqrt(5.0) - 1.0) / 2.0) * cos(x))) + Float64(Float64(Float64(3.0 - sqrt(5.0)) / 2.0) * cos(y)))))
end
↓
function code(x, y)
return Float64(fma(Float64(sqrt(2.0) * Float64(sin(x) + Float64(sin(y) * -0.0625))), Float64(Float64(cos(x) - cos(y)) * Float64(sin(y) + Float64(sin(x) * -0.0625))), 2.0) / Float64(3.0 * Float64(fma(cos(x), Float64(-0.5 + Float64(sqrt(5.0) * 0.5)), 1.0) + Float64(cos(y) / Float64(1.5 + sqrt(1.25))))))
end
code[x_, y_] := N[(N[(2.0 + N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * N[(N[Sin[x], $MachinePrecision] - N[(N[Sin[y], $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[y], $MachinePrecision] - N[(N[Sin[x], $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Cos[x], $MachinePrecision] - N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(3.0 * N[(N[(1.0 + N[(N[(N[(N[Sqrt[5.0], $MachinePrecision] - 1.0), $MachinePrecision] / 2.0), $MachinePrecision] * N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(3.0 - N[Sqrt[5.0], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision] * N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x_, y_] := N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * N[(N[Sin[x], $MachinePrecision] + N[(N[Sin[y], $MachinePrecision] * -0.0625), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Cos[x], $MachinePrecision] - N[Cos[y], $MachinePrecision]), $MachinePrecision] * N[(N[Sin[y], $MachinePrecision] + N[(N[Sin[x], $MachinePrecision] * -0.0625), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision] / N[(3.0 * N[(N[(N[Cos[x], $MachinePrecision] * N[(-0.5 + N[(N[Sqrt[5.0], $MachinePrecision] * 0.5), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + N[(N[Cos[y], $MachinePrecision] / N[(1.5 + N[Sqrt[1.25], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{2 + \left(\left(\sqrt{2} \cdot \left(\sin x - \frac{\sin y}{16}\right)\right) \cdot \left(\sin y - \frac{\sin x}{16}\right)\right) \cdot \left(\cos x - \cos y\right)}{3 \cdot \left(\left(1 + \frac{\sqrt{5} - 1}{2} \cdot \cos x\right) + \frac{3 - \sqrt{5}}{2} \cdot \cos y\right)}
↓
\frac{\mathsf{fma}\left(\sqrt{2} \cdot \left(\sin x + \sin y \cdot -0.0625\right), \left(\cos x - \cos y\right) \cdot \left(\sin y + \sin x \cdot -0.0625\right), 2\right)}{3 \cdot \left(\mathsf{fma}\left(\cos x, -0.5 + \sqrt{5} \cdot 0.5, 1\right) + \frac{\cos y}{1.5 + \sqrt{1.25}}\right)}