Next Article in Journal
Arousing Early Strategic Thinking about SDGs with Real Mathematics Problems
Next Article in Special Issue
Visualizing a Cubic Linkage through the Use of CAS and DGS
Previous Article in Journal
On the Residual Continuity of Global Attractors
Previous Article in Special Issue
Discovering Geometric Inequalities: The Concourse of GeoGebra Discovery, Dynamic Coloring and Maple Tools
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Behavior of a Surface Implicitly Defined

by
Elena Campo-Montalvo
,
Marián Fernández de Sevilla
and
Sonia Pérez-Díaz
*,†
Departamento de Automática, Universidad de Alcalá, E-28871 Madrid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(9), 1445; https://doi.org/10.3390/math10091445
Submission received: 13 April 2022 / Revised: 23 April 2022 / Accepted: 24 April 2022 / Published: 25 April 2022
(This article belongs to the Special Issue Symbolic Computation for Mathematical Visualization)

Abstract

:
In this paper, we introduce the notion of infinity branches and approaching surfaces. We obtain an algorithm that compares the behavior at the infinity of two given algebraic surfaces that are defined by an irreducible polynomial. Furthermore, we show that if two surfaces have the same asymptotic behavior, the Hausdorff distance between them is finite. All these concepts are new and represent a great advance for the study of surfaces and their applications.

1. Introduction

Algebraic curves and surfaces are essential entities for practical applications (see, for example [1,2]). In fact, one may find a lot of literature dealing with different problems related to curves defined by irreducible polynomials (see, e.g., [3,4]). In this paper, we introduce the concept of infinity branches for surfaces. Intuitively speaking, an infinity branch represents the behavior of an algebraic surface at the points with “sufficiently large coordinates”. Informally speaking and generalizing the situation of algebraic curves, an infinity branch is associated with a projective “place” centered at an “infinity point”, and it can be “parametrized” by means of “Puiseux series” (the formal definition of these notions can be seen in Section 2 and Section 3).
Infinity branches are necessary and essential for the study of surfaces since they reveal the behavior at a point at infinity of a real algebraic surface. For instance, the infinity branches of an implicit algebraic plane curve are an important tool to sketch its graph as well as to analyze its topology (see e.g., [5,6,7,8,9]). It should be mentioned that for the case of curves, in [10], the notion of a g-asymptote is introduced. A g-asymptote generalizes the well-known notion of linear asymptotes, and these asymptotes can be computed from the infinity branches. More precisely, we say that a curve C ¯ is a generalized asymptote (or g-asymptote) of an input curve C if C ¯ approaches C at some infinity branch. Furthermore, the input curve C cannot be approached by a new curve having a lower degree. In this paper, we generalize some of these notions introduced previously for curves in [5] to the case of surfaces.
The notion of an infinity branch allows us to define convergent branches and approaching surfaces (see Section 4). More precisely, we say that two infinity branches converge if they get closer as they tend to infinity. Furthermore, an algebraic surface V ¯ approaches V at its infinity branch B if V ¯ has another infinity branch B ¯ such that B ¯ is convergent with B. We obtain important results that characterize whether two algebraic surfaces are approaching.
Using these results, we obtain an algorithm that compares the the behavior of two surfaces at infinity (see Section 5). Finally, it is shown that if two algebraic surfaces defined implicitly have the same asymptotic behavior, then the Hausdorff distance between them is finite.
The results that we obtain in this paper (and also for the case of curves) are essential for some applications in the framework of computer-aided geometric design (CAGD) as, for instance, in the approximate parametrization problem (see [11]). This problem can be stated as follows: we are given a non-rational affine surface V (we assume that it is a perturbation of a rational surface), and we would like to compute, if it exists, a  parametrization of a new rational affine surface denoted by V ¯ , which is near to the input surface. The effectiveness of the algorithm depends on the closeness of V and V ¯ , and then, here, one needs to show that the Hausdorff distance between V and V ¯ is finite.
This paper is structured as follows: in Section 2, we present the notions and preliminaries that will be used throughout the paper. Section 3 introduces the concept of infinity branch, and here, we prove some important properties. Section 4 provides the notions of convergent branches and approaching surfaces. Additionally, we characterize whether or not two algebraic surfaces approach each other. The results presented in this section will be used in Section 5, where an algorithm to compare the asymptotic behavior of two algebraic surfaces is designed. Finally, we prove that if two given algebraic surfaces have the same asymptotic behavior, the Hausdorff distance between them is finite. We finish with a section for conclusions and future work (Section 6).

2. Preliminaries and Terminology

In the following, we present some concepts and notions that will be used throughout the paper. In particular, we introduce some previous results concerning local parametrizations and Puiseux series. For further details, one may see for instance, Section 2.5 in [3], Chapter 4 (Section 2) in [4], [12,13,14,15], etc.
We represent by C [ [ t ] ] the domain of formal power series in the variable t with coefficients in the field of complex numbers C . That is, C [ [ t ] ] is the set of all the sums i = 0 a i t i , a i C . The quotient field of C [ [ t ] ] is the field of a formal Laurent series, which is denoted by C ( ( t ) ) . Every non-zero formal Laurent series A C ( ( t ) ) can be written in the form A ( t ) = t k · ( a 0 + a 1 t + a 2 t 2 + ) , where a 0 0 and k Z . Furthermore, the field C t : = n = 1 C ( ( t 1 / n ) ) is the field of a formal Puiseux series. By Puiseux’s Theorem, one gets that the field K t is algebraically closed. Observe that Puiseux series are power series with fractional exponents. Furthermore, given a Puiseux series, φ , one has a bound for the denominators of exponents with non-vanishing coefficients of φ . This bound is known as the ramification index of φ , and we represent it as ν ( φ ) (see [13]).
The order of a (Puiseux or Laurent) series A is the smallest exponent of a term with a non-vanishing coefficient in A. We represent it by ord ( A ) , and we say that the order of 0 is ∞.
Let Y ( x ) be a Puiseux series solving f ( x , y ) = 0 , ord ( Y ) > 0 , and let n be the least integer for which Y ( x ) C ( ( x 1 n ) ) (i.e., ν ( Y ) = n ). We set x 1 n = t , and then ( t n , Y ( t n ) ) is a local parametrization with center at the origin. The solutions of f ( x , y ) of order 0 are places with a center on the y-axis different from the origin. The solutions of negative order are places at infinity (that is, places with center at an infinity point).
Let Y ( x ) = i r a i x i / n be a Puiseux series with ν ( Y ) = n . The series σ ϵ ( Y ) , ϵ n = 1 are defined as the conjugates of Y, where σ ϵ ( Y ) = i r ϵ i a i x i / n . The set of all the conjugates of Y is called the conjugacy class of Y. The number of different conjugates of Y is ν ( Y ) . Two Puiseux series provide the same place if they belong to the same conjugacy class (see [13,16]).
For the case of Puiseux power series in several variables, one may use the notation introduced in [14]. More precisely, let us fix a vector of variables x = ( x 1 , , x n ) , n > 1 , and an integer d > 0 . We will use the lexicographic order l e x on R n , which can be defined as follows: since we can assign to every monomial x 1 a 1 x n a n , a i Z , i = 1 , , n , the vector ( a 1 , , a n ) Z n of its exponents, we consider the lexicographic order of the group M = { x a } { a Z n } of monomials by writing
x 1 a 1 x n a n l e x x 1 b 1 x n b n ( a 1 , , a n ) l e x ( b 1 , , b n ) .
This extension will also be called the lexocographic order for monomials, and it is a group order. The same argument follows for the group M 1 / d , d Z , d > 0 of monomials with vectors of exponents in ( 1 / d ) · Z .
Let F n , d be the set of all the functions f : ( 1 / d ) · Z C ; then F n , d is an abelian group with respect to the usual addition of functions. If we fix a vector x of variables, we may write every f F n , d as a formal sum f = a Z n f a x a / d where f a = f ( a / d ) K and, if  a = ( a 1 , , a n ) , then x a / d = x 1 a 1 / d x n a n / d . In this case, we set F n , d = F x , d . We call the support of f the set
E ( f ) = { a / d ( 1 / d ) · Z | a Z n , f a 0 } .
Finally, let us denote by K x , d the subgroup K x , d = C ( ( x n 1 / d ) ) C ( ( x 1 1 / d ) ) of F x , d , which is a field constructed by induction (see [14]). Under these conditions, if  0 f F x , d , then f K x , d if E ( f ) is a well-ordered subset of ( 1 / d ) · Z for the lexicographic order. The elements of K x , d will be called generalized Puiseux power series.
In the next definition, we introduce the concept of projective local parametrization for a projective algebraic surface.
Definition 1.
Let V * P 3 ( C ) be a projective algebraic surface defined by the homogeneous polynomial F ( x , y , z , w ) R [ x , y , z , w ] . Let A * , B * , C * , D * be series in C ( ( t 1 , t 2 ) ) such that: (i) F ( A * ( t 1 , t 2 ) : B * ( t 1 , t 2 ) : C * ( t 1 , t 2 ) : D * ( t 1 , t 2 ) ) = 0 (where the three series converge), and (ii) there is no K C ( ( t ) ) \ { 0 } such that K · ( A * , B * , C * , D * ) C 4 . Then P * = ( A * : B * : C * : D * ) P 3 ( C ( ( t 1 , t 2 ) ) ) is called a projective local parametrization of C * .
One can always find such a parametrization such that min { ord ( A * ) , ord ( B * ) , ord ( C * ) , ord ( D * ) } = 0 , and the point P * ( 0 ) V * is called the center of P * .
For a given affine surface, the previous notion can be stated as follows:
Definition 2.
Let V be a real algebraic surface over C implicitly defined by the irreducible polynomial f ( x , y , z ) R [ x , y , z ] . Let A , B , C be series in C ( ( t 1 , t 2 ) ) such that: (i) f ( A ( t 1 , t 2 ) , B ( t 1 , t 2 ) , C ( t 1 , t 2 ) ) = 0 (where the series converge), and (ii) not A and B and C, are constants. Then P = ( A , B , C ) is called an (affine) local parametrization of V . If  ord ( A ) , ord ( B ) , ord ( C ) 0 , and the point P ( 0 ) = ( a , b , c ) V is called the center of P .
In the following, we deal with affine surfaces. The results and notions presented can be adapted easily for projective algebraic surfaces.
For our purposes, we will need a generalization of the above definition. More precisely, we will be in the conditions of Definition 2, but the place P ( 0 , t 2 ) = ( a ( t 2 ) , b ( t 2 ) , c ( t 2 ) ) V , a ( t 2 ) , b ( t 2 ) , c ( t 2 ) C t 2 , will be the center of P .
Definition 3.
An equivalence class of irreducible local parametrizations of the surface V is called a place of V . The common center of the local parametrizations (if it exists) is the center of the place.
Now, we introduce the notion of a branch of a surface.
Definition 4.
Given a local parametrization ( X , Y , Z ) of a surface V , the set of all points ( X ( t 1 , t 2 ) , Y ( t 1 , t 2 ) , Z ( t 1 , t 2 ) ) obtained by allowing t 1 to vary within some neighborhood of 0 where X ( t 1 , t 2 ) and Y ( t 1 , t 2 ) and Z ( t 1 , t 2 ) converge is called a branch of V .
It can be proved that two equivalent local parametrizations provide the same branch. Hence, one gets a branch for each place of the input surface.
Furthermore, the center of a local parametrization of V is a point on V . Reciprocally, from the following theorems, we also get that every point on V is the center of at least one place of V .
Theorem 1.
Let V be a surface defined by f ( x , y , z ) R [ x , y , z ] . To each root ( X ( z ) , Y ( z ) ) C z of f ( x , y , z ) = 0 with ord ( X ) > , ord ( Y ) > 0 , there corresponds a unique place of V with a center at the origin. Conversely, to each place ( X ( t 1 , t 2 ) , Y ( t 1 , t 2 ) , Z ( t 1 , t 2 ) ) of V with a center at the origin, there correspond ord ( Z ) roots of f ( x , y , z ) = 0 , each of order greater than zero.

