Next Article in Journal
Fixed-Point Theorems in Fuzzy Metric Spaces via Fuzzy F-Contraction
Next Article in Special Issue
Geometric Algebra Framework Applied to Symmetrical Balanced Three-Phase Systems for Sinusoidal and Non-Sinusoidal Voltage Supply
Previous Article in Journal
Modeling Recidivism through Bayesian Regression Models and Deep Neural Networks
Previous Article in Special Issue
A Topological View of Reed–Solomon Codes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The μ-Basis of Improper Rational Parametric Surface and Its Application

by
Sonia Pérez-Díaz
1,*,† and
Li-Yong Shen
2,3,†
1
Dpto de Física y Matemáticas, Universidad de Alcalá, E-28871 Madrid, Spain
2
School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
3
Key Laboratory of Big Data Mining and Knowledge Management, CAS, Beijing 100190, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(6), 640; https://doi.org/10.3390/math9060640
Submission received: 12 February 2021 / Revised: 11 March 2021 / Accepted: 12 March 2021 / Published: 17 March 2021
(This article belongs to the Special Issue New Trends in Algebraic Geometry and Its Applications)

Abstract

:
The μ -basis is a newly developed algebraic tool in curve and surface representations and it is used to analyze some essential geometric properties of curves and surfaces. However, the theoretical frame of μ -bases is still developing, especially of surfaces. We study the μ -basis of a rational surface V defined parametrically by P ( t ¯ ) , t ¯ = ( t 1 , t 2 ) not being necessarily proper (or invertible). For applications using the μ -basis, an inversion formula for a given proper parametrization P ( t ¯ ) is obtained. In addition, the degree of the rational map ϕ P associated with any P ( t ¯ ) is computed. If P ( t ¯ ) is improper, we give some partial results in finding a proper reparametrization of V . Finally, the implicitization formula is derived from P (not being necessarily proper). The discussions only need to compute the greatest common divisors and univariate resultants of polynomials constructed from the μ -basis. Examples are given to illustrate the computational processes of the presented results.

1. Introduction

The study of representations of rational curves and surfaces is a fundamental task in computer aided geometric design (CAGD) and computer algebra. There exist two typical problems in the study of representations.
  • Implicitization: for a rational parametric curve or surface, implicitization is to find an algebraic expression of the curve or surface.
  • Proper Reparametrization: for an improper rational parametric curve or surface, proper reparametrization is to find a proper parametric expression of the curve or surface.
