arxiv: v1 [math.na] 26 Sep 2016

Transcription

1 Pythagorean-Hodograph B-Spline Curves Gudrun Albrecht a Carolina Vittoria Beccari b Jean-Charles Canonne a Lucia Romani c a Univ Lille Nord de France UVHC LAMAV FR CNRS 96 F-9 Valenciennes France. b Department of Mathematics University of Bologna P.zza Porta San Donato 407 Bologna Italy c Department of Mathematics and Applications University of Milano-Bicocca Via R. Cozzi 0 Milano Italy arxiv: v [math.na] 6 Sep 06 Abstract We introduce the new class of planar Pythagorean-Hodograph PH B-Spline curves. They can be seen as a generalization of the well-nown class of planar Pythagorean-Hodograph PH Bézier curves presented by R. Faroui and T. Saalis in 990 including the latter ones as special cases. Pythagorean-Hodograph B-Spline curves are nonuniform parametric B-Spline curves whose arc-length is a B-Spline function as well. An important consequence of this special property is that the offsets of Pythagorean-Hodograph B-Spline curves are non-uniform rational B-Spline NURBS curves. Thus although Pythagorean-Hodograph B-Spline curves have fewer degrees of freedom than general B-Spline curves of the same degree they offer unique advantages for computer-aided design and manufacturing robotics motion control path planning computer graphics animation and related fields. After providing a general definition for this new class of planar parametric curves we present useful formulae for their construction discuss their remarable attractive properties and give some examples of their practical use. Keywords: Plane curve; Non-uniform B-Spline; Pythagorean-Hodograph; Arc-length; Offset; G /C Hermite Interpolation. Introduction The purpose of the present article is to introduce the general concept of Pythagorean Hodograph PH B Spline curves. On the one hand B Spline curves since their introduction by Schoenberg [] in 946 have become the standard for curve representation in all areas where curve design is an issue see e.g. [ ]. On the other hand the concept of polynomial PH curves has widely been studied since its introduction by Faroui and Saalis in [0]. The essential characteristic of these curves is that the Euclidean norm of their hodograph is also polynomial thus yielding the useful properties of admitting a closed form polynomial representation of their arc length as well as exact rational parameterizations of their offset curves. These polynomial curves are defined over the space of polynomials using its Bernstein basis thus yielding a control point or so called Bézier representation for them. Rational and spatial counterparts of polynomial PH curves have as well been proposed and most recently an algebraic trigonometric counterpart so called Algebraic Trigonometric Pythagorean Hodograph ATPH curves have been introduced in [0]. So far a general theory for B Spline curves having the PH property is missing. To the best of the authors nowledge the only partial attempt in this direction has been made in [7] where the problem of determining a B Spline form of a C PH quintic spline curve interpolating given points is addressed. Prior to this based on [ 8] in [9] a relation between a planar C PH quintic spline curve and the control polygon of a related C cubic B-Spline curve is presented. The present article shows how to construct a general PH B Spline curve of arbitrary degree over an arbitrary not sequence. To this end we start by defining the complex variable model of a B Spline curve zt of degree n defined Corresponding author. addresses: gudrun.albrecht@univ-valenciennes.fr Gudrun Albrecht carolina.beccari@unibo.it Carolina Vittoria Beccari jean-charles.canonne@univ-valenciennes.fr Jean-Charles Canonne lucia.romani@unimib.it Lucia Romani Preprint submitted to Elsevier September 7 06

