  There are three main types of numerical computations for the Bessel function
of the second kind: series expansion, continued fraction, and asymptotic
expansion. In addition, they are combined in the appropriate domain for each.
However, there are some regions where the combination of these types requires
sufficient computation time to achieve sufficient accuracy, however, efficiency
is significantly reduced when parallelized. In the proposed method, we adopt a
simple numerical integration concept of integral representation. We coarsely
refine the integration range beforehand, and stabilize the computation time by
performing the integration calculation at a fixed number of intervals.
Experiments demonstrate that the proposed method can achieve the same level of
accuracy as existing methods in less than half the computation time.