The parametric expression of a curve or surface is widely used in geometric modeling, such as NURBS representations. The algebraic equation, which is also called implicit equation, is another important representation, and this is much better than the parametric expression in determining whether or not a point is on the curve or surface. Hence the implicitization problem is classical in CAGD and there are implicitization algorithms for rational curves and surfaces proposed over the past several decades [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. Among all of these techniques, the Gröbner bases [2] is well-known, since it is theoretically complete. However, this method has exponential computational complexity and, thus, it is inefficient. This is the reason that people can not apply the Gröbner basis method for practical implicitization in application. Alternatively, in computational application, people prefer to find the implicit equations from certain implicit matrices. The implicit matrices can be constructed as resultant matrices or matrices of moving curves/surfaces. The implicit matrix of the curve is much simpler than that of the surface, since the curve only introduces one variable. Actually, the construction of bivariate resultant is not uniform and it is still a developing technique in computer algebra  [15,16,17,18,19,20,21]. The implicit matrix of moving curves/surfaces was introduced in [10] and developed by more researchers [5,11,22,23,24,25]. The μ -basis of a curve or surface was later defined by moving lines or moving planes with certain properties [5,23].
The implicit equation of the surface is generally included as a factor in its constructed resultants, but a constructed resultant may have extraneous factors that are not easy to identify and remove. For implicit matrices, some works attempt to construct the matrix whose determinant is exactly the implicit equation [10,24,26,27], but the ways to construct such implicit matrices are not complete for general surfaces. In the implicit matrix, there is more information than the implicit equation, such as singularities with multiplicity counting. Accordingly, people sometimes construct the implicit matrix by simple way and then the determinant of this matrix may have extraneous factors other than the implicit equation [22]. For the matrix constructed from the μ -basis of arbitrary three linearly independent syzygies, the extraneous factors are identified completely based on the analysis of base points or to parameters at infinity of tensor product surfaces [11,25], but some computations need the Gröbner bases of a zero dimensional algebraic variety.
When considering the proper reparametrization problem, an essential question is to decide whether a rational parametrization is proper. If a given rational parametrization is not proper, also called improper, a generic point lying on the variety corresponds to more than one parameter. On the other side, if a rational parametrization is improper, we ask whether it can be reparameterized, such that we can get a new proper parametrization. For algebraic curves, is well-known that the existence of a proper reparametrization for a given improper rational parametrization is certified by Lüroth’s Theorem [28]. One can have a look, for instance, at a previous bibliography, as, for instance, [29,30,31], where some efficient methods are proposed to find a proper reparametrization for an improper parametrization of an algebraic curve. For a given algebraic surface, Castelnuovo’s Theorem states that unirationality and rationality are equivalent over algebraically closed fields, but only some partial algorithmic methods approaching the problem are known (see [30,32]).
The μ -basis was first used in [17]. Here, the authors provided a representation for the implicit equation of a given curve defined parametrically. The μ -basis developing as a new algebraic tool can be used to obtain the parametric equation of a rational curve or a rational surface, in order to compute the implicit equation defining these varieties, and to study singularities and intersections [33]. There are several methods to compute the μ -basis for rational curves by computing two moving lines that satisfy the required properties [17], based on Gröbner basis [34] or based on vector elimination [23]. The μ -basis has also been generalized to rational surfaces [5], although the case of rational surfaces is different; for instance, the degrees of the μ -basis elements can be different. An algorithm to compute a μ -basis of a rational surface is designed that is based on polynomial matrix factorization [35]. Another possible way is to compute a basis of the syzygy module of the surface with the application of Quillen–Suslin Theorem [36,37]. In order to avoid the extraneous factor in implicitization, people tried to find the strong μ -bases of surfaces that have the very similar properties of the μ -bases of curves. However, the surfaces with strong μ -bases are relatively rare [25,38].
The μ -basis has shown different advantages by assuming that the given parametrization is always proper. A recent result attempts to find inverse formula, proper reparametrization, and algebraic equation for an improper parametrization of an algebraic curve by using μ -basis [39]. In this paper, we pay attention to the μ -basis of an improper parametrization of an algebraic surface and then apply the μ -basis in the problems of proper reparametrization and implicitization. There are intrinsical differences in the discussions between the surface and the curve; hence, some results of the curve can not be extended to the surface straightforward. After we give the definition of μ -basis of a rational parametric surface defined parametrically by P ( t ¯ ) , we find the inversion formula (if P is proper) and the degree of the rational map that is induced by P while using the μ -basis. Although the proper reparametrization problem is still opening, we address the problem of proper reparametrization partially based on the latest results and the properties of the μ -basis. As an important application, we derive the implicitization formula from the μ -basis from a given parametrization not being necessarily proper. Starting from the μ -basis, the computations only involve greatest common divisors (gcds) and univariate resultants of some polynomials constructed from the μ -basis. While the surface implicitization form μ -basis involved the computation of Gröbner bases in [5].
We have structured the article, as indicated below. In Section 2, we present some important definitions and properties for the μ -basis of rational surfaces. In Section 3, given a rational parametrization P , we study the inversion computation (if P is proper) and the degree of the induced rational map. In Section 4, we focus on the proper reparametrization problem while using μ -basis. In Section 5, we come to the implicitization problem from a given rational parametrization not necessarily proper using μ -basis. We finish the paper with Section 6, where we present a brief summary of our work.

2. μ -Basis for Rational Surfaces: Definition and Previous Results

Let R denote the polynomial ring K [ t 1 , t 2 ] over an algebraically closed field K of characteristic zero and R m denote the set of m dimensional row vectors with entries in the polynomial ring R. A submodule M of R m is a subset of R m , for which this condition holds: for any f 1 , f 2 M and h 1 , h 2 R , we have h 1 f 1 + h 2 f 2 M . A set of elements f i M , for i = 1 , , k , is called a generating set of M if for any m M , there exist h i R , for i = 1 , , k satisfying that
m = h 1 f 1 + + h k f k .
The Hilbert Basis Theorem states that every submodule M R m has a finite generating set. If, for any m M , the above expression is unique, then { f 1 , , f k } is called a basis of the module M. If a module has a basis, then it is called a free module. For any ( f 1 , , f k ) R k , the set
syz ( f 1 , , f k ) : = { ( h 1 , , h k ) R k | h 1 f 1 + + h k f k = 0 }
is a module over R, called a syzygy module [40]. An important result regarding syzygy modules is that if a , b , c , d R [ t 1 , t 2 ] are four relatively prime polynomials then, the syzygy module syz ( a , b , c , d ) is a free module of rank 3 (see [5]).
Let V a be a rational affine irreducible surface, and let
P a ( t ¯ ) = 1 ( t ¯ ) 4 ( t ¯ ) , 2 ( t ¯ ) 4 ( t ¯ ) , 3 ( t ¯ ) 4 ( t ¯ ) K ( t ¯ ) 3 , t ¯ = ( t 1 , t 2 )
be a rational affine parametrization of V a , where gcd ( 1 , 2 , 3 , 4 ) = 1 . Sometimes, we write the parametrization in the homogenous coordinate form P ( t ¯ ) = ( 1 ( t ¯ ) : 2 ( t ¯ ) : 3 ( t ¯ ) : 4 ( t ¯ ) ) and, in this case, we denote the surface in the projective space as V .
A moving surface of degree l is a family of algebraic surfaces with parameter pairs  ( t 1 , t 2 )
S ( x ¯ , t ¯ ) = i = 1 σ f i ( x ¯ ) b i ( t ¯ )
, where f i ( x ¯ ) , i = 1 , , σ are polynomials of degree l, and b i ( t ¯ ) R [ t ¯ ] , i = 1 , , σ (called blending functions) are linearly independent. We say that a moving surface follows the rational surface P if
4 ( t ¯ ) l S P a ( t ¯ ) , t ¯ = 0 .
We observe that the implicit equation of a given rational surface V is a moving surface of P . A moving plane is a moving surface of degree one. We denote the next moving plane
A ( t ¯ ) x 1 + B ( t ¯ ) x 2 + C ( t ¯ ) x 3 + D ( t ¯ ) x 4
by L ( t ¯ ) : = ( A ( t ¯ ) , B ( t ¯ ) , C ( t ¯ ) , D ( t ¯ ) ) R [ t ¯ ] 4 . In the following, we denote, by L t ¯ , the set of the moving planes that follow the rational surface that is parametrized by P . Thus, L t ¯ is exactly the syzygy module syz ( 1 , 2 , 3 , 4 ) . Now, we define the μ -basis of the rational surface P .
Definition 1.
Let p , q , r L t ¯ be three moving planes satisfying that [ p , q , r ] = k P ( t ¯ ) , where k is a nonzero constant. Subsequently, it is said that p , q , r form a μ-basis of the rational surface P . In the following, [ p , q , r ] denotes the outer product of p , q , r .
Geometrically, the above definition means that the point of rational surface P can be represented as the intersection of three moving planes p , q , r . This definition is generalized from the moving lines of a rational curves [23]. Notice that the result in the curve case was proposed twenty years ago, but the surface case has been a mystery for a long time. The μ -bases surfaces that have the very similar properties of the μ -bases of curves is called strong μ -basis, but the strong μ -bases are relatively rare [25]. Therefore, we have to study the μ -basis of the surface from initial definition and, here, we review some basic theorems in [5].
Theorem 1.
For any rational surface P , there always exist three moving planes p , q , r , such that [ p , q , r ] = k · P ( t ¯ ) holds. In fact, any basis p , q , r of syz ( 1 ( t ¯ ) , 2 ( t ¯ ) , 3 ( t ¯ ) , 4 ( t ¯ ) ) satisfies the above equality.
Theorem 2.
Let p , q , r be a μ-basis of the rational surface P . Thus, p , q , r provide a basis for the module L t ¯ (hence, L t ¯ is a free module). That is, for any l ( t ¯ ) L t ¯ , there exist some polynomials h i ( t ¯ ) , i = 1 , 2 , 3 , satisfying that
l ( t ¯ ) = h 1 p + h 2 q + h 3 r .
In addition, the above expression is unique.
An immediate consequence of the above theorems is that if p , q , r form a μ -basis if and only if p , q , r are a basis of syz ( 1 ( t ¯ ) , 2 ( t ¯ ) , 3 ( t ¯ ) , 4 ( t ¯ ) )

3. Inversion and Degree Using μ -Basis

Let K be an algebraically closed field of characteristic zero. We denote, by f ( x 1 , x 2 , x 3 ) K [ x 1 , x 2 , x 3 ] , the defining polynomial of a rational affine irreducible surface V a defined by the rational affine parametrization
P a ( t ¯ ) = 1 ( t ¯ ) 4 ( t ¯ ) , 2 ( t ¯ ) 4 ( t ¯ ) , 3 ( t ¯ ) 4 ( t ¯ ) K ( t ¯ ) 3 , t ¯ = ( t 1 , t 2 ) .
The homogeneous implicit polynomial defining the corresponding the projective rational surface V will be denoted as F ( x 1 , x 2 , x 3 , x 4 ) K [ x 1 , x 2 , x 3 , x 4 ] , where F ( x 1 , x 2 , x 3 , x 4 ) = x 4 deg ( f ) f ( x 1 / x 4 , x 2 / x 4 , x 3 / x 4 ) , and the parametrization in the homogenous coordinate form is given as P ( t ¯ ) = ( 1 ( t ¯ ) : 2 ( t ¯ ) : 3 ( t ¯ ) : 4 ( t ¯ ) ) .
Besides implicitization, other applications of μ -basis include, as in the case of algebraic curves, point inversion and, in general, the computation of the fiber. The point inversion problem can be stated, as follows: given a point Q on the space, decide whether the point is on a rational surface V defined parametrically by P ( t ¯ ) or not, and, in the affirmative case, compute the corresponding parameter values t 1 , t 2 , such that P ( t 1 , t 2 ) = Q . In this section, we recall some efficient algorithms that allow for computing the point inversion and, in general, the computation of the degree of the rational map that is induced by P . For this purpose, we will use μ -basis.
In order to deal with these problems, we first recall that, associated with the parametrization P a ( t ¯ ) , we consider the induced rational map ϕ P : K V a K 3 ; t ¯ P a ( t ¯ ) . We denote, by deg ( ϕ P ) , the degree of the induced rational map ϕ P (see [41] p. 143, and [42] p. 80). Observe that the birationality of ϕ P , which is the properness of the input parametrization, is characterized by deg ( ϕ P ) = 1 (see [41,42]). We additionally remind that deg ( ϕ P ) determines the cardinality of the fiber of a generic element (see Theorem 7, p. 76 in [41]). The degree measures the number of times the parametrization traces the curve when the parameter takes values in K 2 . Finally, let F P ( Q ) be the fiber of a point Q V a ; that is F P ( Q ) = P a 1 ( Q ) = { t ¯ K 2 | P a ( t ¯ ) = Q } .
In the following, given the projective parametrization P ( t ¯ ) of a surface V and
p ( t ¯ ) = ( p 1 , p 2 , p 3 , p 4 ) , q ( t ¯ ) = ( q 1 , q 2 , q 3 , q 4 ) , r ( t ¯ ) = ( r 1 , r 2 , r 3 , r 4 )
a μ -basis for P ( t ¯ ) , we consider a generic point Q = ( x 1 , x 2 , x 3 , x 4 ) on the surface, and the polynomials
p P ( t ¯ , x ¯ ) = p ( t ¯ ) · Q , q P ( t ¯ , x ¯ ) = q ( t ¯ ) · Q , r P ( t ¯ , x ¯ ) = r ( t ¯ ) · Q , x ¯ = ( x 1 , x 2 , x 3 , x 4 ) .
Remind that p P ( t ¯ , P ( t ¯ ) ) = q P ( t ¯ , P ( t ¯ ) ) = r P ( t ¯ , P ( t ¯ ) ) = 0 ( p , q , r is a μ -basis for P ( t ¯ ) ). We denote, by V 1 , V 2 , V 3 , the auxiliary curves over K ( V ) , defined, respectively, by the polynomials p P ( t ¯ , x ¯ ) , q P ( t ¯ , x ¯ ) , r P ( t ¯ , x ¯ ) K ( V ) [ t ¯ ] , where K ( V ) is the field of rational functions of the given surface.
Finally, let
S P ( t 1 , x ¯ ) = pp x ¯ ( Cont Z ( Res t 2 ( p P , q P + Z r P ) ) ) K ( V ) [ t 1 ] ,
where Cont x ( h ) returns the content of a polynomial h w.r.t. some variable x, pp x ( h ) returns the primitive part of a polynomial h with respect to a variable x and Res x ( h 1 , h 2 ) returns the resultant of two polynomials h 1 and h 2 w.r.t. some variable x. Similarly, one also considers the polynomial
T P ( t 2 , x ¯ ) = pp x ¯ ( Cont Z ( Res t 1 ( p P , q P + Z r P ) ) ) K ( V ) [ t 2 ] .
The computation of S P , T P can be done in two different ways. First, we consider that the implicit equation defining the input surface is known. In this case, we carry out the arithmetic over K ( V ) while using this implicit equation. We observe that, since I ( V ) = F ( I ( V ) represents the ideal of V ), the basic arithmetic on K [ V ] can be carried out by computing polynomial remainders. Thus, we conclude that the quotient field K ( V ) is computable. Furthermore, we note that we calculate the resultants of polynomials in K ( V ) [ t ¯ ] , which is a unique factorization domain, and we compute gcds of univariate polynomials over K ( V ) and, thus, in an Euclidean domain. In the second way, we avoid the requirement on the implicit equation. More precisely, the elements are represented (not uniquely) as function of polynomials in the variables x 1 , x 2 , x 3 , x 4 . We check the zero equality while using the input rational parametrization. This way could be too time consuming. In order to avoid this problem, one may test zero–equality by substituting a random point on the surface. The result of this test is correct with probability almost one. Additionally, one may also test the correctness of the computation of the inverse by checking it on a randomly chosen point on the given surface. In this way, we avoid the computation of the implicit polynomial.
In the following theorem, we provide the technique for computing the components of the inverse of a given rational proper parametrization P ( t ¯ ) . Additionally, we characterize the properness of P ( t ¯ ) . We illustrate this result in Example 1.
Theorem 3.
The rational parametrization P ( t ¯ ) is proper if and only if for a generic point Q = ( x 1 , x 2 , x 3 , x 4 ) on the surface, it holds that deg t 1 ( S P ) = 1 . In this case, the t 1 -coordinate of the inverse of P ( t ¯ ) is given by solving S P ( t 1 , x ¯ ) = 0 w.r.t the variable t 1 .
Proof. 
Using the results shown in [43] (see Proposition 1), we deduce that the non-constant t 1 -coordinates of the intersections points in V i , i = 1 , 2 , 3 are given by the roots of the polynomial S P ( t 1 , x ¯ ) . Thus, we only have to prove that P ( t ¯ ) is proper if and only there exists one unique point A = ( A 1 , A 2 ) ( V 1 V 2 V 3 ) ( K K ) 2 (K denotes the algebraic closure of the field K ( V ) ; i.e K = K ( V ) ¯ ). Indeed: first, let M = ( M 1 ( x ¯ ) , M 2 ( x ¯ ) ) be the inverse of the rational proper parametrization P ( t ¯ ) . Subsequently, M ( P ) = t ¯ and, thus, p P ( M ( P ) , P ) = q P ( M ( P ) , P ) = r P ( M ( P ) , P ) = 0 , which implies that
M V 1 V 2 V 3 K 2 .
In addition, since M is the inverse of P ( t ¯ ) , one has that M ( K K ) 2 . Hence, M ( V 1 V 2 V 3 ) ( K K ) 2 . Now, let us see that M is unique. Let M ( V 1 V 2 V 3 ) ( K K ) 2 . The equalities p P ( M , x ¯ ) = q P ( M , x ¯ ) = r P ( M , x ¯ ) = 0 imply that
p P ( R ( t ¯ ) , P ) = q P ( R ( t ¯ ) , P ) = r P ( R ( t ¯ ) , P ) = 0 , R ( t ¯ ) = M ( P ) .
Afterwards, by the properties of resultants and by Lemma 1, we get that P ( t ¯ ) = k P ( R ( t ¯ ) ) and since P is proper we deduce that R ( t ¯ ) = M ( P ( t ¯ ) ) = t ¯ . Thus, left composing by P 1 , we get that M = P 1 = M .
Reciprocally, because there exists a unique point in ( V 1 V 2 V 3 ) ( K K ) 2 , A is fixed under the action of the Galois group and, thus, A K ( V ) 2 . Reasoning similarly, as we did for the uniqueness in the above implication, one gets that A P = t ¯ and, then, we conclude that A is the inverse of P .    □
Remark 1.
Theorem 3 can be stated similarly for T P ( t 2 , x ¯ ) . More precisely, P ( t ¯ ) is proper if and only if, for a generic point Q = ( x 1 , x 2 , x 3 , x 4 ) on the surface, it holds that deg t 2 ( T P ) = 1 . In this case, the t 2 -coordinate of the inverse of P ( t ¯ ) is given by solving T P ( t 2 , x ¯ ) = 0 w.r.t the variable t 2 .
We also note that we may work over the affine space (i.e., x 4 = 1 ) and the obtained results are the same, but over the affine space. If x 4 = 1 , then the computations are more efficient.
Example 1.
Let V be the rational surface that is defined by the parametrization
P ( t ¯ ) = ( t 2 t 1 : t 2 + t 1 : t 2 t 1 : t 2 2 + t 1 2 + 2 ) .
First, we compute the polynomials
p P ( t ¯ , x ¯ ) = p ( t ¯ ) · x ¯ = 2 x 1 + t 1 x 2 + t 1 x 3 , q P ( t ¯ , x ¯ ) = q ( t ¯ ) · x ¯ = 2 t 1 x 1 + ( t 1 2 + 1 + t 2 t 1 ) x 2 x 3 t 1 x 4 , r P ( t ¯ , x ¯ ) = r ( t ¯ ) · x ¯ = ( 2 t 2 + 2 t 1 ) x 1 + ( t 2 2 t 1 2 ) x 2 + 2 x 3 + ( t 2 + t 1 ) x 4 ,
where the μ-basis is given as
p ( t ¯ ) = ( 2 , t 1 , t 1 , 0 ) , q ( t ¯ ) = ( 2 t 1 , t 1 2 + 1 + t 2 t 1 , 1 , t 1 ) , r ( t ¯ ) = ( 2 t 2 + 2 t 1 , t 2 2 t 1 2 , 2 , t 2 + t 1 ) .
Now, we determine S P ( t 1 , x ¯ ) and T P ( t 2 , x ¯ ) . We obtain
S P ( t 1 , x ¯ ) = 2 x 1 + t 1 x 2 + t 1 x 3 ,
T P ( t 2 , x ¯ ) = 2 x 1 x 4 x 2 2 x 1 x 4 x 3 + 2 x 1 x 2 2 t 2 + 2 x 1 x 2 t 2 x 3 4 x 1 2 x 3 x 2 x 3 2 x 3 3 + x 2 3 + x 3 x 2 2 .
Because deg t 1 ( S P ) = 1 , we conclude that P is proper and the first coordinate of the inverse is given as
I 1 = 2 x 1 x 2 + x 3 .
Reasoning similarly with T P , we obtain the second coordinate of the inverse, which is given as
I 2 = 2 x 1 x 4 x 2 + 2 x 1 x 4 x 3 + x 3 3 x 2 3 + 4 x 1 2 x 3 + x 2 x 3 2 x 3 x 2 2 2 x 1 x 2 ( x 2 + x 3 ) .
Based upon the above theorem, we may compute F P ( Q ) for a generic point Q and, thus, to obtain the degree of the rational map induced by deg ( ϕ P ) . For this purpose, we consider Q = P a ( s ¯ ) , where s ¯ = ( s 1 , s 2 ) are new variables, and the polynomials
S P ( t 1 , s ¯ ) = pp s ¯ ( Cont Z ( Res t 2 ( p P , q P + Z r P ) ) ) K [ t 1 , s ¯ ] ,
T P ( t 2 , s ¯ ) = pp s ¯ ( Cont Z ( Res t 1 ( p P , q P + Z r P ) ) ) K [ t 2 , s ¯ ] ,
where
p P ( t ¯ , s ¯ ) = p ( t ¯ ) · P ( s ¯ ) , q P ( t ¯ , s ¯ ) = q ( t ¯ ) · P ( s ¯ ) , r P ( t ¯ , s ¯ ) = r ( t ¯ ) · P ( s ¯ ) .
Remind that p P ( t ¯ , t ¯ ) = q P ( t ¯ , t ¯ ) = r P ( t ¯ , t ¯ ) = 0 ( p , q , r is a μ -basis for P ( t ¯ ) ). We denote, by V 1 , V 2 , V 3 , the auxiliary curves over K ( s ¯ ) defined, respectively, by the polynomials p P ( t ¯ , s ¯ ) , q P ( t ¯ , s ¯ ) , r P ( t ¯ , s ¯ ) K [ t ¯ , s ¯ ] . Thus, one obtains the following proposition.
Theorem 4.
For a generic point Q = P a ( s ¯ ) , it holds that
deg ( ϕ P ) = Card ( F P ( Q ) ) = deg t 1 ( S P ( t 1 , s ¯ ) ) = deg t 2 ( T P ( t 2 , s ¯ ) ) .
Proof. 
First, we use Proposition 1 in [43], and we deduce that the non-constant t 1 -coordinates of the intersections points in V i , i = 1 , 2 , 3 are given by the roots of the polynomial S P ( t 1 , s ¯ ) . Thus, we only have to prove that M F P ( Q ) if and only if M V 1 V 2 V 3 ( K K ) 2 , where K = K ( s ¯ ) ¯ is the algebraic closure of the field. Indeed, if M V 1 V 2 V 3 ( K K ) 2 thus p P ( M , s ¯ ) = q P ( M , s ¯ ) = r P ( M , s ¯ ) = 0 , which implies that
p ( M ) · P ( s ¯ ) = q ( M ) · P ( s ¯ ) = r ( M ) · P ( s ¯ ) = 0 .
Because
p ( M ) · P ( M ) = q ( M ) · P ( M ) = r ( M ) · P ( M ) = 0
, we get that P ( M ) = k P ( s ¯ ) with k 0 (since M K 2 ), which implies that P a ( M ) = P a ( s ¯ ) . Hence, M F P ( Q ) .
Reciprocally, let M F P ( Q ) . Subsequently, P a ( M ) = P a ( s ¯ ) which implies that P ( M ) = k P ( s ¯ ) with k 0 . Because
p ( M ) · P ( M ) = q ( M ) · P ( M ) = r ( M ) · P ( M ) = 0 ,
we get that
p ( M ) · P ( s ¯ ) = q ( M ) · P ( s ¯ ) = r ( M ) · P ( s ¯ ) = 0 .
Thus, p P ( M , s ¯ ) = q P ( M , s ¯ ) = r P ( M , s ¯ ) = 0 and, hence, M V 1 V 2 V 3 . Furthermore, since M F P ( Q ) , we also get that M ( K K ) 2 .    □
Clearly, Theorem 4 can be also stated for a generic point Q = ( x 1 , x 2 , x 3 , x 4 ) on the surface, and the polynomials
p P ( t ¯ , x ¯ ) = p ( t ¯ ) · Q , q P ( t ¯ , x ¯ ) = q ( t ¯ ) · Q , r P ( t ¯ , x ¯ ) = r ( t ¯ ) · Q , x ¯ = ( x 1 , x 2 , x 3 , x 4 ) .
(remind that p P ( t ¯ , P ( t ¯ ) ) = q P ( t ¯ , P ( t ¯ ) ) = r P ( t ¯ , P ( t ¯ ) ) = 0 ( [ p , q , r ] = k P ( t ¯ ) ). For this purpose, one considers the polynomials
S P ( t 1 , x ¯ ) = pp x ¯ ( Cont Z ( Res t 2 ( p P , q P + Z r P ) ) ) K ( V ) [ t 1 ] ,
and
T P ( t 2 , x ¯ ) = pp x ¯ ( Cont Z ( Res t 1 ( p P , q P + Z r P ) ) ) K ( V ) [ t 2 ] .
where computation can be done, as we described in the paragraph before Theorem 3, i.e., over the field of rational functions K ( V ) . Thus, one has the following corollary.
Corollary 1.
For a generic point Q = ( x 1 , x 2 , x 3 , x 4 ) V , it holds that
deg ( ϕ P ) = Card ( F P ( Q ) ) = deg t 1 ( S P ( t 1 , x ¯ ) ) = deg t 2 ( T P ( t 2 , x ¯ ) ) .
Remark 2.
From the proof of Theorem 4, we deduce that F P ( Q ) = V 1 V 2 V 3 ( K K ) 2 .
Example 2.
Let V be the rational surface that is defined by the parametrization
P ( t ¯ ) = ( t 2 2 t 1 2 t 1 4 : t 2 + t 2 3 + t 2 t 1 4 : t 1 2 + t 2 t 1 2 + t 2 2 t 1 2 t 1 4 : t 1 2 + t 2 2 t 1 2 + t 1 4 ) .
We determine the polynomials p P ( t ¯ , s ¯ ) = p ( t ¯ ) · P ( s ¯ ) , q P ( t ¯ , s ¯ ) = q ( t ¯ ) · P ( s ¯ ) , r P ( t ¯ , s ¯ ) = r ( t ¯ ) · P ( s ¯ ) , where the μ-basis is given by
p ( t ¯ ) = ( t 2 t 1 2 + 2 t 1 2 + 4 t 2 4 + 6 t 2 3 4 t 2 4 , 2 t 1 2 , 2 t 2 t 1 2 t 1 2 4 t 2 4 4 t 2 3 + 4 t 2 + 2 , t 2 t 1 2 + t 1 2 + 2 t 2 3 2 ) q ( t ¯ ) = ( 2 t 2 t 1 2 3 t 1 2 + t 2 + 1 , 0 , 2 t 2 t 1 2 + 2 t 1 2 t 2 1 , t 1 2 + 1 ) r ( t ¯ ) = ( 2 t 2 2 t 2 + 2 , 0 , 2 t 2 2 1 , t 2 + 1 ) .
Now, we determine S P ( t 1 , s ¯ ) (similarly, if we compute T P ( t 2 , s ¯ ) ), and we obtain
S P ( t 1 , s ¯ ) = pp s ¯ ( Cont Z ( Res t 2 ( p P , q P + Z r P ) ) ) = t 1 2 s 1 2 K [ t 1 , s ¯ ] .
Therefore, applying Theorem 4, we conclude that P is not proper and in fact deg ( ϕ P ) = 2 . Furthermore,
F P ( P a ( s ¯ ) ) = { ( s 1 , s 2 ) , ( s 1 , s 2 ) } .

4. On the Problem of the Reparametrization Using μ -Basis

In this section, we consider the problem of computing a rational proper reparametrization of a given algebraic surface defined by an improper parametrization. That is, given an algebraically closed field K , and P ( t ¯ ) , t ¯ = ( t 1 , t 2 ) , a rational parametrization of a surface V over K , we want to compute a proper parametrization of V , Q ( t ¯ ) , and R ( t ¯ ) ( K ( t ¯ ) K ) 2 , such that
P ( t ¯ ) = Q ( R ( t ¯ ) ) .
Notice that we consider Q ( R ( t ¯ ) ) , with R ( t ¯ ) = ( r 1 ( t ¯ ) / r ( t ¯ ) , r 2 ( t ¯ ) / r ( t ¯ ) ) ( K ( t ¯ ) K ) 2 , in homogenous form, i.e., P ( t ¯ ) = Q ( R ( t ¯ ) ) means that P ( t ¯ ) = Q r 1 ( t ¯ ) r ( t ¯ ) , r 2 ( t ¯ ) r ( t ¯ ) r ( t ¯ ) deg ( Q ) , which is a polynomial vector in homogenous form.
In this section, although we do not provide a solution to the general reparametrization problem, we show how the μ -basis can be used to provide some information concerning R ( t ¯ ) ( K ( t ¯ ) K ) 2 . We address the problem partially and the idea presented is based in the results in [32], but we trust that we could develop deeply these new approaches in future works, and more results concerning this topic allow us to get more advances.
The approach that is presented in this section is based on the computation of polynomial gcds and univariate resultants. These techniques always work and the time performance is very effective. The algorithm presented follows directly from the algorithm that was developed in [39], which solves the problem for the case of curves. Accordingly, we first outline this approach and illustrate it with an example.
Example 3.
Let C be the rational curve that is defined by the parametrization
P ( t ) = ( t 4 + t 3 + 2 t 2 + 2 t + 1 : ( t + 1 ) ( 3 t 2 + 2 t + 2 ) : 2 t 4 + 3 t 3 + 5 t 2 + 4 t + 2 ) .
In Step 2 of the algorithm, we determine the polynomials p P ( t , s ) and q P ( t , s ) , where the μ-basis is
p ( t ) = 40 t 2 + 40 t + 40 , 30 t 2 35 t 35 , 20 t 2 + 15 t + 15 T
q ( t ) = 30 t 2 + 20 t + 20 , 25 t 2 20 t 20 , 15 t 2 + 10 t + 10 T .
Now, we compute G P ( t , s ) ,
G P ( t , s ) = C 0 ( t ) + C 1 ( t ) s + C 2 ( t ) s 2 ,
where C 0 ( t ) = 10 t 2 ,    C 1 ( t ) = 10 t 2 , and C 2 ( t ) = 10 + 10 t . Because m : = deg t ( G P ) > 1 , we go to Step 5 of Algorithm 1, and we consider
R ( t ) = C 2 ( t ) C 0 ( t ) = 1 t t 2 .
Note that gcd ( C 0 , C 2 ) = 1 . Now, we determine the polynomials
L 1 ( s , x 1 ) = Res t ( x 1 3 ( t ) 1 ( t ) , s C 0 ( t ) C 2 ( t ) ) = ( x 3 + s x 3 + 2 x 1 s 2 x 3 3 s x 1 + 2 s 2 x 1 ) 2 ,
L 2 ( s , x 2 ) = Res t ( x 2 3 ( t ) 2 ( t ) , s C 0 ( t ) C 2 ( t ) ) = ( 3 s x 3 + 2 x 2 3 s x 2 + 2 s 2 x 2 2 s 2 x 3 ) 2 .
Finally, the algorithm outputs the proper parametrization Q ( t ) , and the rational function R ( t ) (see Step 7)
Q ( t ) = t 2 1 t : t ( 3 + 2 t ) : 2 3 t + 2 t 2 , R ( t ) = 1 t t 2 .
Algorithm 1 Proper Reparametrization for Curves using μ -Basis
Input a rational parametrization P ( t ) = 1 ( t ) : 2 ( t ) : 3 ( t ) , of a plane algebraic curve  C .
Output a rational proper parametrization Q ( t ) of C , and a rational function R ( t ) such that P ( t ) = Q ( R ( t ) ) .
Steps
1.
Compute a μ -basis of P . Let p ( t ) , q ( t ) be this μ -basis.
2.
Compute p P ( t , s ) = p ( t ) · P ( s ) , q P ( t , s ) = q ( t ) · P ( s ) .
3.
Compute G P ( t , s ) = gcd ( p P ( t , s ) , q P ( t , s ) ) = C m ( t ) s m + + C 0 ( t ) . Let m : = deg t ( G P ( t , s ) ) .
4.
If m = 1 , return Q ( t ) = P ( t ) , and R ( t ) = t . Otherwise go to Step 5.
5.
Consider R ( t ) = C i ( t ) C j ( t ) K ( t ) , such that C j ( t ) , C i ( t ) are not associated polynomials (i.e., C j ( t ) k C i ( t ) , k K ).
6.
For i = 1 , 2 , compute
L i ( s , x i ) = Res t ( x i 3 ( t ) i ( t ) , s C j ( t ) C i ( t ) ) = ( q i 2 ( s ) x i q i 1 ( s ) ) deg ( R ) .
7.
Return Q a ( t ) = q 11 ( t ) / q 12 ( t ) , q 21 ( t ) / q 22 ( t ) or the equivalent projective parametrization Q ( t ) , and R ( t ) = C i ( t ) / C j ( t ) .
The main idea of the result that we develop in this paper consists in computing a reparametrization of two auxiliary parametrizations (defining two space curves), P 1 and P 2 , directly defined from a given rational parametrization of the surface P (see Definition 2). Moreover, using that the degree of a rational map is multiplicative under composition, we get some results that relate the degree of the rational map that is induced by P with the degree of Q , and the degree of the rational maps induced by P 1 and P 2 . In addition, we also prove the relation with the degree w.r.t. the variables t 1 , t 2 of R ( t ¯ ) K ( t ¯ ) 2 .
To start with, we first provide the following lemma that analyzes the behavior of the μ -basis under reparametrizations.
Lemma 1.
Let p ˜ ( t ¯ ) , q ˜ ( t ¯ ) and r ˜ ( t ¯ ) be a μ-basis for a parametrization Q ( t ¯ ) of a surface V . Let R ( t ¯ ) ( K ( t ¯ ) K ) 2 . Subsequently, p ( t ¯ ) = p ˜ ( R ( t ¯ ) ) , q ( t ¯ ) = q ˜ ( R ( t ¯ ) ) , r ( t ¯ ) = r ˜ ( R ( t ¯ ) ) is a μ-basis for the reparametrization P ( t ¯ ) = Q ( R ( t ¯ ) ) .
Proof. 
Taking into account that p ˜ ( t ¯ ) , q ˜ ( t ¯ ) and r ˜ ( t ¯ ) is a μ -basis for Q ( t ¯ ) , from Theorem 1, it follows that [ p ˜ , q ˜ , r ˜ ] = k Q ( t ¯ ) for some non-zero constant k. Therefore, we easily get that [ p , q , r ] = k P ( t ¯ ) for some non-zero constant k. Hence, from Theorem 1, we conclude that p ( t ¯ ) , q ( t ¯ ) , r ( t ¯ ) is a μ -basis for P ( t ¯ ) .    □
In the next proposition, we assume that we know p , q and  r and R ( t ¯ ) = ( r 1 ( t ¯ ) , r 2 ( t ¯ ) ) ( K ( t ¯ ) K ) 2 , and we present a method for computing p ˜ , q ˜ and  r ˜ from p , q and  r , respectively. We state Proposition 1 for p = ( p 1 , p 2 , p 3 , p 4 ) and p ˜ = ( p ˜ 1 , p ˜ 2 , p ˜ 3 , p ˜ 4 ) . One reasons, similarly, to obtain q ˜ from q and r ˜ from r .
We assume w.l.o.g that p 4 0 m which implies that p ˜ 4 0 (otherwise, we consider another non-zero component of p ). Thus, we may write
( p 1 / p 4 , p 2 / p 4 , p 3 / p 4 ) = ( p ˜ 1 / p ˜ 4 , p ˜ 2 / p ˜ 4 , p ˜ 3 / p ˜ 4 ) ( R ( t ¯ ) ) .
Let us assume that gcd ( p 1 , p 4 ) = gcd ( p ˜ 1 , p ˜ 4 ) = 1 (otherwise, we simplify these rational functions). In addition, we note that, if p 1 = 0 , then we easily get that p ˜ 1 = 0 . For the case of p 1 0 , we consider ( r 1 , r 2 , p 1 / p 4 ) that can be seen as an affine parametrization of the surface defined by the irreducible polynomial p ˜ 1 ( x 1 , x 2 ) x 3 p ˜ 4 ( x 1 , x 2 ) K [ x 1 , x 2 , x 3 ] (note that gcd ( p ˜ 1 , p ˜ 4 ) = 1 and p ˜ 4 p ˜ 1 0 ). Hence, we only have to compute the implicit equation of that surface by applying, for instance, the method that is presented in [44].
Reasoning, similarly, ( r 1 , r 2 , p i / p 4 ) can be seen as a parametrization of the surface defined by the irreducible polynomial p ˜ i ( x 1 , x 2 ) x 3 p ˜ 4 ( x 1 , x 2 ) K [ x 1 , x 2 , x 3 ] , for i = 2 , 3 . Summarizing, we have the following proposition.
Proposition 1.
Under the conditions that are stated above, it holds that the implicit equation of the parametrization ( r 1 , r 2 , p i / p 4 ) is given as p ˜ i ( x 1 , x 2 ) x 3 p ˜ 4 ( x 1 , x 2 ) , for i = 1 , 2 , 3 .
In Remark 3, we apply the same idea that is stated in Proposition 1, but for the particular case of curves.
Given a rational projective parametrization N ( t ¯ ) of a surface over K, in Definition 2 we introduce some auxiliary parametrizations over K ( t i ) that are defined from N .
Definition 2.
Let N ( t ¯ ) be a parametrization with coefficients in K. We define thepartial parametrizations associated to N as the parametrizations N i ( t j ) : = N ( t ¯ ) with coefficients in K [ t i ] (i.e., N i is defined over K [ t i ] ), for i , j { 1 , 2 } and i j .
We note that the partial parametrization N i ( t j ) ( i j ) determines a space curve in K ( t i ) ¯ 3 , where K ( t i ) ¯ is the algebraic closure of K [ t i ] . In addition, we note that Definition 2 can also be stated for any N ( t ¯ ) K ( t ¯ ) 2 . That is, given N ( t ¯ ) K ( t ¯ ) 2 , one may consider N i ( t j ) : = N ( t ¯ ) ( K [ t i ] ) ( t j ) 2 (i.e., N is seen defined over K [ t i ] and in the variable t j ), for i , j { 1 , 2 } and i j . Similarly, one also may adapt Definition 2 for any polynomial n ( t ¯ ) K ( t ¯ ) .
The properness of the input parametrization P of a surface V can be characterized by means of the properness of its partial parametrizations. In particular, it is proved that P is birational if and only if its associated partial parametrizations, P i , are proper and P i 1 K ( x ¯ ) K ( t i ) , where x ¯ = ( x 1 , x 2 , x 3 , x 4 ) (see [32]).
In the following, given a rational affine parametrization P ( t ¯ ) of a surface V , we develop an algorithm that computes a rational parametrization Q ( t ¯ ) of V , and R ( t ¯ ) K ( t ¯ ) 2 , such that P ( t ¯ ) = Q ( R ( t ¯ ) ) . The algorithm is based on the computation of polynomial gcds and univariate resultants whose computing time performance is very satisfactory.
We prove that the partial parametrizations that correspond to the output parametrization, Q ( t ¯ ) , are proper (see Theorem 5), and we get properties relating the degree of ϕ P with the degree of the rational map ϕ Q , and the degree of R ( t ¯ ) (see Theorem 6). More precisely, we prove that
deg ( ϕ P ) = deg ( ϕ Q ) deg t 1 ( S ) deg t 2 ( T )
where
R ( t ¯ ) = ( S ( t ¯ ) , T ( S ( t ¯ ) , t 2 ) ) , S , T K ( t ¯ ) .
In Corollaries 2 and 3, we analyze in which conditions deg ( ϕ Q ) = 1 or, otherwise, the degree of the rational map induced by Q ( t ¯ ) is lower than the degree that is induced by the input parametrization P ( t ¯ ) .
In Theorem 5, we have that the partial parametrizations associated to the output parametrization, Q ( t ¯ ) , are proper (see [32]) but we can not ensure that Q is proper.
Theorem 5.
The partial parametrizations Q 1 ( t 2 ) and Q 2 ( t 1 ) associated to the parametrization Q computed by Algorithm 2 are proper.
Algorithm 2 Proper Reparametrization for Surfaces using μ -Basis
Input a rational parametrization P ( t ¯ ) = 1 ( t ¯ ) : 2 ( t ¯ ) : 3 ( t ¯ ) : 4 ( t ¯ ) of an algebraic surface V .
Output a rational parametrization Q ( t ¯ ) of V , and R ( t ¯ ) ( K ( t ¯ ) K ) 2 such that P ( t ¯ ) = Q ( R ( t ¯ ) ) .
Steps
1.
Compute a μ -basis of P . Let p ( t ¯ ) , q ( t ¯ ) , r ( t ¯ ) be this μ -basis.
2.
Apply Algorithm 1 to P 2 ( t 1 ) . If P 2 is not proper, then the algorithm returns the proper parametrization M 2 ( t 1 ) , and S 2 ( t 1 ) ( K [ t 2 ] ) ( t 1 ) ( S 2 ( t 1 ) = S ( t 1 , t 2 ) seen with coefficients in K [ t 2 ] ), such that P 2 ( t 1 ) = M 2 ( S 2 ( t 1 ) ) . Otherwise, the algorithm returns M ( t ¯ ) = P ( t ¯ ) (i.e M 2 ( t 1 ) = P 2 ( t 1 ) , and S 2 ( t 1 ) = t 1 ).
3.
Apply Algorithm 1 to M 1 ( t 2 ) . If M 1 is not proper, the algorithm returns the proper parametrization Q 1 ( t 2 ) , and T 1 ( t 2 ) ( K [ t 1 ] ) ( t 2 ) ( T 1 ( t 2 ) = T ( t 1 , t 2 ) seen with coefficients in K [ t 2 ] ) such that M 1 ( t 2 ) = Q 1 ( T 1 ( t 2 ) ) . Otherwise, the algorithm returns Q ( t ¯ ) = M ( t ¯ ) (i.e Q 1 ( t 2 ) = M 1 ( t 2 ) , and T 1 ( t 2 ) = t 2 ). Then,
P ( t ¯ ) = M ( S ( t ¯ ) , t 2 ) = Q ( t 1 , T ( t ¯ ) ) ( S ( t ¯ ) , t 2 ) = Q ( S ( t ¯ ) , T ( S ( t ¯ ) , t 2 ) ) .
4.
Return the rational parametrization Q ( t ¯ ) of the surface V , and
R ( t ¯ ) = ( S ( t ¯ ) , T ( S ( t ¯ ) , t 2 ) ) K ( t ¯ ) 2 .
From Algorithm 2, and while using that the degree of a rational map is multiplicative under composition, we deduce some properties that relate the degree the rational map ϕ P with the degree of the rational maps ϕ Q , ϕ M , ϕ P i , i = 1 , 2 , and with deg ( R ) , where R ( t ¯ ) = ( S ( t ¯ ) , T ( S ( t ¯ ) , t 2 ) ) . One reasons, similarly as in [32].
Theorem 6.
It holds that
deg ( ϕ P ) = deg ( ϕ Q ) deg t 1 ( S ( t ¯ ) ) deg t 2 ( T ( t ¯ ) ) , a n d
deg ( ϕ P 2 ) = deg t 1 ( S ( t ¯ ) ) , deg ( ϕ M 1 ) = deg t 2 ( T ( t ¯ ) ) .
In addition,
deg ( ϕ P ) = deg ( ϕ M ) deg t 1 ( S ( t ¯ ) ) , deg ( ϕ M ) = deg ( ϕ Q ) deg t 2 ( T ( t ¯ ) ) .
Corollary 2.
The following statements are equivalent:
(1)
Q is proper.
(2)
deg ( ϕ M ) = deg t 2 ( T ) .
(3)
deg ( ϕ P ) = deg t 1 ( S ) deg t 2 ( T ) .
Finally, in Corollary 3, we show in which conditions Algorithm 2 does not return a better reparametrization than the input one (in the sense of the degree of the rational map that is induced by the rational parametrization).
Corollary 3.
It holds that deg ( ϕ Q ) = deg ( ϕ P ) if and only if deg ( ϕ P 2 ) = deg ( ϕ M 1 ) = 1 . In particular, if deg ( ϕ P i ) = 1 for i = 1 , 2 , then deg ( ϕ Q ) = deg ( ϕ P ) .
We observe that, while using previous results, one may easily analyze whether some families of surfaces can be properly reparametrized using the approach presented in this paper. For instance, if deg ( ϕ P i ) 1 for some i = 1 , 2 , and deg ( ϕ P ) = n , where n is a prime number, then deg ( ϕ Q ) = 1 .
To finish this section, we illustrate Algorithm 2 with one example. The times of our implementation performance is similar to the times that were obtained in [32].
Example 4.
Let V be the rational surface that is defined by the parametrization
P ( t ¯ ) = ( 4 ( t ¯ ) ) : 2 ( t ¯ ) : 3 ( t ¯ ) : 4 ( t ¯ ) =
t 2 4 t 1 2 t 1 4 : t 2 2 + t 2 6 + t 2 2 t 1 4 : t 1 2 + t 2 2 t 1 2 + t 2 4 t 1 2 t 1 4 : t 1 2 + t 2 4 t 1 2 + t 1 4 .
For this purpose, in Step 1 of Algorithm 2, we compute a μ-basis of P and we get that is given by
p ( t ¯ ) = ( t 2 2 t 1 2 + 2 t 1 2 + 4 t 2 8 + 6 t 2 6 4 t 2 2 4 , 2 t 1 2 , 2 t 2 2 t 1 2 t 1 2 4 t 2 8 4 t 2 6 + 4 t 2 2 + 2 , t 2 2 t 1 2 + t 1 2 + 2 t 2 6 2 )
q ( t ¯ ) = ( 2 t 2 2 t 1 2 3 t 1 2 + t 2 2 + 1 , 0 , 2 t 2 2 t 1 2 + 2 t 1 2 t 2 2 1 , t 1 2 + 1 )
r ( t ¯ ) = ( 2 t 2 4 t 2 2 + 2 , 0 , 2 t 2 4 1 , t 2 2 + 1 ) .
Using Theorem 4, one gets that deg ( ϕ P ) = 4 . Now, we apply Algorithm 1 to P 2 ( t 1 ) . We obtain that
G P 2 ( t 1 , s 1 ) = s 1 2 t 1 2 ( K [ t 2 ] ) [ t 1 , s 1 ] ,
and S 2 ( t 1 ) = t 1 2 ( K [ t 2 ] ) [ t 1 ] (remind that S 2 ( t 1 ) = S ( t 1 , t 2 ) is seen as a polynomial in the variable t 1 and with coefficients in K [ t 2 ] ). Subsequently, we determine the polynomials
L i ( s 1 , x i ) = Res t 1 ( x i 4 ( t ¯ ) i ( t ¯ ) , s 1 S 2 ( t 1 ) ) = ( m i , 2 ( s 1 ) x i m i , 1 ( s 1 ) ) deg t 1 ( S )
for i = 1 , 2 , 3 . We obtain that
M ( t ¯ ) = ( ( t 2 4 t 1 ) t 1 : t 2 2 ( t 2 4 1 + t 1 2 ) : ( 1 t 2 2 t 2 4 t 1 ) t 1 : t 1 ( 1 t 2 4 + t 1 ) ) .
Now, in Step 3 of the algorithm, we apply Algorithm 1 to M 1 ( t 2 ) . For this purpose, we first compute a μ-basis of M and we get that it is given by (see Remark 3)
p M ( t ¯ ) = ( 4 t 2 2 + 4 4 t 2 8 6 t 2 6 + t 2 2 t 1 + 2 t 1 , 2 t 1 , 4 t 2 2 2 + 4 t 2 8 + 4 t 2 6 2 t 2 2 t 1 t 1 , 2 t 2 6 + 2 + t 2 2 t 1 + t 1 )
q M ( t ¯ ) = ( t 2 2 + 1 + 2 t 2 2 t 1 + 3 t 1 , 0 , t 2 2 1 2 t 2 2 t 1 2 t 1 , 1 + t 1 )
r M ( t ¯ ) = ( 2 + 2 t 2 4 + t 2 2 , 0 , 2 t 2 4 + 1 , 1 + t 2 2 ) .
We obtain that
G M 1 ( t 2 , s 2 ) = s 2 2 t 2 2 ( K [ t 1 ] ) [ t 2 , s 2 ]
that is, M 1 is not proper. Afterwards, we compute T 1 ( t 2 ) = t 2 2 ( K [ t 1 ] ) [ t 2 ] , and the polynomials
L i ( s 2 , x i ) = Res t 2 ( x i m 4 ( t ¯ ) m i ( t ¯ ) , s 2 T 1 ( t 2 ) ) = ( q i , 2 ( s 2 ) x i q i , 1 ( s 2 ) ) deg t 2 ( T )
for i = 1 , 2 , 3 (remind that T 1 ( t 2 ) = T ( t 1 , t 2 ) is seen as a polynomial in the variable t 2 and with coefficients in K [ t 1 ] ). We obtain that
Q ( t ¯ ) = ( t 1 t 2 2 ) t 1 : t 2 ( t 2 2 1 + t 1 2 ) : ( t 1 + 1 + t 2 t 2 2 ) t 1 : t 1 ( 1 + t 1 t 2 2 ) .
In Step 4, the algorithm returns the parametrization Q ( t ¯ ) , and R ( t ¯ ) = ( S ( t ¯ ) , T ( S ( t ¯ ) , t 2 ) ) = ( t 1 2 , t 2 2 ) . We observe that
deg t 2 ( T ) = 2 , a n d deg t 1 ( S ) = 2 .
Thus, since deg ( ϕ P ) = 4 , by Theorem 6, we conclude that deg ( ϕ Q ) = 1 and, hence, Q is proper.
Remark 3.
Using Lemma 1, we may compute a μ-basis, p M , q M , r M , of M from the μ-basis, p , q , r , of P . Remind that P ( t ¯ ) = M ( S 2 ( t 1 ) , t 2 ) , and S 2 ( t 1 ) = S ( t 1 , t 2 ) is seen as a polynomial in the variable t 1 with coefficients in K [ t 2 ] . Thus, we have a particular case (case of curves) of the reasoning that is presented in Proposition 1 (which is stated for surfaces). More precisely, we write p = ( p 1 , p 2 , p 3 , p 4 ) and p M = ( p m 1 , p m 2 , p m 3 , p m 4 ) . Observe that the implicit equation of the parametrization ( S 2 ( t 1 ) , p i / p 4 ) (seen with coefficients in K ( t 2 ) and in the variable t 1 ) is given by the polynomial p m i ( x 1 , t 2 ) x 2 p m 4 ( x 1 , t 2 ) ( K ( t 2 ) ) [ x 1 , x 2 ] for i = 1 , 2 , 3 (i.e., the coefficients of the the implicit equation are in K ( t 2 ) ). In order to compute this implicit equation, we may use that
Res t 1 ( x i p 4 ( t ¯ ) p i ( t ¯ ) , x 1 S 2 ( t 1 ) ) = ( p m 4 ( x 1 , t 2 ) x i p m i ( x 1 , t 2 ) ) deg t 1 ( S ) , i = 1 , 2 , 3
(see, e.g., [44]). Similarly one reasons to get q M from q and r M from r . Observe that this is a particular case of the result presented in Proposition 1 (we apply the same idea stated in Proposition 1, but for the particular case of curves).

5. Implicitization Using μ -Basis

In the following, we assume that we are in the affine space (i.e., x 4 = 1 ; this simplifies the time on the computations), and we consider the polynomials
G 1 P ( t ¯ , x ¯ ) : = p P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) p 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 2 P ( t ¯ , x ¯ ) : = q P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) q 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 3 P ( t ¯ , x ¯ ) : = r P ( t ¯ , x ¯ ) K [ t ¯ , x ¯ ] , x ¯ = ( x 1 , x 2 , x 3 ) .
Note that, then, we may write
G 1 P ( t ¯ , x ¯ ) : = x 1 ( p 1 r 3 p 3 r 1 ) + x 2 ( p 2 r 3 p 3 r 2 ) + ( p 4 r 3 p 3 r 4 )
G 2 P ( t ¯ , x ¯ ) : = x 1 ( q 1 r 3 q 3 r 1 ) + x 2 ( q 2 r 3 q 3 r 2 ) + ( q 4 r 3 q 3 r 4 )
G 3 P ( t ¯ , x ¯ ) : = x 1 r 1 + x 2 r 2 + x 3 r 3 + r 4 .
In addition, let
S 12 P ( t 1 , x ¯ ) = pp x ¯ ( Res t 2 ( G 1 P , G 2 P ) ) K [ t 1 , x 1 , x 2 ] ,
T 12 P ( t 2 , x ¯ ) = pp x ¯ ( Res t 1 ( G 1 P , G 2 P ) ) K [ t 2 , x 1 , x 2 ] .
Finally, F P 12 ( x 1 , x 2 ) denotes the fiber of a point Q 12 : = π 12 ( Q ) = ( x 1 , x 2 ) , where Q = ( x 1 , x 2 , x 3 ) V a and π 12 ( V a ) is the ( 1 , 2 ) -projection of V a . That is
F P 12 ( Q 12 ) = P 12 1 ( Q 12 ) = { t ¯ K 2 | P 12 ( t ¯ ) = Q 12 } ,
where P 12 : = ( 1 / 4 , 2 / 4 ) : = π 12 ( P a ) .
Lemma 2.
It holds that
deg t 1 ( S 12 P ) = deg t 2 ( T 12 P ) = Card ( F P 12 ( x 1 , x 2 ) ) .
Proof. 
Because p , q , r is a μ -basis of P ( t ¯ ) , we have
p 1 1 + p 2 2 + p 3 3 + p 4 4 = 0 , q 1 1 + q 2 2 + q 3 3 + q 4 4 = 0 , r 1 1 + r 2 2 + r 3 3 + r 4 4 = 0 .
Consider a generic point Q = ( x 1 , x 2 , x 3 ) on the variety generated by ( 1 : 2 : 4 ) and the associated polynomials
p P ( t ¯ , x ¯ ) = p ( t ¯ ) · Q , a n d q P ( t ¯ , x ¯ ) = q ( t ¯ ) · Q ,
where p ( t ¯ ) = ( p 1 r 3 p 3 r 1 , p 2 r 3 p 3 r 2 , p 4 r 3 p 3 r 4 ) , and q ( t ¯ ) = ( q 1 r 3 q 3 r 1 , q 2 r 3 q 3 r 2 , q 4 r 3 q 3 r 4 ) . It holds that p P ( t ¯ , P ( t ¯ ) ) = q P ( t ¯ , P ( t ¯ ) ) = 0 . In fact, p P ( t ¯ , P ( t ¯ ) ) = ( p 1 r 3 p 3 r 1 ) 1 + ( p 2 r 3 p 3 r 2 ) 2 + ( p 4 r 3 p 3 r 4 ) 4 = 0 is derived by eliminating 3 from the first and third equations in (1). Similarly, to find q P ( t ¯ , P ( t ¯ ) ) = 0 from the last two equations in (1).
Thus, one may reason as in Theorem 4 and Corollary 1 (also see Remark 2) to get that
deg t 1 ( S 12 P ( t 1 , x ¯ ) ) = deg t 2 ( T 12 P ( t 2 , x ¯ ) ) = Card ( F P 12 ( x 1 , x 2 ) )
(remind that P 12 : = ( 1 / 4 , 2 / 4 ) = π 12 ( P a ) ). □
Theorem 7.
Let p ( t ¯ ) , q ( t ¯ ) and r ( t ¯ ) a μ-basis for P ( t ¯ ) . It holds that
pp x 3 ( h ( x ¯ ) ) = f ( x ¯ ) deg ( ϕ P )
where
h ( x ¯ ) = Cont { Z , W } ( Res t 2 ( T 12 P ( t 2 , x ¯ ) , K ( t 2 , Z , W , x ¯ ) ) ) K [ x ¯ ] ,
K ( t 2 , Z , W , x ¯ ) = Res t 1 ( S 12 P ( t 1 , x ¯ ) , H P ( t ¯ , Z , W , x ¯ ) ) K [ t 2 , Z , W , x ¯ ] ,
and
H P ( t ¯ , Z , W , x ¯ ) = G 3 P ( t ¯ , x ¯ ) + Z G 1 P ( t ¯ , x ¯ ) + W G 2 P ( t ¯ , x ¯ ) K [ t ¯ , Z , W , x ¯ ] .
Proof. 
First, we recall that
deg t 1 ( S 12 P ) = deg t 2 ( T 12 P ) = Card ( F P 12 ( x ¯ ) ) .
Let d 12 be this quantity. Clearly, d 12 1 . In addition, let m = deg t 1 ( H P ) and K = deg t 2 ( H P ) . Regarding S 12 P and H P as polynomials in K ( t 2 , Z , W , x ¯ ) [ t 1 ] , and using that the resultant of two univariate polynomials is the product of the evaluations of one of them in the roots of the other, we get
K ( t 2 , Z , W , x ¯ ) = Res t 1 ( S 12 P , H P ) = A ( x ¯ ) m i = 1 d 12 H P ( α i , t 2 , Z , W , x ¯ ) ,
where A is the leading coefficient of S 12 P w.r.t. t 1 , and where { α 1 , , α d 12 } are the roots of S 12 P in the algebraic closure K ( x 1 , x 2 ) ¯ of K ( x 1 , x 2 ) (we regard S 12 P as an univariate polynomial in t 1 ). Similarly,
Res t 2 ( T 12 P , K ) = B ( x ¯ ) k j = 1 d 12 K ( β j , Z , W , x ¯ ) ,
where B is the leading coefficient of T 12 P w.r.t. t 2 , and { β 1 , , β d 12 } are the roots of T 12 P in K ( x 1 , x 2 ) ¯ (we regard T 12 P as a univariate polynomial in t 2 ). Therefore,
Res t 2 ( T 12 P , K ) = B k A m · d 12 i = 1 d 12 j = 1 d 12 H P ( α i , β j , Z , W , x ¯ ) .
By Lemma 2, there exist d 12 pairs of points ( α i , β j ) being in F P 12 ( x 1 , x 2 ) , and for each U ( x 1 , x 2 ) F P 12 ( x 1 , x 2 ) it holds that G 1 P ( U , x ¯ ) = G 2 P ( U , x ¯ ) = 0 . Thus,
Res t 2 ( T 12 P , K ) = B k A m · d 12 Q ( x ¯ , Z , W ) U F P 12 ( x 1 , x 2 ) G 3 P ( U , x ¯ ) ,
where
Q ( x ¯ , Z , W ) = 1 i , j d 12 ( α i , β j ) F P 12 ( x 1 , x 2 ) H ( α i , β j , Z , W , x ¯ ) .
Note that for each root α i there exists a unique b j satisfying that the pair ( α i , β j ) is in the fiber. Furthermore, for ( α i , β j ) F P 12 ( x 1 , x 2 ) , either G 1 P ( α i , β j , x ¯ ) 0 or G 2 P ( α i , β j , x ¯ ) 0 (see Lemma 2). Hence, Q ( x ¯ , Z , W ) depends on Z or W. In addition, each H P ( α i , β j , Z , W , x ¯ ) does depend on Z or W.
Next, we show that Q ( x ¯ , Z , W ) , regarded as polynomial in K [ x ¯ ] ¯ [ x ¯ , Z , W ] , is primitive w.r.t. the variables { Z , W } . For this purpose, we denote, by N ( x 3 ) K [ x 1 , x 2 ] ¯ [ x 3 ] , the content of Q w.r.t. { Z , W } . Thus, there exists ( α i , β j ) F P 12 ( x 1 , x 2 ) satisfying that the polynomial N divides H ( α i , β j , Z , W , x ¯ ) ; that is, N ( x 3 ) divides G 3 P ( α i , β i , x ¯ ) + Z G 1 P ( α i , β j , x ¯ ) + W G 2 P ( α i , β j , x ¯ ) and, then, N ( x 3 ) divides G 1 P ( α i , β j , x 1 , x 2 ) and G 2 P ( α i , β j , x 1 , x 2 ) . Taking into account that at least one of them is not zero, we get that N K [ x 1 , x 2 ] ¯ and, thus, Q is primitive w.r.t. { Z , W } . Now, using that
h ( x ¯ ) = Cont { Z , W } ( Res t 2 ( T 12 P , K ) ) ,
we obtain that
h ( x ¯ ) = B k A m · d 12 · N ( x 1 , x 2 ) · U F P 12 ( x 1 , x 2 ) G 3 P ( U , x ¯ ) ,
where N K [ x 1 , x 2 ] ¯ . Thus,
pp x 3 ( h ( x ¯ ) ) = pp x 3 U F P 12 ( x 1 , x 2 ) G 3 P ( U ( x ¯ ) , x ¯ ) .
Under these conditions, it holds that deg x 3 ( pp x 3 ( h ( x ¯ ) ) ) = d 12 . Indeed, clearly one has deg x 3 ( pp x 3 ( h ( x ¯ ) ) ) d 12 . If deg x 3 ( pp x 3 ( h ( x ¯ ) ) ) < d 12 , thus, there exists U F P 12 ( x 1 , x 2 ) , such that r 3 ( U ) = 0 . Moreover, U K 2 , which is impossible since U F P 12 ( x 1 , x 2 ) .
Now, we prove that
pp x 3 U F P 12 ( x 1 , x 2 ) G 3 P ( U ( x ¯ ) , x ¯ ) = f ( x ¯ ) r .
Indeed, clearly one has that
pp x 3 U F P 12 ( x 1 , x 2 ) G 3 P ( U ( x ¯ ) , x ¯ ) = f ( x ¯ ) r g ( x ¯ ) .
Furthermore, r deg ( ϕ P ) , since G 3 P ( U ( x ¯ ) , x ¯ ) = G 3 P ( V ( x ¯ ) , x ¯ ) for U , V F P ( x ¯ ) (observe that F P ( x ¯ ) F P 12 ( x 1 , x 2 ) ) ). Thus, since deg x 3 ( f ) = d 12 / deg ( ϕ P ) (see [44]) and deg x 3 ( pp x 3 ( h ( x ¯ ) ) ) = d 12 , we get that
d 12 / deg ( ϕ P ) · r + deg ( g ) = d 12
, which implies that d 12 ( 1 r / deg ( ϕ P ) ) = deg ( g ) and, hence, r deg ( ϕ P ) . Because r deg ( ϕ P ) , we conclude that deg ( g ) = 0 and r = deg ( ϕ P ) . □
In the following examples, we illustrate the above theorem. These examples are taken from [5].
Example 5.
Let V be the rational surface that is defined by the parametrization
P ( t ¯ ) = ( t 2 2 t 1 t 1 2 : t 2 + t 2 3 + t 2 t 1 2 : t 1 + t 2 t 1 + t 2 2 t 1 t 1 2 : t 1 + t 2 2 t 1 + t 1 2 ) .
We determine the polynomials p P ( t ¯ , x ¯ ) = p ( t ¯ ) · x ¯ , q P ( t ¯ , x ¯ ) = q ( t ¯ ) · x ¯ , r P ( t ¯ , x ¯ ) = r ( t ¯ ) · x ¯ , where the μ-basis is given by
p ( t ¯ ) = ( t 2 t 1 + 2 t 1 + 4 t 2 4 + 6 t 2 3 4 t 2 4 , 2 t 1 , 2 t 2 t 1 t 1 4 t 2 4 4 t 2 3 + 4 t 2 + 2 , t 2 t 1 + t 1 + 2 t 2 3 2 )
q ( t ¯ ) = ( 2 t 2 t 1 3 t 1 + t 2 + 1 , 0 , 2 t 2 t 1 + 2 t 1 t 2 1 , t 1 + 1 )
r ( t ¯ ) = ( 2 t 2 2 t 2 + 2 , 0 , 2 t 2 2 1 , t 2 + 1 ) .
Now, we determine
G 1 P ( t ¯ , x ¯ ) : = p P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) p 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 2 P ( t ¯ , x ¯ ) : = q P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) q 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 3 P ( t ¯ , x ¯ ) : = r P ( t ¯ , x ¯ ) K [ t ¯ , x ¯ ] , x ¯ = ( x 1 , x 2 , x 3 )
and we compute
S 12 P ( t 1 , x ¯ ) = pp x ¯ ( Res ( G 1 P , G 2 P , t 2 ) ) = t 1 + 2 t 1 2 + x 1 + 2 x 1 2 x 2 2 t 1 2 x 1 2 + x 2 2 t 1 2 x 1 4 x 2 2 t 1 3 x 1 + 4 x 2 2 t 1 4 x 1 + 4 x 2 2 t 1 3 x 1 2 x 2 2 t 1 2 x 1 3 4 x 2 2 t 1 4 x 1 3 + 4 x 2 2 t 1 3 x 1 3 4 x 2 2 t 1 4 x 1 2 + 4 x 2 2 t 1 4 4 x 2 2 t 1 3 + x 2 2 t 1 2 5 x 1 t 1 4 t 1 5 x 1 2 x 1 4 t 1 + 2 x 1 3 t 1 4 + 2 x 1 3 t 1 5 + 2 x 1 2 t 1 5 + x 1 5 t 1 2 + 3 x 1 5 t 1 4 x 1 4 t 1 5 3 x 1 5 t 1 3 7 x 1 3 t 1 + x 1 t 1 3 x 1 5 t 1 5 + 16 x 1 3 t 1 2 4 x 1 2 t 1 4 6 x 1 2 t 1 3 + 8 x 1 4 t 1 2 5 x 1 t 1 + 14 x 1 2 t 1 2 9 x 1 2 t 1 + 7 x 1 t 1 2 + 6 x 1 4 t 1 4 11 x 1 4 t 1 3 14 x 1 3 t 1 3 t 1 5 + x 1 3 2 t 1 4 + t 1 3 ,
T 12 P ( t 2 , x ¯ ) = pp x ¯ ( Res ( G 1 P , G 2 P , t 1 ) ) = 2 t 2 5 x 1 + 4 x 1 t 2 3 2 t 2 4 x 2 + t 2 2 x 2 + t 2 5 + x 1 2 t 2 5 x 1 2 t 2 3 + 2 x 1 t 2 4 x 2 3 x 1 t 2 2 x 2 + x 1 x 2 t 2 2 x 1 t 2 + t 2 3 .
Also, let
H P ( t ¯ , Z , W , x ¯ ) = G 3 P ( t ¯ , x ¯ ) + Z G 1 P ( t ¯ , x ¯ ) + W G 2 P ( t ¯ , x ¯ ) K [ t ¯ , Z , W , x ¯ ]
and
K ( t 2 , Z , W , x ¯ ) = Res t 1 ( S 12 P ( t 1 , x ¯ ) , H P ( t ¯ , Z , W , x ¯ ) ) K [ t 2 , Z , W , x ¯ ] .
Finally, we compute
h ( x ¯ ) = Cont { Z , W } ( Res t 2 ( T 12 P ( t 2 , x ¯ ) , K ( t 2 , Z , W , x ¯ ) ) ) K [ x ¯ ] ,
and
pp x 3 ( h ( x ¯ ) ) = f ( x ¯ ) deg ( ϕ P ) = 1 46 x 1 2 x 2 x 3 + 38 x 1 x 2 x 3 2 8 x 1 x 2 x 3 + 4 x 2 2 x 1 2 x 3 12 x 2 2 x 1 x 3 + 10 x 1 3 x 2 x 3 + x 1 2 x 2 x 3 2 10 x 1 x 2 x 3 3 + 4 x 1 x 2 2 x 3 2 10 x 2 x 3 3 + 8 x 2 2 x 3 2 4 x 2 2 x 3 + 5 x 2 x 3 2 5 x 1 4 x 2 + 19 x 1 3 x 2 4 x 2 2 x 3 3 + 4 x 2 x 3 4 + 2 x 1 2 x 2 2 x 3 2 x 1 3 + 47 x 3 2 x 1 2 12 x 1 4 x 3 + 14 x 1 x 3 2 + 11 x 1 2 x 3 3 6 x 1 x 3 4 32 x 1 x 3 3 22 x 3 x 1 3 12 x 3 x 1 2 + 4 x 3 x 1 4 x 1 3 x 2 2 + 4 x 1 2 x 2 2 + 4 x 1 x 2 2 x 1 x 2 + 7 x 3 4 5 x 3 3 5 x 3 2 + 8 x 1 5 + 2 x 1 2 + 2 x 1 3 + x 3 5 + 3 x 3 2 x 1 + x 2 .
Observe that we may conclude that deg ( ϕ P ) = 1 and, thus, P ( t ¯ ) is a proper rational parametrization.
We have implemented this method while using Maple 2016 on a Lenovo ThinkPad Intel(R) Core(TM) i7-7500U CPU @ 2.70 GHz 2.90 GHz and 16 GB of RAM, OS-Windows 10 Pro. The time, in CPU seconds, for this example is 10.907 and using Gröbner basis, we get 0.187 .
Example 6.
Let V be the rational surface defined by the parametrization
P ( t ¯ ) = ( t 1 2 + t 2 t 1 + 2 t 2 2 2 t 2 2 t 1 : t 1 2 + 2 t 2 t 1 + t 2 t 1 2 + 2 t 2 2 t 2 2 t 1 + 2 t 2 2 t 1 2 : t 1 2 + t 2 t 1 + 2 t 2 t 1 2 + 2 t 2 2 t 2 2 t 1 2 t 2 2 t 1 2 : 2 t 2 t 1 2 t 2 t 1 2 2 t 2 2 t 1 t 2 2 t 1 2 ) .
We determine the polynomials p P ( t ¯ , x ¯ ) = p ( t ¯ ) · x ¯ , q P ( t ¯ , x ¯ ) = q ( t ¯ ) · x ¯ , r P ( t ¯ , x ¯ ) = r ( t ¯ ) · x ¯ , where the μ-basis is given by
p ( t ¯ ) = ( 1344390 t 2 2 t 1 + 34075368 t 2 t 1 22657890 t 1 5710808 t 2 3 181563 t 2 2 23392736 t 2 4984080 , 1344390 t 2 2 t 1 25195836 t 2 t 1 + 10711400 t 1 + 3569255 t 2 3 + 1074194 t 2 2 + 18408656 t 2 + 4984080 , 1344390 t 2 2 t 1 17483628 t 2 t 1 11946490 t 1 + 2141553 t 2 3 892631 t 2 2 + 4984080 t 2 , 9704572 t 2 11391246 t 2 t 1 + 6590790 t 1 + 2855404 t 2 3 + 6075203 t 2 2 2492040 )
q ( t ¯ ) = ( 229530 t 2 t 1 50278 t 1 + 139288 t 2 2 174717 t 2 + 194136 , 131160 t 2 t 1 + 155206 t 1 87055 t 2 2 + 85766 t 2 194136 , 65580 t 2 t 1 + 104928 t 1 52233 t 2 2 + 88951 t 2 , 131160 t 2 t 1 + 100556 t 1 69644 t 2 2 58603 t 2 + 97068 )
r ( t ¯ ) = ( 8 t 2 3 + 11 t 2 2 4 t 2 + 4 , 5 t 2 3 6 t 2 2 + 8 t 2 4 , 3 t 2 3 5 t 2 2 4 t 2 , 4 t 2 3 + t 2 2 + 2 ) .
Now, we determine
G 1 P ( t ¯ , x ¯ ) : = p P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) p 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 2 P ( t ¯ , x ¯ ) : = q P ( t ¯ , x ¯ ) r 3 ( t ¯ ) r P ( t ¯ , x ¯ ) q 3 ( t ¯ ) K [ t ¯ , x 1 , x 2 ]
G 3 P ( t ¯ , x ¯ ) : = r P ( t ¯ , x ¯ ) K [ t ¯ , x ¯ ] , x ¯ = ( x 1 , x 2 , x 3 )
and we compute S 12 P ( t 1 , x ¯ ) and T 12 P ( t 2 , x ¯ ) . Additionally, let
H P ( t ¯ , Z , W , x ¯ ) = G 3 P ( t ¯ , x ¯ ) + Z G 1 P ( t ¯ , x ¯ ) + W G 2 P ( t ¯ , x ¯ ) K [ t ¯ , Z , W , x ¯ ]
and
K ( t 2 , Z , W , x ¯ ) = Res t 1 ( S 12 P ( t 1 , x ¯ ) , H P ( t ¯ , Z , W , x ¯ ) ) K [ t 2 , Z , W , x ¯ ] .
Finally, we compute
h ( x ¯ ) = Cont { Z , W } ( Res t 2 ( T 12 P ( t 2 , x ¯ ) , K ( t 2 , Z , W , x ¯ ) ) ) K [ x ¯ ] ,
and from pp x 3 ( h ( x ¯ ) ) , we get that
-4.6cm0cm f ( x ¯ ) = 449792 + 51270879 x 3 x 1 x 2 2 + 13092929 x 3 2 x 1 x 2 3482416 x 3 x 1 2 x 2 + 22904376 x 3 x 1 x 2 + 675054 x 1 x 3 6 862596 x 2 x 3 6 29225028 x 1 2 x 2 x 3 3 + 11830146 x 1 x 2 x 3 3 + 32231373 x 1 x 2 2 x 3 3 + 110760512 x 1 5 x 2 x 3 90099948 x 1 4 x 2 2 x 3 129717124 x 1 4 x 2 x 3 2 + 40844546 x 1 3 x 2 3 x 3 + 124810810 x 1 3 x 2 2 x 3 2 + 74702726 x 1 3 x 2 x 3 3 19905824 x 1 2 x 2 4 x 3 64077866 x 1 2 x 2 3 x 3 2 75736662 x 1 2 x 2 2 x 3 3 + 18645124 x 1 x 2 4 x 3 2 + 9481980 x 1 x 2 5 x 3 + 33323142 x 1 x 2 3 x 3 3 40980736 x 1 6 x 2 50645760 x 1 6 x 3 + 34182816 x 1 5 x 2 2 + 54875936 x 1 5 x 3 2 24633612 x 1 4 x 2 3 26906244 x 1 4 x 3 3 + 30171008 x 1 3 x 2 4 25238190 x 1 2 x 2 5 + 9961900 x 1 x 2 6 1918400 x 2 6 x 3 + 30499314 x 1 x 2 4 x 3 28601154 x 3 2 x 1 2 x 2 2 + 22849453 x 3 2 x 1 x 2 3 + 27085291 x 3 2 x 1 3 x 2 2534290 x 2 5 x 3 2 5374808 x 2 4 x 3 3 + 3235863 x 1 3 x 3 3 9477510 x 2 3 x 3 3 + 2815992 x 1 2 x 3 3 8545002 x 2 2 x 3 3 + 3831717 x 1 x 3 3 5284450 x 2 x 3 3 15996836 x 3 2 x 1 3 + 119044501 x 1 x 2 4 25566784 x 1 5 x 2 + 9070776 x 1 4 x 2 2 + 80443033 x 1 3 x 2 3 + 9622080 x 1 5 x 3 109224707 x 1 2 x 2 4 + 53021640 x 1 x 2 5 9557440 x 2 5 x 3 19184968 x 3 2 x 1 4 6794107 x 3 2 x 2 4 + 16282880 x 1 4 x 2 + 119771644 x 1 3 x 2 2 205506824 x 1 2 x 2 3 + 49430528 x 3 x 1 4 + 41274000 x 1 4 x 3 x 2 63080603 x 1 3 x 3 x 2 2 3378276 x 1 2 x 3 x 2 3 96074840 x 3 x 1 3 x 2 + 9101736 x 3 x 1 2 x 2 2 + 53646214 x 3 x 1 x 2 3 + 10909032 x 3 2 x 1 2 x 2 21171376 x 3 x 2 4 10777540 x 3 2 x 2 3 + 7054782 x 3 4 x 1 3 6762492 x 2 2 x 3 4 + 2808045 x 1 2 x 3 4 + 16635132 x 3 2 x 1 x 2 2 1041561 x 1 x 3 4 2701863 x 2 x 3 4 5920710 x 3 4 x 2 3 3201822 x 3 5 x 2 2 367011 x 3 5 x 1 1776762 x 3 5 x 2 2007234 x 3 5 x 1 2 + 10217583 x 2 x 1 x 3 4 26748144 x 3 4 x 2 x 1 2 + 23463306 x 3 4 x 2 2 x 1 + 6779862 x 3 5 x 2 x 1 + 4466880 x 1 3813952 x 2 520768 x 3 16392064 x 1 2 13928144 x 2 2 857264 x 3 2 + 23714880 x 1 3 28078300 x 2 3 573568 x 1 4 33492106 x 2 4 24141568 x 1 5 23550085 x 2 5 196096 x 1 6 9046200 x 2 6 + 39690 x 3 5 638436 x 3 4 1440780 x 3 3 + 3591536 x 3 x 1 2849856 x 3 x 1 2 + 4138196 x 3 2 x 1 + 31582352 x 1 x 2 4613120 x 3 x 2 90566848 x 1 2 x 2 + 92152884 x 1 x 2 2 15211892 x 3 x 2 2 + 88688752 x 1 3 x 2 195019560 x 1 2 x 2 2 + 140970299 x 1 x 2 3 25290352 x 3 x 1 3 24626550 x 3 x 2 3 5042468 x 3 2 x 2 + 1087832 x 3 2 x 1 2 10534970 x 3 2 x 2 2 + 18348032 x 1 7 1475500 x 2 7 119313 x 3 6 96066 x 3 7 ,
and deg ( ϕ P ) = 1 . That is, P ( t ¯ ) is a proper rational parametrization.
The time, in CPU seconds, for this example is 71.703 , and using Gröbner basis, we get a time that is > 5000 .
Remark 4.
In order to improve the time of computations, one may compute the polynomial
h ( x ¯ ) = Cont { Z , W } ( Res t 2 ( T 12 P ( t 2 , x ¯ ) , K ( t 2 , Z , W , x ¯ ) ) )
as gcd ( R 1 , R 2 ) , where
R i = Res t 2 ( T 12 P ( t 2 , x ¯ ) , K ( t 2 , a i , b i , x ¯ ) ) , i = 1 , 2
and a i , b i K are random constants. The answer is correct with a probability of almost one, since, taking into account the behavior of the gcd under specializations, this property holds in an open Zariski subset (see e.g., Lemmas 7 and 8 in [45]).

