FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE
PARAMETERIZATIONS

arXiv:1310.4915v1 [] 18 Oct 2013

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN
Abstract. Given a birational parameterization φ of an algebraic surface S
in the projective space P3 , the purpose of this paper is to investigate the
sets of points on S whose preimage consists in k or more points, counting
multiplicities. They are described explicitly in terms of Fitting ideals of some
graded parts of the symmetric algebra associated to the parameterization φ.

1. Introduction
Parameterized algebraic surfaces are ubiquitous in geometric modeling because
they are used to describe the boundary of 3-dimensional shapes. To manipulate
them, it is very useful to have an implicit representation in addition to their given
parametric representation. Indeed, a parametric representation is for instance well
adapted for drawing or sampling whereas an implicit representation allows significant improvements in intersection problems that are fundamental operations
appearing in geometric processing for visualization, analysis and manufacturing of
geometric models. Thus, there exists a rich literature on the change of representation from a parametric to an implicit representation under the classical form of a
polynomial implicit equation. Although this problem can always be solved in principle, for instance via Gröbner basis computations, its practical implementation is
not enough efficient to be useful in practical applications in geometric modeling for
general parameterized surfaces.
In order to overcome this difficulty, alternative implicit representations of parameterized surfaces under the form of a matrix have been considered. The first
family of such representations comes from the resultant theory that produces a
non-singular matrix whose determinant yields an implicit equation from a given
surface parameterization. But the main advantage is also the main drawback of
these resultant matrices: since they are universal with respect to the coefficients
of the given surface parameterization, they are very easy to build in practice, but
they are also very sensitive to the presence of base points. As a consequence, a
particular resultant matrix has to be designed for each given particular class of
parameterized surfaces. A second family of implicit matrix representations is based
on the syzygies of the coordinates of a surface parameterization. Initiated by the
geometric modeling community [SC95], these matrices have been deeply explored
in a series of papers (see [BJ03, BC05, BDD09, Bot11] and the references therein).
Compared to the resultant matrices, they are still very easy to build, although
not universal, but their sensitivity to the presence of base points is much weaker.
Date: October 21, 2013.
Nicolás Botbol was partially supported by Marie-Curie Network “SAGA” FP7 contract PITNGA-2008-214584, EU.
1

2

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

However, these matrices are in general singular matrices and the recovering of the
classical implicit equation from them is more involved. Therefore, to be useful these
matrices have to be seen as implicit representations on their own, without relying
on the more classical implicit equation. In this spirit, the use of these singular
matrix representations has recently been explored for addressing the curve/surface
and surface/surface intersection problems [BLB12, LBBM09].
As a continuation, the purpose of this paper is to investigate the self-intersection
locus of a surface parameterization through its matrix representations. More precisely, let φ be a parameterization from P2 to P3 of a rational algebraic surface S
and let M (φ) be one of its matrix representations. The main result of this paper is
that the drop of rank of M (φ) at a given point P ∈ P3 is in relation with the fiber
of the graph of φ over P . Thus, the Fitting ideals attached to M (φ) provide a filtration of the surface which is in correspondence with the degree and the dimension
of the fibers of the graph of the parameterization φ.
It turns out that this kind of results have already been investigated for different purposes in the field of intersection theory under the name of multiplepoint formulas. Given a finite map of schemes ϕ : X → Y of codimension
one, its source and target double-point cycles have been extensively studied (see
e.g. [KLU96, KLU92, Pie78, Tei77]). Moreover, in [MP89] Fitting ideals are used
to give a scheme structure to the multiple-point loci. Therefore, in the particular
case where φ : P2 → P3 has no base point (hence is a finite map), our results are
partially contained in the above literature. Thus, the main contribution of this
paper is to extend these previous works to the case where φ is not necessarily a
finite map under ubiquitous conditions. The theory in this article is thus applicable
to most of the parameterizations that appear in applications.
In what follows, we will first briefly overview matrix representations and define
some Fitting ideals attached to them in Section 2. The main result of this paper
is then proved in Section 3. Section 4 is devoted to the computational aspects of
our results. In Section 5 we will discuss on the link with multiple-points formulas developed in the field of intersection theory. Finally, we will treat the case of
parameterizations φ whose source is P1 × P1 instead of P2 . This type of parameterizations being widely used in geometric modeling (see for instance [SSV12] for a
detailed study of a special case).
2. Fitting ideals associated to surface parameterizations
Let k be a field and φ : P2k 99K P3k be a rational map given by four homogeneous
polynomials f0 , f1 , f2 , f3 of degree d. We will denote by S the closure of its image.
Set s := (s0 , s1 , s2 ), x := (x0 , x1 , x2 , x3 ), and let S := k[s] and R := k[x] be
the polynomial rings defining P2k and P3k respectively. Finally, define the ideal
I = (f0 , f1 , f2 , f3 ) ⊂ Sd . In the sequel, we will assume that:

− S is a surface in P3k ,

− φ is a birational map onto S,
(H)

− V (I) is finite and locally a complete intersection.
2.1. Matrix representations of φ. The approximation complex of cycles associated to the sequence of polynomials f0 , f1 , f2 , f3 is of the form (see e.g. [BJ03, §4]

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

3

for an introduction to these complexes in this context)
Z• :

M (φ)

0 → Z3 → Z2 → Z1 −→ Z0 → 0.
(1)
P3
The map M (φ) sends a 4-tuple (g0 , g1 , g2 , g3 ) to the form
i=0 gi xi , with the
P3
condition that (g0 , g1 , g2 , g3 ) ∈ Z1 if and only if i=0 gi fi = 0. The complex Z•
inherits a grading from the canonical grading of the polynomial ring S. Thus, the
notation M (φ)ν stands for the matrix in (Z• )ν that corresponds to the degree ν
part of M (φ).
Under our hypothesis, the symmetric algebra SI := SymS (I) of the ideal I ⊂ S
is projectively isomorphic to the Rees algebra RI := ReesS (I). In other words,
the graph Γ ⊂ P2k × P3k of the parameterization φ is scheme defined by SI . As
a consequence, the approximation complex of cycles Z• can be used to determine
an implicit representation of the closed image of φ, i.e. the surface S (recall that
H0 (Z• ) = SymS (I)).
Definition 1. Under the assumptions (H), for all integer ν > ν0 the matrix M (φ)ν
of the graded map of free R-modules (Z1 )ν → (Z0 )ν is called a matrix representation
of φ.
For the moment, the integer ν0 can be taken equal to 2d − 2, but we will come
back to its definition in the next section. Here are the main features of the collection
of matrix representations M (φ)ν of φ, for ν > ν0 :
• Its entries are
 linear forms in R = k[x0 , . . . , x3 ],
• it has ν+2
rows and at least as many columns as rows,
2

• it has maximal rank ν+2
at the generic point of P3k , that is to say
2