2 over a not partition µ. We then square zt by using results for the product of normalized B Spline basis functions from [ 8]. Here the determination of the required coefficients involves the solution of linear systems of equations. Finally the result is integrated in order to obtain the general expression of the PH B Spline curve i.e. its B Spline control points and its not partitionρ. General formulae are derived also for the parametric speed the arc length and the offsets of the resulting curves. The interesting subclasses of clamped and closed PH B spline curves are discussed in great detail. When the degree is and explicit expressions of their control points are given together with the B-Spline representation of the associated arc-length and the rational B-Spline representation of their offsets. Finally clamped quintic PH B Spline curves are used to solve a second order Hermite interpolation problem. The remainder of the paper is organized as follows. In section we recall the basic definition of B Spline curves as well as the Pythagorean Hodograph PH concept thus defining the notion of a PH B Spline curve. In section the general construction of PH B Spline curves is developed section. and then adapted to the important particular cases of clamped and closed PH B Spline curves section.. In section 4 general formulae for their parametric speed arc length and offsets are given. Section is devoted to presenting the explicit expressions regarding clamped and closed PH B Spline curves of degree and. Finally in section 6 we solve a second order Hermite interpolation problem by clamped PH B Spline curves of degree. Conclusions are drawn in section 7.. Preliminary notions and notation While a Bézier curve is univocally identified by its degree a B-Spline curve involves more information namely an arbitrary number of control points a not vector and a degree which are related by the formula number of nots - number of control points = degree. For readers not familiar with B-Spline curves we first recall the definition of normalized B-Spline basis functions and successively the one of planar B-Spline curve see e.g. []. Definition. Letµ={t i R t i t i } i Z be a sequence of non-decreasing real numbers called nots and let n N. The i-th normalized B-spline basis function of degree n defined over the not partitionµ is the function Niµ n t having support [t i t in ] and defined recursively as where and 0 0 = 0. N n iµ t= t t i t in t i N n iµ t t in t t in t i N n iµ t N 0 iµ t= { if t [ti t i 0 otherwise Definition. Let m n N with m nµ={t i }...mn be a finite not partition and s 0... s m R. Then the planar parametric curve m st= s i Niµ n t t [t n t m ] is called a planar B Spline curve of degree n associated with the not partitionµ with de Boor points or control points s 0... s m. Remar. If M i denotes the multiplicity of the not t i then st is of continuity class C n max im i t n t m. If the not vectorµ does not have any particular structure the B-Spline curve st will not pass through the first and last control points neither will be tangent to the first and last legs of the control polygon. In this case st is simply called open B-spline curve. In order to clamp st so that it is tangent to the first and the last legs at the first and last control points respectively as a Bézier curve does the multiplicity of the first and the last not must be adapted. For later use the precise conditions we use for identifying a clamped B-Spline curve are the following. Remar. If t 0 = t =...=t n and t m =...=t mn then the B-Spline curve st given in Definition verifies st n =s 0 st m =s m

3 as well as s t n = n s s 0 s n t m = s m s m t n t t mn t m i.e. st is a clamped B Spline curve. Moreover if all the nots t n... t m are simple then st C n t n t m. On the other hand to mae the B-Spline curve st closed some not intervals and control points must be repeated such that the start and the end of the generated curve join together forming a closed loop. The precise conditions to be satisfied by not intervals and control points in order to get a closed B-Spline curve are recalled in the following. Remar. If in Definition we replace m by mn consider the not partitionµ={t i }...mn with t m t m = t t for =... n assume s 0...s m R to be distinct control points and s m = s 0... s mn = s n then the B-Spline curve mn st= s i Niµ n t t [t n t mn ] has the additional property st n =st mn i.e. st is a closed B Spline curve. Moreover if all the nots t n... t mn are simple then st C n [t n t mn ]. At this point we have all the required preliminary notions to generalize the definition of Pythagorean-Hodograph Bézier curves see [0] to Pythagorean-Hodograph B-Spline curves. Definition. For p n N p n let ut vt and wt be non-zero degree-n spline functions over the not partition µ={t i }...pn i.e. let p p p ut= u i Niµ n t vt= v i Niµ n t wt= w i Niµ n t t [t n t p ] with u i v i w i R for all... p such that ut and vt are non-constant and do not have a non constant spline function over the partition µ as common factor. Then the planar parametric curve xt yt whose coordinate components have first derivatives of the form x t=wt u t v t and y t=wtutvt is called a planar Pythagorean-Hodograph B-Spline curve or a planar PH B-Spline curve of degree n. Indeed as in the case of PH polynomial Bézier curves [9 0] the parametric speed of the plane curve xt yt is given by σt := x t y t = wt u tv t and its unit tangent unit normal and signed curvature are given respectively by t= u v uv n= uv v u κ= uv u v u v u v wu v where for conciseness in the parameter t is omitted. In the following we will restrict our attention to the so-called primitive case wt=. Since in this case equation simplifies asσt=u tv t the primitive case coincides with the regular case. In this case the representation may be obtained by squaring the complex function zt=utivt yielding z t=u t v tiutvt. The coordinate components x t y t of the hodograph r t of the parametric curve rt=xt yt are thus given by the real and imaginary part of z t respectively. In the remainder of the paper we will exclusively use this complex notation and we will thus write r t= x tiy t=u t v tiutvt=z t 4 as also previously done for planar PH quintics [9 0]. Since by construction r t is a degree-n B-spline curve then the PH B-Spline curve rt= r tdt has degree n. In the next section we will construct the corresponding not vector and thus now the continuity class of the resulting PH B-Spline curve rt.

4 . Construction of Pythagorean Hodograph B Spline curves.. The general approach We start with a not partition of the form µ={t i }...pn over which a degree-n B-Spline curve zt=utivt is defined for t [t n t p ]. Thus according to Definition the planar parametric curve zt can be written as zt= p z i Niµ n t t [t n t p ] 6 where z i = u i iv i... p. To express the product z t as a B-Spline curve according to [ 8] we have to augment the multiplicity of each single not t i to n. We thus obtain the not partition ν={s i }...pnn ={< t i > n }...pn 7 where< t i > denotes a not t i of multiplicity. The product z t is thus a degree-n B-Spline curve over the not partitionν which can be written in the form pt=z t= p p q z i z j Niµ n tnn jµ t= p Nν n t 8 j=0 =0 with q=npn according to Definition. Our goal is thus to obtain the explicit expressions of the coefficients p for = 0... q. To this end we set f i j t := Niµ n tnn jµ t and loo for the unnown coefficients χ i j := χ i j 0 χi j...χi q j T i j=0... p such that f i j t= q =0 χ i j Nn ν t. 9 For accomplishing this we apply the method from [] as follows. Let< > be an inner product of the linear space of B-splines of degree n with not vectorν. According to [] for any pair of functions at bt defined over the interval [t 0 t pn ] we use< at bt>= t pn atbt dt to construct the q q linear equation system t 0 with and Aχ i j = b i j 0 tpn A=a l l=0...q a l := Nν n Nn lν = Nν n t Nn lν t dt t 0 tpn b i j = b i j l l=0...q b i j l := f i j Nlν n = t 0 f i j t Nlν n t dt. Since{N n ν t} =0...q are linearly independent the matrix A is a Gramian and therefore nonsingular. This allows us to wor out the unnown coefficients{χ i j } =0...q by solving the linear system in 0. Remar 4. Since b i j = b ji for all i j=0... p thenχ i j =χ ji for all i j=0... p. Therefore the unnown vectors to be obtained from 0 are indeedχ i j... p j=0... i. Being A a non-singular Gramian matrix all the corresponding linear systems always have a unique solution. Moreover since a h = 0 if h >n A is not only symmetric and positive definite but also of band form. Thus by applying the Cholesy decomposition algorithm one 4

5 can compute the factorization A=LL T where L is a lower triangular matrix of the same band form of A i.e. such that l h = 0 if h >n. The non-zero elements of L may be determined row by row by the formulas l h l hh = a h h s=h n l hs l s /l =h n... h χ i j h = a hh h s=h n l hs with the convention that l rc = 0 if c 0 or c>r. Hence the solution of each linear system in 0 can be easily obtained by solving the two triangular linear systems Ly i j = b i j and L T χ i j = y i j via the formulas y i j h = b h h =h n l h y i j /lhh h=0... q = y i j h hn =h l hχ i j /lhh h=0... q where a similar convention as above is adopted with respect to suffices outside the permitted ranges see [7]. From the computed expressions ofχ i j = 0... q 0 i j p we thus get q p p p p pt=z t= =0 j=0 χ i j z iz j N n ν t and p = The resulting PH B-Spline curve rt is now obtained by integrating pt as: j=0 χ i j z iz j =0... q. q rt= ptdt = r i Niρ n t t [t n t p ] whereρ={s i }...pnn with s i = s i for i=... pnn t = s 0 s and s pnn s pnn = t pn i.e. ρ={t {< t > n } =0...pn t pn } with the additional nots t t pn as well as for... q and arbitrary r 0. r i = r i s in s i n p i = r i s in s i n p i 4 Remar. Note that by construction s n = s n = t n as well as s q = s q = t p namely the B-spline curves zt pt and rt are defined on the same domain. If the not partitionµcontains simple inner nots t n... t p then the degree-n spline zt C n t n t p. As a consequence the degree-n spline r t C n t n t p and the degree-n spline rt C n t n t p... Construction of clamped and closed PH B-Spline curves We now consider the conditions for obtaining a clamped respectively closed PH B-Spline curve rt. Proposition. Let rt be the PH B-Spline curve in defined over the not partitionρin where for a clamped respectively closed PH B-Spline curve rt we assume p = m respectively p = m n. a For rt to be clamped i.e. satisfying the following conditions have to be fulfilled rt n =r 0 rt m =r q r t n = n s r r 0 r n t m = n s s qn s q r q r q n r n n B n α=r 0 =0 and n p n n B n α=p 0 with α= t n t n 6 t n t n =0

6 n r mn B n β=r q and =0 n p mn B n β=p q with β= t m t m 7 t m t m =0 where B n t= n i t i t n i =0... n denote the Bernstein polynomials of degree n. b For rt to be closed and of continuity class C n at the junction point rt n =rt mn we require the fulfillment of the conditions mnn s jn s j p j = 0 for =0... n 8 and j=nn t m t m = t t for =n n. 9 Proof. According to [] for every degree-n B-Spline curve xu = q x in n iµ u Rd over the not partition µ={t i }...nq there exists a unique multi-affine symmetric application or blossom X : R n R d u... u n Xu...u n such that xu= Xu... u= X< u> n. Its control points are x i = Xt i... t ni... q. Thus denoting Pu...u n the blossom of the curve pu from 8 over the not partitionνfrom 7 the control points p i may be written as: p = P< t 0 > n < t > min{n n} < t > max{ 0} for =... n p jn = P< t j > n < t j > min{n n} < t j > max{ 0} for j=... pn and =0... n p pnn = P< t pn > n < t pn > min{n n} for = 0. Analogously denoting Ru... u n the blossom of the curve ru from over the not partitionρfrom the control points r i may be written as r jn = R< t j > n < t j > min{nn} < t j > max{ 0} for j=0... pn and =0... n r pnn = R< t pn > n < t pn > min{nn} for =0. Recalling de Boor s algorithm and the properties of blossoms the control points involved for calculating a point rt j =R< t j > n are the following: R< t j > n < t j > n =r jn n... R< t j > n < t j > n =r jn 0 a We wish to obtain a clamped curve satisfying conditions. In order to satisfy the positional constraints we thus apply de Boor s algorithm for calculating rt n respectively rt m to the control points 0 for j=n respectively j=m. For j=nwe obtain R< t n > n l < t n > n < t n > l = α R< t n > n l < t n > n < t n > l α R< t n > n l < t n > n < t n > l for l=... n. This yields the following condition for α which results thus to be independent of the indices l: α t n α t n = t n. De Boor s algorithm thus degenerates to de Casteljau s algorithm yielding the first equation of condition 6. The first equation of condition 7 is obtained analogously. In order to satisfy the tangential constraints of we first note that they are equivalent to the following positional constraints for pt: pt n =p 0 pt m =p q. We thus apply the same reasoning as above to pt and obtain the second equations in 6 and 7. b We wish to obtain a closed curve rt with rt n =rt mn. Recalling de Boor s algorithm and the properties of blossoms the control points involved for calculating a point rt j =R< t j > n are the following: R< t j > n < t j > n =r jn n... R< t j > n < t j > n =r jn 6

7 In order for condition to hold the following points and their corresponding not intervals thus have to coincide: r nn = r mnn for =0... n as well as t m t m = t t for =n n. 4 Setting f j = s jn s j n for j = 0... q and considering condition 4 condition is equivalent to or equivalently nn j=0 f j p j = mnn j=nn mnn j=0 f j p j for =0... n f j p j = 0 for =0... n. These conditions also guarantee the maximum possible continuity class at the junction point. In the clamped case of the above proposition we notice that if t n = t n thenα=0 which yields rt n =r nn n and pt n =p nn n and if t m = t m thenβ= which yields rt m =r mn and pt m =p mn. If in the not partitionµ from we have t 0 =...=t n and t m =...=t mn i.e. if zt from 6 is a clamped curve itself the first respectively last n n control points of rt and pt coincide i.e. as well as r 0 =...=r nn n and r mn =...=r mnn p 0 =...=p nn n and p mn =...=p mnn. 6 In this case condition from Proposition a is automatically satisfied yielding a more intuitive way of obtaining a clamped PH B-Spline curve. In order to simplify the notation we remove redundant nots in the not partition ν from 7 together with the control point multiplicities from and summarize the result in the following Corollary. Corollary. Let zt= m z i N n iµ t t [t n t m ] be a clamped B-Spline curve over the not partition µ={< t n > n {t i } i=n...m < t m > n } 7 as in Remar. Then where q=nnm n and pt=z t= q =0 p Nν n t 8 as well as ν={s i }...4nnm n ={< t n > n {< t i > n } i=n...m < t m > n } q rt= ptdt = r i Niρ n t t [t n t m ] whereρ={s i }...4nnm n with s i= s i for i=... 4nnm n s 0 = s and s 4nnm n = s 4nnm n i.e. ρ={< t n > n {< t i > n } i=n...m< t m> n } 9 and the control points r i satisfy 4. 7

8 Proof. The result is obtained by removing n n of the multiple control points from together with n of the multiple nots ofνfrom 7 at the beginning and at the end. The same result is obtained by proceeding with the general construction of the PH B-Spline curve starting with a clamped B-Spline curve over the not partitionµ from 7. In this way in both cases the clamped and the closed one the resulting PH B-Spline curve is of degree n and of continuity class C n. For example for n= we obtain PH B-Spline curves of degree and continuity class C while for n= we have PH B-Spline curves of degree and continuity class C. 4. Parametric speed arc length and offsets According to the parametric speed of the regular PH curve rt = xt iyt is given by σt= r t = z t =zt zt. Exploiting 6 we thus obtain σt= p p q z i z j Niµ n tnn jµ t= σ Nν n t 0 j=0 =0 where in analogy to. σ = p p j=0 χ i j z i z j 4.. Arc-length where The arc length of the PH B-Spline curve is thus obtained as l i = l i s in s i n with l 0 = 0. The cumulative arc length is given by lξ= and the curve s total arc length thus is L=lt p = q σtdt = l i Niρ n t t [t n t p ] ξ t n σ i = l i s in s i n q σtdt= l i N n iρ ξ Nn iρ t n tp t n σ i q σtdt= l i N n iρ t p Niρ n t n. 4 This general formula simplifies in the clamped and closed cases as follows. In the clamped case for the not partition ρ from 9 we notice that { if Niρ n t n= 0 else and { if i=nnm n Niρ n t m= 0 else. 8

9 Considering l 0 = 0 and recalling Corollary in the clamped case the total arc length L in 4 thus becomes L=l nnm n. Due to the structure of the not partitionρin we notice that { Niρ n 0 if n n<i<nn t n = 0 else and In this case the total arc length L from 4 thus becomes L = = pn i=pn n =0 { 0 if pn<i<pn Niρ n t p = 0 else. l i N n iρ t p nn i=n n l i N n iρ t n lpn N n pnρ t p l n n N n n nρ t n. 6 In the case of a closed curve from Proposition b we have p=mn and conditions 9. On the not partitionρ from the normalized B Spline basis functions having as support [t n t n ] are N n t for i=n n... nn and those having as support [t mn t mn ] are Niρ n t for i=mnn... mnn. With conditions 9 this means N n n nρ t n=n n mnnρ t mn for =0... n. In the case of a closed curve its total arc length L from 6 thus reads L= n =0 which by taing into account becomes n L= =0 lpn l n n N n n nρ t n 7 mnn j=n n s jn s j n iρ σ j Nn n nρ t n Offsets The offset curve r h t at signed distance h of a PH B-Spline curve rt is the locus defined by where nt= r h t=rth nt y t x t t = iz t. x t y t = ir σt σt Note that since we are dealing with the regular caseσt 0and the offset curve is always well defined. Thus Herein the product σtrt reads as r h t= σtrt i h z t σt. σtrt = q q j=0 σ j r i N n iρ tnn jν t. 9

10 Again according to [ 8] we can write with the not partition and w Niρ n tnn jν t= ζ i j N4n τ =0 t 9 τ={< t > n {< t > n } =0...pn < t pn > n } 40 Remar 6. In the clamped case from Corollary we obtain w=npn. 4 τ={< t n > 4n {< t > n } =n...m < t m > 4n } 4 and Differently in the closed case we have w=4nm nn. 4 τ={< t > n {< t > n } =0...mn < t mn > n } 44 and w=nmn. 4 To wor out the unnown coefficientsζ i j := ζ i j 0 ζi j...ζi j w T in 9 we solve the linear system Cζ i j = e i j 46 with and tpn C=c h h=0...w c h := Nτ 4n N4n hτ = Nτ 4n t N4n hτ t 0 t dt tpn e i j = e i j h h=0...w e i j h := g i j Nhτ 4n = g i j t Nhτ 4n t dt where g i jt := Niρ n t 0 t Nn Lie in the previous case C is a banded Gramian and thus nonsingular. This guarantees that each of the linear systems in 46 has a unique solution that can be efficiently computed by means of the Cholesy decomposition algorithm for symmetric positive definite band matrices. The computed expressions ofζ i j =0... w... q j=0... q thus yield σtrt = w q q =0 j=0 ζ i j σ jr i N 4n τ t. jν t. By writing and we thus obtain pt=pt = q q p j N n jν t Niρ n t j=0 j=0 σt=σt = q q σ j N n jν t Niρ n t pt= w q q =0 j=0 0 ζ i j p jn 4n τ t

11 The offset curve r h t finally has the form w q q σt= ζ i j σ jnτ 4n t. =0 j=0 where for =0... w and r h t= q = q γ = w =0 q Nτ 4n t w =0 γ 47 Nτ 4n t q j=0 q σ j r i i h p j ζ i j 48 q j=0 σ j ζ i j. 49. General explicit formulas for the cubic and quintic case.. Clamped cubic PH B-Splines n = Let m N m. For a general not vectorµ={ 0 < t <...<t m < t m } satisfying the constraints t 0 = t = 0 and t m = t m see Figure first row by applying the above method we construct the not partitions ν={ 0 < t <...< t m < t m } ρ={ 0 4 < t <...< t m < t m 4 } τ={ 0 6 < t <...< t m < t m 6 } illustrated in Figure. Then by solving the linear systems 0 we calculate the coefficientsχ i j 0 i j m 0 m. All of them turn out to be zero with the exception of = = 0... m χ =χ = =0... m. 0 χ In addition we compute the coefficientsζ i j 0 i m 0 j m 0 m as the solutions to the linear systems 46. All of them turn out to be zero with the exception of ζ = d D ζ = d D =0... m ζ = d D ζ = ζ = d D =0... m ζ = d 0D ζ = ζ = d 0D ζ = 0 =0... m ζ = 0 ζ = d 0D ζ = ζ = d 0D = 0... m ζ 4 = d D ζ 4 = ζ 4 = d D = 0... m where D := d d = 0... m and d 0 = d m := 0. By means of the computed coefficients{χ i j j m }0 i 0 m we can thus shortly write the control points of r t as p = z =0... m p = z z =0... m and the coefficients of the parametric speedσt as σ = z z =0... m σ = z z z z = 0... m.

12 d d d m d m 0 t t... t m t m t m d d d m d m 0 t t... t m t m t m d d d m d m 0 t t... t m t m t m d d d m d m 0 t t... t m t m t m Figure : Knot partitions for the clamped case n =. From top to bottom: µ ν ρ τ. Thus according to the clamped cubic PH B-Spline curve defined over the not partition ρ is given by rt= m r i N iρ t t [t t m ] t = 0 with control points and arbitrary r 0. r = r 0 d z 0 r i = r i d i z iz i... m r i = r i d id i z i... m r m = r m d m z m Remar 7. Note that when m= and t = t = the expressions of the control points coincide with those of Faroui s PH Bézier cubic from [4 0]. According to the total arc length of the clamped PH B-Spline curve of degree is given by L=l m where l 0 = 0 l = l 0 d z 0 z 0 l i = l i d i 6 zi z i z i z i... m l i = l i d id i z i z i... m l m = l m d m z m z m. Over the not partitionτ the offset curve r h t has the rational B-Spline form r h t= m =0 q N τ t m =0 γ N τ t t [t t m ] where by exploiting the explicit expressions of the coefficients{ζ i j 0 j m }0 i m 0 m from weights and control points can easily be obtained; they are reported in the Appendix. Some examples of clamped cubic PH B-Spline curves are shown in Figure and their offsets are displayed in Figure... Clamped quintic PH B-Splines n = Let m N m. For a general not vectorµ={ 0 < t <...<t m < t m } satisfying the constraints t 0 = t = t = 0 and t m = t m = t m see Figure 4 first row by applying the above method we construct the not

13 r 0 n= m= n= m= n= m= r r r r 0 r r r 4 r 0 r r 4 r r r r r 6 r 7 r Figure : Clamped cubic PH B-Spline curves with: m = left m = center and m = right. n= m= n= m= n= m= Figure : Offsets of clamped cubic PH B-Spline curves from Figure with. row and without. row control polygon where: m= left column m=center column and m= right column. partitions ν={ 0 < t <...< t m < t m } ρ={ 0 6 < t <...< t m < t m 6 } τ={ 0 0 < t 8 <...< t m 8 < t m 0 }

14 illustrated in Figure 4. Then by solving the linear systems 0 we calculate the coefficientsχ i j 0 i j m 0 m. All of them turn out to be zero with the exception of χ 00 0 = χ 0 =χ 0 = χ =χ = 6 χ χ = d d d d d d =χ = 6 d d d d d d χ =χ = 6 χ = d d d =... m χ =χ = χ =χ = =... m d d d d d =... m =... m d d d d d d =... m d d d =... m d d d =... m χ = d d d =... m χ m m m χ mm m = =χmm m = where d = t t =... m and d 0 = d m = 0. In addition we compute the coefficientsζ i j 0 i m 0 j m 0 8m 7 as the solutions to the linear systems 46. All of them turn out to be zero with the exception of ζ 8 = d ζ D 8 = d d ζ 9D 8 = 4d ζ 9D 8 = d d ζ 9D 8 = d = 0... m 9D ζ 8 = d ζ 9D 8 = d d ζ 9D 8 = d d ζ 9D 8 = d ζ D 8 = 4d =0... m 9D ζ 8 = d ζ 6D 8 = d d ζ D 8 = d ζ 6D 8 = d 9D ζ 8 = d 8D ζ 8 = d 8D = 0... m ζ 8 = d ζ D 8 = d d ζ D 8 = d ζ D 8 = d 4D ζ 8 = d 4D ζ 8 = 0 ζ 8 = d 4D ζ 8 = d 4D = 0... m ζ 4 d 84 = ζ 4 6D D 84 = d d ζ 4 6D D 84 = d d ζ 6D D 84 = d d 6D D ζ 84 = d d d 0d 6D D 6D ζ 84 = d d 0d 6D D 6D ζ 84 = 0 ζ 84 = d d 6D D 0 6 ζ 4 84 = d d 6D D ζ 84 = d d 6D D ζ 4 84 = d 6D D =0... m ζ 4 8 = d 6D D ζ 4 8 = d d 6D D ζ 8 = d d 6D D ζ 8 = d d 6D D 0 6 ζ 8 = 0 ζ 8 = d d 0d 6D D 6D ζ 4 8 = d d d 0d 6D D 6D ζ 8 = d d ζ 6D D 8 = d d ζ 4 6D D 8 = d d ζ 6D D 8 = = 0... m 6D D ζ 4 86 = d 4D ζ 86 = d 4D ζ 86 = 0 ζ4 86 = d 4D ζ 86 = d 4D ζ 86 = d ζ 4 D 86 = d d ζ D 86 = d =0... m D ζ 4 87 = d 8D ζ 87 = d 8D ζ 4 87 = d 9D ζ 87 = d 6D ζ 4 87 = d d ζ D 87 = d =0... m 6D d 4 where D := d d = 0... m with d 0 = d m := 0. 4

15 d d d m d m 0 t t 4... t m t m t m d d d m d m 0 t t 4... t m t m t m d d d m d m 0 t t 4... t m t m t m d d d m d m 0 t t 4... t m t m t m Figure 4: Knot partitions for the clamped case n =. From top to bottom: µ ν ρ τ. By means of the computed coefficients{χ i j j m }0 i 0 m we can thus shortly write the control points of r t as p 0 = z 0 p = z 0 z p = z d z d z d z d z d d d d =... m d p = z z d z d d =... m d p = z z d z d d =... m p m = z m z m p m = z m and the coefficients of the parametric speedσt as σ 0 = z 0 z 0 σ = z 0 z z z 0 d d σ = 6 d d d d z z z z 6 d d d d z z z z σ = d d d z z σ = σ m = z m z m z m z m σ m = z m z m d d d d d d d z z d d 6 d d d d z z z z =... m d d d z z z z =... m d d d z z z z d d d z z =... m with d 0 := 0 and d m := 0. Thus according to the clamped quintic PH B-Spline curve defined over the not partition ρ is given by rt= m r i N iρ t t [t t m ] t = 0

16 with control points and arbitrary r 0. r = r 0 d z 0 r = r d z 0 z r i = r i d i z i di z i d i z i di z i d i z i d i d i d i d i i=... m r i = r i z i d iz i d i z i i=... m r i = r i z i d i z i d i z i i=... m r m = r m d m z m z m r m = r m d m z m Remar 8. Note that when m= and t = t 4 = t = the expressions of the control points coincide with those of Faroui s PH Bézier quintic from [4 0]. According to the total arc length of the clamped PH B-Spline curve of degree is given by L=l m where l 0 = 0 l = l 0 d z 0 z 0 l = l d 0 z 0 z z z 0 l i = l i d i z i z i d id i d i d i d i d i d i di d i z i z i d i z i z i z i z i z i z i z i z i d i d i z i z i z i z i i=... m l i = l i di z i z i d i z i z i z i z i i=... m l i = l i di z i z i d i z i z i z i z i i=... m l m = l m d m 0 z m z m z m z m l m = l m d m z m z m. Over the not partitionτ the offset curve r h t has the rational B-Spline form 8m 7 =0 q N 9 τ t where by exploiting the explicit expressions of the coefficients{ζ i j 0 8m 7 points are easily obtained; the explicit formulae are reported in the Appendix. r h t= t [t t 8m 7 m] 6 γ Nτ 9 t =0 }0 i m 0 j m from 4 weights and control n= m= n= m= r 8 n= m=4 r 0 r r 0 r r r r r 7 r r r 6 r 4 r 0 r r r 0 r 9 r r r 6 r 8 r 4 r 4 r r r 7 Figure : Clamped quintic PH B-Spline curves with: m=left m= center and m=4 right. Some examples of clamped quintic PH B-Spline curves are shown in Figure and their offsets are displayed in Figure 6. 6

17 n= m= n= m= n= m=4 Figure 6: Offsets of clamped quintic PH B-Spline curves from Figure with. row and without. row control polygon where: m= left column m=center column and m=4 right column... Closed cubic PH B-Splines n = Let m N m. For a general not vectorµ={0=t 0 < t <...<t m } see Figure 7 first row by applying the above method we construct the not partitions ν={ 0 < t <...< t m } ρ={t < 0 < t <...< t m < t m4 } τ={ t < 0 < t <...< t m < t m4 } illustrated in Figure 7. Then by solving the linear systems 0 we calculate the coefficientsχ i j 0 i j m 0 m 4. All of them turn out to be zero with the exception of χ = =0... m χ =χ = =0... m. In addition we compute the coefficientsζ i j 0 i m 0 j m4 0 m9 as the solutions to the linear systems 46. All of them turn out to be zero with the exception of ζ 00 = d 0 D 0 ζ 00 4 = ζ0 4 = d 0 0D 0 ζ 0 = 0 ζ0 = ζ0 6 = d D ζ 6 = ζ0 6 = d D ζ 7 = d D ζ 7 = d D =0... m ζ 8 = d D ζ 8 = ζ 8 = d D =0... m ζ 9 = d 0D ζ 9 = ζ 9 = d 0D ζ 9 = 0 =0... m ζ 0 = 0 ζ 0 = d 0D ζ 0 = ζ4 0 = d 0D =0... m ζ = d D ζ = ζ4 = d D =0... m ζ mm4 m = d m D m ζ m4m m = ζm4m4 m = d m D m ζ m4m4 m4 = ζmm m4 = 0 ζ mm m = d m4 0D m ζ mm4 m = ζmm4 m6 = d m4 D m 7 7

18 d d d m d m d m 0 t t... t m t m t m t m d d d m d m d m 0 t t... t m t m t m t m d 0 d d d m d m d m d m4 t 0 t t... t m t m t m t m t m4 d 0 d d d m d m d m d m4 t 0 t t... t m t m t m t m t m4 Figure 7: Knot partitions for the closed case n =. From top to bottom: µ ν ρ τ. where D := d d = 0... m. By means of the computed coefficients{χ i j j m }0 i 0 m4 we can thus shortly write the control points of r t as and the coefficients of the parametric speedσt as p 0 = 0 p = z =0... m p = z z =0... m p m4 = 0 σ 0 = 0 σ = z z =0... m σ = z z z z =0... m σ m4 = 0. Thus according to the closed cubic PH B-Spline curve defined over the not partition ρ is given by with control points rt= m r i N iρ t t [t t m ] t 0 = 0 r = r 0 r i = r i d id i z i... m r i = r i d i z iz i... m r m = r m4 and arbitrary r 0. Note that due to condition 9 d m = d d m = d. Moreover z m and z m must be suitably fixed in order to satisfy condition 8. According to 7 and using the partition of unity N ρ t N ρ t = the total arc length of the closed PH B-Spline curve of degree is given by where L=l m l m4 l m l N ρ t l = l 0 = 0 l i = l i d id i z i z i... m l i = l i d i 6 zi z i z i z i... m l m = l m4. Over the not partitionτthe offset curve r h t has the rational B-Spline form m9 =0 q Nτ r h t= t m9 =0 γ Nτ t t [t t m ] 8 8

19 m zt open/closed Conditions on zt in order to satisfy 8 Illustration closed z = i z 0 z = z 0 Figure 8 first column closed z = d z 0 d z r d d z = z 0 Figure 8 second column closed z = d z 0 d 4 z R d d 4 z 4 = z 0 Figure 8 third column open z = d d d d d d i d d z 0 z = z 0 Figure 9 first column open z = d z 0 d z r d d z = z 0 Figure 9 second column open z = d z 0 d 4 z R d d 4 z 4 = z 0 Figure 9 third column Table : Data for the examples in the closed case for n=. where by exploiting the explicit expressions of the coefficients{ζ i j 0 m9 from 7 weights and control points can easily be obtained; their explicit formulae are reported in the Appendix. We now illustrate some closed degree PH B-Spline curves. To this end we define }0 i m 0 j m4 r ± = 4d d 4d d 4d d d z 0 4d d ± d d 4d d z 0 z 4d d 4d d 4d d d z R ± = 4d d 4d d 4 4d d 4 d z 0 4d d 4d d 4 z 0 z ± d d 4 z 0 z 4d d 4d d 4d d 4 4d d 4 z 4d d 4d d 4 z z 4d d 4d d 4 4d d 4 d 4 z. Figures 8 and 9 contain examples of cubic PH B-Spline curves obtained for the data in Table. Figure 0 shows the offsets of the PH B-Spline curves of Figure 9. Deg- spline zt- m= Deg- spline zt- m= z z 0 z Deg- spline zt- m= z z z 0 z 4 z 0 z z z z r Deg- PH-spline rt - m= r 7 r r 0 r 6 Deg- PH-spline rt - m= Deg- PH-spline rt - m= r 7 r 9 r 8 r r r 7 r 6 r r 0 r r r 6 r 4 r r 0 r r r 4 r r 9 r r 0 r 8 r 4 r Figure 8: Closed cubic PH B-Spline curves with m = left m = center m = right originated from closed degree- splines zt..4. Closed quintic PH B-Splines n = Let m N m. For a general not vectorµ={0=t 0 < t <...<t m } see Figure first row by applying the above method we construct the not partitions ν={ 0 < t <...< t m } ρ={t < 0 < t <...< t m < t m6 } τ={ t < 0 8 < t 8 <...< t m 8 < t m6 } 9

20 Deg- spline zt- m= z 0 z Deg- spline zt- m= Deg- spline zt- m= z 0 z z z z z 4 z z z 0 z Deg- PH-spline rt - m= Deg- PH-spline rt - m= Deg- PH-spline rt - m= r 7 r 6 r r 0 r r r 6 r r 6 r 7 r r 0 r r 9 r 8 r r r 0 r r r 4 r r 7 r r 9 r r 0 r 4 r 8 r 4 Figure 9: Closed cubic PH B-Spline curves with m = left m = center m = right originated from open degree- splines zt. n= m= n= m= n= m= Figure 0: Offsets of closed cubic PH B-Spline curves with: m = left m = center m = right. illustrated in Figure. Then by solving the linear systems 0 we calculate the coefficientsχ i j 0 i j m 0 m. All of them turn out to be zero with the exception of χ 00 = d d d χ χ =χ = 6 =χ = 6 χ = χ =χ = 6 d d d d d d d d d d d d d d d d d =... m =... m =... m d d d d d d =... m 0

21 where d = t t =... m. χ = d d d =... m χ =χ = d d d =... m χ =χ = d d d =... m χ = d d d =... m χ mm d m0 = m d m4 d m d d m d m d m4 d m 0 t... t m t m t m t m4 t m d d m d m d m4 d m 0 t... t m t m t m t m4 t m d 0 d d m d m d m4 d m d m6 t 0 t... t m t m t m t m4 t m t m6 d 0 d d m d m d m4 d m d m6 t 0 t... t m t m t m t m4 t m t m6 Figure : Knot partitions for the closed case n =. From top to bottom: µ ν ρ τ. In addition we compute the coefficientsζ i j 0 i m 0 j m 0 8m47 as the solutions to the linear systems 46. All of them turn out to be zero with the exception of ζ 00 = d 0 ζ 00 6D 6 = d 0 4D 0 0 ζ 0 6 = d 0 D 0 ζ 00 7 = 0 ζ0 7 = d 0 d 0d 0 6D 0 D 6D 0 ζ 0 7 = d 0 d 6D 0 D ζ 0 8 = d 0d 6D 0 D 0 6 ζ0 8 = d 0d 6D 0 D ζ 0 8 = 0 ζ 0 9 = d 4D ζ 0 9 = d 4D ζ 0 9 = d 4D ζ 9 = 0 ζ0 9 = d 4D ζ 0 0 = d ζ 6D 0 = d 8D ζ 0 = d 8D ζ 0 0 = d d ζ D 0 = d 9D ζ 0 0 = d 6D ζ 8 = d ζ D 8 = d d ζ 9D 8 = d d ζ 9D 8 = d ζ 9D 8 = 4d =0... m 9D ζ 8 = 4d ζ 9D 8 = d ζ 9D 8 = d d ζ 9D 8 = d d ζ 9D 8 = d =0... m D ζ 8 = d ζ 6D 8 = d 9D ζ 8 = d d ζ D 8 = d 8D ζ 8 = d 8D ζ 8 = d =0... m 6D ζ 4 84 = d ζ D 84 = d 4D ζ 4 84 = d d ζ D 84 = 0 ζ 84 = d 4D ζ 4 84 = d ζ 4 D 84 = d 4D ζ 4 84 = d 4D =0... m ζ 4 8 = d d ζ d 6D D 8 = ζ 4 6D D 8 = d d d 0d 6D D 6D ζ 8 = d d 6D D ζ 8 = 0 ζ4 8 = d d 0d 6D D 6D ζ 8 = d d ζ 4 6D D 8 = d d 6D D ζ 4 8 = d d 6D D 0 6 ζ d 8 = 6D D ζ 8 = d d 6D D =0... m 9

22 ζ 4 86 = d d 6D D ζ d 86 = 6D D ζ 4 86 = d d 6D D 0 6 ζ 86 = d d 6D D ζ 4 86 = d d ζ 4 6D D 86 = d d 0d 6D D 6D ζ 4 86 = 0 ζ 86 = d d ζ 6D D 86 = d d d 0d 6D D 6D ζ 6 86 = d 6D D ζ 4 87 = d 4D ζ 87 = d 4D ζ 4 87 = d ζ 4 D 87 = d 4D ζ = 0 ζ 87 = d d ζ D 87 = d 4D ζ 6 87 = d =0... m D ζ 4 88 = d ζ 44 6D 88 = d 8D ζ 4 88 = d 8D ζ 88 = d d D ζ 4 88 = d 9D ζ 6 88 = d =0... m 6D ζ m0m 8m7 = d m ζ mm 6D 8m7 = d m 9D m4 m4 ζ mm 8m7 = d m4d m D m4 ζ 6 86 = d d =0... m 6D D ζ mm0 8m7 = d m 8D m4 ζ mm 8m7 = d m4 8D m4 ζ mm 8m7 = d m4 6D m4 ζ mm 8m8 = d m 4D m4 ζ mm 8m8 = 0 ζmm 8m8 = d m4 4D m4 ζ mm0 8m8 = d m 4D m4 ζ mm 8m8 = d m4 4D m4 ζ mm 8m9 = 0 ζmm0 8m9 = d md m6 6D m4 D m ζ mm 8m9 = d m4d m6 6D m4 D m 0 6 ζ mm0 8m40 = d mdm6 ζ mm 6D m4 D 8m40 = d m4dm6 0d m6 6D m m4 D 6D m m ζ mm 8m40 = 0 ζ mm 8m4 = d m6 ζ mm D 8m4 = d m6 4D m m ζ mm 8m4 = d m6 6D m where D := d d = 0... m. By means of the computed coefficients{χ i j j m }0 i 0 m we can thus shortly write the control points of r t as p 0 = p = 0 p = d d d z 0 p = z d z d z d z d z d d d d =... m d p = z z d z d d =... m d p = z z d z d d =... m d p m0 = m d m4 d m z m p m = p m = 0 and the coefficients of the parametric speedσt as σ 0 =σ = 0 σ = d d d z 0 z 0 d d σ = 6 d d d d z z z z 6 d d d d d d z z 6 d d d d d z z z z d d d d d d z z z z =... m d d d z z z z =... m σ = d d d z z σ = d d d z z z z d d d z z =... m d σ m0 = m d m4 d m z m z m σ m =σ m = 0.

23 Deg- spline zt- m= Deg- spline zt- m= z 0 z n= m= z z 0 z z 4 z z z 4 z Deg- PH-spline rt - m= r 8 r 9 r r r 0 r r r Deg- PH-spline rt - m= r 9 r 0 r r 0 r r r 0 r 8 r r 9 r 8 r 7 r 4 r r r 7 r 9 r 8 r 7 r 4 r r 6 r 7 r 6 r r 6 r r 4 r r 6 r r 4 Figure : Closed quintic PH B-Spline curve with m= originated from a closed/open degree- spline zt first/second column and offsets of the non self-intersecting curve with and without control polygon third column. Thus according to the closed quintic PH B-Spline curve defined over the not partitionρis given by with control points r = r = r 0 r = r d z 0 rt= m z i r i N iρ t t [t t m ] t 0 = 0 r i = r i d i di z i d i z i di z i d i z i d i d i d i d i i=... m r i = r i z i di z i d i z i i=... m r i = r i z i di z i d i z i i=... m r m = r m0 d m z m r m = r m = r m and arbitrary r 0. Note that due to condition 9 d m = d d m4 = d. Moreover z m z m z m must be suitably fixed in order to satisfy condition 8. According to 7 the total arc length of the closed PH B-Spline curve of degree is given by L=l m7 l 4 N 4ρ t l m8 l N ρ t l m9 l 6 N 6ρ t

24 where N 4ρ t N ρ t N 6ρ t = and l = l = l 0 = 0 l = l d z 0 z 0 l i = l i d i z i z i d i d i d i d i d i di d i z i z i d i zi z i z i z i d id i zi z i z i z i d i d i zi z i z i z i i=... m l i = l i di z i z i d i zi z i z i z i i=... m l i = l i di z i z i d i zi z i z i z i i=... m l m = l m0 d m z m z m l m = l m = l m. Over the not partitionτthe offset curve r h t has the rational B-Spline form r h t= 8m47 =0 q N 9 τ t 8m47 =0 γ N 9 τ t t [t t m ] 60 where by exploiting the explicit expressions of the coefficients{ζ i j 0 j m }0 i m 0 8m47 from 9 weights and control points can easily be obtained; their explicit formulae are reported in the Appendix. Some examples of closed PH B-Spline curves of degree are shown together with offsets in Figure. 6. A practical interpolation problem We are concerned with C continuous quintic PH B-Spline curves defined by setting n= and m=. According to Corollary starting from zt= z i Niµ t over the partitionµ={0 0 0 a } we obtain the PH B-Spline curve 8 rt= r i Niρ t t [0 ] 6 over the not vectorρ={ a a a } from. We will now use these curves in order to solve the following practical interpolation problem several variants of which have been of interest to the PH community up to date see e.g. [ 6 0 ]. Problem. Given arbitrary points p 0 p tangents d 0 d and curvature valuesκ 0 κ to be interpolated by a PH quintic B-Spline curve from 6 such that r0=p 0 r=p r 0=d 0 r =d κ0=κ 0 κ=κ we loo for the control points r 0 r... r 8 of the curve 6. According to this means that we have to determine z 0... z. The positional interpolation constraints clearly imply r 0 = p 0 and r 8=p. According to 4 8 the tangential interpolation constraints yield the following conditions. r 0=z 0=p 0 = z 0 = d 0 r =z =p 7 = z = d. 6 By applying de Moivre s theorem to the two equations in 6 we obtain the following solutions for z 0 and z : z 0 = ± d 0 exp i ω 0 =± d 0 ω0 ω0 cos i sin z = ± d exp i ω =± d ω ω cos i sin 6 whereω = argd for =0. 4

25 Note that we can limit ourselves to considering the combinations of z 0 z with signs and since the others will give rise to the same solutions. Next recalling the general formula for the curvature of a PH curve see e.g. [0] κt= Im zt z t zt 4 64 we can express the curvature constraints at t=0 and t= as κ0= 4 a I z 0z z 0 4 = 4 a u0v u v 0 u 0 v 0 =κ 0 κ= 4 a I z z z 4 = 4 a uv u v u v =κ. 6 If u 0 0 and u 0 we can derive v and v from 6 obtaining If v = u v 0 a u 0 4 κ 0u 0 v 0 v = u v a κ u u 4 v. 66 u 0 = 0 respectively u = 0 67 i.e. if d 0 respectively d has the form const. 0 we can not proceed as in 66. In this case by exploiting the affine invariance of B-Spline curves we propose to rotate the initial data in order to avoid the case 67. We can thus always assume u 0 0 and u 0. To determine the remaining unnowns u and u in 66 we consider the equation r 7 r = 6 r i with r i = r i r i. 68 i= Since r 0= a r r 0 =d 0 r = a r 8 r 7 =d 69 the left hand side of equation 68 is completely determined by the given Hermite data as: r 7 r = r 8 r 0 a d a d 0= p p 0 a d a d By the right hand side of 68 reads a z az z z a z z a4 az 0z a a z z a z 0 z. 7 By separating the real and imaginary parts of this equation we obtain the following system of two real quadratic equations in the two real unnowns u and u u u A u = 0 u u B u = 0 7 u u with real symmetric matrices A = a i j 0 i j and B = b i j 0 i j whose coefficients are reported in the Appendix. The two equations in 7 represent two conic sections in the real Euclidean plane. Finding the solutions u and u of these equations is equivalent to determining the intersection points of the two conic sections. In order to obtain them we consider the pencil of conic sections defined by the given conics as u u A λb u u = 0 λ R. 7 Among the one parameter set of conics there are three degenerate conics given by theλ-values obtained as solutions of the cubic equation inλ: deta λb= λ detbλ detb detb detb λdeta deta deta deta=0 74

26 where the matrices A for = are obtained by replacing the -th column of A by the -th column of B and the matrices B for = are obtained by replacing the -th column of B by the -th column of A. These degenerate conics are pairs of intersecting lines. The intersection points of the conics of the pencil 7 are then easily obtained in the following way. For a solution of the cubic equation 74 by inserting the correspondingλ-value into the equation 7 a quadratic equation easily decomposable into two linear factors is obtained. These two linear equations in u and u can be solved for e.g. u in dependency of u which are then inserted into another conic of the pencil for example one of the given conics or another degenerate one yielding a quadratic equation in u. We thus obtain the coordinates of the intersection points of the pencil conics. The control points of the corresponding PH B-Spline curves are obtained by inserting the obtained u i v i... into the expression for r i in. Letδ respectivelyδ be the number of real intersection points of the conics from 7 for the sign choice respectively in 6. Thus the total number of solutions of our interpolation problem isδ=δ δ. Among all solutions as it is classical seee.g. [ 9] the good curve is represented by the curve with minimum absolute rotation index and bending energy. The following figures show the cubic curve interpolating the assigned values and derivatives and the different quintic PH B-spline curves determined as the solution of the Hermite interpolation problem where the endpoint curvaturesκ 0 andκ are sampled from the cubic curve. The corresponding conics with coefficient matrices A and B from 7 and a degenerate conic of the pencil corresponding to a pair of lines are also displayed. 6.. Existence of solutions In the following we discuss the problem of studying the existence of a solution for Problem. Assigned a set of values tangents and curvatures there are two situations where a solution may not be found. The first corresponds to the case where either A or B or both in 7 yield imaginary conics. The second is the case where the two conics in 7 are real but have no intersection. The latter situation has never been encountered in our experiments. Thus we will focus on discussing the former problematic case. In particular we will consider the following problem. Suppose that points p 0 p and tangents d 0 d are given how shall we prescribe the curvaturesκ 0 κ in order to guarantee that A and B represent real conics? Type of conic I [M] I [M] I [M] Ellipse 0 > 0 I [M] I [M] < 0 Hyperbola 0 < 0 Parabola 0 = 0 Non-degenerate imaginary conic 0 > 0 I [M] I [M] > 0 Pair of intersecting real lines = 0 < 0 Pair of intersecting imaginary lines = 0 > 0 Pair of parallel real lines = 0 rani [M] = = 0 < 0 Pair of parallel imaginary lines = 0 rani [M] = = 0 > 0 Double line pair of coinciding real lines = 0 rani [M] = = 0 Table : Classification of conic sections. For a generic conicc M : m u m u u m u m 0u m 0 u m 00 = 0 with symmetric matrix M := m i j 0 i j m i j = m ji we suppose m 00 0 and we denote by M := m i j i j the matrix of the associated quadratic form. We then consider the following invariants: I [M] := trace M I [M] := det M I [M] := detm. 7 Then according to Table C M represents an imaginary conic in the following cases: P I [M] 0 I [M] > 0 I [M] I [M] > 0 C M is a non-degenerate imaginary conic; P I [M] = 0 I [M] > 0 C M is a pair of intersecting imaginary lines; P I [M] = 0 rani [M] = I [M] = 0 I [M] > 0. C M is a pair of parallel imaginary lines. 6

27 When choosingκ 0 κ we shall then avoid thatc A orc B fall into one of the points above. To this aim it is useful to observe that for M {A B} I [M] and I [M] depend on the location of not a only whereas I [M] depends on aκ 0 κ see the formulae of the coefficients of A and B in the Appendix. For a fixed value of a the sign of I [M] and I [M] is thus determined whereas I [M] = 0 is a quadratic equation in the two unnownsκ 0 κ. The general idea would thus be to first choose the value of a and based on the corresponding signs of I [M] and I [M] fixκ 0 κ in order to prevent that M falls into one of the items above. We will illustrate the idea by means of the four test cases contained in Examples to 4. This is not an exhaustive discussion of all the possible forms of the two conicsc A andc B but it suffices to illustrate the general criterion for choosing suitable values of aκ 0 andκ. The general procedure is summarized in Algorithm. Algorithm Algorithm for the solution of an Hermite interpolation problem with given interpolation points p 0 p and tangents at these points d 0 d. Input: interpolation points p 0 p and tangents at these points d 0 d. From 6 and 6 determine all possible values z 0 and z i.e. u 0 v 0 u v and for each combination with signs or in 6 perform the following steps:. compute the matrices A and B of the two conics in 7. These matrices will depend on the parameters a κ 0 andκ ;. choose the value of the not a 0 and compute I [M] I [M] for M = A B. Since these quantities depend on a only their sign is determined;. Chooseκ 0 andκ in such a way that bothc A andc B are real. In particularc A orc B or both may represent an imaginary conic when I [A] > 0 or I [B] > 0 or both. I [M] = 0 M=A B is a quadratic function ofκ 0 andκ and represents a conic in the Euclidean plane. Plotting the graph of I [M] = 0 we can chooseκ 0 andκ in the region where the relative signs of I [M] I [M] and I [M] yield a real conic. 4. Solve the system 7 for u and u.. For each solution of the system 7 compute v and v from 66. The obtained values u i v i... completely determine a set of complex coefficients z i i=... and thus identify one PH B-spline curve. Output: PH B-spline curves corresponding to the given initial data. Example. Let us consider the initial data p 0 = 0 p = 0. d 0= d = 0.. Using these data from 6 we can compute z 0 and z namely u 0 v 0 u v. For any combination of z 0 and z to be considered i.e. taing signs and in 6 we should proceed to computing the intersection of the two conics in equation 7. Each intersection point ofc A andc B will give rise to a solution for u v u v. The first couple of values corresponds to signs. In this case any choice of a in 0 yields I [A] > 0 I [A] < 0 I [B] > 0 I [B] < 0. Hence bothc A andc B are real conics independently on the choice of aκ 0 andκ. In particular we will set a=0. and sample the curvature values from the cubic polynomial interpolating p 0 p and d 0 d which results inκ 0 =.0409 andκ = The two conicsc A andc B are displayed in the left part of Figure while the center and right figures show the PH B-spline curves corresponding to the two intersection points of the conics. The second couple of values corresponds to signs. In this setting I [M] I [M] M {A B} have the same signs as for the previous case. Choosing the same values for aκ 0 andκ C A andc B turn out to be the two hyperbolas plotted in the left part of Figure 4 whereas the PH B-spline curves corresponding to each of their intersection points are shown in the center and right part of the figure. Based on the absolute rotation index and bending energy displayed in Figures and 4 the best obtained PH B-spline curve turns out to be the one in the center of Figure. Example. Let us consider the initial data p 0 = 6 p = 0 d 0= 0 d = 0. Using these data from 6 we can compute z 0 and z namely u 0 v 0 u v. For any combination of z 0 and z to be considered 7