The construction of C2 Pythagorean-hodograph (PH) quintic spline curves that interpolate a sequence of points p0, . . . , pN and satisfy prescribed end conditions incurs a “tridiagonal” system of N quadratic equations in N complex unknowns. Albrecht and Farouki [1] invoke the homotopy method to compute all 2N+k solutions to this system, among which there is a unique “good” PH spline that is free of undesired loops and extreme curvature variations (k ∈ {−1, 0,+1} depends on the adopted end conditions). However, the homotopy method becomes prohibitively expensive when N 10, and efficient methods to construct the “good” spline only are desirable. The use of iterative solution methods is described herein, with starting approximations derived from “ordinary” C2 cubic splines. The system Jacobian satisfies a global Lipschitz condition in CN, yielding a simple closed-form expression of the Kantorovich condition for convergence of Newton–Raphson iterations, that can be evaluated with O(N2) cost. These methods are also generalized to the case of non-uniform knots.

Efficient solution of the complex quadratic tridiagonal system for C^2 PH quintic splines / FAROUKI R.; KUSPA B.; MANNI C.; A. SESTINI. - In: NUMERICAL ALGORITHMS. - ISSN 1017-1398. - STAMPA. - 27:(2001), pp. 35-60. [10.1023/A:1016621116240]

Efficient solution of the complex quadratic tridiagonal system for C^2 PH quintic splines

SESTINI, ALESSANDRA
2001

Abstract

The construction of C2 Pythagorean-hodograph (PH) quintic spline curves that interpolate a sequence of points p0, . . . , pN and satisfy prescribed end conditions incurs a “tridiagonal” system of N quadratic equations in N complex unknowns. Albrecht and Farouki [1] invoke the homotopy method to compute all 2N+k solutions to this system, among which there is a unique “good” PH spline that is free of undesired loops and extreme curvature variations (k ∈ {−1, 0,+1} depends on the adopted end conditions). However, the homotopy method becomes prohibitively expensive when N 10, and efficient methods to construct the “good” spline only are desirable. The use of iterative solution methods is described herein, with starting approximations derived from “ordinary” C2 cubic splines. The system Jacobian satisfies a global Lipschitz condition in CN, yielding a simple closed-form expression of the Kantorovich condition for convergence of Newton–Raphson iterations, that can be evaluated with O(N2) cost. These methods are also generalized to the case of non-uniform knots.
2001
27
35
60
FAROUKI R.; KUSPA B.; MANNI C.; A. SESTINI
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificatore per citare o creare un link a questa risorsa: https://hdl.handle.net/2158/223701
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 50
  • ???jsp.display-item.citation.isi??? 46
social impact