ν+2
rank M (φ)ν ⊗R Frac(R) =
,
2
• when specializing M (φ)ν at a point P ∈ P3 , its rank drops if and only if
P ∈ S.
2.2. Fitting ideals associated to φ. From the properties of matrix representations of φ, for all ν > ν0 the support of the Fitting ideal F0 ((SI )ν ) is S and hence
provides a scheme structure for the closure of the image of φ. Following [Tei77]
(and [EH00, V.1.3]), we call it the Fitting image of φ.
Remark 2. Observe that by definition, the ideals F0 ((SI )ν ) depend on the integer
ν (it is generated in degree ν) whereas annR ((SI )ν ) = ker(φ] ) for all ν > ν0 (see
[BJ03, Proposition 5.1] and [BC05, Theorem 4.1]), where φ] : R → S is the map of
rings induced by φ (it sends each xi to fi (s) for all i = 0, . . . , 3) and ker(φ] ) is the
usual ideal definition of the image of φ.
In this paper, we will push further the study of F0 ((SI )ν ) by looking at the other
Fitting ideals Fi ((SI )ν ), i > 0, since they provide a natural stratification :
ν+2
F0 ((SI )ν ) ⊂ F1 ((SI )ν ) ⊂ F2 ((SI )ν ) ⊂ · · · ⊂ F( 2 ) ((SI )ν ) = R.

As we will see, these Fitting ideals are closely related to the geometric properties of
the parameterization φ. For simplicity, the Fitting ideals Fi ((SI )ν ) will be denoted

Fiν (φ). We recall that Fiν (φ) ⊂ R is generated by all the minors of size ν+2
− i of
2
any matrix representation M (φ)ν .

4

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

Example 3. Consider the following parameterization of the sphere
φ : P2C

99K

P3C

(s0 : s1 : s2 )

7→

(s20 + s21 + s22 : 2s0 s2 : 2s0 s1 : s20 − s21 − s22 ).

Its matrix representations M (φ)ν have the expected properties for all ν > 1 (see
the next section). The computation of the smallest such matrix yields


0
X1
X2
−X0 + X3
.
0
−X0 − X3
X2
M (φ)1 =  X1
−X2 −X0 − X3
0
X1
Here is a Macaulay2 [GS] code computing this matrix (and others by tuning its
inputs: the parameterization φ and the integer ν):
>A=QQ[s0,s1,s2];
>f0=s0^2+s1^2+s2^2; f1=2*s0*s2; f2=2*s0*s1; f3=s0^2-s1^2-s2^2;
>F=matrix{{f0,f1,f2,f3}};
>Z1=kernel koszul(1,F);
>R=A[X0,X1,X2,X3];
>d=(degree f0)_0; nu=2*(d-1)-1;
>Z1nu=super basis(nu+d,Z1);
>Xnu=matrix{{X0,X1,X2,X3}}*substitute(Z1nu,R);
>Bnu=substitute(basis(nu,A),R);
>(m,M)=coefficients(Xnu,Variables=>{s0_R,s1_R,s2_R},Monomials=>Bnu);
>M -- this is the matrix representation in degree nu

A primary decomposition of the 3 × 3 minors of M (φ)1 , i.e. F01 (φ), returns
(X02 − X12 − X22 − X32 ) ∩ (X2 , X1 , X02 + 2X0 X3 + X32 )
which corresponds to the implicit equation of the sphere plus one embedded double
point (1 : 0 : 0 : −1). Now, a primary decomposition of the 2 × 2 minors of M (φ)1 ,
i.e. F11 (φ), returns
(X2 , X1 , X0 + X3 ) ∩ (X3 , X22 , X1 X2 , X0 X2 , X12 , X0 X1 , X02 )
which corresponds to the same embedded point (1 : 0 : 0 : −1), now with multiplicity one, plus an additional component supported at the origin. Finally, the
ideal of 1-minors of M (φ)1 , i.e. F21 (φ), is supported at the origin (i.e. is empty as a
subscheme of P2 ).
The point (1 : 0 : 0 : −1) is actually a singular point of the parameterization
φ (but not of the sphere itself). Indeed, the line L = (0 : s1 : s2 ) is a P1 that is
mapped to (s21 + s22 : 0 : 0 : −(s21 + s22 )). In particular, the base points of φ, (0 : 1 : i)
and (0 : 1 : −i), are lying on this line, and the rest of the points are mapped to
the point (1 : 0 : 0 : −1). Outside L at the source and P at the target, φ is an
isomorphism.
2.3. Regularity of the symmetric algebra. Hereafter, we give an upper bound
for the regularity of the symmetric algebra SI in our setting. We will use in the
course of this paper. We begin with some classical notation.
Given a finitely generated graded S-module N , we will denote by HFN its Hilbert
function and by HPN its Hilbert polynomial. Recall that for every µ,
HFN (µ) := dimk (Nµ ),

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

5

and
HPN (µ) := HFN (µ) −

X
i (N ) (µ)
(−1)i HFHm

(2)

i

i (N ) (µ)
where m denotes the ideal m = (s0 , s1 , s2 ) ⊂ S. It is well-known that HFHm
is finite for any i and µ, and that HPN is a polynomial function.
We also recall the definition of three classical invariants attached to a graded
module N . We will denote

i
ai (N ) := sup{µ | Hm
(N )µ 6= 0}, a∗ (N ) := max{ai (N )},
i

and
reg(N ) := max{ai (N ) + i}
i

for the ai -invariant, the a∗ -invariant and the Castelnuovo-Mumford regularity of N ,
respectively. Finally, the notation µ(J) for a given ideal J stands for its minimal
number of generators.
The following lemma will be useful in some places (see [Sza08, Corollary A12]
for a self contained elementary proof) :
Lemma 4. Let J be a graded ideal in S, generated in degree d such that dim(S/J) 6
1 and µ(Jp ) 6 3 for every prime ideal m ) p ⊃ J. Then,
reg(S/J sat ) 6 2d − 3

and
reg(S/J) 6 3d − 3 − indeg(J sat )

unless µ(J) = 2, in which case J = J sat and reg(S/J) = 2d − 2.
Proof. We may assume that k is infinite. Let f and g be two general forms of degree
d in J. They form a complete intersection. Hence, by liaison (see e.g. [CU02, 4.1
(a)])
reg(S/J sat ) = reg(S/(f, g)) − indeg(((f, g) : J sat )/(f, g)).
Hence reg(S/J sat ) 6 2d − 3, unless (f, g) ⊇ J sat in which case J = (f, g).
Let J 0 ⊆ J be an ideal generated by 3 general forms of degree d. As µ(Jp ) 6 3 for
every prime ideal m ) p, J 0 and J have the same saturation. Therefore, it suffices
0
(S/J 0 )µ+d = 0. Observe that by Koszul duality (see for instance
to show that Hm
0
0
[Cha04, Lemma 5.8]) Hm
(S/J 0 )µ+d = 0 is equivalent to Hm
(S/J 0 )2d−3−µ = 0. But
since 2d − 3 − µ < d we indeed have
0
sat
Hm
(S/J 0 )2d−3−µ = (J 0 )sat
2d−3−µ = J2d−3−µ = 0

since 2d − µ − 3 < indeg(J sat ).



Proposition 5. Assume that dim(S/I) 6 1 and µ(Ip ) 6 3 for every prime ideal
m ) p ⊃ I, then Z• is acyclic, a∗ (SI ) + 1 6 ν0 := 2d − 2 − indeg(I sat ), and
reg(SI ) 6 ν0 ,
unless I is a complete intersection of two forms of degree d∗, in which case a∗ (SI ) =
d − 3 and reg(SI ) = d − 1 = ν0 + 1.
∗

Notice that in this case the closed image of φ cannot be a surface.

6

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

i
Proof. The acyclicity of Z• is proved in [BC05, Lemma 4.2]. Furthermore, Hm
(SI ) =
i
0 for i ∈
/ {0, 1, 2, 3}, so we only need to examine Hm (SI ) for i ∈ {0, 1, 2, 3}. The
case i = 0 is treated in [BC05, Theorem 4.1]. The proof uses the spectral sequence
j
i
Hm
(Zi+j )µ+(i+j)d ⊗k R ⇒ Hm
(SI )µ

(3)

from which we will also deduce the vanishing for i > 0.

2
3
For i = 1, let µ > ν0 . We will show that Hm
(Z1 )µ+d = Hm
(Z2 )µ+2d = 0 as by
1
1
1
(3) it implies that Hm (SI )µ = 0 as claimed (recall that Hm (Z0 ) = Hm
(S) = 0).
3
∗
By [BC05, Lemma 4.1], Hm (Z2 ) ' (Z1 ) [3 − 4d] and we deduce that
3
Hm
(Z2 )µ+2d ' ((Z1 )∗ [3 − 4d])µ+2d ' (Z1 )∗µ−2d+3 .

But 2d − µ − 3 < indeg(I sat ) 6 d as µ > ν0 . Hence (Z1 )2d−µ−3 = 0 since Z1 ⊂
3
S(−d)4 , proving that Hm
(Z2 )µ+2d = 0.
The exact sequence 0 → Z1 → S(−d)4 → S → S/I → 0 provides a graded
0
2
(S/I). By Lemma 4, reg(S/I) 6 3d − 3 − indeg(I sat ) =
(Z1 ) ' Hm
isomorphism Hm
2
d + ν0 − 1, hence Hm (Z1 )µ+d = 0.
For i = 2, consider the canonical exact sequence defining K
0 → K → SI → RI → 0.

As K is the direct sum of the S-modules ker(SymtS (I) → I t ) that are supported
i
2
2
on V (I), we have Hm
(K) = 0 for i > dim S/I. It follows that Hm
(SI ) ' Hm
(RI ).
2
2 t
1
t
Now, a (RI ) = maxt>1 {a (I ) − td} = maxt>1 {a (S/I ) − td}. The surjective map
t+3
(S/I)(−td)( 3 ) → I t /I t+1 provides the inequality a1 (I t /I t+1 ) 6 a1 (S/I) + td, and
it remains to estimate a1 (S/I). By Lemma 4,
a2 (RI ) 6 a1 (S/I) − d 6 d − 4,

unless µ(I) = 2. As indeg(I sat ) 6 d it follows that a2 (RI ) 6 ν0 − 2, unless
I = (f, g), in which case a2 (RI ) = ν0 − 1.

3
3
(S)µ = 0 for all µ > −3, hence a3 (SI ) 6 −3. Finally,
(Z0 )µ = Hm
For i = 3, Hm
3

a (SI ) = ⊕t>0 Hm (S)(td), hence a3 (SI ) = −3.
3

Remark 6. If V (I) = ∅ then a2 (SI ) = −∞ and a1 (SI ) = a1 (I) − d = reg(S/I) − d.
Also, when dim(S/I) = 1 and V (I) is locally a complete intersection, one has
a2 (SI ) = a2 (I) − d as above and a1 (SI ) = max{a1 (I) − d, a1 (I 2 ) − 2d}. This shows
in particular that a1 (SI ) and a2 (SI ) can be (very quickly) computed from I, using
a dedicated software.
3. Fibers of the canonical projection of the graph of φ onto S

Let Γ ⊂ P2P3 be the graph of φ and R+ = (X0 , . . . , X3 ) the graded maximal ideal
of R. We have the following diagram

/ P2 3
Γ
(4)
P
k

π1


P2k

π2



P3k

where all the maps are canonical. For any p ∈ Proj(R), we will denote by κ(p) its
residue field Rp /pRp . The fiber of π2 at p ∈ Proj(R) is the subscheme
π2−1 (p) := Proj(SI ⊗R κ(p)) ⊂ P2κ(p) .

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

7

For simplicity, we set Nνp := dimκ(p) (SI ⊗R κ(p))ν for the Hilbert function of the
fiber π2−1 (p) in degree ν. It turns out that this quantity governs the support of the
Fitting ideals Fiν (φ) we are interested in.
Proposition 7. For any p ∈ Proj(R)
(
Fiν (φ) ⊗R κ(p) =

0
κ(p)

if 0 6 i < Nνp
if Nνp 6 i

Proof. First, as a consequence of the stability under base change of the Fitting
ideals, one has
Fi ((SI )ν ) ⊗R κ(p) = Fi ((SI )ν ⊗R κ(p)).
Now, for any ν > 0 the κ(p)-vector space
Nνp .

(SI )ν ⊗R κ(p) = (SI ⊗R κ(p))ν

The result then follows from the fact that the i-th Fitting
is of dimension
ideals of a κ(p)-vector space of dimension N is 0 for i < N and κ(p) for i > N . 
This leads us to focus on the behavior of the Hilbert functions Nνp . By (2), it
suffices to control the local cohomology of the fiber SI ⊗ κ(p) with respect to the
ideal m. Part of this control is already known,
Lemma 8 ([Cha12, Prop. 6.3]). Let p ∈ Proj(R) and vp := dim(SI ⊗R κ(p)). Then,
avp (SI ⊗R κ(p)) = avp (SI ⊗R Rp ) 6 avp (SI )

and
avp −1 (SI ⊗R κ(p)) 6 max{avp −1 (SI ⊗R Rp ), avp (SI ⊗R Rp )} 6 max{avp −1 (SI ), avp (SI )}.
We are now ready to describe the support of the ideals Fiν (φ) with i > 0 and
ν > ν0 . We begin with the points on S whose fiber is 0-dimensional.
Theorem 9 (0-dimensional fibers). Let p ∈ Proj(R). If dim π2−1 (p) = 0 then for
all ν > ν0
p ∈ V (Fiν (φ)) ⇔ deg(π2−1 (p)) > i + 1.
(5)
i
Proof. By Proposition 5 and Lemma 8, we have Hm
(SI ⊗R κ(p))ν = 0 for all ν > ν0
and i > 0. Then, the conclusion follows from Proposition 7 and the equality (2)
since the Hilbert polynomial of the fiber π2−1 (p) is a constant polynomial which is
equal to the degree of this fiber.


We now turn to points on S whose fiber is 1-dimensional.
Lemma 10. Assume that the fi ’s are linearly independent, the fiber over a closed
point p of coordinates (p0 : p1 : p2 : p3 ) is of dimension 1, and its unmixed component is defined by hp ∈ S. Let `p be a linear form with `p (p0 , . . . , p3 ) = 1 and set
`i (x) := xi − pi `p (x). Then, hp = gcd(`0 (f ), . . . , `3 (f )) and
I = (`p (f )) + hp (g0 , . . . , g3 )
with `i (f ) = hp gi and `p (g0 , . . . , g3 ) = 0. In particular
(`p (f )) + hp (g0 , . . . , g3 )sat ⊆ I sat ⊆ (`p (f )) + (hp ).

8

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

P
P
Proof. A syzygy L = i ai xi = i ai (pi `p (x) + `i (x)) provides an equation for the
P
fiber : L = i ai pi . Recall that hp = gcd(L1 , . . . , Lt ) where the Lj are generators
of the syzygies of I.
The particular syzygy L(i) := `i (f )`p (x) − `p (f )`i (x) satisfies L(i) = `i (f ). It
follows that hp divides h := gcd(`0 (f ), . . . , `3 (f )). Set `i (f ) = hgi for some gi ∈ S.
The `i ’s span the linear forms vanishing at p, and are related by the equation
`p (`0 , . . . , `3 ) = 0. It follows that
I = (`p (f ), `0 (f ), . . . , `3 (f )) = (`p (f )) + h(g0 , . . . , g3 ),
P
and `p (g0 , . . . , g3 ) = 0. Now if one has a relation a`p (f ) + i bi hgi = 0, then h
divides a, as the fi ’s have no common factor. This in turn shows that h divides hp
and completes the proof.

Notice that if pi 6= 0, gi is a linear combination of the gj ’s for j 6= i as
`p (g0 , . . . , g3 ) = 0. For instance if p0 6= 0, I = (`p (f )) + h(g1 , . . . , g3 ).
Remark 11. The above lemma shows that fibers of dimension 1 can only occur
when V (I) 6= ∅ as V (I) ⊇ V ((`p (f ), hp )). It also shows that
d deg(hp ) 6 deg(I)
if there is a fiber of dimension 1 of unmixed part given by hp . Furthermore, any
element F in I sat of degree < d has to be a multiple of hp for any fiber of dimension
1. As hp and hp0 cannot have a common factor for p 6= p0 , the product of these
forms divide F . This gives a simple method to compute 1 dimensional fibers if such
an F exists.
Theorem 12 (1-dimensional fibers). Let p ∈ Proj(R) be a k-rational closed point.
If dim π2−1 (p) = 1 then, setting δ := deg(π2−1 (p))),
a0 (SI ⊗ κ(p)) 6 ν0 − δ − 1
and
reg(SI ⊗ κ(p)) 6 max{ν0 − δ − 1, reg(I) − d} 6 ν0 − 1.

As a consequence, a∗ (SI ⊗ κ(p)) 6 ν0 − 2.

Proof. By Lemma 10, I = (f ) + h(g1 , g2 , g3 ) with h the equation of the unmixed
part of the fiber. Let I 0 := (g1 , g2 , g3 ).
We may assume for simplicity that p = (1 : 0 : 0 : 0) (i.e. f = f0 and fi = hgi
for i = 1, 2, 3), and we then have an exact sequence
0 → Z10 (d) → Z1 (d) → S → SI ⊗ κ(p) → 0,

with Z10 the syzygy module of (f1 , f2 , f3 ), the first map sending (b1 , b2 , b3 ) to
(0, b1 , b2 , b3 ) and the second (a0 , a1 , a2 , a3 ) to a0 .
It follows that reg(SI ⊗ κ(p)) 6 max{reg(S), reg(Z1 ) − d − 1, reg(Z10 ) − d − 2}
0
1
(S) = Hm
(Z1 ) = 0.
and a0 (SI ⊗ κ(p)) 6 a2 (Z10 ) − d 6 reg(Z10 ) − d − 2 as Hm
0
By Lemma 10, h = gcd(f1 , f2 , f3 ) is of degree δ. Set fi = hfi for all i = 1, 2, 3.
The syzygies of (f1 , f2 , f3 ) and I 0 := (f10 , f20 , f30 ) are equal, up to a degree shift by
δ, and the exact sequence 0 → Z10 (δ) → S(−d + δ)3 → I 0 → 0 shows that :
reg(Z10 ) = reg(I 0 ) + 1 + δ.

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

9

Notice that hI 0 ⊂ I, hence h(I 0 )sat ⊆ I sat and indeg(I sat ) 6 δ + indeg((I 0 )sat ). By
Lemma 4 we obtain
reg(I 0 ) 6 3(d − δ) − 2 − indeg((I 0 )sat ) 6 3d − 2δ − indeg(I sat ) − 2
hence
a0 (SI ⊗ κ(p)) 6 2d − δ − 3 − indeg(I sat ) = ν0 − δ − 1.
Finally, reg(Z1 ) − d − 1 = reg(I) − d 6 ν0 − 1 by Lemma 4.



To conclude the above case discussion, let us state a simple corollary :
Corollary 13. Let p ∈ Proj(R).
(i) If π2−1 (p) = ∅ then a∗ (SI ⊗ κ(p)) = reg(SI ⊗ κ(p)) 6 ν0 − 1.
(ii) If dim(π2−1 (p)) = 0 then a∗ (SI ⊗ κ(p)) 6 ν0 − 1 and reg(SI ⊗ κ(p)) 6 ν0 .
(iii) If dim(π2−1 (p)) = 1 then a∗ (SI ⊗ κ(p)) 6 ν0 − 2 and reg(SI ⊗ κ(p)) 6 ν0 − 1.
Let us now remark that there is no fiber of dimension two (i.e. equal to P2 ) over
a point p ∈ Proj(R). More generally :
Lemma 14. Let φ : Pn−1
99K Pm
k be a rational map and Γ be the closure of the
k
graph of φ. The following are equivalent : (i) π : Γ → Pm
k has a fiber Fx equal to
,
(ii)
φ
is
the
restriction
to
the
complement
of
a
hypersurface
of the constant
Pn−1
x
map sending any point to the zero dimensional and k-rational point x.
Proof. For (i)⇒(ii), notice that Pkn−1 × {x} is a subscheme of the absolutely irreducible scheme Γ that has same dimension n − 1. Hence Fx = Γ showing (ii). 
Theorem 15 (2-dimensional fibers). For all ν > ν0 ,



ν+2
=∅
p ∈ Proj(R) : HPπ−1 (p) (ν) = HPP2 (ν) =
2
2
Proof. It follows from Proposition 5 that the defining ideal of the fiber is generated
in degree at most ν0 . As this ideal
 is non trivial for any point in Proj(R) by Lemma