6. Conclusions

The μ -basis has shown as a bridge tool between the parametric form and the implicit form of curves and surfaces. Moreover, the μ -basis has also been introduced into applications in singularities analysis and collision detections. The μ -basis theory of curves are more complete than that of surfaces, but surfaces would certainly attract more attention, although the discussion is more difficult. We study the μ -basis further for improper rational surfaces. The results are essential to the theoretical completeness of the μ -basis of surface.
We show how the μ -basis allows for computing the inversion of a given proper parametrization P ( t ¯ ) of an algebraic surface. If P ( t ¯ ) is not proper, we show how the degree of the rational map that is induced by P ( t ¯ ) can be computed as well as the elements of the fiber. Furthermore and directly from P ( t ¯ ) , we propose a method to find a μ -basis for a proper reparametrization Q ( t ¯ ) with some assumptions. If P ( t ¯ ) is improper, we give some partial results in finding a proper reparametrization of V . Finally, we show how the μ -basis of a given not being necessarily proper parametrization also allows for computing the implicit equation of a given surface bysubsequence univariate resultants.
As the further work, the numerical consideration could be an interesting extension of the μ -basis theory. One possible way would consist in generalizing the symbolic computation to numerical situation using the ideas and techniques that have already been implemented in some other problems, such as the numerical proper reparametrization of surfaces [46].

