A Regularized B-Spline-Heaviside Collocation Method for Cauchy Singular Integral Equations with Piecewise Hölder Solutions

Abstract

We develop a B-spline-Heaviside collocation method for Cauchy singular integral equations on a smooth closed $C^2$ contour when the exact solution is piecewise Hölder continuous with finitely many prescribed jumps. Since the Cauchy singular integral of a discontinuous function generally has logarithmic terms at the jump points, we study $M=cI+dS+K:X_α\to Y_α$, $X_α=PH^α(Γ,D)$, where $Y_α$ is a logarithmically enlarged piecewise Hölder space. The discontinuous component is represented by a nonredundant system of normalized relative Heaviside functions adapted to the closed contour. Collocation uses point evaluations at spline nodes separated from the jump set together with logarithmic-coefficient functionals at the jumps. Assuming continuous stability of $M:X_β\to Y_β$, mesh-uniform scaled discrete stability of the regularized collocation operators, and a scaled consistency estimate for exact-jump approximants, we prove existence and uniqueness for sufficiently fine meshes and the error bound \[ \|φ-φ^H_{n_B}\|_{X_β} \le C h_B^{α-β}\|φ\|_{X_α}, \qquad 0<β<α<1. \] We give a matrix realization of the regularized scheme, including the logarithmic and point-collocation blocks, singularity-subtracted evaluation of the Cauchy action on splines, principal-value-safe arc formulas for the Heaviside terms, and an implementation algorithm. An abstract perturbation result shows that the same rate is preserved under sufficiently accurate quadrature. Numerical experiments with arcwise errors and finite-dimensional stability and consistency indicators support the theoretical assumptions over the tested meshes.