14, we deduce that HPνp < ν+2
for any prime in Proj(R) if ν > ν0 .

2
It is worth mentioning two facts that are direct consequences of the above results.
First, the embedded components of F0ν (φ) are exactly supported on V (F1ν (φ)) ⊂
V (F0ν (φ)) = S . Second, the set-theoretical support of each Fiν (φ), i fixed, stabilizes
for ν  0; its set-theoretical support then correspond to those points on S whose
fiber is either a curve and or a finite scheme of degree greater or equal to i + 1.
4. Computational aspects
In this section, we detail the consequences of the previous results for giving a
computational description of the singular locus of the parameterization φ of the
surface S.

10

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

96
95
94
93
92
91
90
89
88
87
86
85
84
83
82
81
80
79
78
77
76
75
74
73
72
71
70
69
68
67
66
65
64
63
62
61
60
59
58
57
56
55
54
53
52
51
50
49
48
47
46
45
44
43
42
41
40
39
38
37
36
35
34
33
32
31
30
29
28
27
26
25
24
23
22
21
20
19
18
17
16
15
14
13
12
119
10
8
7
6
5
4
3
2
1
0

HPS (ν)

HPSI ⊗R κ(p) (ν)
ν0 − 1 ν0

Figure 1. Graph of Nνp as a function of ν for a given p ∈ Proj(R).
4.1. Description of the fiber of a given point on S . Given a point p on S ,
we summarize the results we obtained with Figure 1. The blue curve corresponds
to a finite fiber, i.e. a point p such that dim π2−1 (p) = 0, the red curve corresponds a
one dimensional fiber, i.e. a point p such that dim π2−1 (p) = 1, and the black curve
represents the Hilbert polynomial of S, the coordinate ring of P2k .
The fact that p belongs or not to the support of Fiν (φ) can be checked by the
computation of the rank of a matrix representation Mν (φ) evaluated at the point
p. More precisely, we get the following properties:
• If dim π2−1 (p) 6 0 then for all ν > ν0

corank(Mν (φ)(p)) = deg(π2−1 (p)).

• If dim π2−1 (p) = 1 then for all ν > ν0 − 1

corank(Mν (φ)(p)) = deg(π2−1 (p))ν + c

where c ∈ Z is such that HPπ−1 (p) (ν) = deg(π2−1 (p))ν + c.
2

From here, we get an algorithm to determine the Hilbert polynomial of a fiber
π2−1 (p) of a given point p by comparing the rank variation of the matrix Mν (φ)(p)
for two consecutive integers ν and ν + 1, assuming ν > ν0 . Before giving a more
precise description in the form of an algorithm, we refine the distinction between
finite and non-finite fibers.
Proposition 16. Let p be a point on S and choose an integer ν such that
ν0 6 ν 6 2d − 2.

If corank(Mν (φ)(p)) 6 ν then

π2−1 (p)

is necessarily a finite fiber.

Proof. We know that the Hilbert function of the fiber π2−1 (p) coincides with its
Hilbert polynomial for all ν > ν0 . Therefore, to prove our claim it is enough to
show that if the fiber π2−1 (p) of p is of dimension 1, then
corank(Mν (φ)(p)) = HPπ−1 (p) (ν) > ν + 1.
2

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

11

So, assume that π2−1 (p) is 1-dimensional and set δ := deg(π2−1 (p)) so that


δ−1
HPπ−1 (p) (ν) = δν + 1 −
+N
2
2

where δ−1
is the arithmetic genus of the unmixed dimension one part of the fiber
2
π2−1 (p) and N > 0 is the degree of its remaining finite part. Observe that δ 6 d.
Now, a straightforward computation yields
δ(2ν + 3 − δ)
+ N.
2
From here, we see that HPπ−1 (p) (ν) is an increasing function of δ on the interval
2
[−∞; ν + 23 ] which takes value ν + 1 + N for δ = 1. Observe moreover that ν + 32 > δ
as soon as ν > d − 1. And if ν = d − 2, then HPπ−1 (p) (d − 2) takes the same value
2
for δ = d − 1 and δ = d.

HPπ−1 (p) (ν) =
2

An interesting consequence of the above proposition is that for all integer 0 6
n 6 2d − 2
corank(M2d−2 (φ)(p)) = n ⇔ π2−1 (p) is finite of degree n.

Therefore, as soon as d > 2 we have
corank(M2d−2 (φ)(p)) = 1 ⇔ p has a unique preimage.
Observe in addition that once we know that a point p has a unique preimage,
this preimage can be computed from the kernel of the transpose of M2d−2 (φ)(p).
Indeed, this kernel is generated by a single vector. But it is clear by the definition of
M2d−2 (φ) that the evaluation of the basis of Sν chosen to build M2d−2 (φ) evaluated
at the point (s0 : s1 : s2 ) is also a nonzero vector in the kernel of the transpose
of M2d−2 (φ)(p). Hence, from this property it is straightforward to compute the
pre-image of p if it is unique.
Algorithm : Determination of the characters of a fiber π2−1 (p)
Input: a parameterization φ of a surface S satisfying (H) and a point p ∈ P3 .
Output: the hilbert polynomial of the fiber π2−1 (p).
1.
2.
3.
4.
5.
6.
7.
8.
9.

Pick an integer ν > ν0 (e.g. ν = 2d − 2 which is always valid).
Compute a matrix representation Mν (φ).
Compute rν := corank(Mν (φ))(p).
If rν = 0 then p does not belong to S (and stop).
If 0 < rν 6 ν then π2−1 (p) is a finite set of rν points, counted with multiplicity.
If rν > ν then compute Mν+1 (φ).
Compute rν+1 := corank(Mν+1 (φ))(p).
If rν = rν+1 then π2−1 (p) is a finite set of rν points, counted with multiplicity.
If rν < rν+1 then dim(π2−1 (p)) = 1 and
HPπ−1 (p) (ν) = (rν+1 − rν )ν + rν − (rν+1 − rν )ν.
2

12

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

In other words, π2−1 (p) is made of a curve of degree rν+1 − rν and of a finite
set of


rν+1 − rν − 1
+ rν − (rν+1 − rν )ν − 1.
2
4.2. Pull-back to the parameter space. So far, we have seen that the fibers of
π2 can be described from the Fitting ideals Fiν (φ) which are homogeneous ideals in
the ring R of implicit variables. However, for applications in geometric modeling,
it is also interesting to get such a description of the singular locus of φ in the source
space P2 in place of the target space P3 . From a computational point of view, this
also allows to reduce the number of ambient variables by one. For that purpose,
we just have to pull-back our Fitting ideals through φ. For all integers i > 0 and
ν > ν0 we define the ideals of S
φ−1 (Fiν ) := φ] (Fiν ).S
It is clear that for all i we have φ−1 (Fiν ) ⊂ I, that is to say that the base points
of φ are always contained in the support of φ−1 (Fiν ) for all i > 0 and all ν > ν0 . By
Lemma 15, this support is exactly the base points of φ if i = ν+2
− 1. In other
2
words we have
 ν+2 