Author Contributions

Writing—original draft, S.P.-D. and L.-Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially supported by Beijing Natural Science Foundation under Grant Z190004, NSFC under Grant 61872332, the Fundamental Research Funds for the Central University and FEDER/Ministerio de Ciencia, Innovación y Universidades-Agencia Estatal de Investigación/MTM2017-88796-P (Symbolic Computation: new challenges in Algebra and Geometry together with its applications).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

First author belongs to the Research Group ASYNACS (Ref. CT-CE2019/683).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adkins, W.A.; Hoffman, J.W.; Wang, H.H. Equations of parametric surfaces with base points via syzygies. J. Symb. Comput. 2005, 39, 73–101. [Google Scholar] [CrossRef] [Green Version]
  2. Buchberger, B. Groebner-Bases: An Algorithmic Method in Polynomial Ideal Theory. In Multidimensional Systems Theory-Progress, Directions and Open Problems in Multidimensional Systems; Bose, N., Ed.; Reidel Publishing Company: Dordrecht, The Netherlands, 1985; Chapter 6; pp. 184–232. [Google Scholar]
  3. Busé, L.; Chardin, M.; Jouanolou, J.P. Torsion of the Symmetric Algebra and Implicitization. Proc. Am. Math. Soc. 2009, 137, 1855–1865. [Google Scholar] [CrossRef] [Green Version]
  4. Busé, L.; Cox, D.A.; D’andrea, C. Implicitization of surfaces in P3 in the presence of base points. J. Algebra Its Appl. 2003, 2, 189–214. [Google Scholar] [CrossRef]
  5. Chen, F.; Cox, D.A.; Liu, Y. The μ-basis and implicitization of a rational parametric surface. J. Symb. Comput. 2005, 39, 689–706. [Google Scholar] [CrossRef] [Green Version]
  6. Cox, D.A.; Little, J.B.; O’Shea, D. Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, 4th ed.; Springer: New York, NY, USA, 2015. [Google Scholar]
  7. D’Andrea, C.; Emiris, I.Z. Sparse Resultant Perturbations. In Algebra, Geometry and Software Systems; Joswig, M., Takayama, N., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; pp. 93–107. [Google Scholar] [CrossRef]
  8. Emiris, I.Z.; Kotsireas, I. Implicitization exploiting sparseness. In Geometric and Algorithmic Aspects of Computer-Aided Design and Manufacturing; Dutta, D., Smid, M., Janardan, R., Eds.; DIMACS Series in Discrete Mathematics and Theoretical Computer Science 67; American Mathematical Society, Providence RI: Berlin/ Heidelberg, Germany, 2005; pp. 281–298. [Google Scholar]
  9. Manocha, D.; Canny, J.F. Algorithm for implicitizing rational parametric surfaces. Comput. Aided Geom. Des. 1992, 9, 25–50. [Google Scholar] [CrossRef]
  10. Sederberg, T.W.; Chen, F. Implicitization Using Moving Curves and Surfaces. In Proceedings of the 22nd Annual Conference on Computer Graphics and Interactive Techniques, SIGGRAPH ’95, Los Angeles, CA, USA, 6–11 August 1995; ACM: New York, NY, USA, 1995; pp. 301–308. [Google Scholar] [CrossRef]
  11. Shen, L.Y.; Goldman, R. Implicitizing Rational Tensor Product Surfaces Using the Resultant of Three Moving Planes. ACM Trans. Graph. 2017, 36, 167:1–167:14. [Google Scholar] [CrossRef]
  12. Shen, L.Y.; Goldman, R. Combining complementary methods for implicitizing rational tensor product surfaces. Comput.-Aided Des. 2018, 104, 100–112. [Google Scholar] [CrossRef]
  13. Shen, L.Y.; Pérez-Díaz, S.; Goldman, R.; Feng, Y. Representing rational curve segments and surface patches using semi-algebraic sets. Comput. Aided Geom. Des. 2019, 74, 101770. [Google Scholar] [CrossRef]
  14. Shen, L.Y.; Yuan, C.M. Implicitization using univariate resultants. J. Syst. Sci. Complex. 2010, 23, 804–814. [Google Scholar] [CrossRef]
  15. Busé, L.; Mantzaflaris, A.; Tsigaridas, E. Matrix formulae for resultants and discriminants of bivariate tensor-product polynomials. J. Symb. Comput. 2020, 98, 65–83. [Google Scholar] [CrossRef] [Green Version]
  16. Chionh, E.W.; Zhang, M.; Goldman, R. Fast Computation of the Bezout and Dixon Resultant Matrices. J. Symb. Comput. 2002, 33, 13–29. [Google Scholar] [CrossRef] [Green Version]
  17. Cox, D.A.; Sederberg, T.W.; Chen, F. The moving line ideal basis of planar rational curves. Comput. Aided Geom. Des. 1998, 15, 803–827. [Google Scholar] [CrossRef]
  18. Dickenstein, A.; Emiris, I.Z. Multihomogeneous resultant formulae by means of complexes. J. Symb. Comput. 2003, 36, 317–342. [Google Scholar] [CrossRef] [Green Version]
  19. Emiris, I.Z.; Mantzaflaris, A. Multihomogeneous resultant formulae for systems with scaled support. J. Symb. Comput. 2012, 47, 820–842. [Google Scholar] [CrossRef] [Green Version]
  20. Shi, X.; Jia, X.; Goldman, R. Using a bihomogeneous resultant to find the singularities of rational space curves. J. Symb. Comput. 2013, 53, 1–25. [Google Scholar] [CrossRef]
  21. Shi, X.; Wang, X.; Goldman, R. Using μ-bases to implicitize rational surfaces with a pair of orthogonal directrices. Comput. Aided Geom. Des. 2012, 29, 541–554. [Google Scholar] [CrossRef]
  22. Busé, L. Implicit matrix representations of rational Bézier curves and surfaces. Comput.-Aided Des. 2014, 46, 14–24. [Google Scholar] [CrossRef] [Green Version]
  23. Chen, F.; Wang, W. The μ-basis of a planar rational curve—Properties and computation. Graph. Model. 2002, 64, 368–381. [Google Scholar] [CrossRef]
  24. Lai, Y.; Chen, F.; Shi, X. Implicitizing rational surfaces without base points by moving planes and moving quadrics. Comput. Aided Geom. Des. 2019, 70, 1–15. [Google Scholar] [CrossRef]
  25. Shen, L.Y.; Goldman, R. Strong μ-bases for Rational Tensor Product Surfaces and Extraneous Factors Associated to Bad Base Points and Anomalies at Infinity. SIAM J. Appl. Algebra Geom. 2017, 1, 328–351. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, M.; Goldman, R.; Chionh, E.W. Efficient Implicitization of Rational Surfaces by Moving Planes. In Proceedings of the ASCM; Gao, X.-S., Wang, D., Eds.; World Scientific: Singapore, 2000; pp. 142–151. [Google Scholar]
  27. Zheng, J.; Sederberg, T.W.; Chionh, E.W.; Cox, D.A. Implicitizing rational surfaces with base points using the method of moving surfaces. In Topics in Algebraic Geometry and Geometric Modeling; American Mathematical Society: Providence, RI, USA, 2003; pp. 151–168. [Google Scholar] [CrossRef]
  28. Walker, R.J. Algebraic Curves; Princeton University Press: Princeton, NJ, USA, 1950. [Google Scholar]
  29. Gao, X.S.; Chou, S.C. On the Parameterization of Algebraic Curves. Appl. Algebra Eng. Commun. Comput. 1992, 3, 27–38. [Google Scholar] [CrossRef] [Green Version]
  30. Pérez-Díaz, S. On the problem of proper reparametrization for rational curves and surfaces. Comput. Aided Geom. Des. 2006, 23, 307–323. [Google Scholar] [CrossRef]
  31. Sederberg, T.W. Improperly parametrized rational curves. Comput. Aided Geom. Des. 1986, 3, 67–75. [Google Scholar] [CrossRef]
  32. Pérez-Díaz, S. A partial solution to the problem of proper reparametrization for rational surfaces. Comput. Aided Geom. Des. 2013, 30, 743–759. [Google Scholar] [CrossRef]
  33. Jia, X.; Shi, X.; Chen, F. Survey on the theory and applications of μ-bases for rational curves and surfaces. J. Comput. Appl. Math. 2018, 329, 2–23. [Google Scholar] [CrossRef]
  34. Zheng, J.; Sederberg, T.W. A Direct Approach to Computing the μ-basis of Planar Rational Curves. J. Symb. Comput. 2001, 31, 619–629. [Google Scholar] [CrossRef] [Green Version]
  35. Deng, J.; Chen, F.; Shen, L. Computing μ-bases of Rational Curves and Surfaces Using Polynomial Matrix Factorization. In Proceedings of the 2005 International Symposium on Symbolic and Algebraic Computation, ISSAC ’05, Beijing, China, 24–27 July 2005; ACM: New York, NY, USA, 2005; pp. 132–139. [Google Scholar] [CrossRef]
  36. Fabiańska, A.; Quadrat, A. Applications of the Quillen-Suslin theorem to multidimensional systems theory. In Gröbner Bases in Control Theory and Signal Processing; Park, H., Regensburger, G., Eds.; Radon Series on Computation and Applied Mathematics 3; de Gruyter Publisher: Berlin, Germany, 2007; pp. 23–106. [Google Scholar]
  37. Fabiańska, A.; Quadrat, A. QuillenSuslin: A Maple Implementation of a Constructive Version of the Quillen-Suslin Theorem. 2007. Available online: https://who.rocq.inria.fr/Alban.Quadrat/QuillenSuslin/ (accessed on 22 December 2020).
  38. Shen, L.Y.; Goldman, R. Algorithms for computing strong μ-bases for rational tensor product surfaces. Comput. Aided Geom. Des. 2017, 52–53, 48–62. [Google Scholar] [CrossRef]
  39. Pérez-Díaz, S.; Shen, L.Y. Inversion, degree, reparametrization and implicitization of improperly parametrized planar curves using μ-basis. Comput. Aided Geom. Des. 2021, 84, 101957. [Google Scholar] [CrossRef]
  40. Cox, D.A.; Little, J.B.; O’Shea, D. Using Algebraic Geometry; Graduate Texts in Mathematics; Springer: New York, NY, USA, 1998. [Google Scholar]
  41. Shafarevich, I.R. Basic Algebraic Geometry Schemes; 1 Varieties in Projective Space; Springer: Berlin/Heidelberg, Germany, 1994. [Google Scholar]
  42. Harris, J. Algebraic Geometry: A First Course; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  43. Pérez-Díaz, S.; Schicho, J.; Sendra, J. Properness and Inversion of Rational Parametrizations of Surfaces. Appl. Algebra Eng. Commun. Comput. 2002, 13, 29–51. [Google Scholar] [CrossRef]
  44. Pérez-Díaz, S.; Sendra, J. A univariate resultant-based implicitization algorithm for surfaces. J. Symb. Comput. 2008, 43, 118–139. [Google Scholar] [CrossRef] [Green Version]
  45. Pérez-Díaz, S.; Sendra, J. Computation of the degree of rational surface parametrizations. J. Pure Appl. Algebra 2004, 193, 99–121. [Google Scholar] [CrossRef] [Green Version]
  46. Shen, L.Y.; Pérez-Díaz, S.; Yang, Z. Numerical Proper Reparametrization of Space Curves and Surfaces. Comput.-Aided Des. 2019, 116, 102732. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pérez-Díaz, S.; Shen, L.-Y. The μ-Basis of Improper Rational Parametric Surface and Its Application. Mathematics 2021, 9, 640. https://doi.org/10.3390/math9060640

AMA Style

Pérez-Díaz S, Shen L-Y. The μ-Basis of Improper Rational Parametric Surface and Its Application. Mathematics. 2021; 9(6):640. https://doi.org/10.3390/math9060640

Chicago/Turabian Style

Pérez-Díaz, Sonia, and Li-Yong Shen. 2021. "The μ-Basis of Improper Rational Parametric Surface and Its Application" Mathematics 9, no. 6: 640. https://doi.org/10.3390/math9060640

APA Style

Pérez-Díaz, S., & Shen, L. -Y. (2021). The μ-Basis of Improper Rational Parametric Surface and Its Application. Mathematics, 9(6), 640. https://doi.org/10.3390/math9060640

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop