Arxiv:1506.00967V2 [Cs.GR] 6 Jan 2016 N Lsiyn T Hsslto Ssoni 1]T Eone Be to [11] Elsevier in to Submitted Shown Preprint Is Solution This Quadric It

Home , Oval (projective plane)

arXiv:1506.00967v2 [cs.GR] 6 Jan 2016 n lsiyn t hsslto ssoni 1]t eone be to [11] Elsevier in to submitted shown Preprint is solution This quadric it. a classifying to check- belongs and for Bézier triangle rational algorithm rational a an whether of provides ing classification [10] thorough triangles. a classifies ge- quadratic provides their [8] [9] studies and ometry. to quadrics. approach surfaces on algebraic parametrised surfaces an quadratically and use [7] curves Bézier con- and construct quadratic [6] a polyhedron. construct- on for trol grounded method patches a quadric describes [5] ing view. are of points patches many quadric cases some in [4]. only obtained quartic and general hand, [3] in other surfaces are the Bézier triangles as On quadratic represented rational Bézier triangle. be hand, quadratic can one rational patch On a triangular surfaces. quadric to every go not we when situation different topologies singularities. different producing to without for adapting have quadrilaterals, and to surfaces latter compared constructing to the triangles, due of but theory flexibility element common, the finite and most animation in the applications are The ones [2]. Bézierformer triangles and patches product tensor CAD: engi- [1]. in architecture used in commonly and spheres, are neering as quadrics such paraboloids of surfaces cylinders, quadric inclusion cones, since the fashion, polyno- is exact piecewise ones an in or rational to polynomial surfaces from mial CAD in used faces Introduction 1. framework t this relevant within Keywords: also characterised is are coefficient inv revolution This coefficient, of one quadrics. quadrics in ruled and encoded from that is oval quadrics quadric apart of the tell geom pencil about information algebraic a the considering projective of fashion, from resu coordinate-free techniques The o a using axis patch. in on the quadric rely and the results vertex for weights the Our infinity, the ele at and geometric conic vertices for the formulas control center, closed the derive and as classify (such we paper this In Abstract udi raglrpthshv ensuidfrom studied been have patches triangular Quadric a have we conics, are curves quadratic rational Whereas in surfaces implementing of ways two mainly are There sur- of framework the extending for reasons the of One emti lmnsadcasfiaino udisi ration in quadrics of classification and elements Geometric leri rjciegoer,rtoa ´ze ace,quadr Bézier patches, rational geometry, projective Algebraic .Cnń .Frńe-abia .Rsd aá ..Vá Mar´ıa, M.J. Rosado Fernández-Jambrina, Cantón, E. L. A. nvria oićiad Madrid Politécnica de Universidad -84-ard Spain E-28040-Madrid, Matemática Aplicada hc sncsayfriscasfiaini opie na in comprised is classification its surface the for fact, of necessary In information is geometric which relevant of step. the use intermediate of geo- the an most only avoiding as of patch, expressions use the coordinate the about is originality information formulation major metric The our patch. of the advantage spanned of and plane vertices the corner and of the closed terms planes by in tangent a quadrics their find these weights, of to their equation in- implicit us for the allows for of, form instead coordinates, cartesian ) . . stance, conics. Using boundary 17]. planes, containing (tangent 15, meaning planes [14, geometric clear in a with used forms of been in linear pencils also of have form use which the tangential from quadrics, as such borrowed and geometry, techniques projective point algebraic using in fashion, both coordinate-free quadrics, designing for the help the great of of One be [16]. can them. geometric in quadric with for conics a formulas closed for of that done elements is weights this is the doing it and for as net reasons control patch, the elements the just geometric of using their form, calculate closed patches. and in product them tensor classify We and triangles Bézier quadric nal vectors normal the the [15] and In them. points at patches. three surface are param- quadric parameters shape design shape as to used [14] are of ratio- in weights eters axes Three constructing three the and for quadrics. points obtain non-degenerate tool corner on to a patches used provides quadratic is [13] nal algorithm quadric. the [12] the condition- In numerical different ing. with possibilities, fifteen of t r xeddas oqarctno rdc patches. product tensor quadric to also extended are lts ihti oli id edrv iierfrsfor forms bilinear derive we mind, in goal this With ratio- from information geometric draw we paper this In aaood n h rnia lns,uigjs the just using planes), principal the and paraboloids f r agn otegvnqarcaogacnc Most conic. a along quadric given the to tangent are . et fqarc nrtoa ´ze raglrform Bézier triangular rational in quadrics of ments eemn h ffietp fteqarc Spheres quadric. the of type affine the determine o tyt n utbeblna om o h quadric the for forms bilinear suitable find to etry ligtewihso h ac,wihalw sto us allows which patch, the of weights the olving c,Senrsurfaces. Steiner ics, zquez-Gallo lBézier form al ue3,2018 30, June single coefficient, the parameter of the pencil of quadrics, such that u + v + w = 1, u,v,w [0, 1], where the coeffi- ∈ which depends on the weights of the patch. cients cijk are the control points for the surface and ωijk The coordinate-free bilinear forms for the quadrics en- are their corrresponding weights. able the derivation of closed formulas for several geomet- The surface patch is bounded by three rational Bézier ric elements (center of the quadrics, bilinear forms for the curves of degree n, which are obtained by fixing u = 0, conic at infinity, diametral planes) calculated in terms of v = 0, w = 0. Their respective control points and weights weights and vertices. To our knowledge, such closed for- are respectively the ones with index i = 0, j = 0, k = 0. mulas have not been produced before. Similarly, linear We are interested in the special case of quadratic sur- forms for principal planes are obtained up to solving a faces. In this case the control net and the matrix of weights cubic equation. The degeneracy of the solutions of this are equation allows us to identify quadrics of revolution and spheres. c002 c011 c020 ω002 ω011 ω020 Finally, the previous results for triangular quadric c c , ω ω , (2)  101 110   101 110  patches can be extended to rational biquadratic quadric c200 ω200 patches using the tangent planes at three corners of the     boundary of the tensor product patch. and the surface patch is completely determined by the This paper is organised as follows: We revisit ratio- three boundary conic curves: nal Bézier triangular patches and the characterisation of The conic at u = 0 has control points c002,c011,c020 quadric patches in Section 2. In Section 3 we construct the and weights ω ,ω ,ω , the one at{v = 0 has con-} { 002 011 020} pencil of quadrics, in point and in tangential form, which trol points c002,c101,c200 and weights ω002,ω101,ω200 , are tangent to our quadric at the conic through the three and the one{ at w = 0 has} control points{ c ,c ,c } { 020 110 200} corner vertices of the patch and fix the only free parameter and weights ω020,ω110,ω200 . We just consider the case in terms of the weights of the quadric. This provides us bi- of non-degenerate{ boundary} curves. We name the planes linear forms for the quadrics and hence implicit equations. where such conics are located respectively as u,v,w. A coordinate-free expression for the center of the quadric Quadratic rational Bézier triangular surfaces comprise is obtained in Section 4 as a barycentric combination of the quadrics as a subcase, but in general they are quartic corner vertices of the patch and the intersection of their surfaces named Steiner surfaces [3]. On the contrary, tangent planes. It is used to tell paraboloids from centered not every quadric patch bounded by three arcs can be quadrics. In Section 5 we calculate a bilinear form for the parametrised as a quadratic rational Bézier triangle. Ac- conic at infinity, which is useful for determining the affine cording to [4]: type of the quadric. Section 6 is devoted to degenerate quadrics. With this information we provide in Section 7 a If the Steiner surface is a non-degenerate quadric, the way to classify quadric patches using the signature of the • three conic boundary curves meet at a point S and bilinear form, the center and the bounding conic arcs. We their respective tangent vectors at S define a plane. introduce the scalar product in Section 8 in order to calculate Euclidean geometric elements such as principal planes If the three conic boundary curves meet at a point and axes of quadrics and vertices of paraboloids. This al- • S and their respective tangent vectors at S define a lows us to characterise spheres and quadrics of revolution. plane, the Steiner surface is a quadric patch. In Section 9 we show several examples of application of our results. We show how to extend our results to quadric If there are three alligned points, one on each bound- tensor product patches in Section 10. A final section of • ary conic curve, and the tangent plane to the surface conclusions is included. is the same for all three points, the Steiner surface is a degenerate quadric patch. 2. Quadratic rational Bézier triangular patches The last case just implies that the quadric (cylinder or cone) is ruled and the three points lie on the same ruling. Rational Bézier triangular patches of degree n are defined using trivariate Bernstein polynomials of degree n in We are considering Steiner surfaces with such a point parameters u,v,w, S. This case comprises non-degenerate quadrics. Since S belongs to the three boundary conics, there are values n! us, vs, ws of the parameters such that Bn (u,v,w)= uivj wk, i + j + k = n, ijk i!j!k! 2 2 2 2 n ω0j2 j c0j2 j Bj (vs) ωi2 i0ci2 i0Bi (us) ωijkcijk Bijk(u,v,w) − − − − j=0 i=0 i+j+k=n S = X = X c(u,v,w)= X , (1) 2 2 n 2 2 ωijkBijk(u,v,w) ω0j2 j Bj (vs) ωi2 i0Bi (us) − − i+j+k=n j=0 i=0 X X X 2 C U

p Q u t P q w v W Figure 1: Tangent cone to a quadric along a conic on a plane t V R r

2 Figure 2: Conic circumscribed by a triangle 2 ω2 k0kc2 k0kBk(ws) − − k=0 = , c110 X 2 200 2 c ω2 k0kBk(ws) w=0 − r k=0 X c101 q c020 where we have introduced quadratic Bernstein polynomials, v=0 p u=0 c011 B2(t)=(1 t)2, B2(t)=2t(1 t), B2(t)= t2. 0 − 1 − 2 However, since the set of weights of rational Bézier c002 curves is unique up to Möbius transformations [18] of the Figure 3: Tangent planes p,q,r parameter, we can use this freedom to set us, vs, ws at infinity. It is a simple exercise to check that this is possible if and only if the tangent vectors at S are coplanar [19]. For instance, we can take the cone C tangent to our Hence, this choice is not a restriction. quadric along the conic on t and the double plane t (see From now on we assume that the set of weights for the Fig. 1). Hence, the bilinear form for this pencil of quadrics surface fulfills this condition and therefore we may write is just C + λt2, using the same letter for a surface and its S in three different forms as form. Finally, we can determine our quadric just using the ω002c002 2ω011c011 + ω020c020 existence of point S, S = − ωu ω c 2ω c + ω c C(S)+ λt(S)2 =0. = 002 002 − 101 101 200 200 ω v We just have then to provide a bilinear form for the cone ω c 2ω c + ω c = 200 200 − 110 110 020 020 , (3) C. Since λ = C(S)/t(S)2, we notice that the sign of λ ω w depends on whether− S lies in or out of the cone C. where ωu,ωv,ωw guarantee that these are barycentric com- It will be useful to start with the bilinear form for the binations, conic on t (see Fig. 2). In this figure u,v,w are the straight lines which are respectively the intersections of the planes ωu = ω 2ω + ω , ωv = ω 2ω + ω , 002 − 011 020 002 − 101 200 u,v,w (where the boundary conics lie) with the plane t. ωw = ω200 2ω110 + ω020, We use the same letter for a point (uppercase) and its − polar line (lower case). If a point A lies on the conic, if S is a proper point. If S is point at infinity, the coeffi- its polar line is the tangent a to the conic at such point. cients ωu,ωv,ωw are chosen in order to have the same vec- Hence, p, q, r are respectively the tangent lines to the conic tor as representantive for S on the three boundary curves. at P , Q, R (see Fig. 3). We shall use the same letters p,q,r also for the tangent planes at P,Q,R. 3. Pencils of Steiner quadrics On the other hand, the polar line to a point A not on the conic is the line a linking the tangency points of the In order to obtain a bilinear form for the quadric, we tangent lines to the conic drawn from A. Hence, u, v, describe the pencil of quadrics which are tangent to our w are the polar lines to the vertices of the circumscribed quadric at the conic on the plane t defined by the points triangle U, V , W . at the corners of the surface patch, P = c002, Q = c020, The pencil of conics circumscribed by the triangle de- R = c200. This pencil of quadrics has been used also in fined by the lines p, q, r is easily described in tangential [14, 15, 17]. form. The bilinear form for this pencil is just To this aim, we just need two quadrics belonging to the pencil [20], which can be degenerate. αUV + βUW + γV W, 3 with coefficients α, β, γ. It is easy to check that this form Since we require that points P , Q, R lie on the conic, vanishes on p, q, r, since βq(P )= γr(P ), αp(Q)= γr(Q), αp(R)= βq(R), p(U)=0= p(V ), q(U)=0= q(W ), r(V )=0= r(W ). we identify the coefficients of the bilinear form, But, for our purposes, we require the point bilinear form 1 1 1 for the conic. Referred to the lines r, q, p, which form a α = , β = , γ = , dual reference to the one formed by U, V , W , provided ΩW ΩV ΩU that their linear forms satisfy up to a common factor. p(W )=1, q(V )=1, r(U)=1, (4) Hence, a bilinear form for the conic on t is the matrix of this point bilinear form is the inverse of the p2 q2 r2 2pq 2pr 2qr 2 + 2 + 2 (8) one for the tangential bilinear form, ΩW ΩV ΩU − ΩW ΩV − ΩW ΩU − ΩV ΩU 1 0 α β − γ2 βγ αγ and, in tangential form, 1 α 0 γ = βγ− β2 −αβ , 2αβγ  β γ 0  −  −αγ αβ− α2  ΩU ΩV UV +ΩU ΩW UW +ΩV ΩW V W. (9) − − and hence such a bilinear form is  The terms of the bilinear form can be factored as 2 2 2 α2p2 2αβpq 2αγpr + β2q2 2βγqr + γ2r2. p q r q r q r − − − + + , Ω − Ω − Ω Ω − Ω − Ω Ω In order to determine the coefficients α,β,γ, it is useful W V U V U V U to take into account that (3) provides relations between and thereby it has signature (+ + ) regardless of the values of the linear forms for p,q,r, values of the coefficients. − ω020 ω200 We move now back from plane to space: If p, q, r are the p(Q)= p(S)= p(R), ωu ωv tangent planes to the quadric at P , Q, R, the bilinear form ω200 ω002 (8) describes a degenerate quadric of signature (+ + 0) q(R)= q(S)= q(P ), − ωw ωu which contains the conic on t and it is tangent to our ω ω 002 r(P )= r(S)= 020 r(Q), quadric along it. Since the rank of the bilinear form is ωv ωw three, it is a cone (or a cylinder, if the intersection point which can be used to write U,V,W in terms of P,Q,R, of p, q, r goes to infinity), since these are the only non- plane degenerate quadrics. Hence, it is the cone we are ω ωwP + ω ωvQ ω ωuR U = 002 020 − 200 , looking for. ΩU The bilinear form for the pencil of quadrics referred to ω ωwP ω ωvQ + ω ωuR V = 002 − 020 200 , the planes r,q,p,t, ΩV p2 q2 r2 2pq 2pr 2qr ω002ωwP + ω020ωvQ + ω200ωuR 2 2 + 2 + 2 + λt , W = − , (5) Ω Ω Ω − ΩW ΩV − ΩW ΩU − ΩV ΩU ΩW W V U or, conversely, produces a bilinear form for our quadric if λ = C(S)/t(S)2. ΩU U +ΩV V ΩW W +ΩU U ΩV V +ΩW W P = ,Q = , R = , − We see that the sign of λ determines the signature of the 2ω ω 2ω ω 2ω ω 002 w 020 v 200 u bilinear form, defined as the difference between the number where the normalisation denominators are of its positive and negative eigenvalues. This is relevant since ruled quadrics have null signature, (+ + ) and ΩU = ω002ωw + ω020ωv ω200ωu, − oval quadrics have signature two (+ + +). This−− is useful ΩV = ω002ωw ω020ωv + ω200ωu, − − to classify the quadric. ΩW = ω ωw + ω ωv + ω ωu, (6) − 002 020 200 In tangential form referred to the points U,V,W,T , the provided that U,V,W are proper points. If one of these bilinear form for the quadric is points goes to infinity, the respective denominator is taken T 2 to be one. The interpretation of these coefficients as ΩU ΩV UV +ΩU ΩW UW +ΩV ΩW V W , − λ weights for U,V,W on the conic on t is discussed in Sec- tion 6. if t(T ) = 1, where T is the intersection point of the planes Finally, the normalisation condition (4) allows us to fix p,q,r, that is, the vertex of the cone C. the linear forms for p,q,r since With this information we can write T as

2ω ω 2ω ω ω002ωw(ωu + ωv ωw)P + ω020ωv(ωu ωv + ωw)Q p(W )= u v p(S), q(V )= u w q(S), T = − − Ω Ω ǫU ΩU ωu + ǫV ΩV ωv + ǫW ΩW ωw 2ǫSωuωvωw W V − 2ωvωw ω ωu( ωu + ωv + ωw)R 2ωuωvωwS r(U)= r(S). (7) + 200 − − (10), ΩU ǫU ΩU ωu + ǫV ΩV ωv + ǫW ΩW ωw 2ǫSωuωvωw − 4 or alternatively, 4. Center of a quadric Ω Ω Ω S = U U + V V + W W + 2ωvωw 2ωuωw 2ωuωv The bilinear form for the quadric provides a way to ob- 2ǫSωuωvωw ǫU ΩU ωu ǫV ΩV ωv ǫW ΩW ωw tain its center Z as the pole of the plane z at infinity. Since + − − − T, 2ωuωvωw the plane at infinity is formed by vectors, we may write its elements as barycentric combinations where ǫA takes the value one if the corresponding point A is proper or zero if it is a point at infinity. aU + bV + cW + dT, a + b + c + d =0. Now we can compute λ and determine the quadric: Theorem 1. A rational triangular quadratic patch with We have to consider the possibility of having any of the control net c ,c ,c ,c ,c ,c and weights points of the reference at infinity. In such case, the null { 002 011 020 101 110 200} ω002,ω011,ω020,ω101,ω110,ω200 , such that the three sum is restricted to proper points, {boundary conics meet at a point} S, which is written as in (3), is a quadric with bilinear form aU +bV +cW +dT, ǫU a+ǫV b+ǫW c+ǫT d =0.(14) 2 2 2 p q r 2pq 2pr 2qr 2 Hence, a linear form for the plane at infinity in this = 2 + 2 + 2 +λt , (11) Q ΩW ΩV ΩU −ΩW ΩV −ΩW ΩU −ΩV ΩU reference is just where t is the linear form of the plane containing c002,c020,c200, p is the linear form of the plane containing z = ǫW p + ǫV q + ǫU r + ǫT t. (15) c002, c011, c101, q is the linear form of the plane containing c020, c011, c110 and r is the linear form of the plane The pole of this plane, containing c200, c101, c110 which satisfy ˜(z) = ǫW ΩW (ΩU U +ΩV V )+ ǫV ΩV (ΩW W +ΩU U) p(W )=1, q(V )=1, r(U)=1, t(T )=1, Q 2ǫT + ǫU ΩU (ΩV V +ΩW W ) T, (16) at the points U (intersection of the planes p,q,t), V (in- − λ tersection of the planes p,r,t) and W (intersection of the can be written in a simpler way in terms of P,Q,R,T planes q,r,t) given by (5) and the point T (intersection of in order to produce an expression for the center of the the planes p,q,r), given by (10). quadric, The coefficients are given in (6) except for λ, 2 2 2 ǫ Ω ω ω P + ǫ Ω ω ω Q + ǫ Ω ω ω R ωu + ωv + ωw 2ωuωv 2ωvωw 2ωwωu Z = W W w 002 V V v 020 U U u 200 λ = − − − 2 , (12) −(2ǫSωuωvωw ǫU ΩU ωu ǫV ΩV ωv ǫW ΩW ωw) ΩZ − − − ǫT and the tangential bilinear form for the quadric is T, (17) − λΩZ T 2 ˜ =ΩU ΩV UV +ΩU ΩW UW +ΩV ΩW V W , (13) where the denominator Ω is one if the center is a point Q − λ Z at infinity or, if it is a proper point, Furthermore, for λ> 0 the quadric is oval and for λ< 0 the quadric is ruled. ǫT ΩZ = Ω˜ W ωwω + Ω˜ V ωvω + Ω˜ U ωuω , (18) 002 020 200 − λ The expression we have obtained for the bilinear form of the quadric has the advantage of encoding most of the and we have introduced for simplicity, information about the surface in just the coefficient λ for 2 t . Ω˜ U = ǫU ΩU , Ω˜ V = ǫV ΩV , Ω˜ W = ǫW ΩW . (19) This result provides a procedure for computing a bilinear form, and hence the implicit equation, for a non- Since the center of a paraboloid is a point at infinity, we degenerate Steiner quadric patch in a coordinate-free fash- have a simple characterisation: ion using just the vertices of the control net and their respective weights: Corollary 1. A rational triangular quadratic patch with 1. Obtain S as intersection of the planes u,v,w. control net c002,c011,c020,c101,c110,c200 and weights ω ,ω ,ω{ ,ω ,ω ,ω , such that} the three 2. Compute the normalised linear forms for the planes { 002 011 020 101 110 200} t,p,q,r. boundary conics meet at a point S, which is written as 3. Obtain an equivalent list of weights fulfilling (3). in (3), is a paraboloid if the quadric is non-degenerate and

4. Use Theorem 1 to obtain the bilinear form for the ǫT Q Ω˜ W ωwω + Ω˜ V ωvω + Ω˜ U ωuω =0, quadric patch. 002 020 200 − λ 5. The implicit equation for the quadric patch is then (X,X) = 0. where the coefficients are given in (6), (12) and (19). Q 5 5. The conic at infinity an ellipsoid then, since it does not intersect the plane at infinity. The conic at infinity is the intersection of the quadric On the other hand, if the determinant (21) is negative, with the plane z at infinity and it is formed by its asymp- the signature of the bilinear form is (+ + ) and we have totic directions. It is useful for classifying quadric patches, a proper conic. The quadric is a hyperboloid− in this case: as it is done in [10]. Since points X on z satisfy Corollary 2. A rational triangular quadratic patch with z(X)= ǫ p(X)+ ǫ q(X)+ ǫ r(X)+ ǫ t(X)=0, W V U T control net c ,c ,c ,c ,c ,c and weights { 002 011 020 101 110 200} we can use as bilinear form for the conic at infinity on z ω002,ω011,ω020,ω101,ω110,ω200 , such that the three {boundary conics meet at a point}S, which is written as in 2 2 2 p q r 2pq 2pr 2qr (3) and with positive λ, given by (12), is an oval quadric = 2 + 2 + 2 Z ΩW ΩV ΩU − ΩW ΩV − ΩW ΩU − ΩV ΩU and: + λ(ǫ p + ǫ q + ǫ r)2, (20) W V U An ellipsoid if det > 0. • Z except when T is a point at infinity. An elliptic paraboloid if det =0. In order to draw information about the conic at infinity, • Z we may factor its bilinear form, A two-sheeted hyperboloid if det < 0. • Z λ = (Ω˜ 2 λ + 1)σ2 + σ2 where the expression for det is given in (21). W 1 ˜ 2 2 Z Z ΩW λ +1 1 2 6. Conics on planes t,u,v,w and degenerate + Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U σ , − λ 3 quadrics where we have introduced three linear forms, σ1, σ2, σ3, in In order to draw more information about the surface order to diagonalise , Z patch, we take a look at the conic arcs on planes t,u,v,w. A conic with weights ω ,ω ,ω can be classified [18] us- p Ω˜ V Ω˜ W λ 1 q Ω˜ U Ω˜ W λ 1 r 0 1 2 σ1 := + − + − , ing the canonical weight ω = ω1/√ω0ω2: ΩW Ω˜ 2 λ +1 ΩV Ω˜ 2 λ +1 ΩU W W If ω < 1, it is an ellipse; if ω = 1, it is a parabola and if ω > 1, it is a hyperbola. q σ2 := (Ω˜ V + Ω˜ W ) The canonical weights for the conics on planes u,v,w ΩV are respectively 2 Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U Ω˜ 2/λ r + − W − , ω011 ω101 ω110 ˜ ˜ , , . ΩV + ΩW ΩU √ω002ω020 √ω002ω200 √ω020ω200 2 r A set of weights for an arc of the conic on t is readily ob- σ3 := . Ω˜ V + Ω˜ W ΩU tained. For instance, we can use P,U,Q as control polygon and use the same kind of construction of (3) to assign an The case of T at infinity is simpler, as points at infinity infinite parameter to the point R, so that its coordinates satisfy referred to P,U,Q, z(X)= ǫ p(X)+ ǫ q(X)+ ǫ r(X)=0, W V U ω ωwP ΩU U + ω ωvQ R = 002 − 020 , but it can be handled similarly. ω200ωu Combining both cases, we obtain a general expression provides us a set of weights ω002ωw, ΩU /2, ω020ωv for the for the determinant of in this reference, Z conic arc. Similarly, we obtain weights ω002ωw, ΩV /2, ω200ωu for det = λ Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U ǫT . (21) Z − the arc with control polygon P,V,R and ω020ωv, ΩW /2, ω ω for the arc with control polygon Q, W, R. Hence The conic at infinity of a paraboloid is degenerate. 200 u the canonical weights for these arcs are Hence, det vanishes for these quadrics. This condition Z is equivalent to the one obtained in the previous section. ΩU ΩV ωP UQ = , ωPVR = , The classification of the conic at infinity allows us to 2√ω002ω020ωvωw 2√ω002ω200ωuωw finish the classification of quadrics. Since paraboloids are ΩW non-centered, one-sheeted hyperboloids are ruled and cen- ωQW R = . (22) 2√ω ω ωuωv tered, we just have to tell ellipsoids from two-sheeted hy- 020 200 perboloids, since they are both oval and centered: We can use either of these to classify the conic on t. Since λ is positive in this case, we notice that if the Furthermore, this result furnishes an interpretation of determinant (21) is positive, the signature of the bilinear ΩU , ΩV , ΩW as weights for the points U,V,W if ω002ωw, form is (+++) and hence the conic is imaginary. We have ω020ωv and ω200ωu are respectively the weights for P,Q,R. 6 We have calculated this set of weights resorting to the point S, but it is clear that it can be obtained indepen- dently from the barycentric combinations v

ΩQQ ΩW W +ΩRR ΩP P ΩV V +ΩRR P = − ,Q = − , ΩQ ΩW +ΩR ΩP ΩV +ΩR − − ΩP P ΩU U +ΩQQ R = − , (23) ΩP ΩU +ΩQ − up to a multiplicative factor. This is useful for degenerate quadrics (cones and cylinders), since for their triangular quadric patches the bound- Figure 4: The pole of a diametral plane is the direction of the tangent ary conics do not meet in general at a point S, but we can cylinder to the quadric along its intersection. still use (8) as the bilinear form for the tangent cone to the quadric along the conic on t: – Centered: One-sheeted hyperboloids. Theorem 2. A rational triangular quadratic patch with – Non-centered: Hyperbolic paraboloids. control net c002,c011,c020,c101,c110,c200 and weights ω ,ω ,ω{ ,ω ,ω ,ω , on a degenerate} quadric We may tell ellipsoids from two-sheeted hyperboloids in { 002 011 020 101 110 200} has a bilinear form with coefficients ΩU , ΩV , ΩW satisfying other ways. For instance, if the conics at planes t,u,v,w (23), are all ellipses, the quadric is an ellipsoid. We can use (22) for this. p2 q2 r2 2pq 2pr 2qr = 2 + 2 + 2 , (24) Q ΩW ΩV ΩU −ΩW ΩV −ΩW ΩU −ΩV ΩU 8. Diametral planes and axes with coefficients given by (6) and where p is the linear form of the plane containing c002, c011, c101, q is the linear form If a plane contains the center of the quadric, it is called of the plane containing c020, c011, c110 and r is the linear diametral. As the center of the quadric is the pole of the form of the plane containing c200, c101, c110 which satisfy plane at infinity, polar planes of points ~v at infinity are p(W )=1, q(V )=1, r(U)=1, diametral. That is, the tangent cone to the quadric along its intersection with a diametral plane degenerates to a at the points U (intersection of the planes p,q,t), V (in- cylinder (see Fig. 4). The direction of the cylinder is given tersection of the planes p,r,t) and W (intersection of the by the pole ~v of the diametral plane. planes q,r,t) given by (5). We choose a basis of vectors U,~ V,~ W~ , where If the planes p, q, r meet at a proper point T , the quadric { } is a cone with vertex T . If T is a point at infinity, the U~ = −→TU, V~ = −→TV, W~ = −−→TW, (25) quadric is a cylinder and the direction of its axis is given by T . if U,V,W are proper points. If one of them is a point at infinity, we take it as vector of the basis. For instance, if U is a point at infinity, we take U~ = U. 7. Classification of quadrics The case of an improper point T can be handled simi- ~ ~ ~ ~ The classification of quadric patches is refined now that larly. We use in this case a basis U, V, T with U = −−→WU, ~ ~ { } we know whether the quadric has a center or not: V = −−→W V , T = T , if W is a proper point. Otherwise, we choose U or V as origin and use W as one of the vectors λ> 0: Oval quadrics: • of the basis. – Centered: For a direction Ellipsoids: det > 0. ~v = αU~ + βV~ + γW,~ (26) ∗ Z Two-sheeted hyperboloids: det < 0. ∗ Z the polar plane is a diametral plane with linear form given – Non-centered: Elliptic paraboloids. by λ = 0: Degenerate quadrics: γ α β p β α γ q • ~v = + – Cones: if the vertex T is a proper point. Q ΩW − ΩU − ΩV ΩW ΩV − ΩU − ΩW ΩV – Cylinders: if T is a point at infinity. The type of α β γ r + λ(ǫU α + ǫV β + ǫW γ)t.(27) cylinder is determined classifying any of its conic ΩU − ΩV − ΩW ΩU − sections [16]. Before going on, we need information about the normal λ< 0: Ruled quadrics: of a plane given by a linear form: • 7 Lemma 1. A plane l with linear form ap + bq + cr + dt, In the case of improper T and, for instance, proper W , where p,q,r,t are linear forms satisfying the equations for the coordinates of ~v = αU~ +βV~ +δT~ are

p(W )= q(V )= r(U)= t(T )=1, b ǫV a = µ(gV U α + gV V β + gV T δ), − contains the vectors c ǫU a = µ(gUU α + gUV β + gUT δ), − ~v = (b ǫV d)U~ + (ǫU d c)V,~ d = µ(gT U α + gT V β + gT T δ), (31) 1 − − ~v = (a ǫW d)V~ + (ǫV d b)W,~ where 2 − − ~ ~ ǫU α + ǫV β α β 1 ~v3 = (ǫW d a)U + (c ǫU d)W, a = , − − − ΩW − ΩU − ΩV ΩW and hence its normal vector is given by β α ǫU α + ǫV β 1 ~ ~ ~ ~ ~ ~ b = + , ~n = (a ǫW d)U V + (c ǫU d)V W + (b ǫV d)W U. Ω − Ω Ω Ω − × − × − × V U W V The proof is simple, since a vector ~v = αU~ + βV~ + γW~ α β ǫU α + ǫV β 1 c = + , belongs to the plane if and only if Ω − Ω Ω Ω U V W U 0= l(~v) = (c ǫU d)α + (b ǫV d)β + (a ǫW d)γ. d = λδ. (32) − − − A diametral plane is called a principal plane or a plane These conditions can be seen as arising from an alter- of symmetry if it is orthogonal to its pole ~v. The principal native definition of principal axes as lines with direction axis of the quadric are the lines which are intersections of given by eigenvectors of the bilinear form of the conic at two principal planes. That is, the poles of the principal infinity. The values of the coefficient µ are the correspond- planes are the directions of the axes. ing eigenvalues, which are obtained by imposing that the Since for this definition we need to include a scalar prod- system (29) has non-trivial solutions for α,β,γ. Hence uct, it is necessary to provide another symmetric bilinear µ has to satisfy a cubic equation and there are in gen- form G, such that ~v ~w = G(~v, ~w), denoting by a dot the eral three principal planes and axes, except for quadrics of · scalar product. revolution and spheres: The matrix of such form is usually called Gram matrix and it is Corollary 3. A rational triangular quadratic patch with control net c ,c ,c ,c ,c ,c and weights g g g 002 011 020 101 110 200 UU UV UW ω ,ω ,ω{ ,ω ,ω ,ω , such that} the three G = g g g where g := A~ B,~ (28) 002 011 020 101 110 200 V U V V V W AB {boundary conics meet at a point}S, which is written as in  g g g  · W U W V W W (3), has diametral planes with linear forms given by (27)   in the basis U,~ V,~ W~ . with a pole ~v = αU~ + βV~ + γW~ , or ~v = αU~ + βV~ + γT~, { } In order to derive conditions for principal planes, we if T is a point at infinity, in a vector basis (25), defined have to impose that the pole ~v be orthogonal to a basis with points (5) and (10). of vectors of the diametral plane, which can be two of the The poles of principal planes are the directions of princi- ones which have been calculated in Lemma 1, pal axes and, if T is a proper point, they satisfy the linear system (29) for values of µ for which the determinant 0 = ~v1 ~v = (b ǫV d)(gUU α + gUV β + gUW γ) · − 1 1 1 W U W U W V W V W W W + (ǫ d c)(g α + g β + g γ), λ Ω Ω µg λ Ω Ω µg λ + Ω2 µg U V U V V V W − W U − − W V − W − − 1 1 1 λV U µgV U λV + µgV V λV W µgV W 0 = ~v ~v = (a ǫ d)(g α + g β + g γ) ΩV ΩU Ω2 ΩV ΩW 2 W V U V V V W − − V − − − · − 1 1 1 U UU UV UV UW UW λ + Ω2 µg λ Ω Ω µg λ Ω Ω µg + (ǫV d b)(gW U α + gW V β + gW W γ), U − − U V − − U W − − and these equations can be easily solved up to a propor- vanishes, with λ := ǫ ǫ λ and coefficients given by (6) AB A B tionality factor µ, and (28).

a ǫW d = µ(gW U α + gW V β + gW W γ), If T is a point at inﬁnity, the coordinates of the pole − satisfy the linear system (31) for values of µ for which the b ǫV d = µ(gV U α + gV V β + gV W γ), − determinant c ǫU d = µ(gUU α + gUV β + gUW γ), (29) − AUU µgUU AUV µgUV µgUT where − − − AUV µgV U AV V µgV V µgV T − − − γ α β 1 µgT U µgT V λ µgT T a = , − − − ΩW − ΩU − ΩV ΩW vanishes, with β α γ 1 b = , 2 2 ΩV − ΩU − ΩW ΩV ǫU 1 ǫV 1 AUU := + , AV V := + , α β γ 1 ΩW ΩU ΩW ΩV c = , Ω − Ω − Ω Ω U V W U ǫU ǫV ǫV ǫU 1 AUV := + + . d = λ(ǫU α + ǫV β + ǫW γ). (30) Ω2 Ω Ω Ω Ω Ω Ω − W W U W V − U V 8 Z γ =ΩW (Ω˜ U + Ω˜ V ),

we get the diﬀerences A = a ǫW d, B = b ǫV d, C = − − c ǫU d for the planes which are orthogonal to the center. That− is, we have the linear form for the tangent plane at the vertex of the paraboloid, except for the coeﬃcient d. Since tangent planes are solutions of the implicit equation O of the quadric in tangential form (13), o 0=ΩU ΩV (C + ǫU d)(B + ǫV d)

Figure 5: The vertex O of a paraboloid is the pole of a tangent plane + ΩV ΩW (B + ǫV d)(A + ǫW d) with the center Z as normal. + ΩW ΩU (A + ǫW d)(C + ǫU d) 1

In general there are three different values for µ. If there − Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U is a double non-null solution, the quadric is a surface of = ΩU ΩV BC +ΩV ΩW AB +ΩW ΩU AD revolution. If there is a triple solution, the surface is a + Ω Ω (ǫ B + ǫ C)+Ω Ω (ǫ A + ǫ B) sphere. U V U V V W V W + Ω W ΩU (ǫU A + ǫW C) d, 3 2 The discriminant of a cubic equation a3µ +a2µ +a1µ+ a0 = 0, the coefficients of the plane are readily obtained 2 2 3 3 2 2 ∆= 27a a + 18a a a a 4a a 4a a + a a , ΩU ΩV BC +ΩV ΩW AB +ΩW ΩU AC − 3 0 0 1 2 3 − 3 1 − 2 0 2 1 d = , − ΩU ΩV (ǫU B+ǫV C)+ΩV ΩW (ǫV A+ǫW B)+ΩW ΩU (ǫU A+ǫU C) provides a simple way of checking whether a quadric is a ˜ ˜ ˜ ˜ surface of revolution: A vanishing discriminant ∆ is equiv- a = ǫW d + gW U ΩU (ΩV + ΩW )+ gW V ΩV (ΩW + ΩU ) alent to having a double root. + gW W ΩW (Ω˜ U + Ω˜ V ), Besides, a vanishing second derivative of the equation, b = ǫV d + gV U ΩU (Ω˜ V + Ω˜ W )+ gV V ΩV (Ω˜ W + Ω˜ U ) 3a µ+2a = 0, implies a triple root and hence the quadric 3 2 + g Ω (Ω˜ + Ω˜ ), would be a sphere. V W W U V ˜ ˜ ˜ ˜ It is easily checked that the eigenvalue µ = 0 just ap- c = ǫU d + gUU ΩU (ΩV + ΩW )+ gUV ΩV (ΩW + ΩU ) 1 pears for λ− = Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U in the case of + gUW ΩW (Ω˜ U + Ω˜ V ), (33) proper T . That is, for paraboloids. For these surfaces the plane at infinity is a principal plane: The plane at infinity where a,b,c,d are given by (30). is diametral, as it comprises the center, and it is principal, The vertex O is just the pole of the plane o = ap + bq + since it fulfills (29). Besides, it is the polar plane of the cr + dt: center, since it is tangent to the paraboloid at the center. In the case of improper T , the eigenvalue µ = 0 appears Corollary 4. A rational triangular quadratic patch for a only if λ = 0. Hence, cylinders have a null eigenvalue, paraboloid with control net c002,c011,c020,c101,c110,c200 and weights ω ,ω ,ω { ,ω ,ω ,ω , such that} corresponding to a pole T , the direction of the axis, which { 002 011 020 101 110 200} has no polar plane. For parabolic cylinders there is another the three boundary conics meet at a point S, which is writ- null eigenvalue, since the plane at infinity is an improper ten as in (3), has a vertex given by principal plane. (bΩ + aΩ )Ω U + (cΩ + aΩ )Ω V There are then just two proper principal planes for O = V W U U W V ΩU ΩV (b+c 2ǫU ǫV d)+ΩV ΩW (a+b 2ǫV ǫW d)+ΩW ΩU (a+c 2ǫW ǫU d) paraboloids, except for paraboloids of revolution. The in- − − − tersection of these principal planes is the only proper axis (cΩU + bΩV )ΩW W 2d(Ω˜ U Ω˜ V + Ω˜ V Ω˜ W + Ω˜ W Ω˜ U )T + − , of the paraboloid, a line with direction given by the center ΩU ΩV (b+c 2ǫU ǫV d)+ΩV ΩW (a+b 2ǫV ǫW d)+ΩW ΩU (a+c 2ǫW ǫU d) − − − of the paraboloid. The axis meets the paraboloid at the where A = a ǫW d, B = b ǫV d, C = c ǫU d and (33). center and at a single proper point named vertex. − − − We may calculate the vertex solving a quadratic equa- The coefficients are given in (6) and (19) and the points tion, but there is a simpler way, taking into account that U,V,W,T are defined in (5) and (10). the tangent plane o at the vertex O is orthogonal to the Finally, the axis of a cylinder is easily determined, since center Z (see Fig. 5). We can use this property to compute it has the direction of T and contains the center of every the vertex: conic section. For instance, we can use the center of the A plane with linear form ap+bq+cr+dt is orthogonal to conic at t, which according to (17) is given by a vector (26) if their coordinates fulfill (29). Hence reading the coordinates of Z from (16) with, for instance, µ = 1, ǫW ΩW ωwω002P + ǫV ΩV ωvω020Q + ǫU ΩU ωuω200R Zt = . α =ΩU (Ω˜ V + Ω˜ W ), β =ΩV (Ω˜ W + Ω˜ U ), ΩZ 9 c101 c101

c200 c200

c110 c110

c002 c011 c002 c011 c020 c020

Figure 6: Ellipsoid Figure 7: Two-sheeted hyperboloid

9. Examples The three eigenvalues calculated according to Corol- lary 3 are different, µ = 1.78, 0.90, 0.52, and so this el- Now we apply our results to several quadric patches: lipsoid is not a surface of revolution. (0, 0, 0) (1, 0, 1) (2, 0, 0) The three principal planes are Example 1. Net: ( 2, 2, 3) (1, 1, 1) and 0.61p +0.37q 1.32r 2.04t,  − 4 2  − − (0, 3 , 3 ) 1 1 1   0.90p + 15.05q 4.49r +5.04t, 2 − − weights: 1 2 (Fig. 6):  3 3  10.46p + 17.24q 0.57r 19.12t, 1 − − −   The normalised linear forms for the tangent planes are with respective implicit equations in cartesian coordinates 10x + 25y 10z 2 x z 12 3x 9y 0.39x +1.67y 1.13z 0.90=0, p = − , q = − − , r = − − , − − 32 8 16 1.32x +0.72y +0.60z +0.39=0, and meet at a point ( 8/7, 12/7, 22/7). − − The three boundary conics meet at a point, 5.32x 3.67y 7.25z +3.88=0, − − − S = (1, 0, 1) = c c + c − 002 − 011 020 and obviously meet at the center. c 2c /3+ c c 4c /3+ c = 002 − 101 200 = 200 − 110 020 , 4/3 2/3 (0, 0, 0) (1, 0, 1) (2, 0, 0) Example 2. Net: ( 1, 1, 2) (1, 1 , 1) and  − 2  and hence ωu = 1, ωv =4/3, ωw =2/3. (0, 2, 0) The normalised linear form for the plane through the 12 3   corners of the net is weights: 1 2 (Fig. 7):   7y + 14z 1 t = − . 32   The normalised linear forms for the tangent planes are The quadric is oval, since the coefficient λ = 207/49 for x +3y z x + z 2 12 4x 6y 5z this quadric patch is positive. p = − , q = − , r = − − − , It is not a paraboloid, since the center is the proper point 4 4 8 (1/2, 2/5, 1/30). Since the boundary curves and the and meet at a point (2/5, 2/5, 8/5). − conic on t are ellipses, the quadric is an ellipsoid. One ar- The three boundary conics meet at a point at infinity, rives to the same conclusion checking that det = 480/49 Z S = (2, 0, 4) = c 4c +3c is positive. − 002 − 011 020 The bilinear form that we get for this surface is = c002 2c101 + c200 = c200 4c110 +3c020, − − 18x2 + 45xy + 81y2 81yz + 54z2 36x 90y + 36z and hence ω = ω = ω = 1. = − − − , u v w Q 64 The normalised linear form for the plane through the and the implicit equation, in cartesian coordinates is corners of the net is 5z 0=18x2 + 45xy + 81y2 81yz + 54z2 36x 90y + 36z. t = . − − − 8 10 c110 The three boundary conics meet at a point at infinity,

S = (0, 0, 2) = c002 2c011 + c020 c200 − c020 c 2c + c c 2c + c = 002 − 101 200 = 200 − 110 020 , 1/4 3/4 c011 − c101 and hence ωu = 1, ωv = 1/4, ωw =3/4. The normalised linear− form for the plane through the c002 corners of the net is 4x 2y 4z Figure 8: Hyperbolic paraboloid t = − − . 3 The quadric is ruled, since the coefficient λ = 1 for The quadric is oval, since the coefficient λ = 3/25 for this quadric patch is negative. − this quadric patch is positive. It is a hyperbolic paraboloid, since the center is a point It is not a paraboloid, since the center is the proper point at infinity (0, 0, 3/64), which is also the direction of the (7/4, 1/2, 5/2). As the boundary curves are not ellipses, − − axis. but two hyperbolas and one parabola, the quadric is not The bilinear form that we get for this surface is an ellipsoid, but a two-sheeted hyperboloid. Accordingly, det = 16/25 is negative. 8z 4x2 +4y2 Z − = − , The bilinear form that we get for this surface is Q 3 4x2 + 12xy + 12y2 +6yz z2 8x 24y +8z and the implicit equation, in cartesian coordinates is = − − − , Q 48 0=2z x2 + y2. − and the implicit equation, in cartesian coordinates is The eigenvalues calculated according to Corollary 3 are 0=4x2 + 12xy + 12y2 +6yz z2 8x 24y +8z. different, µ = 0, 4/3, 4/3, and one of them is null, as it − − − is expected for a hyperbolic− paraboloid. The eigenvalues for the normal directions of the princi- The three principal planes are pal planes are different, µ = 0.33, 0.44, 0.30, and hence the hyperboloid is not a surface of− revolution. q + r + t, The principal planes are 2p + q +4r + t, 1.58p +0.45q 0.91r +0.08t, 3 3 − p + q +3r +3t. 0.36p +0.05q +0.30r +0.22t, −2 4 The first principal plane is the plane at infinity and the 0.79p +4.65q +1.21r 0.96t, − other ones have y = 0, x = 0 as implicit equations in cartesian coordinates. They all meet at the center. and have the respective equations in cartesian coordinates The vertex is the point (0, 0, 0), as it is clear from the 0.97x +1.87y +0.34z 1.60=0, form of the implicit equation. − 0.05x +0.05y 0.13z +0.43=0, (0, 1, 0) (2, 0, 0) (0, 1, 0) − − Example 4. Net: −(1, 1 , 1 ) (1, 1 , 1) and  − 2 2 2  0.75x 0.32y 0.39z 0.51=0, (1, 0, 1) − − − 11 1   and meet at the center as expected. weights: 1 1 (Fig. 9):  1  (0, 0, 0) (1, 0, 0) (2, 0, 2)   Example 3. Net: (0, 1 , 0) (1, 1 , 0) and The normalised linear forms for the tangent planes are  2 2  (0, 1, 1 ) − 2 x x 1 11   p =1 + y, q =1 y, r = x 1, − 2 − 2 − − weights: 1 1 (Fig. 8)  1  and meet at a point at infinity (0, 0, 1).   The three boundary conics are parabolas and do not The normalised linear forms for the tangent planes are meet at any point. Hence the patch does not belong to a non-degenerate quadric. If it is a degenerate quadric, it 4 4x +2z 1+2y +2z p = z, q = − , r = − , is then a parabolic cylinder with direction (0, 0, 1). It is − 3 3 easy to check, for instance, that their respective centers and meet at a point (1, 1/2, 0). are aligned and hence it is a degenerate quadric. 11 c200 c110 c110

c101

c002 c020 c200 c020

c101 c011 c011 c002

Figure 9: Parabolic cylinder Figure 10: Cone

The normalised linear form for the plane through the The normalised linear form for the plane through the corners of the net is corners of the net is 2+2x 3z t = x + z. t = − . − 2 Using (23) we find a set of weights ΩU = 1/2, ΩV = − The quadric is degenerate, since the coefficient for this 1=ΩW and hence the bilinear form that we get for this quadric patch is λ = 0. It is a cone, since the intersection surface is of the tangent planes is a proper point, which is the vertex = 4(x + y2 1), (0, 0, 0). Q − The bilinear form that we get for this surface is and the implicit equation, in cartesian coordinates is 16 = (x2 + y2 z2), x =1 y2. Q 3 − − The eigenvalues for the poles of the principal planes are and the implicit equation, in cartesian coordinates is µ =0, 0, 4, as it is expected for a parabolic cylinder. 0= x2 + y2 z2. The principal planes are − The eigenvalues for the poles of the principal planes are p + q + r, p q. µ = 16/3, 16/3, 16/3, and so the surface is a cone of − − The first one is the plane at infinity and the second one revolution. is the only proper principal plane of a parabolic cylinder, The principal planes are with equation in cartesian coordinates given by y = 0. p q r, − − (2, 0, 2) ( 4 , 4√3 , 4 ) ( 1 , √3 , 1) 3γp + (1 2γ)q (1 + γ)r, 3 9 3 2 2 − − Example 5. Net: ( 4 , 4√3 , 4 ) (2, 0, 1)  and 3 − 9 3 for every value of γ and their respective equations in carte-  ( 1 , √3 , 1)  sian coordinates are  2 − 2  3 9 1  8 16 z =0, weights: 9 1 (Fig. 10):  16 2  1 3γx + (2 γ)√3y =0.   − The normalised linear forms for the tangent planes are 10. Tensor product quadric patches z x x √3y +2z x + √3y +2z p = − , q = − − , r = − , 2 4 4 Tensor product patches are the most common way to model surfaces in CAD. In particular, in some cases and meet at a point (0, 0, 0). quadrics can be parametrised by biquadratic rational The three boundary conics meet at a point Bézier patches, 3c /8 9c /8+ c S = ( 1, 0, 1) = 002 011 020 2 2 − 2 2 − 1/4 ωij cij Bi (u)Bj (v) 3c /8 9c /8+ c i=0 j=0 = 002 − 101 200 = c c + c , c(u, v)= X X , u,v [0, 1], 1/4 200 − 110 020 2 2 ∈ 2 2 ωij Bi (u)Bj (v) and hence ωu =1/4= ωv, ωw = 1. i=0 j=0 X X 12 c22 S c12 c21 c01 c11 c02 c00 c02 q r c20 01 c10 c c12 c10 c11 u p v c00 S c20 c22 c21

Figure 11: Planes on a biquadratic tensor product patch Figure 12: Spherical rational tensor product patch

(a, 0,a) (a,a,a) (0,a,a) for a control net c00,c01,c02,c10,c11,c12,c20,c21,c22 and their respective weights.{ } Example 6. Net: (2a, 0, 0) (2a, 2a, 0) (0, 2a, 0)  (a, 0, a) (a, a, a) (0, a, a)  The patch is bounded by four conic arcs with control − − −  1 1 2  polygons c00,c01,c02 , c20,c21,c22 , c00,c10,c20 and { } { } { } with a =1/√2 and weights: 1/2 1/2 1 (Fig. 12): c02,c12,c22 , meeting two by two at the four corner ver- { }  1/2 1/2 1  tices c00, c02, c20, c22. Not every rational biquadratic patch is a quadric patch   We use a triangular patch through P = c , Q = c , [21, 22], but we can apply our knowledge about quadric 00 02 R = c with the following control net and weights triangular patches to them. 20

For instance, we can take P = c00, Q = c02, R = c20 (a, 0,a) (a,a,a) (0,a,a) 1 1 2 and define a triangular patch with these three corners, (2a, 0, 0) c110 , 1/2 ω110 , as we know that the tangent planes are defined by the  (a, 0, a)   1/2  − neighbouring vertices: p contains c00,c01,c10, q contains     and notice that the conic at u and the conic at v meet at c02,c01,c12 and r contains c20, c10, c21 (see Fig. 11). the point S = ( √2/2, 0, √2/2) and We already know the conic at u, defined by the control − points c ,c ,c and their weights, and the conic at { 00 01 02} c00 c10 + c20/2 v, defined by c ,c ,c and their weights. For the S = c00 2c01 +2c02 = − . { 00 10 20} − 1/2 conic at w we have the control points c02 and c20 and their weights, but we lack the intermediate control point There is no need to perform Möbius transformations, c110 and the weight ω110. since the weights already satisfy In order to have a triangular quadric patch, we can use 1 the other point S where the conics at u and v meet, besides ωu = ω 2ω + ω =1, ωv = ω 2ω + ω = , 00 − 01 02 00 − 10 20 2 P . We may reparametrise both conics as we did in (3) so that their weights satisfy but the other denominator is not determined,

ω c 2ω c + ω c ω c 2ω c + ω c 5 00 00 01 01 02 02 00 00 10 10 20 20 ωw = ω 2ω + ω = 2ω . S = − = − , 20 − 110 02 2 − 110 ωu ωv The control points and S deﬁne the planes p,q,r,t and where denominators, if S is a proper point, are their intersections,

ωu = ω00 2ω01 + ω02, ωv = ω00 2ω10 + ω20. √2 √2 3√2 √2 − − U = , , , V = √2, , 0 , 4 4 4 ! − 2 ! Now we can define the plane w as the one containing S,c ,c and complete Fig. 2 by computing the intersec- 02 20 T = (1, √2, √2, 0), tion points U,V,W on plane t. The barycentric combinations for U,V,W provide us the except for W , which is a point at infinity with direction value of ωw and hence of ω110 and c110. (1, 1, 1). − If the biquadratic patch is in fact part of a quadric sur- This means that ΩW = 1, but the normalisation term face, Theorem 1 provides its bilinear forms and we can for W calculate its geometric elements. We see it with an example: ω ωw + ω ωv + ω ωu =2ω 1 − 002 020 200 110 − 13 must vanish and hence ω110 = 1/2, ωw = 3/2 and the References representative for W is the vector [1] H. Pottmann, A. Asperl, M. Hofer, A. Kilian., Architectural geometry, Bentley Institute Press, Exton, 2007. √2 [2] G. Farin, Triangular Bernstein-Bézier patches., Comput. Aided W = ω ωwP + ω ωvQ + ω ωu = ( 1, 1, 1) . − 00 02 21 2 − − Geom. Design 3 (2) (1986) 83–127, ISSN 0167-8396. [3] T. Sederberg, D. Anderson, Steiner Surface Patches, IEEE The coefficients Computer Graphics and Applications 5 (1985) 23–36. [4] W. Boehm, D. Hansford, Bézier Patches on Quadrics, in: G. Farin (Ed.), NURBS for Curves and Surface Design, SIAM, ΩU =2, ΩV =1, ΩW =1, λ =2, 1–14, 1991. [5] S. Lodha, J. Warren, Bézier representation for quadric surface yield the expression for the bilinear form for the quadric, patches, Computer-Aided Design 22 (9) (1990) 574 – 579. [6] R. Dietz, J. Hoschek, B. Jüttler, An algebraic approach to r2 curves and surfaces on the sphere and on other quadrics, Com- = p2 + q2 + 2pq pr qr +2t2. Q 4 − − − puter Aided Geometric Design 10 (3-4) (1993) 211 – 229. [7] R. Dietz, J. Hoschek, B. Jüttler, Rational patches on quadric The normalised forms for the planes are surfaces, Computer-Aided Design 27 (1) (1995) 27 – 40. [8] A. Coffman, A. J. Schwartz, C. Stanton, The algebra and geometry of Steiner and other quadratically parametrizable surfaces, √2 √2 p = (x + z √2), q = (√2 y z), Computer Aided Geometric Design 13 (3) (1996) 257 – 286. − 2 − 3 − − [9] W. Degen, The Types of Triangular Bézier Surfaces, in: G. Mullineux (Ed.), The Mathematics of Surfaces VI, Claren- don Press, 153–170, 1996. √2 √2 √2 r = (√2 x + z), t = x + y , [10] G. Albrecht, Determination and classification of triangular 3 − 3 − 2 ! quadric patches, Computer Aided Geometric Design 15 (7) (1998) 675 – 697. [11] G. Albrecht, Rational quadratic Bezier triangles on quadrics, and so the implicit equation for the surface in cartesian in: F.-E. Wolter, N. Patrikalakis (Eds.), Computer Graphics coordinates is International, 1998. Proceedings, IEEE, Los Alamitos, CA, 34– 40, 1998. x2 + y2 + z2 1=0. [12] G. Albrecht, Geometric invariants of parametric triangular − quadric patches, International Electronic Journal of Geometry 4 (2) (2011) 63 – 84. [13] J. Sánchez-Reyes, M. Paluszny, Weighted radial displacement: 11. Conclusions A geometric look at Bézier conics and quadrics, Computer Aided Geometric Design 17 (3) (2000) 267 – 289. [14] G. Albrecht, An Algorithm for Parametric Quadric Patch Con- We have derived closed formulas in terms of construction, Computing 72 (1-2) (2004) 1–12. trol points and weights for several geometric elements of [15] G. Albrecht, M. Paluszny, M. Lentini, An intuitive way for con- quadrics in rational Bézier form, both in triangular and structing parametric quadric triangles, Computational and Ap- plied Mathematics (2015) 1–23. tensor product representation. To our knowledge, these [16] A. Cantón, L. Fernández-Jambrina, E. Rosado-Mar´ıa, Geomet- formulas have not been produced before. The main dif- ric characteristics of conics in Bézier form, Computer-Aided De- ference with other procedures for drawing geometric infor- sign 43 (11) (2011) 1413 – 1421. mation from rational triangular patches [3] is the use of [17] G. Albrecht, Rational Triangular Bézier Surfaces - Theory and Applications, Habilitationschrift, Fakultät für Mathematik, TU geometric entities such as tangent planes to the quadric München, Shaker-Verlag, Aachen, 1999. and their intersections as ingredients for obtaining bilin- [18] G. Farin, Curves and surfaces for CAGD: a practical guide, ear forms, and hence, implicit equations, for the surface. Morgan Kaufmann Publishers Inc., San Francisco, CA, USA, There are many ways of implicitising a parametric surface 5th edn., ISBN 1-55860-737-4, 2002. [19] A. Cantón, L. Fernández-Jambrina, E. Rosado Mar´ıa, [23], but the use of linear forms with geometrical meaning M. Váquez-Gallo, Implicit Equations of Non-degenerate Ratio- instead of cartesian coordinates simplifies this problem for nal Bézier Quadric Triangles, in: J.-D. Boissonnat, A. Cohen, quadric patches. Besides, these geometric entities appear O. Gibaru, C. Gout, T. Lyche, M.-L. Mazure, L. L. Schumaker naturally in the formulas for geometric elements because (Eds.), Curves and Surfaces, vol. 9213 of Lecture Notes in Com- puter Science, Springer International Publishing, ISBN 978-3- they are already present in the expressions for the bilin- 319-22803-7, 70–79, 2015. ear forms for the quadric. The use of projective algebraic [20] J. G. Semple, G. T. Kneebone, Algebraic projective geometry, geometry allows us to perform calculations in a synthetic Oxford University Press, London, 1952. [21] W. Boehm, Some remarks on quadrics, Computer Aided Geo- fashion, instead of resorting to cartesian coordinates. metric Design 10 (3–4) (1993) 231 – 236. Additionally we classify affine quadrics using one coeffi- [22] U. Fink, Biquadratische Bézier-Flächenstücke auf Quadriken, cient involving the weights of the patch. This can be done Master’s thesis, Fakultät Mathematik der Universität Stuttgart, without implicitising the quadric patch [10, 12], but the Mathematisches Institut der Universität Stuttgart, 1992. [23] T. Sederberg, D. Anderson, R. Goldman, Implicit represen- closed form for the implicit equations is what enables us tation of parametric curves and surfaces, Computer Vision, to derive closed formulas for geometric elements. Graphics, and Image Processing 28 (1) (1984) 72 – 84. The results are obtained initially for Bézier triangles, but are also extended to quadric patches in tensor product form. 14