Fix dpnp.interp returning nan at exact knot when fp contains inf

[antonwolfy]

dpnp.interp returned NaN when a query point landed exactly on an interior knot of xp and the adjacent fp value was Inf. Interpolating at such a point should return the knot value directly.

The root cause is IEEE-754 arithmetic in the linear formula: when fp[x_idx+1] is inf, slope becomes Inf, and since x_val - xp[x_idx] == 0 at the knot, the result is inf * 0 = nan. The existing NaN-recovery fallback also yields nan in this case.

import dpnp as np

x = np.array([2.0])
xp = np.array([1.0, 2.0, 3.0, 4.0])
fp = np.array([1.0, 2.0, np.inf, 4.0])
np.interp(x, xp, fp)
# before: array([nan])
# after:  array([2.])   (matches numpy.interp)

• Have you provided a meaningful PR description?
• Have you added a test, reproducer or referred to an issue with a reproducer?
• Have you tested your changes locally for CPU and GPU devices?
• Have you made sure that new changes do not introduce compiler warnings?
• Have you checked performance impact of proposed changes?
• Have you added documentation for your changes, if necessary?
• Have you added your changes to the changelog?

[github-actions]

View rendered docs @ https://intelpython.github.io/dpnp/index.html

[coveralls]

No base build to compare — fix-interp-nan-at-knots into master

[github-actions]

Array API standard conformance tests for dpnp=0.21.0dev3=py313h509198e_6 ran successfully.
Passed: 1373
Failed: 2
Skipped: 5

[vlad-perevezentsev]

LGTM
Thank you @antonwolfy