( )−1
= I, ∀ ν > ν0 .
φ−1 Fν 2
For all ν > ν0 , we clearly have φ−1 (F0ν ) = (0) ⊂ S. Notice that this behavior
is similar to substituting the parameterization φ into an implicit equation of S
which returns 0. However, φ−1 (F1ν ) is nonzero since we have assumed φ birational.
Actually, as a consequence of the results of Section 3, the sequence of ideals
 ν+2 
( )−1
−1
1
−1
2
−1
=I
φ (Fν ) ⊂ φ (Fν ) ⊂ · · · ⊂ φ
Fν 2
yields a filtration of P2 which is in correspondence with the singularities of φ. Here
is an illustrative example where we focus on the singularities of the parameterization
φ in a situation where S is smooth and very simple.
Example 17. Take four general linear forms l0 , l1 , l2 , l3 in S1 and consider the parameterization
φ : P2C → P3C : (s0 : s1 : s2 ) → (l1 l2 : l0 l3 : l0 l1 : l0 l1 ).

The closed image of φ is the plane of equation S : X2 − X3 = 0. There are three
base points that are all local complete intersections:
P01 = {l0 = 0} ∩ {l1 = 0}, P02 = {l0 = 0} ∩ {l2 = 0}, P13 = {l1 = 0} ∩ {l3 = 0}.

The fiber of each base point by π1 (see the diagram (4)) gives a line on S , the projection of the associated exceptional divisor, say D01 , D02 and D13 . The equations
of the graph Γ are l0 X0 − l2 X2 , l1 X1 − l3 X2 , (l3 − l2 )X2 , X2 − X3 so that we have
D01 : {X2 = 0} ∩ S , D02 : {l1 (P02 )X1 − l3 (P02 )X2 = 0} ∩ S ,

D13 : {l0 (P13 )X0 − l2 (P13 )X2 = 0} ∩ S .
Now, let P be a point on S . If P ∈
/ D01 ∪ D02 ∪ D13 then π1 (π2−1 (P )) does not
contain any base point and is necessarily zero-dimensional. Moreover, π1 π2−1 (D01 ∩
D13 ) = {l1 = 0}, so that we have the one-dimensional fiber over D01 ∩ D13 = (0 :

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

13

1 : 0 : 0) that goes through only two of the three base points. A similar property
appears for π1 π2−1 (D01 ∩ D02 ) = {l0 = 0}.
Notice that the three lines D01 , D02 , D13 are not in the image of φ except for the
two points: (0 : 1 : 0 : 0) and (1 : 0 : 0 : 0). Moreover, the fibers over these two
points are lines and they are the only points with a one-dimensional fiber.
5. Link with multiple-point schemes of finite maps
Given a finite map of schemes ϕ : X → Y of codimension one, its source and
target double-point cycles have been extensively studied in the field of intersection
theory (see e.g. [KLU96, KLU92, Pie78, Tei77]). In [MP89], Fitting ideals are used
to give a scheme structure to the multiple-point loci.
Under the hypothesis that φ : P2 → P3 is a finite map, our results are partially
contained in the above literature. Recall that if φ has no base point then φ is finite
(see Remark 11). The contribution of our paper is hence mainly an extension to
the case where φ is not finite. It is known that it is always possible to remove the
base points by blowing-up, and this is exactly the role of the symmetric algebra in
our approach of matrix representations, but still some one-dimensional fibers might
remain.
In loc. cit. the basic ingredient is the direct image of the structural sheaf of P2 :
ϕ∗ (OP2 ). Since ϕ is assumed to be finite, then ϕ∗ (OP2 ) is a coherent sheaf and
hence one can consider its Fitting ideal sheaves [MP89]. In our approach, we did
not considered this sheaf, but the sheafification of the graded part of degree ν of
the symmetric algebra SI (ν > ν0 ). It turns out that they are isomorphic as soon
as φ is finite. Let us be more precise.
Assume that φ is finite (e.g. φ has no base point). In the following diagram
Γ ⊂ P2k × P3k
π1


P2k

π2
φ

%/

P3k

the morphism π1|Γ : Γ → P2k is an isomorphism, so that π1] : OP2 ' π1∗ (OΓ ) by
definition. In particular, we have the sheaf isomorphisms
π2∗ (OΓ ) = φ∗ π1∗ (OΓ ) ' φ∗ (OP2 ).

(6)

Given an integer ν, we can consider two sheaves of OP3 -modules:
π2∗ (OΓ ⊗OP2 ×P3 π1∗ (OP2 (ν)))

and the sheaf associated to the graded R-module (SI )ν that we will denote by
^
(S
I )ν .
Lemma 18. For all ν > ν0 we have
^
φ∗ (OP2 ) ' (S
I )ν
Proof. The exact sequence of R-modules
0
1
0 → Hm
(SI )ν → (SI )ν → H 0 (P2Spec(R) , OΓ ⊗ π1∗ (OP2 (ν))) → Hm
(SI )ν → 0

shows that for all ν > ν0 we have (extending π2 to P2 × Spec(R) → Spec(R))

(SI )ν = H 0 (P2Spec(R) , OΓ ⊗π1∗ (OP2 (ν))) = H 0 (Spec(R), π2∗ (OΓ ⊗π1∗ (OP2 (ν)))) (7)

14

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

where the last equality follows from the definition of π2∗ (−) and the fact that
π2−1 (Spec(R)) = P2 × Spec(R). Since (7) is an equality of graded R-modules, we
deduce that
∗
^
(S
I )ν ' π2∗ (OΓ ⊗OP2 ×P3 π1 (OP2 (ν))).

Comparing this with (6) we obtain that for all ν > ν0
^
φ∗ (OP2 ) ' (S
I )ν

(because π1∗ (OP2 (ν)) is an invertible sheaf).



In particular, we get that
0
^
^
F0P3 (π1∗ (OP2 (ν))) ' F0P3 ((S
I )ν ) ' Fν (φ)

where the last isomorphism follows from the stability of Fitting ideals under base
change and the classical properties of the sheaves associated to graded modules.
6. The case of bi-graded parameterizations of surfaces
In this paper, we have treated rational surfaces that are parameterized by a
birational map from P2 to P3 . In the field of geometric modeling such surfaces
are called rational triangular Bézier surfaces (or patches). In this section, we will
(briefly) show that we can treat similarly rational surfaces that are parameterized
by a birational map from P1 × P1 to P3 ; such surfaces are called rational tensor
product surfaces in the field of geometric modeling where they are widely used.
Suppose given a rational map φ : P1k × P1k 99K P3k where f0 , f1 , f2 , f3 are four
bi-homogeneous polynomials of bi-degree d := (d1 , d2 ), that is to say that fi ∈
S := k[s0 , s1 , t0 , t1 ] is homogeneous of degree d1 with respect to the variables s0 , s1
and is homogeneous of degree d2 with respect to the variables t0 , t1 . As in Section
2, denote by S the closure of the image of φ, set x = (x0 , x1 , x2 , x3 ), R := k[x],
I = (f0 , f1 , f2 , f3 ) ⊂ Sd and assume that conditions (H) hold. In this setting,
the approximation complex of cycles can still be considered and it inherits a bigraduation from the ring S with respect to the two sets of variables s0 , s1 and t0 , t1 .
It is proved in [Bot11] that matrix representations of φ still exist in this setting.
To be more precise, define the following region (see Figure 2)
Rd := {(n1 , n2 ) ∈ Z2 such that (n1 6 d1 − 2), or (n2 6 d2 − 2), or

(n1 6 2d1 − 2 and n2 6 2d2 − 2)} ⊂ Z2 .

