Improve handling of large values by inverse hyperbolic and trigonometric functions (#2928)

The `dpnp.tensor.acosh` function was returning incorrect results
(infinity) for complex numbers with large negative real parts (see
#2924, #2896).

The issue was caused by the threshold used `1/epsilon` to decide when to
switch from acos in SYCL's experimental complex extension
(`exprm_ns::acos()`) to a log-based formula.
The SYCL `exprm_ns::acos()` returns incorrect results for large negative
real parts, due to used formula `z + sqrt(z² - 1)`, which leads to
`log(0) = -infinity`.

The PR proposes to change the threshold from `1/epsilon` to
`sqrt(1/epsilon)/2` for all inverse trig/hyperbolic functions. This is
the precision loss point where `sqrt(z² ± 1)` calculations become
unstable.

Also the PR reduces code duplication by creating shared `casinh()`,
`casin()`, `catanh()`, `catan()` functions in new `complex_math.hpp`
header. 9154e64

Author: github-actions[bot]

reference/comparison.html

@@ -1040,7 +1040,7 @@
 <h1>Comparison Table NumPy/ DPNP/ CuPy<a class="headerlink" href="#comparison-table-numpy-dpnp-cupy" title="Link to this heading">¶</a></h1>
 <p>Here is a list of NumPy and CuPy APIs and its corresponding DPNP implementations.</p>
 <p><code class="docutils literal notranslate"><span class="pre">-</span></code> in DPNP column means that DPNP implementation is not provided yet.</p>
-<p><strong>NumPy(v2.5.0) / DPNP(v0.21.0dev2+8.g8588c5abe1b) / CuPy(v14.1.1) APIs</strong></p>
+<p><strong>NumPy(v2.5.0) / DPNP(v0.21.0dev2+9.g9154e642b41) / CuPy(v14.1.1) APIs</strong></p>
 <section id="summary">
 <h2>Summary<a class="headerlink" href="#summary" title="Link to this heading">¶</a></h2>
 <div class="table-wrapper docutils container">