3. Infinity Branches

In the following, we define an infinity branch (see Definition 5), and we get some important properties which will be essential for the results we will obtain.
For this purpose, let V be an algebraic affine surface over C implicitly defined by the irreducible polynomial f ( x , y , z ) R [ x , y , z ] . Let V * be its corresponding projective surface defined by the homogeneous polynomial F ( x , y , z , w ) R [ x , y , z , w ] . In addition, let p = ( m 1 ( t 2 ) , m 2 ( t 2 ) ) , where m 1 ( t 2 ) , m 2 ( t 2 ) C t 2 , be a local parametrization of an infinity curve of V * implicitly defined by the irreducible polynomial g ( y , z , 0 ) that divides F ( 1 : y : z : 0 ) (see Definition 2). Observe that F ( 1 : m 1 ( t 2 ) : m 2 ( t 2 ) : 0 ) = 0 , and by abuse of notation, we refer to P = ( 1 : m 1 ( t 2 ) : m 2 ( t 2 ) : 0 ) as an infinity point of the input surface  V * .
We compute the series expansion for the solutions of g ( y , z , t 1 ) = 0 w.r.t ( y , z ) in some neighborhood of t 1 = 0 . We obtain solutions given by different Puiseux series that can be grouped into conjugacy classes. Let one of these solutions be given by the Puiseux series φ ( t 1 , t 2 ) = ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) ) = :
m 1 ( t 2 ) + N N 2 h 1 , N ξ N / N , m 2 ( t 2 ) + N N 2 h 2 , N ξ N / N C t 1 , t 2 2 ,
where, for  j = 1 , 2 , N = ( N j 1 , N j 2 ) , ξ N / N = t 1 N j 1 / N 1 t 2 N j 2 / N 2 , ν t i ( φ ) = N i N , i = 1 , 2 , N N 2 , and the lexicographic order for monomials is considered. We have that g ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , t 1 ) = 0 in some neighborhood of t 1 = 0 where φ ( t 1 , t 2 ) converges. Then, there exists some M R + such that
F ( 1 : φ 1 ( t 1 , t 2 ) : φ 2 ( t 1 , t 2 ) : t 1 ) = g ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , t 1 ) = 0 ,
where ( t 1 , t 2 ) C 2   and  | t 1 | < M , which implies that
F ( t 1 1 : t 1 1 φ 1 ( t 1 , t 2 ) : t 1 1 φ 2 ( t 1 , t 2 ) : 1 ) = f ( t 1 1 , t 1 1 φ 1 ( t 1 , t 2 ) , t 1 1 φ 2 ( t 1 , t 2 ) ) = 0 ,
for ( t 1 , t 2 ) C 2 and 0 < | t 1 | < M . We set t 1 1 t 1 , and we obtain that
f ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) = 0 , ( t 1 , t 2 ) C 2 a n d | t 1 | > M 1 ,
where
r j ( t 1 , t 2 ) = t 1 φ j ( t 1 1 , t 2 ) = m j ( t 2 ) t 1 + N N 2 h j , N χ N / N C t 1 , t 2 ,
N = ( N j 1 , N j 2 ) , χ N / N = t 1 1 N j 1 / N 1 t 2 N j 2 / N 2 for j = 1 , 2 .
Since ν t 1 ( φ ) = N 1 , we get that there are N 1 different series in its conjugacy class (with respect to the variable t 1 ). Let φ 11 , , φ 1 N 1 be these series, and  r 1 i ( t 1 , t 2 ) = t 1 φ 1 i ( t 1 1 , t 2 ) =
m 1 ( t 2 ) t 1 + a 1 ( t 2 ) c i N 11 t 1 1 N 11 / N 1 + a 2 ( t 2 ) c i N 21 t 1 1 N 21 / N 1 + a 3 ( t 2 ) c i N 31 t 1 1 N 31 / N 1 +
where c 1 , , c N 1 are the N 1 complex roots of x N 1 = 1 . Similarly, one reasons for r 2 i ( t 1 , t 2 ) . Now, we introduce the notion of an infinity branch.
Definition 5.
The set B = i = 1 N 1 L i where
L i = { ( t 1 , r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M i }
is called an infinity branch of the affine surface V . The subsets L 1 , , L N 1 are called the leaves of the infinity branch B.
Remark 1.
  • Note that, up to conjugation, an infinity branch is uniquely determined from one leaf. That is, if  B = i = 1 N 1 L i , where L i = { ( t 1 , r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M i } , and  r 1 i ( t 1 , t 2 ) = t 1 φ 1 i ( t 1 1 , t 2 ) =
    m 1 ( t 2 ) t 1 + a 1 ( t 2 ) t 1 1 N 11 / N 1 + a 2 ( t 2 ) t 1 1 N 21 / N 1 + a 3 ( t 2 ) t 1 1 N 31 / N 1 +
    then r 1 j = r 1 i , j = 1 , , N 1 , up to conjugation; i.e.,  r 1 j ( t 1 , t 2 ) = t 1 φ 1 j ( t 1 1 , t 2 ) =
    m 1 ( t 2 ) t 1 + a 1 ( t 2 ) c j N 11 t 1 1 N 11 / N 1 + a 2 ( t 2 ) c j N 21 t 1 1 N 21 / N 1 + a 3 ( t 2 ) c j N 31 t 1 1 N 31 / N 1 +
    where c j N 1 = 1 , j = 1 , , N 1 , and  N 1 , N i 1 N . Similarly, one reasons for r 2 j .
  • Let M : = max { M 1 , , M N } . In the following, we consider L i = { ( t 1 , r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } .
Let φ j i ( t 1 , t 2 ) = m j ( t 2 ) + N N 2 h j , i , N ξ N / N , j = 1 , 2 , be a series expansion for a solution of g ( φ 1 i ( t 1 , t 2 ) , φ 2 i ( t 1 , t 2 ) , t 1 ) = 0 . We consider ψ j i ( t 1 , t 2 ) : = φ j i ( t 1 N 1 , t 2 N 2 ) , and we observe that ( 1 : φ 1 i ( t 1 N 1 , t 2 N 2 ) , φ 2 i ( t 1 N 1 , t 2 N 2 ) : t 1 N 1 ) is a local projective parametrization, with a center at P = ( 1 : m 1 ( t 2 ) : m 2 ( t 2 ) : 0 ) , of the projective surface V * .
Thus, from  ψ j i ( t 1 , t 2 ) : = φ j i ( t 1 N 1 , t 2 N 2 ) , i = 1 , , N 1 , j = 1 , 2 ( φ j i are the N 1 different series in the conjugacy class of φ j i ). We obtain N 1 equivalent local projective parametrizations, ( 1 : ψ 1 i ( t 1 , t 2 ) : ψ 2 i ( t 1 , t 2 ) : t 1 N 1 ) (note that they are equivalent since φ j 1 , φ j N 1 belong to the same conjugacy class). Therefore, the leaves of B are all associated to a unique infinity place.
Conversely, from a given infinity place defined by a local projective parametrization ( 1 : φ 1 ( t 1 N 1 , t 2 N 2 ) , φ 2 ( t 1 N 1 , t 2 N 2 ) : t 1 N 1 ) , we obtain N 1 Puiseux series, φ j i ( t 1 , t 2 ) = ψ ( c i t 1 1 / N 1 t 1 1 / N 1 ) , c i N 1 = 1 , that provide different expressions r j i ( t 1 , t 2 ) = t 1 φ j i ( t 1 1 , t 2 ) , i = 1 , , N 1 , j = 1 , 2 . Hence, the infinity branch B is defined by the leaves
L j = { ( t 1 , r 1 j ( t 1 , t 2 ) , r 2 j ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , j = 1 , , N 1 .
From the previous disquisitions, we conclude that there exists a one-to-one relation between infinity places and infinity branches, and we may say that each infinity branch is associated with a unique infinity point (which is given by the center of the corresponding infinity place). Reciprocally, from the previous construction, we obtain that every infinity point has associated, at least, one infinity branch. Thus, every algebraic surface has, at least, one infinity branch. Furthermore, every algebraic surface has a finite number of branches.
The above process can be applied to an infinity point of the form ( a ( t 2 ) : b ( t 2 ) : c ( t 2 ) : 0 ) , a ( t 2 ) = 1 that provides a local parametrization ( b ( t 2 ) , c ( t 2 ) ) of an infinity curve. For the case that b ( t 2 ) = 1 , we may reason similarly by considering the surface implicitly defined by the polynomial F ( x , 1 , z , w ) . Observe that g ( p , 0 ) = 0 , where p = ( m 1 ( t 2 ) , m 2 ( t 2 ) ) , m 1 ( t 2 ) , m 2 ( t 2 ) C t 2 , is a local parametrization of an infinity curve of V * implicitly defined by the irreducible polynomial g ( x , z , 0 ) that divides F ( x : 1 : z : 0 ) . In this case, F ( m 1 ( t 2 ) : 1 : m 2 ( t 2 ) : 0 ) = 0 , and by abuse of notation, we refer to P = ( m 1 ( t 2 ) : 1 : m 2 ( t 2 ) : 0 ) as an infinity point of the input surface. In this situation, we get that there exists M R + such that
F ( φ 1 ( t 1 , t 2 ) , 1 , φ 2 ( t 1 , t 2 ) , t 1 ) = h ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , t 1 ) = 0 ,
for  ( t 1 , t 2 ) C 2   and  | t 1 | < M , where for j = 1 , 2 ,
φ j ( t 1 , t 2 ) = m j ( t 2 ) + N N 2 h j , N ξ N / N C t 1 , t 2 ,
N = ( N j 1 , N j 2 ) , ξ N / N = t 1 N j 1 / N 1 t 2 N j 2 / N 2 , ν t i ( φ j ) = N i N , N j i N , i = 1 , (and the lexicographic order for monomials is considered) is a series expansion for a solution of h ( x , z , t 1 ) = 0 with respect to ( x , z ) in some neighborhood of t 1 = 0 . We set t 1 t 1 1 , and we get that
f ( r 1 ( t 1 , t 2 ) , t 1 , r 2 ( t 1 , t 2 ) ) = 0 , ( t 1 , t 2 ) C 2 and | t 1 | > M 1 ,
where r j ( t 1 , t 2 ) = t 1 φ j ( t 1 1 , t 2 ) for j = 1 , 2 .
Thus, we obtain an infinity branch B = i = 1 N 1 L i whose leaves have the form:
L i = { ( r 1 i ( t 1 , t 2 ) , t 1 , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } .
Observe that we may apply this construction to any infinity point of the form ( a ( t 2 ) : b ( t 2 ) : c ( t 2 ) : 0 ) , b ( t 2 ) 0 .
Finally, for the case c ( t 2 ) = 1 , we may consider the surface implicitly defined by the polynomial F ( x , y , 1 , w ) . Observe that g ( p , 0 ) = 0 , where p = ( m 1 ( t 2 ) , m 2 ( t 2 ) ) , m 1 ( t 2 ) , m 2 ( t 2 ) C t 2 , is a local parametrization of an infinity curve of V * implicitly defined by the irreducible polynomial g ( x , y , 0 ) that divides F ( x : y : 1 : 0 ) . In this case, F ( m 1 ( t 2 ) : m 2 ( t 2 ) : 1 : 0 ) = 0 , and by abuse of notation, we refer to P = ( m 1 ( t 2 ) : m 2 ( t 2 ) : 1 : 0 ) as an infinity point of the input surface. In this situation, we get that there exists M R + such that
F ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , 1 , t 1 ) = h ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , t 1 ) = 0 ,
for  ( t 1 , t 2 ) C 2   and  | t 1 | < M , where for j = 1 , 2 ,
φ j ( t 1 , t 2 ) = m j ( t 2 ) + N N 2 h j , N ξ N / N C t 1 , t 2 ,
N = ( N j 1 , N j 2 ) , ξ N / N = t 1 N j 1 / N 1 t 2 N j 2 / N 2 , ν t i ( φ j ) = N i N , N j i N , i = 1 , (and the lexicographic order for monomials is considered) is a series expansion for a solution of h ( x , y , w ) = 0 . We set t 1 t 1 1 and we get that
f ( r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) , t 1 ) = 0 , ( t 1 , t 2 ) C 2 a n d | t 1 | > M 1 , w h e r e
r j ( t 1 , t 2 ) = t 1 φ j ( t 1 1 , t 2 ) for j = 1 , 2 .
Thus, we obtain an infinity branch B = i = 1 N 1 L i whose leaves have the form:
L i = { ( r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) , t 1 ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } .
Observe that we may apply this construction to any infinity point of the form ( a ( t 2 ) : b ( t 2 ) : c ( t 2 ) : 0 ) , c ( t 2 ) 0 .
Definition 6.
Let V be an affine surface over C defined by an irreducible polynomial f ( x , y , z ) R [ x , y , z ] .
  • An infinity branch of V of type 1 associated with the infinity point P = ( 1 : m 1 : m 2 : 0 ) , m 1 , m 2 C t 2 is a set B = i = 1 N 1 L i , where L i = { ( t 1 , r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } i = 1 , , N 1 , M R + , and  r 11 , , r 1 N 1 are the conjugates of
    r 1 ( t 1 , t 2 ) = m 1 ( t 2 ) t 1 + N N 2 h N χ N / N C t 1 , t 2 ,
    N = ( N 11 , N 12 ) , χ N / N = t 1 1 N 11 / N 1 t 2 N 12 / N 2 (similarly for r 2 );
  • An infinity branch of V of type 2 associated with the infinity point P = ( m 1 : 1 : m 1 : 0 ) , m 1 , m 2 C t 2 is a set B = i = 1 N 1 L i , where L i = { ( r 1 i ( t 1 , t 2 ) , t 1 , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } i = 1 , , N 1 , M R + , and  r 11 , , r 1 N 1 are the conjugates of
    r 1 ( t 1 , t 2 ) = m 1 ( t 2 ) t 1 + N N 2 h N χ N / N C t 1 , t 2 ,
    N = ( N 11 , N 12 ) , χ N / N = t 1 1 N 11 / N 1 t 2 N 12 / N 2 (similarly for r 2 );
  • An infinity branch of V of type 3 associated with the infinity point P = ( m 1 : m 2 : 1 : 0 ) , m 1 , m 2 C t 2 is a set B = i = 1 N 1 L i , where L i = { ( r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) , t 1 ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } i = 1 , , N 1 , M R + , and  r 11 , , r 1 N 1 are the conjugates of
    r 1 ( t 1 , t 2 ) = m 1 ( t 2 ) t 1 + N N 2 h N χ N / N C t 1 , t 2 ,
    N = ( N 11 , N 12 ) , χ N / N = t 1 1 N 11 / N 1 t 2 N 12 / N 2 (similarly for r 2 ).
Remark 2.
  • In the following, we work with the type 1 infinity branches of a given algebraic surface V . Similarly, one can reason for the other infinity branches;
  • We will say that N 1 is the ramification index of the branch B with respect to t 1 , and we will write it as ν t 1 ( B ) = N 1 . Note that B has ν t 1 ( B ) leaves.
In the following examples, we compute the infinity branches for some given surfaces.
Example 1.
Let V be the a surface implicitly defined by the irreducible polynomial
f ( x , y , z ) = x 2 + y 2 z 2 1 R [ x , y , z ] .
The corresponding projective surface V * is defined by
F ( x : y : z : w ) = x 2 + y 2 z 2 w 2 R [ x , y , z , w ] .
Note that P = ( 1 : 1 / 2 ( t 2 1 ) ( t 2 + 1 ) / t 2 : 1 / 2 ( 1 + t 2 2 ) / t 2 : 0 ) is an infinity point of  V * .
We determine the infinity branches associated with P. For this purpose, we consider the curve defined by the irreducible polynomial g ( y , z , w ) = F ( 1 : y : z : w ) , and we observe that g ( p , 0 ) = 0 , where p = ( 1 / 2 ( t 2 1 ) ( t 2 + 1 ) / t 2 , 1 / 2 ( 1 + t 2 2 ) / t 2 ) is a rational parametrization of the curve defined implicitly by g ( y , z , 0 ) = F ( 1 : y : z : 0 ) = 1 + y 2 z 2 . Note that in this case, we have more than a local parametrization of g ( y , z , 0 ) but a rational parametrization of g ( y , z , 0 ) .
Now, we compute the series expansion for the solutions of g ( y , z , t 1 ) = 0 with respect to ( y , z ) around t 1 = 0 . In this case, since g ( y , z , t 1 ) is rational over C ( t 1 ) , we compute a parametrization and we get that:
φ 1 ( t 1 , t 2 ) = 1 / 2 ( 1 t 1 2 t 2 2 ) / t 2 C ( t 1 , t 2 ) C t 1 , t 2 , a n d
φ 2 ( t 1 , t 2 ) = 1 / 2 ( 1 t 1 2 + t 2 2 ) / t 2 C ( t 1 , t 2 ) C t 1 , t 2 .
Observe that p = ( φ 1 ( 0 , t 2 ) , φ 2 ( 0 , t 2 ) ) , and  g ( φ 1 ( t 1 , t 2 ) , φ 2 ( t 1 , t 2 ) , t 1 ) = 0 . Note that ν t 1 ( φ j ) = 1 , j = 1 , 2 , which implies that we only have one Puiseux series in the conjugacy class of φ j , j = 1 , 2 . Thus, we obtain one infinity branch:
B 1 = L 1 = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
r 1 ( t 1 , t 2 ) = t 1 φ 1 ( t 1 1 , t 2 ) = 1 / 2 t 1 t 2 1 1 / 2 t 1 1 t 2 1 1 / 2 t 1 t 2 ,
r 2 ( t 1 , t 2 ) = t 1 φ 2 ( t 1 1 , t 2 ) = 1 / 2 t 1 t 2 1 1 / 2 t 1 1 t 2 1 + 1 / 2 t 1 t 2 .
Example 2.
Let V be the surface implicitly defined by the irreducible polynomial
f ( x , y , z ) = x 2 + z 2 x 2 + z y 3 R [ x , y , z ] .
The corresponding projective surface V * is defined by
F ( x : y : z : w ) = x 2 w 2 + z 2 x 2 + z y 3 R [ x , y , z , w ] .
Note that P 1 = ( 1 : t 2 : t 2 3 : 0 ) and P 2 = ( 1 : t 2 : 0 : 0 ) are the two infinity points of V * . We determine the infinity branches associated with P 1 and P 2 . For this purpose, we consider the curve defined by g ( y , z , w ) = F ( 1 : y : z : w ) , and we observe that g ( p 1 , 0 ) = 0 , where p 1 = ( t 2 , t 2 3 ) . Note that in this case, we have more than a local parametrization of g ( y , z , 0 ) but a rational parametrization of g ( y , z , 0 ) .
Now, we compute the series expansion for the solutions of g ( y , z , t 1 ) = 0 with respect to ( y , z ) around t 1 = 0 . In this case, g ( y , z , t 1 ) is not rational over C ( t 1 ) ; thus, we compute a local parametrization, and we get two different solutions:
  • First, we get that φ 1 ( t 1 , t 2 ) = ( φ 11 ( t 1 , t 2 ) , φ 12 ( t 1 , t 2 ) ) , where
    φ 11 ( t 1 , t 2 ) = t 2 C ( t 1 , t 2 ) C t 1 , t 2 , a n d
    φ 12 ( t 1 , t 2 ) = t 2 3 + t 1 2 / t 2 3 + t 1 4 / t 2 9 + C ( t 1 , t 2 ) C t 1 , t 2 .
    Observe that φ 1 ( 0 , t 2 ) = p 1 , and  g ( φ 1 ( t 1 , t 2 ) , t 1 ) = 0 . Note that ν t 1 ( φ 1 j ) = 1 , j = 1 , 2 , which implies that we only have one Puiseux series in the conjugacy class of φ j , j = 1 , 2 . Thus, we obtain one infinity branch:
    B 1 = L 1 = { ( t 1 , r 11 ( t 1 , t 2 ) , r 12 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r 11 ( t 1 , t 2 ) = t 1 φ 11 ( t 1 1 , t 2 ) = t 1 t 2 ,
    r 12 ( t 1 , t 2 ) = t 1 φ 12 ( t 1 1 , t 2 ) = t 1 t 2 3 + t 1 1 t 2 3 + t 1 3 t 2 9 + .
    In Figure 1, we plot the surface V and a surface, V 1 , constructed from the infinity branch B 1 that approach the input surface (see Section 4);
  • We also get that φ 2 ( t 1 , t 2 ) = ( φ 21 ( t 1 , t 2 ) , φ 22 ( t 1 , t 2 ) ) , where
    φ 21 ( t 1 , t 2 ) = t 2 C ( t 1 , t 2 ) C t 1 , t 2 , a n d
    φ 22 ( t 1 , t 2 ) = t 1 2 / t 2 3 t 1 4 / t 2 9 + C ( t 1 , t 2 ) C t 1 , t 2 .
    Observe that φ 2 ( 0 , t 2 ) = p 2 , and  g ( φ 2 ( t 1 , t 2 ) , t 1 ) = 0 . Note that ν t 1 ( φ 2 j ) = 1 , j = 1 , 2 , which implies that we only have one Puiseux series in the conjugacy class of φ j , j = 1 , 2 . Thus, we obtain one infinity branch:
    B 2 = L 2 = { ( t 1 , r 21 ( t 1 , t 2 ) , r 22 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r 21 ( t 1 , t 2 ) = t 1 φ 21 ( t 1 1 , t 2 ) = t 1 t 2 ,
    r 22 ( t 1 , t 2 ) = t 1 φ 22 ( t 1 1 , t 2 ) = t 1 1 t 2 3 t 1 3 t 2 9 + .
    In Figure 2, we plot the surface V and a surface, V 2 , constructed from the infinity branch B 2 that approach the input surface (see Section 4). In  Figure 3, we plot the surface V and the surfaces V 1 and V 2 together.
Remark 3.
The computation of φ i ( t 1 , t 2 ) C t 1 , t 2 , i = 1 , 2 , is not an easy question. In some cases, this problem can be easily solved as in Examples 1 and 2. For a more general cases, and also for the case of surfaces parametrically defined, we will deal with this question in a future work.
In the following, we prove that any point of the surface with sufficiently large coordinates belongs to some infinity branch. For this purpose, we recall that if h is a complex-valued function of a complex variable, h : C 3 C , we say that the limit of h ( z 1 , z 2 , z 3 ) as z i approaches is L, written lim z h ( z ) = L . If whenever { z n } n N is a sequence of points with lim n z n = , it holds that lim n h ( z n ) = L (see, e.g., [17] or [18]).
Lemma 1.
Let V be an algebraic surface. There exists K R + such that for every p = ( a , b , c ) V with | a | > K , it holds that p B p , where B p is an infinity branch of C .
Proof. 
Let us assume that the lemma does not hold, and we consider a sequence { K n } n N R + such that lim n K n = . Then, for every n N , there exists a point p n = ( a n , b n ( t 2 ) , c n ( t 2 ) ) V such that | a n | > K n , and  p n does not belong to any infinity branch of V .
Let P n = ( a n : b n ( t 2 ) : c n ( t 2 ) : 1 ) . Since F ( P n ) = f ( p n ) = 0 , then lim n F ( P n ) = 0 . Thus, we distinguish the following different cases:
(a)
If there exist two not-bounded monotone subsequences, { b n l ( t 2 ) / a n l } l N and { b n l ( t 2 ) / c n l ( t 2 ) } l N , we have that
lim l b n l ( t 2 ) / a n l = lim l b n l ( t 2 ) / c n l ( t 2 ) = ,
and then lim l a n l / b n l ( t 2 ) = 0 and lim l c n l ( t 2 ) / b n l ( t 2 ) = 0 . Hence, lim l F ( Q n l ) = F ( 0 : 1 : 0 : 0 ) = 0 , where Q n l = ( a n l / b n l ( t 2 ) : 1 : c n l ( t 2 ) / b n l ( t 2 ) : 1 / b n l ( t 2 ) ) , which implies that P = ( 0 : 1 : 0 : 0 ) is an infinity point of V * ;
(b)
If there exist a not-bounded monotone subsequence { b n l ( t 2 ) / a n l } l N and a bounded monotone subsequence { b n l ( t 2 ) / c n l ( t 2 ) } l N , we have that lim l b n l ( t 2 ) / a n l = and lim l b n l ( t 2 ) / c n l ( t 2 ) = m ( t 2 ) , and then lim l a n l / b n l ( t 2 ) = 0 and lim l c n l ( t 2 ) / b n l ( t 2 ) = 1 / m ( t 2 ) : = m 2 ( t 2 ) (if m ( t 2 ) 0 ). Hence,
lim l F ( Q n l ) = F ( 0 : 1 : m 2 ( t 2 ) : 0 ) = 0 ,
where Q n l = ( a n l / b n l ( t 2 ) : 1 : c n l ( t 2 ) / b n l ( t 2 ) : 1 / b n l ( t 2 ) ) , which implies that P = ( 0 : 1 : m 2 ( t 2 ) : 0 ) is an infinity point of V * .
Note that if m ( t 2 ) = 0 , then we consider lim l a n l / c n l ( t 2 ) = 0 . Hence,
lim l F ( Q n l ) = F ( 0 : 0 : 1 : 0 ) = 0 ,
where Q n l = ( a n l / c n l ( t 2 ) : b n l / c n l ( t 2 ) : 1 : 1 / c n l ( t 2 ) ) , which implies that P = ( 0 : 0 : 1 : 0 ) is an infinity point of V * ;
(c)
If there exist two bounded monotone subsequences { b n l ( t 2 ) / a n l } l N and { b n l ( t 2 ) / c n l ( t 2 ) } l N , we have that lim l b n l ( t 2 ) / a n l = m 1 ( t 2 ) and
lim l c n l ( t 2 ) / a n l = lim l ( c n l ( t 2 ) / b n l ( t 2 ) ) / ( b n l ( t 2 ) / a n l ) = m ( t 2 ) / m 1 ( t 2 ) : = m 2 ( t 2 )
(if m 1 ( t 2 ) 0 ). Thus,
lim l F ( Q n l ) = F ( 1 : m 1 ( t 2 ) : m 2 ( t 2 ) : 0 ) = 0 ,
where Q n l = ( 1 : b n l ( t 2 ) / a n l : c n l ( t 2 ) / a n l : 1 / a n l ) , which implies that P = ( 1 : m 1 ( t 2 ) : m 2 ( t 2 ) : 0 ) is an infinity point of V * .
If m 1 ( t 2 ) = 0 , that is, lim l b n l ( t 2 ) / a n l = 0 , since lim l b n l ( t 2 ) / c n l ( t 2 ) = n ( t 2 ) (we assume that n ( t 2 ) 0 ), we get that lim l c n l ( t 2 ) / a n l = 0 . Thus,
lim l F ( Q n l ) = F ( 1 : 0 : 0 : 0 ) = 0 ,
where Q n l = ( 1 : b n l ( t 2 ) / a n l : c n l ( t 2 ) / a n l : 1 / a n l ) , which implies that P = ( 1 : 0 : 0 : 0 ) is an infinity point of V * .
If m 1 ( t 2 ) = n ( t 2 ) = 0 , that is, lim l b n l ( t 2 ) / a n l = lim l b n l ( t 2 ) / c n l ( t 2 ) = 0 , then lim l c n l ( t 2 ) / a n l = n ( t 2 ) or lim l a n l / c n l ( t 2 ) = m ( t 2 ) . Let us assume that lim l c n l ( t 2 ) / a n l = n ( t 2 ) . Then, since lim l b n l ( t 2 ) / a n l = 0 , we get that
lim l F ( Q n l ) = F ( 1 : 0 : n ( t 2 ) : 0 ) = 0 ,
where Q n l = ( 1 : b n l ( t 2 ) / a n l : c n l ( t 2 ) / a n l : 1 / a n l ) , which implies that P = ( 1 : 0 : n ( t 2 ) : 0 ) is an infinity point of V * .
(d)
If there exist a bounded monotone subsequence { b n l ( t 2 ) / a n l } l N and a not-bounded monotone subsequence { b n l ( t 2 ) / c n l ( t 2 ) } l N , we have that
lim l b n l ( t 2 ) / a n l = m 1 ( t 2 )
and
lim l c n l ( t 2 ) / a n l = lim l ( c n l ( t 2 ) / b n l ( t 2 ) ) / ( b n l ( t 2 ) / a n l ) = 0 .
Thus,
lim l F ( Q n l ) = F ( 1 : m 1 ( t 2 ) : 0 : 0 ) = 0 ,
where Q n l = ( 1 : b n l ( t 2 ) / a n l : c n l ( t 2 ) / a n l : 1 / a n l ) , which implies that P = ( 1 : m 1 ( t 2 ) : 0 : 0 ) is an infinity point of V * .
Note that if lim l b n l ( t 2 ) / a n l = 0 , then since lim l c n l ( t 2 ) / b n l ( t 2 ) = 0 , we get that lim l c n l ( t 2 ) / a n l = 0 . Thus,
lim l F ( Q n l ) = F ( 1 : 0 : 0 : 0 ) = 0 ,
where Q n l = ( 1 : b n l ( t 2 ) / a n l : c n l ( t 2 ) / a n l : 1 / a n l ) , which implies that P = ( 1 : 0 : 0 : 0 ) is an infinity point of V * .
From both situations, we deduce that there exists a sequence { Q n } n N that approaches to an infinity point P as n tends to infinity; i.e., there exists M R + such that Q n P ϵ for  n M . Therefore, one may conclude that { Q n } n N , n M can be determined by a place centered at P. Thus, p n belongs to some infinity branch of V , which contradicts the hypothesis.    □
Remark 4.
Reasoning as in Lemma 1; one gets that there exists K R + satisfying that for every p = ( a , b , c ) V | b | > K , it holds that p B p , where B p is an infinity branch of V . The reasoning is similar if | c | > K .

4. Convergent Branches and Approaching Surfaces

In the following, we define convergent branches and approaching surfaces. We say that two infinity branches converge if they get closer as they tend to infinity. This notion will allow us to analyze whether two surfaces approach each other at the infinity.
The results obtained in this section will be used in Section 5. In Section 5, we present a method that compares the asymptotic behavior of two surfaces implicitly defined.
Definition 7.
Given two leaves, L = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } and L ¯ = { ( t 1 , r ¯ 1 ( t 1 , t 2 ) , r ¯ 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M ¯ } , we say that they are convergent if lim t 1 ( r ¯ i ( t 1 , t 2 ) r i ( t 1 , t 2 ) ) = 0 , for i = 1 , 2 .
Lemma 2.
Two leaves L = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } and L ¯ = { ( t 1 , r ¯ 1 ( t 1 , t 2 ) , r ¯ 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M ¯ } are convergent if and only if the monomials on the variable t 1 that have a non-negative exponent in the series r i ( t 1 , t 2 ) and r ¯ i ( t 1 , t 2 ) are the same (for i = 1 , 2 ).
Proof. 
Let
r 1 ( t 1 , t 2 ) = m 1 ( t 2 ) t 1 + N N 2 h 1 , N χ N / N C t 1 , t 2 ,
N = ( N 11 , N 12 ) , χ N / N = t 1 1 N 11 / N 1 t 2 N 12 / N 2 , and 
r ¯ 1 ( t 1 , t 2 ) = m ¯ 1 ( t 2 ) t 1 + N ¯ N 2 h ¯ 1 , N ¯ χ ¯ N ¯ / N ¯ C t 1 , t 2 ,
N ¯ = ( N ¯ 11 , N ¯ 12 ) , χ ¯ N ¯ / N ¯ = t 1 1 N ¯ 11 / N ¯ 1 t 2 N ¯ 12 / N ¯ 1 . Then,
r 1 ( t 1 , t 2 ) r ¯ 1 ( t 1 , t 2 ) = m 1 t 1 m ¯ 1 t 1 + N N 2 h 1 , N χ N / N N ¯ N 2 h ¯ 1 , N ¯ χ ¯ N ¯ / N ¯ .
Note that lim t 1 ( r ¯ 1 ( t 1 , t 2 ) r 1 ( t 1 , t 2 ) ) = 0 if and only if r ¯ 1 ( t 1 , t 2 ) r 1 ( t 1 , t 2 ) has no monomials having a non-negative exponent. This situation holds if the monomials on the variable t 1 that have a non-negative exponent in both series, r 1 and r ¯ 1 , are the same.
One reasons similarly for r 2 and r ¯ 2 .    □
Remark 5.
  • From Lemma 2, we deduce that m i ( t 2 ) = m ¯ i ( t 2 ) , i = 1 , 2 , and then, L and L ¯ are associated with the same infinity point;
  • Note that the number of monomials with regard to t 1 that have a positive exponent in both series is finite.
Definition 8.
Two infinity branches, B and B ¯ , are convergent if there exist two convergent leaves L B and L ¯ B ¯ .
Remark 6.
Statement 1 in Remark 5 implies that two convergent infinity branches are associated with the same infinity point.
Proposition 1.
Two infinity branches B and B ¯ are convergent if and only if for each leaf L B there exists a leaf L ¯ B ¯ convergent with L, and reciprocally.
Proof. 
The proof follows reasoning similarly as in Proposition 4.6 in [5].    □
Remark 7.
Two convergent infinity branches may have different ramification indexes; that is, they may have different numbers of leaves. However, n 1 N , which is obtained by simplifying the non-negative exponents in the variable t 1 , is the same in both branches. We refer to it as the degree of the infinity branch with respect to t 1 . Note that from the proof of Proposition 1, we get that two convergent infinity branches have the same degree with respect to t 1 .
Two convergent infinity branches may be contained in the same surface or they may belong to different surfaces. In this second case, we will say that those surfaces approach each other. In order to define this concept in a more formal way, we first introduce the following distance:
Definition 9.
Given an algebraic surface V over C and a point p C 3 , we define the distance from p to V as d ( p , V ) = min { d ( p , q ) : q V } .
Remark 8.
We should note that since V is a closed set, this minimum exists.
Definition 10.
Let V be an algebraic surface over C with an infinity branch B. We say that a surface V ¯ approaches V at its infinity branch B if there exists one leaf L = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } B such that lim t 1 d ( ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) , V ¯ ) = 0 .
We will show that this condition is satisfied for one leaf of B if and only if it is satisfied for every leaf of B. It will be derived as a consequence of the following theorem.
Theorem 2.
Let V be an algebraic surface over C with an infinity branch B. An algebraic surface V ¯ approaches V at B if and only if V ¯ has an infinity branch, B ¯ , such that B and B ¯ are convergent.
Proof. 
The proof follows reasoning similarly as in Theorem 4.11 in [5].    □
Remark 9.
  • Theorem 2 implies that “proximity” is a symmetric relation. More precisely, the surface V ¯ approaches the surface V at some infinity branch B if and only if V approaches V ¯ at some infinity branch B ¯ . In the following, we say that V and V ¯ approach each other or that they are approaching surfaces ;
  • From Theorem 2 and Remark 6, we get that two approaching surfaces have a common infinity point;
  • Theorem 2 and Proposition 1 imply that V ¯ approaches V at an infinity branch B if for every leaf
    L = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : t 1 , t 2 , C , | t 1 | > M } B ,
    it holds that lim z d ( ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) , V ¯ ) = 0 .
Corollary 1.
Let V be an algebraic surface with an infinity branch B. Let V ¯ 1 and V ¯ 2 be two different surfaces that approach V at B. Then, V ¯ 1 and V ¯ 2 approach each other.
Proof. 
From Theorem 2, there exist two infinity branches B 1 V ¯ 1 and B 2 V ¯ 2 , convergent with B. Thus, for each leaf L = { ( t 1 , r 1 ( t 1 , t 2 ) , r 2 ( t 1 , t 2 ) ) C 3 : t i C , | t 1 | > M } B , there exist two leaves L 1 = { ( t 1 , l 1 ( t 1 , t 2 ) , l 2 ( t 1 , t 2 ) ) C 2 : t i C , | t 1 | > M 1 } B 1 and L 2 = { ( t 1 , s 1 ( t 1 , t 2 ) , s 2 ( t 1 , t 2 ) ) C 3 : t i C , | t 1 | > M 2 } B 2 such that lim t 1 ( r 1 ( t 1 , t 2 ) l 1 ( t 1 , t 2 ) ) = 0 and lim t 1 ( r 1 ( t 1 , t 2 ) s 1 ( t 1 , t 2 ) ) = 0 . Then
| l 1 ( t 1 , t 2 ) s 1 ( t 1 , t 2 ) | | l 1 ( z ) r 1 ( t 1 , t 2 ) | + | r 1 ( t 1 , t 2 ) s 1 ( t 1 , t 2 ) | t 1 0
(one reasons similarly for r 2 ). Therefore, V ¯ 1 and V ¯ 2 approach each other.   □
In Example 3, we illustrate the above results.
Example 3.
Let V and V ¯ be two surfaces implicitly defined by the polynomials
f ( x , y , z ) = x 2 + z 2 x 2 + z y 3 R [ x , y , z ] a n d
f ¯ ( x , y , z ) = x 2 + z 2 x 2 + z y 3 + x 2 y 9 R [ x , y , z ] ,
respectively. Let us prove that V and V ¯ approach each other at the infinity branch associated with the infinity points
P 1 = ( 1 : t 2 : t 2 3 : 0 ) , P 2 = ( 1 : t 2 : 0 : 0 )
(note that both surfaces have P 1 = ( 1 : t 2 : t 2 3 : 0 ) and P 2 = ( 1 : t 2 : 0 : 0 ) as infinity points). Reasoning as in Example 2, we get that the infinity branch of V associated with P 1 is given by
B 1 = L 1 = { ( t 1 , r 11 ( t 1 , t 2 ) , r 12 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } ,
where
r 11 ( t 1 , t 2 ) = t 1 φ 11 ( t 1 1 , t 2 ) = t 1 t 2 ,
r 12 ( t 1 , t 2 ) = t 1 φ 12 ( t 1 1 , t 2 ) = t 1 t 2 3 + t 1 1 t 2 3 + t 1 3 t 2 9 + .
The infinity branch of V associated with P 2 is given by
B 2 = L 2 = { ( t 1 , r 21 ( t 1 , t 2 ) , r 22 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } ,
where
r 21 ( t 1 , t 2 ) = t 1 φ 21 ( t 1 1 , t 2 ) = t 1 t 2 ,
r 22 ( t 1 , t 2 ) = t 1 φ 22 ( t 1 1 , t 2 ) = t 1 1 t 2 3 t 1 3 t 2 9 + .
On the other hand, the infinity branch of V ¯ associated with P 1 is given by
B ¯ 1 = L ¯ 1 = { ( t 1 , r ¯ 11 ( t 1 , t 2 ) , r ¯ 12 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } ,
where
r ¯ 11 ( t 1 , t 2 ) = t 1 t 2 ,
r ¯ 12 ( t 1 , t 2 ) = t 1 t 2 3 + t 1 1 t 2 3 2 t 1 2 t 2 3 + t 1 3 t 2 9 + .
The infinity branch of V ¯ associated with P 2 is given by
B ¯ 2 = L ¯ 2 = { ( t 1 , r ¯ 21 ( t 1 , t 2 ) , r ¯ 22 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } ,
where
r ¯ 21 ( t 1 , t 2 ) = t 1 t 2 ,
r ¯ 22 ( t 1 , t 2 ) = t 1 1 t 2 3 2 t 1 2 t 2 2 t 1 2 t 2 3 + .
From Lemma 2, we conclude that both branches converge, since the terms with non-negative exponent in both series, r i j and r ¯ i j , are the same.
Remark 10.
In the above example, the surfaces V and V ¯ are approaching surfaces, since V approaches V ¯ at one of its infinity branches reciprocally.
In fact, V approaches V ¯ at all of its infinity branches and reciprocally. In this case, we say that both surfaces have the same asymptotic behavior. We focus on this special relation in the next section.

5. Asymptotic Behavior

From the results obtained previously, we obtain an algorithm that compares the asymptotic behavior of two surfaces implicitly defined. Two surfaces have the same asymptotic behavior if they approach each other at all of the infinity branches. Furthermore, we show that if two algebraic surfaces have the same asymptotic behavior, the Hausdorff distance between them is finite.
The algorithm developed, as well as the results presented in this section, provide essential tools in the frame of practical applications in computer-aided geometric design (CAGD) such as, for instance, the problem of the approximate parametrization (see Section 1).
To start with, we first introduce the following definition.
Definition 11.
We say that two algebraic surfaces, V and V ¯ , have the same asymptotic behavior if every infinity branch of V converges to another branch of V ¯ and does so reciprocally.
Remark 11.
From Theorem 2, we get that V and V ¯ have the same asymptotic behavior if V approaches V ¯ at all its infinity branches and does so reciprocally.
Now, we recall the notion of Hausdorff distance.
Definition 12.
Given a metric space ( E , d ) and two subsets A , B E \ { } , the Hausdorff distance between them is defined as:
d H ( A , B ) = max { sup x A inf y B d ( x , y ) , sup y B inf x A d ( x , y ) } .
If E = C 3 and d is the Euclidean distance, the Hausdorff distance between two surfaces V and V ¯ can be expressed as:
d H ( V , V ¯ ) = max { sup p V d ( p , V ¯ ) , sup p ¯ V ¯ d ( p ¯ , V ) } .
Proposition 2.
Let V and V ¯ be two algebraic surfaces having the same asymptotic behavior.
Then, the Hausdorff distance between them is finite.
Proof. 
The proof follows reasoning similarly as in Proposition 5.4 in [5].    □
The following algorithm allows us to compare the asymptotic behavior of two surfaces V and V ¯ .
We assume that we have prepared V and V ¯ such that by means of a suitable linear change of coordinates (the same change applied to both surfaces), x = 0 is not a curve of infinity of V * and V ¯ * .
In Example 4, we illustrate the performance of Algorithm 1.
Algorithm 1 Asymptotic Behavior.
  • Given two implicit algebraic surfaces V and V ¯ , the algorithm decides whether V and V ¯ have the same asymptotic behavior.
    1.
    Compute the infinity curves of V and V ¯ . If they are not the same, Return the surfaces do not have the same asymptotic behavior (see Remark 6). Otherwise, let P 1 , , P n be the infinity points corresponding to these infinity curves.
    2.
    For each P k : = ( 1 : m 1 k : m 2 k : 0 ) k = 1 , , n do:
    2.1.
    Compute the infinity branches of V associated to P k . Let B 1 , , B n k be these branches. For each i = 1 , , n k , let L i = { ( t 1 , r 1 i ( t 1 , t 2 ) , r 2 i ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M i } be any leaf of B i .
    2.2.
    Compute the infinity branches of V ¯ associated with P k . Let B ¯ 1 , , B ¯ l k be these branches. For each j = 1 , , l k , let L ¯ j = { ( t 1 , r 1 j ( t 1 , t 2 ) , r 2 j ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M j } be any leaf of B ¯ j .
    2.3.
    For each B i V , find B ¯ j V ¯ such that the terms with a non-negative exponent in the variable t 1 in r i and r ¯ j (for = 1 , 2 ) are the same up to conjugation. If there is not such a branch, Return the surfaces do not have the same asymptotic behavior (see Lemma 2).
    2.4.
    For each B ¯ j V ¯ , find B i V such that the terms with non-negative exponents in the variable t 1 in r i and r ¯ j are the same up to conjugation. If there is not such a branch, Return the surfaces do not have the same asymptotic behavior (see Lemma 2).
    3.
    Returnthe surfaces V and V ¯ have the same asymptotic behavior.
Example 4.
Let V and V ¯ be two surfaces implicitly defined by the polynomials
f ( x , y , z ) = x y x z y 4 + y 3 z + y z 2 z 3 + 1 , a n d
f ¯ ( x , y , z ) = x y x z + y z 2 z 3 y 4 + y 3 z + 5 + 2 x 3 y x 2 + 2 y 2 ,
respectively (see Figure 4).
We apply Algorithm 1 to decide whether V and V ¯ have the same asymptotic behavior:
  • Step 1: Compute the infinity points of V and V ¯ . We obtain that V and V ¯ have the same infinity points that correspond to the curves defined implicitly by y = 0 and y z = 0 .
    We start by analyzing the infinity branches associated with y = z :
  • Step 2.1: We reason similarly as in Example 1, and we get that the only infinity branch associated with y = z in V is given by
    B 1 = L 1 = { ( t 1 , r 11 ( t 1 , t 2 ) , r 12 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r 11 ( t 1 , t 2 ) = t 1 t 2 , r 12 ( t 1 , t 2 ) = t 1 t 2 t 1 3 t 2 3 t 1 4 t 2 4 + .
  • Step 2.2: We also have that there exists only one infinity branch associated with these curves of infinity in V ¯ . It is given by
    B ¯ 1 = L ¯ 1 = { ( t 1 , r ¯ 11 ( t 1 , t 2 ) , r ¯ 12 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r ¯ 11 ( t 1 , t 2 ) = t 1 t 2 ,
    r ¯ 12 ( t 1 , t 2 ) = t 1 t 2 + t 1 1 t 2 3 2 t 1 1 t 2 1 + t 1 1 t 2 4 + t 1 2 t 2 2 2 t 1 2 t 2 3 + .
  • Step 2.3 and Step 2.4: r 1 i ( t 1 , t 2 ) and r ¯ 1 i ( t 1 , t 2 ) , i = 1 , 2 , have the same terms with a non-negative exponent with respect to t 1 . Thus, B 1 and B ¯ 1 converge.
    In Figure 5, we plot the surfaces V 1 constructed from the infinity branch B 1 that approach the input surface V (left), we plot the surfaces V ¯ 1 constructed from the infinity branch B ¯ 1 that approach the input surface V ¯ (right).
    Now we analyze the infinity branches associated with y = 0 :
  • Step 2.1: Reasoning as in Example 1, we get that the only infinity branch associated with y = 0 in V is given by
    B 2 = L 2 = { ( t 1 , r 21 ( t 1 , t 2 ) , r 22 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r 21 ( t 1 , t 2 ) = t 1 t 2 ,
    r 22 ( t 1 , t 2 ) = t 2 3 / 2 t 1 3 / 2 1 / 2 t 2 3 / 2 t 1 1 / 2 1 / 8 t 1 5 / 2 t 2 9 / 5 + 1 / 2 t 1 3 t 2 3 + .
  • Step 2.2: The only infinity branch associated with y = 0 in V ¯ is given by
    B ¯ 2 = L ¯ 2 = { ( t 1 , r ¯ 21 ( t 1 , t 2 ) , r ¯ 22 ( t 1 , t 2 ) ) C 3 : ( t 1 , t 2 ) C 2 , | t 1 | > M } , w h e r e
    r ¯ 21 ( t 1 , t 2 ) = t 1 t 2 ,
    r ¯ 22 ( t 1 , t 2 ) = t 2 3 / 2 t 1 3 / 2 1 / 2 t 2 3 / 2 t 1 1 / 2 1 / 2 t 1 1 t 2 3 t 1 1 t 2 1 1 / 2 t 1 3 / 2 t 2 7 / 2 + t 1 3 / 2 t 2 3 / 2 1 / 2 t 1 2 t 2 4 + .
In Figure 6, we plot the surfaces V 2 constructed from the infinity branch B 2 that approach the input surface V (left), and we plot the surfaces V ¯ 2 constructed from the infinity branch B ¯ 2 that approach the input surface V ¯ (right).
  • Step 2.3 and Step 2.4: r 2 i ( t 1 , t 2 ) and r ¯ 2 i ( t 1 , t 2 ) , i = 1 , 2 , have the same terms with non-negative exponent s with regard to t 1 . Thus, B 2 and B ¯ 2 converge.
Since every infinity branch of V converges to another branch of V ¯ , and reciprocally, the algorithm returns that V and V ¯ have the same asymptotic behavior.
In Figure 7, we plot the surfaces V and V i , i = 1 , 2 together (left), and the surfaces V ¯ and V ¯ i , i = 1 , 2 together (right).
Remark 12.
  • Once we have the infinity branches, we can compute surfaces having the same asymptotic behavior as the given surface at each of the infinity points. For this purpose, one simply has to remove the terms with negative exponents in the variable t 1 from r 1 and r 2 . However, this problem, which will be dealt with in a future work, has to be carefully analyzed. Note that if we remove terms with negative exponents in the variable t 1 , we could be removing necessary terms in the variable t 2 ;
  • If we remove the terms with negative exponents in the variable t 1 from the series r 1 and r 2 defining the branch B, we obtain r ˜ 1 , r ˜ 2 ( C t 2 ) [ t 1 ] . In a future work, we will analyze whether one may compute the surface V ˜ i defined by the local parametrization
    ( t 1 N 1 , r ˜ 1 ( t 1 N 1 , t 2 ) , r ˜ 2 ( t 1 N 1 , t 2 ) ) ( C t 2 ) [ t 1 ] .
    In this case, V ˜ and V would have the same asymptotic behavior at B;
  • For the computation of the series φ i , i = 1 , 2 , we are considering y = t 2 and solving z by using Puiseux series (see Examples 1 and 2). This is not the best solution for some surfaces, and thus, this question should be deeply analyzed in a future work (see Remark 3).

6. Conclusions

In this paper, we introduce the notion of infinity branches and approaching surfaces for the case of a given algebraic surface implicitly defined. From these notions and the obtained properties, we present an algorithm that compares the behavior at the infinity of two algebraic surfaces defined implicitly. As a consequence and taking into account the case of curves (see [5]), we prove that if two algebraic surfaces have the same asymptotic behavior, the Hausdorff distance between them is finite.
As in the case of curves, this first paper opens some important questions that should be answered. In particular, the computation of surfaces having the same asymptotic behavior as the given surface at each of the infinity branches, the definition of a perfect surface, and the properties as well as the computation of the generalized asymptotes for implicitly and parametrically defined surfaces are important points that should be analyzed in a future work (see [10]).
Finally, we should remind the reader that, as we state in Remark 12, for the computation of the series φ i , i = 1 , 2 , we are considering y = t 2 and solving z by using Puiseux series (see Examples 1 and 2). This is not the best solution for some surfaces, and thus, this question for the general case (and also for the case of surfaces defined by a rational parametrization) should be deeply analyzed in a future work (see Remark 3).

Author Contributions

E.C.-M., M.F.d.S. and S.P.-D. contributed equally to this work and they worked together through the whole paper. All authors have read and agreed to the published version of the manuscript.

Funding

Ministerio de Ciencia, Innovación y Universidades—Agencia Estatal de Investigación/ PID2020-113192GB-I00 (Mathematical Visualization: Foundations, Algorithms and Applications).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author S. Pérez-Díaz belongs to the Research Group ASYNACS (Ref. CT-CE2019/683).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoffmann, C.M.; Sendra, J.R.; Winkler, F. Parametric Algebraic Curves and Applications. J. Symb. Comput. 1997, 23, 1006. [Google Scholar]
  2. Hoschek, J.; Lasser, D. Fundamentals of Computer Aided Geometric Design; Peters, A.K., Ed.; Wellesley Ltd.: Wellesley, MA, USA, 1993. [Google Scholar]
  3. Sendra, J.R.; Winkler, F.; Perez-Diaz, S. Rational Algebraic Curves: A Computer Algebra Approach; Series: Algorithms and Computation in Mathematics; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  4. Walker, R.J. Algebraic Curves; Princeton University Press: Princeton, NJ, USA, 1950. [Google Scholar]
  5. Blasco, A.; Pérez-Díaz, S. Asymptotic behavior of an implicit algebraic plane curve. Comput. Aided Geom. Des. 2014, 31, 345–357. [Google Scholar] [CrossRef] [Green Version]
  6. Gao, B.; Chen, Y. Finding the Topology of Implicitly defined two Algebraic Plane Curves. J. Syst. Sci. Complex. 2012, 25, 362–374. [Google Scholar] [CrossRef]
  7. González-Vega, L.; Necula, I. Efficient Topology Determination of Implicitly defined Algebraic Plane Curves. Comput. Aided Geom. Des. 2002, 19, 719–743. [Google Scholar] [CrossRef]
  8. Hong, H. An Effective Method for Analyzing the Topology of Plane Real Algebraic Curves. Math. Comput. Simul. 1996, 42, 572–582. [Google Scholar] [CrossRef]
  9. Zeng, G. Computing the Asymptotes for a Real Plane Algebraic Curve. J. Algebra 2007, 316, 680–705. [Google Scholar] [CrossRef] [Green Version]
  10. Blasco, A.; Pérez-Díaz, S. Asymptotes and Perfect Curves. Comput. Aided Geom. Des. 2014, 31, 81–96. [Google Scholar] [CrossRef] [Green Version]
  11. Pérez–Díaz, S.; Sendra, J.; Sendra, J.R. Parametrization of Approximate Algebraic Surfaces by Lines. Comput. Aided Geom. Des. 2005, 22, 147–181. [Google Scholar] [CrossRef]
  12. Adrovic, D.; Verschelde, J. Computing Puiseux Series for Algebraic Surfaces. In Proceedings of the 37th International Symposium on Symbolic and Algebraic Computation (ISSAC 2012), Grenoble, France, 22–25 July 2012; pp. 20–27. [Google Scholar]
  13. Duval, D. Rational Puiseux Expansion. Compos. Math. 1989, 70, 119–154. [Google Scholar]
  14. Soto, M.J.; Vicente, J.L. The Newton Procedure for several variables. Linear Algebra Its Appl. 2011, 435, 255–269. [Google Scholar] [CrossRef] [Green Version]
  15. Stadelmeyer, P. On the Computational Complexity of Resolving Curve Singularities and Related Problems. Ph.D. Thesis, The Johannes Kepler University Linz, Linz, Austria, 2000. [Google Scholar]
  16. Verger-Gaugry, J.-L. Beta-Conjugates of Real Algebraic Numbers as Puiseux Expansions. Integers: Electronic Journal of Combinatorial Number Theory. In Proceedings of the Leiden Numeration Conference, Leiden, The Netherlands, 7–18 June 2010; Volume 11B. [Google Scholar]
  17. Ahlfors, L.V. Complex Analysis, 3rd ed.; McGraw-Hill: London, UK, 1979. [Google Scholar]
  18. Conway, J.B. Functions of One Complex Variable I. Graduate Texts in Mathematics; Springer: New York, NY, USA, 1995. [Google Scholar]
Figure 1. Surface V (left), and surface V 1 (right).
Figure 1. Surface V (left), and surface V 1 (right).
Mathematics 10 01445 g001
Figure 2. Surface V (left), and surface V 2 (right).
Figure 2. Surface V (left), and surface V 2 (right).
Mathematics 10 01445 g002
Figure 3. Surface V (left), surface V and surfaces V 1 and V 2 (right).
Figure 3. Surface V (left), surface V and surfaces V 1 and V 2 (right).
Mathematics 10 01445 g003
Figure 4. V (left) and V ¯ (right).
Figure 4. V (left) and V ¯ (right).
Mathematics 10 01445 g004
Figure 5. Surface V 1 (left), and surface V ¯ 1 (right).
Figure 5. Surface V 1 (left), and surface V ¯ 1 (right).
Mathematics 10 01445 g005
Figure 6. Surface V 2 (left), and surface V ¯ 2 (right).
Figure 6. Surface V 2 (left), and surface V ¯ 2 (right).
Mathematics 10 01445 g006
Figure 7. Surfaces V and V i , i = 1 , 2 (left), and surfaces V ¯ and V ¯ i , i = 1 , 2 (right).
Figure 7. Surfaces V and V i , i = 1 , 2 (left), and surfaces V ¯ and V ¯ i , i = 1 , 2 (right).
Mathematics 10 01445 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campo-Montalvo, E.; Fernández de Sevilla, M.; Pérez-Díaz, S. Asymptotic Behavior of a Surface Implicitly Defined. Mathematics 2022, 10, 1445. https://doi.org/10.3390/math10091445

AMA Style

Campo-Montalvo E, Fernández de Sevilla M, Pérez-Díaz S. Asymptotic Behavior of a Surface Implicitly Defined. Mathematics. 2022; 10(9):1445. https://doi.org/10.3390/math10091445

Chicago/Turabian Style

Campo-Montalvo, Elena, Marián Fernández de Sevilla, and Sonia Pérez-Díaz. 2022. "Asymptotic Behavior of a Surface Implicitly Defined" Mathematics 10, no. 9: 1445. https://doi.org/10.3390/math10091445

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