Then, for all ν = (ν1 , ν2 ) ∈
/ Rd the matrix M (φ)ν , which corresponds to the νgraded part of the matrix M (φ) in (1), is a matrix representation of φ with all the
expected properties, in particular (see [Bot11] for the details):
• its entries are linear forms in R = k[x0 , . . . , x3 ],
• it has (ν1 + 1)(ν2 + 1) rows and at least as many columns as rows,
• it has maximal rank (ν1 + 1)(ν2 + 1) at the generic point of P3k ,
• when specializing M (φ)ν at a point P ∈ P3 , its rank drops if and only if
P ∈ S.
Therefore, one can also consider the Fitting ideals Fiν (φ) ⊂ R associated to M (φ)ν
for which Proposition 7 holds. Notice that in this setting, the projective plane P2

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

(d1 − 2, 2d2 − 2)

15

b

(2d1 − 2, 2d2 − 2)

b

b

Rd

(2d2 − 1, d2 − 2)

Figure 2. The region Rd is represented in grey color.
is replaced by P1 × P1 in (4), so that for any point p ∈ Proj(R) the fiber of π2 at p
is a subscheme of P1κ(p) × P1κ(p) :
π2−1 (p) = BiProj(SI ⊗R κ(p)) ⊂ P1κ(p) × P1κ(p) .
Thus, the Hilbert function of a fiber π2−1 (p) is a function of Z2 . Similarly to Section
3, we set Nνp := HFSI ⊗R κ(p) (ν). It turns out that an equality similar to (2) also
holds:
X
(−1)i HFHBi (SI ⊗R κ(p)) (ν)
HPSI ⊗R κ(p) (ν) := HFSI ⊗R κ(p) (ν) −
i

where the ideal B denotes the product of ideals B := (s0 , s1 )(t0 , t1 ) ⊂ S. The
function HPSI ⊗R κ(p) is a polynomial function called the Hilbert polynomial (see
[CN08] or [BC13, Chapter 4.6] for more details). The following result allows to
control the vanishing of the Hilbert function of the local cohomology modules of
the fibers.
Theorem 19. For all p ∈ Spec(R), all integer i and ν ∈
/ Rd ,
i
HB
(SI ⊗R κ(p))ν = 0

Proof. We will only outline the proof of this result because it follows the main lines
of the proof of Corollary 13.
Assume first that δ := dim(SI ⊗R κ(p)) = 1. Lemma 8 implies that
1
0
HB
(SI ⊗R κ(p))ν = HB
(SI ⊗R κ(p))ν = 0
0
1
if HB
(SI )ν = HB
(SI )ν = 0. As proved in [Bot11, Theorem 4.7], this latter condition holds if ν ∈
/ Rd .
Now, we turn to the case δ = 2 for which we have to control the vanishing of
0
1
2
HB
(SI ⊗R κ(p))ν , HB
(SI ⊗R κ(p))ν and HB
(SI ⊗R κ(p))ν . We proceed as in the
proof of Theorem 12. We may assume for simplicity that p = (1 : 0 : 0 : 0) so that
I = (f0 ) + h(g1 , g2 , g3 ) with h the equation of the unmixed part of dimension 2 of
the fiber at p; it is a bi-homogeneous polynomial of bi-degree e = (e1 , e2 ). Setting
I 0 := (g1 , g2 , g3 ), we have an exact sequence

0 → Z10 (d) → Z1 (d) → S → SI ⊗ κ(p) → 0,

16

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

with Z10 the syzygy module of (f1 , f2 , f3 ), the first map sending (b1 , b2 , b3 ) to
(0, b1 , b2 , b3 ) and the second (a0 , a1 , a2 , a3 ) to a0 . The two spectral sequences associated to the double complex
•
•
•
0 → CB
(Z10 (d)) → CB
(Z1 (d)) → CB
(S)

show that
0
1
2
• HB
(SI ⊗R κ(p))ν = 0 if HB
(Z1 (d))ν = HB
(Z10 (d))ν = 0,
1
2
3
• HB
(SI ⊗R κ(p))ν = 0 if HB
(Z1 (d))ν = HB
(Z10 (d))ν = 0,
2
2
3
• HB (SI ⊗R κ(p))ν = 0 if HB (S)ν = HB (Z1 (d))ν = 0.

i
3
Recall that HB
(S) = 0 if i 6= 2, 3, that HB
(S)ν 6= 0 if and only if ν1 6 −2 and
2
ν2 6 −2, and that HB (S)ν 6= 0 if and only if ν1 6 −2 and ν2 > 0, or ν1 > 0 and
ν2 6 −2 (see [Bot11, Lemma 6.7 and §7]).
To control the local cohomology of Z1 we consider the long exact sequence of
local cohomology associated to the canonical short exact sequence

0 → Z1 → S(−d)4 → I → 0.
1
1
2
2
We get HB
(Z1 ) = 0 and the exact sequences HB
(I) → HB
(Z1 ) → HB
(S(−d)4 )
2
3
3
4
and HB (I) → HB (Z1 ) → HB (S(−d) ). Since we know the local cohomology of S,
it remains to estimate the vanishing of the local cohomology modules H 1 (I) and
H 2 (I). For that purpose, we examine the two spectral sequences associated to the
double complex
•
CB
(K• ((f0 , f1 , f2 , f3 ); S)).
1
2
3
(I)ν = 0 for all ν such that HB
(S(−2d))ν = HB
(S(−3d))ν = 0
We obtain that HB
2
2
3
and that HB (I)ν = 0 for all ν such that HB (S(−d))ν = HB (S(−2d))ν = 0. From
2
3
(Z1 (d))ν = HB
(Z1 (d))ν = 0 for all ν ∈
here, we deduce that HB
/ Rd .
Similarly, to control the local cohomology of Z10 we consider the long exact
sequence of local cohomology associated to the canonical short exact sequence

0 → Z10 (e) → S(−d + e)3 → I 0 → 0.
1
1
2
2
We get HB
(Z10 ) = 0 and the exact sequences HB
(I 0 ) → HB
(Z10 (e)) → HB
(S(−d +
0
0
3
3
2
3
3
e) ) and HB (I ) → HB (Z1 (e)) → HB (S(−d + e) ). As in the previous case, it
remains to estimate the vanishing of the local cohomology modules H 1 (I 0 ) and
H 2 (I 0 ). We proceed similarly by inspecting the two spectral sequences associated
•
2
to the double complex CB
(K• ((g1 , g2 , g3 ); S)) and finally deduce that HB
(Z10 (d))ν =
3
0
/ Rd .

