25 research outputs found
Discrete orthogonal projections on multiple knot periodic splines
AbstractThis paper establishes properties of discrete orthogonal projections on periodic spline spaces of order r, with knots that are equally spaced and of arbitrary multiplicity M⩽r. The discrete orthogonal projection is expressed in terms of a quadrature rule formed by mapping a fixed J-point rule to each sub-interval. The results include stability with respect to discrete and continuous norms, convergence, commutator and superapproximation properties. A key role is played by a novel basis for the spline space of multiplicity M, which reduces to a familiar basis when M=1
An integral method for solving nonlinear eigenvalue problems
We propose a numerical method for computing all eigenvalues (and the
corresponding eigenvectors) of a nonlinear holomorphic eigenvalue problem that
lie within a given contour in the complex plane. The method uses complex
integrals of the resolvent operator, applied to at least column vectors,
where is the number of eigenvalues inside the contour. The theorem of
Keldysh is employed to show that the original nonlinear eigenvalue problem
reduces to a linear eigenvalue problem of dimension .
No initial approximations of eigenvalues and eigenvectors are needed. The
method is particularly suitable for moderately large eigenvalue problems where
is much smaller than the matrix dimension. We also give an extension of the
method to the case where is larger than the matrix dimension. The
quadrature errors caused by the trapezoid sum are discussed for the case of
analytic closed contours. Using well known techniques it is shown that the
error decays exponentially with an exponent given by the product of the number
of quadrature points and the minimal distance of the eigenvalues to the
contour