HB (Z1 (d))ν = 0 for all ν ∈

As a consequence of this result, all the properties we obtained in Section 4 also
applied to the setting of a bi-graded parameterization. In particular, we have the
following properties:
/ Rd
• If dim π2−1 (p) 6 0 then for all ν ∈
corank(Mν (φ)(p)) = deg(π2−1 (p)).
• If dim π2−1 (p) = 1 then for all ν ∈
/ Rd
corank(Mν (φ)(p)) = e1 ν1 + e2 ν2 + c
where the curve component of π2−1 (p) is of bi-degree (e1 , e2 ) and c ∈ Z.

FITTING IDEALS AND MULTIPLE-POINTS OF SURFACE PARAMETERIZATIONS

17

From here, an algorithm similar to the one presented in §4.1 can be derived, for all
ν∈
/ Rd , in order to determine the characters of a fiber π2−1 (p) by comparing matrix
representations with consecutive indexes. Typically, the first matrix representation
to consider is M (φ)(d1 −1,2d2 −1) (or M (φ)(2d1 −1,d2 −1) ) which has 2d1 d2 rows. Then,
informations on the fiber at a given point p can be obtained by comparing the rank
at p of M (φ)(d1 −1,2d2 −1) with the ranks at p of M (φ)(d1 ,2d2 −1) and M (φ)(d1 −1,2d2 ) .
Acknowledgments. We would like to thank Ragni Piene and Bernard Teissier
for enlightening discussions about the singularities of finite maps. Most of this
work was done while the first author was hosted at INRIA Sophia Antipolis with
a financial support of the european Marie-Curie Initial Training Network SAGA
(ShApes, Geometry, Algebra), FP7-PEOPLE contract PITN-GA-2008-214584.

References
[BC05]

Laurent Busé and Marc Chardin. Implicitizing rational hypersurfaces using approximation complexes. J. Symbolic Comput, 40(4-5):1150–1168, 2005.
[BC13]
Nicolás Botbol and Marc Chardin. Castelnuovo mumford regularity with respect to
multigraded ideals. Preprint arXiv:1107.2494v2, 2013.
[BDD09] Nicolás Botbol, Alicia Dickenstein, and Marc Dohm. Matrix representations for toric
parametrizations. Comput. Aided Geom. Design, 26(7):757–771, 2009.
[BJ03]
Laurent Busé and Jean-Pierre Jouanolou. On the closed image of a rational map and
the implicitization problem. J. Algebra, 265(1):312–357, 2003.
[BLB12]
Laurent Busé and Thang Luu Ba. The surface/surface intersection problem by means
of matrix based representations. Comput. Aided Geom. Design, 29(8):579–598, 2012.
[Bot11]
Nicolás Botbol. Implicit equation of multigraded hypersurfaces. J. Algebra, 348(1):381–
401, 2011.
[Cha04]
Marc Chardin. Regularity of ideals and their powers. Insitut de Mathématiques de
Jussieu, Prépublication 364, 2004.
[Cha12]
Marc Chardin. Powers of ideals and the cohomology of stalks and fibers of morphisms.
Algebra and Number Theory, to appear (preprint arXiv:1009.1271), 2012.
[CN08]
Gemma Colomé-Nin. Multigraded structures and the depth of blow-up algebras. PhD
Thesis. Universitat de Barcelona, 2008.
[CU02]
Marc Chardin and Bernd Ulrich. Liaison and Castelnuovo-Mumford regularity. Am.
J. Math., 124(6):1103–1124, 2002.
[EH00]
David Eisenbud and Joe Harris. The geometry of schemes. Graduate Texts in Mathematics. 197. New York, NY: Springer. x, 294 p., 2000.
[GS]
Daniel R Grayson and Michael E Stillman. Macaulay 2, a software system for research
in algebraic geometry. http://www.math.uiuc.edu/Macaulay2/.
[KLU92] Steven Kleiman, Joseph Lipman, and Bernd Ulrich. The source double-point cycle
of a finite map of codimension one. Complex projective geometry, Sel. Pap. Conf.
Proj. Var., Trieste/Italy 1989, and Vector Bundles and Special Proj. Embeddings,
Bergen/Norway 1989, Lond. Math. Soc. Lect. Note Ser. 179, 199-212 (1992)., 1992.
[KLU96] Steven Kleiman, Joseph Lipman, and Bernd Ulrich. The multiple-point schemes of a
finite curvilinear map of codimension one. Ark. Mat., 34(2):285–326, 1996.
[LBBM09] Thang Luu Ba, Laurent Busé, and Bernard Mourrain. Curve/surface intersection problem by means of matrix representations. In H. Kai and H. Sekigawa, editors, SNC,
pages 71–78, Kyoto, Japan, 2009. ACM Press.
[MP89]
David Mond and Ruud Pellikaan. Fitting ideals and multiple points of analytic mappings. In Algebraic geometry and complex analysis (Pátzcuaro, 1987), volume 1414 of
Lecture Notes in Math., pages 107–161. Springer, Berlin, 1989.
[Pie78]
Ragni Piene. Polar classes of singular varieties. Ann. Sci. École Norm. Sup. (4),
11(2):247–276, 1978.
[SC95]
Tom Sederberg and Falai Chen. Implicitization using moving curves and surfaces.
303:301–308, 1995.

18

NICOLÁS BOTBOL, LAURENT BUSÉ, AND MARC CHARDIN

[SSV12]

[Sza08]
[Tei77]

Hal Schenck, Alexandra Seceleanu, and Javid Validashti. Syzygies and singularities
of tensor product surfaces of bidegree (2,1). Preprint arXiv:1211.1648, to appear in
Mathematics of Computations, 2012.
Agnes Szanto. Solving over-determined systems by the subresultant method. J. Symbolic Comput., 43(1):46–74, 2008. With an appendix by Marc Chardin.
Bernard Teissier. The hunting of invariants in the geometry of discriminants. In Real
and complex singularities (Proc. Ninth Nordic Summer School/NAVF Sympos. Math.,
Oslo, 1976), pages 565–678. Sijthoff and Noordhoff, Alphen aan den Rijn, 1977.

Departamento de Matemática, FCEN, Universidad de Buenos Aires, Argentina.
E-mail address: nbotbol@dm.uba.ar
INRIA Sophia Antipolis - Méditerranée, Galaad team, 2004 route des Lucioles, B.P. 93,
F-06902, Sophia Antipolis France.
E-mail address: Laurent.Buse@inria.fr
Institut Mathématique de Jussieu et Université Pierre et Marie Curie, Boı̂te 247, 4
place Jussieu, F-75252 Paris CEDEX 05, France.
E-mail address: chardin@math.jussieu.fr

