Testing zero-dimensionality of varieties at a point

TL;DR

This paper introduces effective algorithms using comprehensive Gröbner systems to test zero-dimensionality of algebraic varieties at a point, analyze deformations of singularities, and compute local dimensions, with applications to parameter space stratification.

Contribution

It presents novel algorithms leveraging comprehensive Gröbner systems for zero-dimensionality testing and local dimension computation at a point, including parameter space decomposition.

Findings

Algorithms successfully determine zero-dimensionality at a point.

Methods enable analysis of deformations of hypersurface singularities.

Parameter space can be decomposed into strata based on local dimension.

Abstract

Effective methods are introduced for testing zero-dimensionality of varieties at a point. The motivation of this paper is to compute and analyze deformations of isolated hypersurface singularities. As an application, methods for computing local dimensions are also described. For the case where a given ideal contains parameters, the proposed algorithms can output in particular a decomposition of a parameter space into strata according to the local dimension at a point of the associated varieties. The key of the proposed algorithms is the use of the notion of comprehensive Gr\"obner systems.

Tables

Table 1. Table 1. Comparison of the three algorithms

Problem	$F$	Algorithm 1	Algorithm 2-1	Algorithm 2-2
1	$\{f_1,\frac{\partial f_1}{\partial x},\frac{\partial f_1}{\partial y},\frac{\partial f_1}{\partial z}\}$	0.0156	0.0781	0.0156
2	$\{f_2,\frac{\partial f_2}{\partial x},\frac{\partial f_2}{\partial y}\}$	1.375	0.0468	0.0156
3	$\{f_3,\frac{\partial f_3}{\partial x},\frac{\partial f_3}{\partial y}\}$	27.94	0.7031	0.0156
4	$\{\frac{\partial f_4}{\partial x},\frac{\partial f_4}{\partial y}\}$	$>3h$	$0.0156$	0.0156
5	$\{f_5,\frac{\partial f_5}{\partial x},\frac{\partial f_5}{\partial y}\}$	2.719	$>3h$	0.0156
6	$\{\frac{\partial f_6}{\partial x},\frac{\partial f_6}{\partial y},\frac{\partial f_6}{\partial z}\}$	0.063	9.969	0.2188
7	$\{\frac{\partial f_7}{\partial x},\frac{\partial f_7}{\partial y},\frac{\partial f_7}{\partial z}\}$	0.0156	$>3h$	226.8
8	$\{f_7,\frac{\partial f_7}{\partial x},\frac{\partial f_7}{\partial y},\frac{\partial f_7}{\partial z}\}$	0.0156	6.922	0.3125
9	$\{\frac{\partial f_8}{\partial x},\frac{\partial f_8}{\partial y}\}$	$>3h$	0.375	0.2669
10	$\{f_9,\frac{\partial f_9}{\partial x},\frac{\partial f_9}{\partial y},\frac{\partial f_9}{\partial z}\}$	$>3h$	$>3h$	1.391
11	$\{f_{10},\frac{\partial f_{10}}{\partial x},\frac{\partial f_{10}}{\partial y},\frac{\partial f_{10}}{\partial z}\}$	$>3h$	$>3h$	6.188

Equations

$$f$$

$$C_{p}(V)=\operatorname{\mathbb{V}}\Bigl{(}f_{p,min}\mid f\in\operatorname{\mathbb{I}}(V)\Bigr{)}$$

$${\dim }_p(V)=\dim (C_p(V)).$$

$$f^h(x_0,x)=\sum _{i=0}^d{f}_i(x)x_0^{d-i}$$

$$I^h=⟨f^h|f\in I⟩\subset C[x_0,x].$$

$$C_{O}(V)=\operatorname{\mathbb{V}}\Bigl{(}\varepsilon(g_{1})_{O,min},\varepsilon(g_{2})_{O,min},\ldots,\varepsilon(g_{r})_{O,min}\Bigr{)}$$

$$I^h=⟨x_1^3{x}_0^9+t_1{x}_1^2{x}_2^4{x}_0^6+x_2^{12},3x_1^2{x}_0^2+2t_1{x}_2^4,4t_1{x}_1^2{x}_2^3{x}_0^6+12x_2^{11}⟩$$

$$\{(C\backslash V(t_1(4t_1^3+27)),\{(4t_1^3+27)x_1{x}_2^{11},(4t_1^3+27)x_2^{15},3x_0^3{x}_1^2+2t_1{x}_1{x}_2^4,$$

$$2t_1^2{x}_0^3{x}_1{x}_2^7-9x_2^{11}\}),$$

$$(V(4t_1^3+27),\{3x_0^3{x}_1^2+2t_1{x}_1{x}_2^4,3x_0^3{x}_1{x}_2^7+2t_1{x}_2^{11}\}),(V(t_1),\{x_2^{11},x_0{x}_1^2\})\}.$$

$$I:J^{\infty }=\{h\in C[x]|h{J}^r\subset I\text{ for some }r>0\}.$$

$$I:J\subseteq I:J^2\subseteq I:J^3\subseteq \cdots$$

$$V(G)=V(⟨F⟩:{\mathfrak{m}}^{\infty })=\overline{V(F)\backslash \{O\}}=V_2\cup \cdots \cup V_{\nu }.$$

$$V(G)=V(⟨F⟩:{\mathfrak{m}}^{\infty })=\overline{V(F)\backslash \{O\}}\text{ and }O\notin V(G),$$

$$\{1\},$$

$$\{x_2^2+x_3,x_1\}.$$

$$\{(C\backslash V(t_1^2+1),\{1\}),(V(t_1^2+1),\{x_1-t_1{x}_3,x_2^3+2t_1{x}_3^2\})\}.$$

$$\Bigl{(}\langle F\rangle:\langle x_{1}^{\alpha_{1}},x_{2}^{\alpha_{2}},\ldots,x_{n}^{\alpha_{n}}\rangle\Bigr{)}:{\mathfrak{m}}^{\infty}=\langle F\rangle:{\mathfrak{m}}^{\infty}.$$

$$⟨F⟩:{\mathfrak{m}}^{\beta -1}\subset ⟨F⟩:{\mathfrak{m}}^{\beta }=⟨F⟩:{\mathfrak{m}}^{\beta +1}=\cdots .$$

$$\Bigl{(}\langle F\rangle:J\Bigr{)}:{\mathfrak{m}}^{\beta}=\langle F\rangle:J\cdot{\mathfrak{m}}^{\beta}.$$

$$⟨F⟩:{\mathfrak{m}}^{\beta }\subseteq ⟨F⟩:J\cdot {\mathfrak{m}}^{\beta }\subseteq ⟨F⟩:{\mathfrak{m}}^N.$$

$$\text{mdeg}_{x_{i}}(f):=\max\Bigl{(}\{\gamma_{i}|\ {\bf r}=(\gamma_{1},\ldots,\gamma_{i},\ldots,\gamma_{n})\in\operatorname{\mathbb{N}}^{n},a_{{\bf r}}\neq 0\}\Bigr{)},$$

$$\alpha=\max\Bigl{(}\{\text{mdeg}_{x_{1}}(F),\text{mdeg}_{x_{2}}(F),\ldots,\text{mdeg}_{x_{n}}(F)\}\Bigr{)}$$

$$W=\Bigl{(}\operatorname{\mathbb{C}}^{2}\backslash\operatorname{\mathbb{V}}((t_{1}^{4}+t_{1}^{2}t_{2}+1)(t_{2}-2)(t_{2}+2)\Bigr{)}\cup\Bigl{(}\operatorname{\mathbb{V}}(t_{2})\backslash\operatorname{\mathbb{V}}(t_{1}^{4}+1,t_{2})\Bigr{)},$$

$$⟨F⟩=Q_0\cap Q_1\cap \dots \cap Q_{\nu }$$

$$⟨F⟩=Q_0\cap S\text{ and }{Q}_0=⟨F⟩:S.$$

$$\{(V(t_1),\{1\}),(C\backslash V(t_1),\{2910897t_1^3{x}_1+16385050t_1{x}_2+14895500,$$

$$-9801t_1^2{x}_1+52855t_1{x}_2^2+48050x_2\})\}.$$

$$\{(V(t_1),Q_0^{\prime })\},$$

$$\displaystyle\begin{array}[]{c}x_{1}+u_{11}x_{\ell+1}+u_{12}x_{\ell+2}+\cdots+u_{1n-\ell}x_{n},\\ x_{2}+u_{21}x_{\ell+1}+u_{22}x_{\ell+2}+\cdots+u_{2n-\ell}x_{n},\\ x_{3}+u_{31}x_{\ell+1}+u_{32}x_{\ell+2}+\cdots+u_{3n-\ell}x_{n},\\ \vdots\\ x_{\ell}+u_{\ell 1}x_{\ell+1}+u_{\ell 2}x_{\ell+2}+\cdots+u_{\ell n-\ell}x_{n},\\ \\ \end{array}$$

$$\{x_1^2+x_2^3\}.$$

$$\{x_2+u_{11}^2,x_1-u_{11}^3\}.$$

Full text

1. Introduction

The local dimension, the dimension of a variety at a point, is one of the most important invariants in algebraic geometry, complex analysis and singularity theory [2, 7, 15, 20]. Thus, a practical tool to compute the dimension or test zero-dimensionality is required for studying local properties of a variety [3, 11, 12, 14].

In this paper, we propose two methods for testing zero-dimensionality of a variety at a point, and we generalize them to the parametric cases. The main tools of our approach are Gröbner bases and comprehensive Gröbner systems. The proposed methods do not utilize primary ideal decompositions, and are free from computation in local rings.

Definition 1.1**.**

Let $V$ be an affine variety in $\operatorname{\mathbb{C}}^{n}$. For $p\in V$, the dimension of $V$ at $p$, denoted $\dim_{p}(V)$, is the maximum dimension of an irreducible component of $V$ containing the point $p$.

In singularity theory, problems that contain parameters are often studied, for instance, deformations of singularities, a family of hyperplane sections of a variety, etc. In such cases, since structures of relevant ideals or varieties may vary as parameters changes, there is a possibility that the local dimension of varieties may also depend on the values of parameters. We need methods to decompose a parameter space into strata according to the local dimensions of a given family of varieties.

In order to state precisely the problem, we give an example. Let $f_{0}=x_{1}^{4}+x_{1}x_{3}^{2}+x_{2}^{4}$ and consider $f=f_{0}+t_{1}x_{2}x_{3}^{2}$, where $t_{1}$ is a parameter. The hypersurface defined by $f_{0}=0$ has an isolated singularity at the origin $O$ in $\operatorname{\mathbb{C}}^{3}$, i.e., $\dim_{O}(\operatorname{\mathbb{V}}(f_{0},\frac{\partial f_{0}}{\partial x_{1}},\frac{\partial f_{0}}{\partial x_{2}},\frac{\partial f_{0}}{\partial x_{3}}))=0$. Since $f$ has the parameter $t_{1}$, there is a possibility that the family of hypersurfaces defined by $f=0$ have non-isolated singularities at the origin for some values of the parameter $t_{1}.$ In fact, if $t_{1}^{4}+1=0$, then $f$ has a non-isolated singularity at $O$ and $\dim_{O}(\operatorname{\mathbb{V}}(f,\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}}))$ $=1$. If $t_{1}^{4}+1\neq 0$, then $f$ has an isolated singularity at $O$. We really would like the condition $t_{1}^{4}+1\neq 0$, or detect the condition $t_{1}^{4}+1=0$ in an algorithmic manner to study local properties of the deformation of an isolated singularity. How do we obtain such conditions?

Basically, the condition can be obtained by testing zero-dimensionality of the variety $\operatorname{\mathbb{V}}(f,$ $\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}})$ at the origin $O$. We show in the present paper that the methods for testing zero-dimensionality of a variety at a point can be constructed by using Gröbner bases. Furthermore, we generalize the methods to parametric cases by utilizing comprehensive Gröbner systems [5, 9, 10, 13, 18]. We give two different kinds of algorithms for testing zero-dimensionality at a point of a family of varieties with parameters.

Note that the resulting algorithms do not involve computation in local rings and efficiently output the necessary and sufficient conditions for zero-dimensionality.

This paper is organized as follows. Section 2 briefly reviews comprehensive Gröbner systems, and give notations that will be used in this paper. Section 3 consider the use of tangent cone and gives the discussion of the first algorithm for testing zero-dimensionality of varieties at a point. Section 4 consider the use of saturation and discusses the second algorithm for testing zero-dimensionality of varieties at a point. Section 5 gives results of the benchmark tests. Appendix A gives an efficient algorithm for computing ideal quotients with parameters, that utilizes a comprehensive Gröbner system of a module.

2. Preliminaries

Let $t=\{t_{1},\ldots,t_{m}\}$ and $x=\{x_{1},\ldots,x_{n}\}$ be variables such that $t\cap x=\emptyset$ and $\operatorname{\mathbb{C}}[t][x]$ be a polynomial ring with coefficients in a polynomial ring $\operatorname{\mathbb{C}}[t]$. For $f_{1},\ldots,f_{s}\in\operatorname{\mathbb{C}}[x]$ (or $\operatorname{\mathbb{C}}[t][x]$), $\langle f_{1},\ldots,f_{s}\rangle=\{\sum_{i=1}^{s}h_{i}f_{i}|h_{1},\ldots,h_{s}\in\operatorname{\mathbb{C}}[x](\text{ or }\operatorname{\mathbb{C}}[t][x])\}$.

A symbol $\operatorname{Term}(x)$ is the set of terms of $x$. Fix a term order $\succ$ on $\operatorname{Term}(x)$. Let $f\in\operatorname{\mathbb{C}}[x]$ (or $f\in\operatorname{\mathbb{C}}[t][x]$), then, $\operatorname{ht}(f),\operatorname{hm}(f)$ and $\operatorname{hc}(f)$ denote the head term, head monomial and head coefficient of $f$ (i.e., $\operatorname{hm}(f)=\operatorname{ht}(f)\cdot\operatorname{ht}(f)$). For $F\subset\operatorname{\mathbb{C}}[x]$ (or $F\subset\operatorname{\mathbb{C}}[t][x]$), $\operatorname{ht}(F)=\{\operatorname{ht}(f)|f\in F\}$.

For $g_{1},\ldots,g_{r}\in\operatorname{\mathbb{C}}[t]$, $\operatorname{\mathbb{V}}(g_{1},\ldots,g_{r})\subseteq\operatorname{\mathbb{C}}^{m}$ denotes the affine variety of $g_{1},\ldots,g_{r}$, i.e., $\operatorname{\mathbb{V}}(g_{1},\ldots,g_{r})=\{\bar{t}\in\operatorname{\mathbb{C}}^{m}|g_{1}(\bar{t})=\cdots=g_{r}(\bar{t})=0\}$. We call an algebraic constructible set of a from $\operatorname{\mathbb{V}}(g_{1},\ldots,g_{r})\backslash\operatorname{\mathbb{V}}(g^{\prime}_{1},\ldots,g^{\prime}_{r^{\prime}})\subseteq\operatorname{\mathbb{C}}^{m}$ with $g_{1},\ldots,g_{r},g^{\prime}_{1},\ldots,g^{\prime}_{r^{\prime}}\in\operatorname{\mathbb{C}}[t]$, a stratum. Notations $\operatorname{\mathbb{A}}_{1},\operatorname{\mathbb{A}}_{2},$ $\ldots,\operatorname{\mathbb{A}}_{\nu}$ are frequently used to represent strata.

For every $\bar{t}\in\operatorname{\mathbb{C}}^{m}$, the canonical specialization homomorphism $\sigma_{\bar{t}}:\operatorname{\mathbb{C}}[t][x]\rightarrow\operatorname{\mathbb{C}}[x]$ (or $\operatorname{\mathbb{C}}[t]\rightarrow\operatorname{\mathbb{C}}$) is defined as the map that substitutes $t$ by $\bar{t}$ in $f(t,x)\in\operatorname{\mathbb{C}}[t][x]$ (i.e., $\sigma_{\bar{t}}(f)=f(\bar{t},x)\in\operatorname{\mathbb{C}}[x]$). The image $\sigma_{\bar{t}}$ of a set $F$ is denoted by $\sigma_{\bar{t}}(F)=\{\sigma_{\bar{t}}(f)|f\in F\}\subset\operatorname{\mathbb{C}}[x]$. In this paper, the set of natural numbers $\operatorname{\mathbb{N}}$ includes zero.

We adopt the following as a definition of a comprehensive Gröbner system.

Definition 2.1** (comprehensive Gröbner system).**

Let $\succ$ be a term order on $\operatorname{Term}(x)$. Let $F$ be a subset of $\operatorname{\mathbb{C}}[t][x]$, $\operatorname{\mathbb{A}}_{1},\operatorname{\mathbb{A}}_{2},\ldots,\operatorname{\mathbb{A}}_{\nu}$ strata in $\operatorname{\mathbb{C}}^{m}$ and $G_{1},G_{2},\ldots,G_{\nu}$ subsets in $\operatorname{\mathbb{C}}[t][x]$. If a finite set ${\mathcal{G}}=\{(\operatorname{\mathbb{A}}_{1},G_{1}),(\operatorname{\mathbb{A}}_{2},G_{2}),\ldots,(\operatorname{\mathbb{A}}_{\nu},$ $G_{\nu})\}$ of pairs satisfies the following conditions

1. (1)

$\operatorname{\mathbb{A}}_{i}\neq\emptyset$ and $\operatorname{\mathbb{A}}_{i}\cap\operatorname{\mathbb{A}}_{j}=\emptyset$ for $1\leq i\neq j\leq\nu$, 2. (2)

for all $\bar{t}\in\operatorname{\mathbb{A}}_{i}$, $\sigma_{\bar{t}}(G_{i})$ is a minimal Gröbner basis of $\langle\sigma_{\bar{t}}(F)\rangle$ w.r.t. $\succ$ in $\operatorname{\mathbb{C}}[x]$, and 3. (3)

for all $\bar{t}\in\operatorname{\mathbb{A}}_{i}$ and $f\in G_{i}$, $\sigma_{\bar{t}}(\operatorname{hc}(f))\neq 0$,

the set ${\mathcal{G}}$ is called a comprehensive Gröbner system on $\operatorname{\mathbb{A}}_{1}\cup\cdots\cup\operatorname{\mathbb{A}}_{\ell}$ for $\langle F\rangle$ w.r.t. $\succ$. We simply say that ${\mathcal{G}}$ is a comprehensive Gröbner system for $\langle F\rangle$ if $\operatorname{\mathbb{A}}_{1}\cup\cdots\cup\operatorname{\mathbb{A}}_{\ell}=\operatorname{\mathbb{C}}^{m}$.

In several papers [5, 9, 10, 13], algorithms and implementations for computing comprehensive Gröbner systems are introduced.

Example 1*.*

Let $F=\{t_{1}x_{1}x_{2}+x_{2}+1,x_{1}^{2}x_{2}+t_{1}x_{1}+3\}$ be a subset in $\operatorname{\mathbb{C}}[t_{1}][x_{1},x_{2}]$ and $\succ$ the lexicographic term order s.t. $x_{1}\succ x_{2}$. We regard $t_{1}$ as a parameter in $\operatorname{\mathbb{C}}$. Then, a comprehensive Gröbner system of $\langle F\rangle$ w.r.t. $\succ$ is

\Bigl{\{}(\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}^{3}-t_{1}),\{x_{2}^{2}+(2t_{1}^{2}+2)x_{2}-t_{1}^{2}+1,(t_{1}^{3}-t_{1})x_{1}+x_{2}+3t_{1}^{2}+1\}),(\operatorname{\mathbb{V}}(t_{1}^{2}-1),\{4x_{1}+3t_{1},x_{2}+4\}),(\operatorname{\mathbb{V}}(t_{1}),\{x_{1}^{2}-3,x_{2}+1\})\Bigl{\}}.

3. Algorithm 1 (Tangent cone approach)

Here, we present an algorithm for testing zero-dimensionality of a variety at a point. This algorithm is based on the method described in section 9 of the famous textbook [1]. We generalize the method to parametric cases.

Before introducing the algorithm, we prepare some notations and basic facts.

Let $p=(p_{1},\ldots,p_{n})\in\operatorname{\mathbb{C}}^{n}$, $\alpha=(\alpha_{1},\ldots,\alpha_{n})\in\operatorname{\mathbb{N}}^{n}$ and $(x-p)^{\alpha}=(x_{1}-p_{1})^{\alpha_{1}}\cdots(x_{n}-p_{n})^{\alpha_{n}}$. Given any polynomial $f\in\operatorname{\mathbb{C}}[x]$ of total degree $d$, $f$ can be written as a polynomial in $x_{i}-p_{i}$, namely,

[TABLE]

where $f_{p,j}$ is a linear combination of $(x-p)^{\alpha}$ for $\alpha_{1}+\cdots+\alpha_{n}=j\leq d$. Note that $f_{p,0}=f(p)$ and $f_{p,1}=\frac{\partial f}{\partial x_{1}}(p)(x_{1}-p_{1})+\cdots+\frac{\partial f}{\partial x_{n}}(p)(x_{n}-p_{n})$.

The next definition is borrowed from [1].

Definition 3.1** (tangent cone).**

Let $V\subset\operatorname{\mathbb{C}}^{n}$ be an affine variety and let $p=(p_{1},\ldots,p_{n})\in V$.

1. (i)

If $f\in\operatorname{\mathbb{C}}[x]$ is a non-zero polynomial, then $f_{p,min}$ is defined to be $f_{p,j}$, where $j$ the smallest integer such that $f_{p,j}\neq 0$ in (3.1). 2. (ii)

The tangent cone of $V$ at $p$, denoted $C_{p}(V)$, is the variety

[TABLE]

where $\operatorname{\mathbb{I}}(V)=\{f\in\operatorname{\mathbb{C}}[x]|f(\bar{x})=0,\text{ for all }\bar{x}\in V\}$.

The details of the tangent cone are described in [19, 20]. In 1965, H. Whitney gave the following theorem.

Theorem 3.2** (H. Whitney [18]).**

Let $V\subset\operatorname{\mathbb{C}}^{n}$ be an affine variety and let $p=(p_{1},\ldots,p_{n})\in V$. Then,

[TABLE]

In order to compute a tangent cone, we need the following definition.

Definition 3.3**.**

1. (i)

Let $f(x)\in\operatorname{\mathbb{C}}[x]$ be a polynomial of total degree $d$. Let $f(x)=\sum_{i=0}^{d}f_{i}(x)$ be the expansion of $f(x)$ as the sum of its homogeneous components where $f_{i}(x)$ has total degree $i$. Then,

[TABLE]

is a homogeneous polynomial of total degree $d$ in $\operatorname{\mathbb{C}}[x_{0},x]$ where $x_{0}$ is a new variable. 2. (ii )

Let $I$ be an ideal in $\operatorname{\mathbb{C}}[x]$. We define the homogenization of $I$ to be the ideal

[TABLE]

From now on, we assume that the point $p$ is the origin $O=(0,\ldots,0)$ in $\operatorname{\mathbb{C}}^{n}$.

Proposition 3.4** (Proposition 4, p.485 [1]).**

Assume that the origin $O$ is a point of $V\subset\operatorname{\mathbb{C}}^{n}$. Let $\succ$ be a block term order such that $x_{0}\gg x$. Let $I$ be an ideal such that $V=\operatorname{\mathbb{V}}(I)$. If $\{g_{1},\ldots,g_{r}\}$ is a Gröbner basis of $I^{h}$ w.r.t. $\succ$, then

[TABLE]

where $\varepsilon(g_{i})$ is the dehomogenization of $g_{i}$ for $1\leq i\leq r$.

There exist several algorithms for computing the dimension of a variety, thus, the dimension of $C_{O}(V)$ can be obtained. The procedure for computing $\dim_{O}(V)$ is the following.

1. Step 1:

Compute $C_{O}(V)$. 2. Step 2:

Compute $\dim(C_{O}(V))$. 3. Return

$\dim(C_{O}(V))$ (as $\dim_{O}(V)=\dim(C_{O}(V))$).

We turn to the parametric cases. Let $I$ be an ideal in $\operatorname{\mathbb{C}}[t][x]$ where we regard $t$ as parameters. After here we simply say that $I$ is a “parametric” ideal.

As described in section 2, there exist algorithms for computing comprehensive Gröbner systems, it is possible to compute a comprehensive Gröbner system of $I^{h}$ w.r.t. $\succ$ in Proposition 3.4. Therefore, Proposition 3.4 and the procedure above can be extended to the case of parametric ideals.

The following algorithm which utilizes a comprehensive Gröbner system outputs a condition of zero-dimensionality of $\operatorname{\mathbb{V}}(F)$ at $O$.

**Algorithm 1

 ** Input: $F=\{f_{1},f_{2},\ldots,f_{s}\}\subset\operatorname{\mathbb{C}}[t][x]$ s.t. $O\in\operatorname{\mathbb{V}}(F)$.

           $\succ$: a block term order s.t. $x_{0}\gg x$ on $\operatorname{Term}(\{x_{0}\}\cup x)$.

Output: $\operatorname{\mathbb{A}}\subset\operatorname{\mathbb{C}}^{m}$: For all $\bar{t}\in\operatorname{\mathbb{A}}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))=0$ (i.e., $\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F))$ has an isolated point at $O$). For all $\bar{t}\in\operatorname{\mathbb{C}}^{m}\backslash\operatorname{\mathbb{A}}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))\neq 0$.

BEGIN

$\operatorname{\mathbb{A}}\leftarrow\emptyset$;

${\mathcal{G}}\leftarrow$ Compute a comprehensive Gröbner system of $\langle f_{1}^{h},f_{2}^{h},\ldots,f_{s}^{h}\rangle$ w.r.t. $\succ$ in $\operatorname{\mathbb{C}}[t][x_{0},x]$;

while ${\mathcal{G}}\neq\emptyset$ **do

**    Select $(\operatorname{\mathbb{A}}^{\prime},G^{\prime})$ from ${\mathcal{G}}$; ${\mathcal{G}}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

    $M\leftarrow\operatorname{ht}(G^{\prime})$ w.r.t. $\succ$;

    $C_{O}\leftarrow\{\varepsilon(h)|h\in M\}$ in $\operatorname{\mathbb{C}}[x]$;

      if $\dim(C_{O})=0$ **then

**          $\operatorname{\mathbb{A}}\leftarrow\operatorname{\mathbb{A}}^{\prime}\cup\operatorname{\mathbb{A}}$;

      **end-if

end-while

return** $\operatorname{\mathbb{A}}$;

END

Since the algorithms [5, 9, 10, 13, 18] for computing comprehensive Gröbner systems always terminate and return a finite set of pairs, Algorithm 1 also terminates. The correctness follows from Theorem 3.2 and Proposition 3.4.

Note that Algorithm 1 contains a part of computing local dimensions. Thus, it can be naturally generalized to a method for decomposing a parameter space into strata according to the local dimensions of a given family of varieties.

We illustrate Algorithm 1 with the following example.

Example 2*.*

Let $f=x_{1}^{3}+t_{1}x_{1}^{2}x_{2}^{4}+x_{2}^{12}\in\operatorname{\mathbb{C}}[t_{1}][x_{1},x_{2}]$, $F=\{f,\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}}\}$ and $I=\langle F\rangle$ where $t_{1}$ is a parameter. Let $\succ$ be the total degree lexicographic term order with $x_{1}\succ x_{2}$.

A comprehensive Gröbner systems of

[TABLE]

w.r.t. the block term order with $x_{0}\gg\{x_{1},x_{2}\}$, in $\operatorname{\mathbb{C}}[t_{1}][x_{0},x_{1},x_{2}]$, is

[TABLE]

Hence,

1. $\cdot$

if $t_{1}$ belongs to $\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}(4t_{1}^{3}+27)$, then $C_{O}(\operatorname{\mathbb{V}}(I))=\operatorname{\mathbb{V}}(x_{1}x_{2}^{11},x_{2}^{15},x_{1}^{2},x_{1}^{2}x_{2}^{7})$ and $\dim_{O}(\operatorname{\mathbb{V}}(I))$ $=0$, 2. $\cdot$

if $t_{1}$ belongs to $\operatorname{\mathbb{V}}(4t_{1}^{3}+27)$, then $C_{O}(\operatorname{\mathbb{V}}(I))=\operatorname{\mathbb{V}}(x_{1}^{2},x_{1}x_{2}^{7})$ and $\dim_{O}(\operatorname{\mathbb{V}}(I))=1$, 3. $\cdot$

if $t_{1}$ belongs to $\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1})$, then $C_{O}(\operatorname{\mathbb{V}}(I))=\operatorname{\mathbb{V}}(x_{1}^{2},x_{2}^{11})$ and $\dim_{O}(\operatorname{\mathbb{V}}(I))=0$.

Therefore, for all $\bar{t}\in\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(4t_{1}^{3}+27)$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(I)))=0$, namely, $f$ has an isolated singularity at the origin $O$.

4. Algorithm 2 (Saturation approach)

in this section, we consider the use of saturation and introduce an alternative method for testing zero-dimensionality of a variety at a point. We present an algorithm for testing zero-dimensionality of a family of varieties at a point, that utilizes a comprehensive Gröbner system of a module. Furthermore, we improve the algorithm in speed and give an efficient algorithm. We also show that, according to the concept of Chevalley dimension, local dimensions of varieties can also be computed by utilizing saturation.

4.1. Saturation approach

Let $I,J$ be ideals in $\operatorname{\mathbb{C}}[x]$. The ideal quotient of $I$ by $J$ is $I:J=\{h\in\operatorname{\mathbb{C}}[x]|hg\in I\text{ for all }g\in J\}.$ The saturation of $I$ with respect to $J$ is the ideal

[TABLE]

The saturation $I:J^{\infty}$ is the ideal at which the chain

[TABLE]

stabilizes.

Now, we give the following main theorem which is utilized to construct the new algorithm for testing the zero-dimensionality of a variety at a point.

Theorem 4.1**.**

Let $F\subset\operatorname{\mathbb{C}}[x]$, ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle\subset\operatorname{\mathbb{C}}[x]$ and $O\in\operatorname{\mathbb{V}}(F)$. Let $G$ be a basis of the ideal $\langle F\rangle:{\mathfrak{m}}^{\infty}$ in $\operatorname{\mathbb{C}}[x]$. Then, the affine variety $\operatorname{\mathbb{V}}(F)$ has an isolated point at the origin $O$ if and only if there exists $g\in G$ such that $g(O)\neq 0$ (i.e., $g$ has non-zero constant term).

Proof.

As $G$ is a basis of $\langle F\rangle:{\mathfrak{m}}^{\infty}$, $\operatorname{\mathbb{V}}(G)=\operatorname{\mathbb{V}}(\langle F\rangle:{\mathfrak{m}}^{\infty})$ is the Zariski closure $\overline{\operatorname{\mathbb{V}}(F)\backslash\{O\}}$. The variety $\operatorname{\mathbb{V}}(F)$ can be written as $\operatorname{\mathbb{V}}(F)=V_{1}\cup V_{2}\cup\cdots\cup V_{\nu}$ (finite union) where $V_{1},V_{2},\ldots,V_{\nu}$ are distinct irreducible varieties.

First, assume that the affine variety $\operatorname{\mathbb{V}}(F)$ has an isolated point at the origin $O$. Then, one of $V_{i}$ must be $\{O\}$ and other varieties does not contain $O$ where $i\in\{1,2,\ldots,\nu\}$. Without loss of generality, set $V_{1}=\{O\}$. Then,

[TABLE]

As $O\notin V_{2}\cup\cdots\cup V_{\nu}=\operatorname{\mathbb{V}}(G)$, there exists $g\in G$ such that $g(O)\neq 0$.

Next, assume that there exists $g\in G$ such that $g(O)\neq 0$. Since $O\in\operatorname{\mathbb{V}}(F)$,

[TABLE]

there exists an irreducible variety $\{O\}$ in $\operatorname{\mathbb{V}}(F)$. Therefore, $\operatorname{\mathbb{V}}(F)$ has an isolated point at the origin $O$. ∎

The following corollary is a direct consequence of Theorem 4.1.

Corollary 4.2**.**

Let $f\in\operatorname{\mathbb{C}}[x]$, ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle\subset\operatorname{\mathbb{C}}[x]$ and $O\in\operatorname{\mathbb{V}}(f,\frac{\partial f}{\partial x_{1}},\ldots,\frac{\partial f}{\partial x_{n}})$. Let $G$ be a basis of $\langle\frac{\partial f}{\partial x_{1}},\ldots,\frac{\partial f}{\partial x_{n}}\rangle:{\mathfrak{m}}^{\infty}$ (or $\langle f,\frac{\partial f}{\partial x_{1}},\ldots,\frac{\partial f}{\partial x_{n}}\rangle:{\mathfrak{m}}^{\infty}$). Then, the hypersurface, defined by $f$, has an isolated singularity at $O$ if and only if there exists $g\in G$ such that $g(O)\neq 0$.

Example 3*.*

Let us consider $f_{1}=x_{1}^{2}x_{3}+x_{2}x_{3}^{2}+x_{2}^{5}+x_{2}^{3}x_{3},f_{2}=x_{1}^{2}x_{3}+x_{2}x_{3}^{2}+x_{2}^{5}+2x_{2}^{3}x_{3}$ and ${\mathfrak{m}}=\langle x_{1},x_{2},x_{3}\rangle$ in $\operatorname{\mathbb{C}}[x_{1},x_{2},x_{3}]$. Let $\succ$ be the total degree lexicographic term order with $x_{1}\succ x_{2}\succ x_{3}$. Then, the reduced Gröbner basis of $\langle\frac{\partial f_{1}}{\partial x_{1}},\frac{\partial f_{1}}{\partial x_{2}},\frac{\partial f_{1}}{\partial x_{3}}\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$ is

[TABLE]

and the reduced Gröbner basis of $\langle\frac{\partial f_{2}}{\partial x_{1}},\frac{\partial f_{2}}{\partial x_{2}},\frac{\partial f_{2}}{\partial x_{3}}\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$ is

[TABLE]

Therefore, $f_{1}$ has an isolated singularity at $O$, and $f_{2}$ does not have an isolated singularity at $O$.

We turn to the parametric cases. There exists an algorithm for computing a comprehensive Gröbner system of the saturation of $\langle F\rangle$ w.r.t. a given parametric ideal. The algorithm is given in Appendix A. Therefore, Theorem 4.1 is generalized to the parametric cases.

**Algorithm 2-1

 ** Input: $F=\{f_{1},f_{2},\ldots,f_{s}\}\subset\operatorname{\mathbb{C}}[t][x]$ s.t. $O\in\operatorname{\mathbb{V}}(F)$, ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle$.

           $\succ$: a term order on $\operatorname{Term}(x)$.

Output: $\operatorname{\mathbb{A}}\subset\operatorname{\mathbb{C}}^{m}$: For all $\bar{t}\in\operatorname{\mathbb{A}}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))=0$ (i.e., $\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F))$ has an isolated point at $O$). For all $\bar{t}\in\operatorname{\mathbb{C}}^{m}\backslash\operatorname{\mathbb{A}}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))\neq 0$.

BEGIN

$\operatorname{\mathbb{A}}\leftarrow\emptyset$;

${\mathcal{G}}\leftarrow$ Compute a comprehensive Gröbner system of $\langle F\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$;

while ${\mathcal{G}}\neq\emptyset$ **do

**    Select $(\operatorname{\mathbb{A}}^{\prime},G^{\prime})$ from ${\mathcal{G}}$; ${\mathcal{G}}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

      if $\exists g\in G^{\prime}$ s.t. $g(O)\neq 0$ **then

**          $\operatorname{\mathbb{A}}\leftarrow\operatorname{\mathbb{A}}^{\prime}\cup\operatorname{\mathbb{A}}$;

      **end-if

end-while

return** $\operatorname{\mathbb{A}}$;

END

The correctness and termination of Algorithm 2-1 follows from Theorem 4.1 and that of the algorithm for computing comprehensive Gröbner systems.

We illustrate Algorithm 1 with the following example.

Example 4*.*

Let $f=x_{1}^{3}+x_{1}x_{3}^{2}+t_{1}x_{1}x_{2}^{3}+x_{2}^{3}x_{3}\in\operatorname{\mathbb{C}}[t_{1}][x_{1},x_{2},x_{3}]$ and $\succ$ the total degree reverse lexicographic term order with the coordinate $x_{1}\succ x_{2}\succ x_{3}$.

A comprehensive Gröbner system of $\langle f,\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}}\rangle:\langle x_{1},x_{2},x_{3}\rangle^{\infty}$ w.r.t. $\succ$ is

[TABLE]

Hence,

1. $\cdot$

if $t_{1}$ belongs to $\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}^{2}+1)$, then \dim_{O}\Bigl{(}\operatorname{\mathbb{V}}(f,\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}})\Bigr{)}=0, 2. $\cdot$

if $t_{1}$ belongs to $\operatorname{\mathbb{V}}(t_{1}^{2}+1)$, then \dim_{O}\Bigl{(}\operatorname{\mathbb{V}}(f,\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}})\Bigr{)}\neq 0.

Therefore, for all $\bar{t}\in\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}^{2}+1)$, $\sigma_{\bar{t}}(f)$ has an isolated singularity at the origin $O$.

4.2. Improvement

We improve Algorithm 2-1 in computation speed. The following lemma alows us to devise an efficient and practical algorithm for computing the saturation $\langle F\rangle:{\mathfrak{m}}^{\infty}$.

Lemma 4.3**.**

Let $F\subset\operatorname{\mathbb{C}}[x]$ and ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle$. For all $\alpha_{1},\alpha_{2},\ldots,\alpha_{n}\in\operatorname{\mathbb{N}}\backslash\{0\}$,

[TABLE]

Proof.

There exists $\beta\in\operatorname{\mathbb{N}}$ such that

[TABLE]

Let $J=\langle x_{1}^{\alpha_{1}},x_{2}^{\alpha_{2}},\ldots,x_{n}^{\alpha_{n}}\rangle$ and \alpha=\max\Bigl{(}\{\alpha_{1},\alpha_{2},\ldots,\alpha_{n}\}\Bigr{)}. Obviously,

[TABLE]

Since $J\cdot{\mathfrak{m}}^{\beta}\subseteq{\mathfrak{m}}^{\beta}$, thus $\langle F\rangle:{\mathfrak{m}}^{\beta}\subseteq\langle F\rangle:J\cdot{\mathfrak{m}}^{\beta}$. Take a sufficiently large number $N$ such that $N>\alpha+\beta$, then ${\mathfrak{m}}^{N}\subseteq J\cdot{\mathfrak{m}}^{\beta}$. Hence,

[TABLE]

As $\langle F\rangle:{\mathfrak{m}}^{\beta}=\langle F\rangle:{\mathfrak{m}}^{N}$, we have $\langle F\rangle:{\mathfrak{m}}^{\beta}=\langle F\rangle:J\cdot{\mathfrak{m}}^{\beta}$. Therefore, \Bigl{(}\langle F\rangle:\langle x_{1}^{\alpha_{1}},\ldots,x_{n}^{\alpha_{n}}\rangle\Bigr{)}:{\mathfrak{m}}^{\infty}=\langle F\rangle:{\mathfrak{m}}^{\infty}. ∎

The lemma above leads the following procedure for computing $\langle F\rangle:{\mathfrak{m}}^{\infty}$.

1. Step 1:

Compute a basis $G$ of $\langle F\rangle:\langle x_{1}^{\alpha_{1}},x_{2}^{\alpha_{2}},\ldots,x_{n}^{\alpha_{n}}\rangle$. 2. Step 2:

Compute a basis $G^{\prime}$ of $\langle G\rangle:{\mathfrak{m}}^{\infty}$. 3. Return

$G^{\prime}$.

Notice that in the procedure above arbitrary positive integers $\alpha_{1},\ldots,\alpha_{n}$ can be used to compute $\langle F\rangle:{\mathfrak{m}}^{\infty}$. Our strategy of choosing the integers is the following.

Let $\displaystyle f=\sum_{{\bf r}\in\operatorname{\mathbb{N}}^{n}}a_{{\bf r}}x^{{\bf r}}$ be a non-zero polynomial in $\operatorname{\mathbb{C}}[x]$ and $F\subset\operatorname{\mathbb{C}}[x]$. We set

[TABLE]

and $\text{mdeg}_{x_{i}}(F):=\max(\{\text{mdeg}_{x_{i}}(g)|g\in F\})$ where $i\in\{1,\ldots,n\}$. In Algorithm 2-2, we take

[TABLE]

as $\alpha_{1},\ldots,\alpha_{n}$ to compute a basis of $\langle F\rangle:{\mathfrak{m}}^{\infty}$, namely, $\alpha_{1}=\alpha_{2}=\cdots=\alpha_{n}=\alpha$.

Lemma 4.3 and the strategy above yield the following improvement.

 

**Algorithm 2-2

 ** Input: $F=\{f_{1},f_{2},\ldots,f_{s}\}\subset\operatorname{\mathbb{C}}[t][x]$ s.t. $O\in\operatorname{\mathbb{V}}(F)$, ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle$.

           $\succ$: a term order on $\operatorname{Term}(x)$.

Output: $\operatorname{\mathbb{A}}^{\ast}\subset\operatorname{\mathbb{C}}^{m}$: For all $\bar{t}\in\operatorname{\mathbb{A}}^{\ast}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))=0$ (i.e., $\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F))$ has an isolated point at $O$). For all $\bar{t}\in\operatorname{\mathbb{C}}^{m}\backslash\operatorname{\mathbb{A}}^{\ast}$, $\dim_{O}(\operatorname{\mathbb{V}}(\sigma_{\bar{t}}(F)))\neq 0$.

BEGIN

$\operatorname{\mathbb{A}}^{\ast}\leftarrow\emptyset$;  \alpha\leftarrow\max\Bigl{(}\{\text{mdeg}_{x_{1}}(F),\text{mdeg}_{x_{2}}(F),\ldots,\text{mdeg}_{x_{n}}(F)\}\Bigr{)};

${\mathcal{G}}\leftarrow$ Compute a comprehensive Gröbner system of $\langle F\rangle:\langle x_{1}^{\alpha},x_{2}^{\alpha},\ldots,x_{n}^{\alpha}\rangle$ w.r.t. $\succ$;

while ${\mathcal{G}}\neq\emptyset$ **do

**  Select $(\operatorname{\mathbb{A}},G)$ from ${\mathcal{G}}$; ${\mathcal{G}}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}},G)\}$;

  ${\mathcal{G}}^{\prime}\leftarrow$ Compute a comprehensive Gröbner system of $\langle G\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$ on $\operatorname{\mathbb{A}}$;

    while ${\mathcal{G}}^{\prime}\neq\emptyset$ **do

**        Select $(\operatorname{\mathbb{A}}^{\prime},G^{\prime})$ from ${\mathcal{G}}^{\prime}$; ${\mathcal{G}}^{\prime}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

            if $\exists g\in G_{2}$ s.t. $g(O)\neq 0$ **then

**               $\operatorname{\mathbb{A}}^{\ast}\leftarrow\operatorname{\mathbb{A}}^{\prime}\cup\operatorname{\mathbb{A}}^{\ast}$;

            **end-if

**    **end-while

end-while

return** $\operatorname{\mathbb{A}}^{\ast}$;

END

We have implemented Algorithm 2-2 in the computer algebra system Risa/Asir [16]. We give some outputs of our implementation in the following examples.

Example 5*.*

Let $f=x_{1}^{3}x_{2}+t_{1}x_{1}^{2}x_{2}^{4}+x_{2}^{10}+t_{2}x_{2}^{11}\in\operatorname{\mathbb{C}}[t_{1},t_{2}][x_{1},x_{2}]$ and $V=\operatorname{\mathbb{V}}(\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}})$. Our implementation outputs the following.

1. $\cdot$

If $(t_{1},t_{2})$ belongs to $W=\operatorname{\mathbb{C}}^{2}\backslash\operatorname{\mathbb{V}}(4t_{1}^{3}+27,t_{2}),$ then $\dim_{O}(V)=0$, namely, the hypersurface $S$, defined by $f$, has an isolated singularity at $O$. 2. $\cdot$

If $(t_{1},t_{2})$ does not belong to $W$, then $\dim_{O}(V)\neq 0$. The hypersurface $S$ does not have an isolated singularity at $O$.

Example 6*.*

Let $f=x_{1}x_{3}^{2}+x_{1}^{4}+x_{2}^{4}+t_{1}x_{2}x_{3}^{2}+t_{2}x_{1}^{2}x_{2}^{2}\in\operatorname{\mathbb{C}}[t_{1},t_{2}][x_{1},x_{2},x_{3}]$ and $V=\operatorname{\mathbb{V}}(\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}})$. Our implementation outputs the following.

1. $\cdot$

If $(t_{1},t_{2})$ belongs to

[TABLE]

then $\dim_{O}(V)=0$, namely, the hypersurface $S$, defined by $f$, has an isolated singularity at $O$. 2. $\cdot$

If $(t_{1},t_{2})$ does not belong to $W$, then $\dim_{O}(V)\neq 0$. The hypersurface $S$ does not have an isolated singularity at $O$.

We will see in section 5 how the use of $\langle F\rangle:\langle x_{1}^{\alpha},x_{2}^{\alpha},\ldots,x_{n}^{\alpha}\rangle$ reduces the cost of computation of the saturation $\langle F\rangle:{\mathfrak{m}}^{\infty}$ drastically.

4.3. Primary ideal component at $O$

Let $\langle F\rangle$ be the ideal generated by $F$ s.t. $O\in\operatorname{\mathbb{V}}(F)$. Let

[TABLE]

be the primary ideal decomposition of the ideal $\langle F\rangle$. Let $S$ denote the saturation $\langle F\rangle:{\mathfrak{m}}^{\infty}$ where ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle$ is the maximal ideal in $\operatorname{\mathbb{C}}[x]$.

Assume that $O\notin\operatorname{\mathbb{V}}(S)$. Then, we have

[TABLE]

Therefore, the primary component $Q_{0}$ at $O$ of $\langle F\rangle$ such that $\operatorname{\mathbb{V}}(Q_{0})=\{O\}$ can also be computed by using the saturation $S=\langle F\rangle:{\mathfrak{m}}^{\infty}$.

The method above works for parametric cases, too.

Example 7*.*

Let $f=x_{1}^{3}x_{2}+x_{1}^{2}x_{2}^{4}+t_{1}x_{2}^{10}\in\operatorname{\mathbb{C}}[t_{1}][x_{1},x_{2}]$ and $F=\{\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}}\}$. Let $\succ$ be the total degree lexicographic term order with $x_{1}\succ x_{2}$. A comprehensive Gröbner system of $\langle F\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$ is

[TABLE]

1. (i)

A comprehensive Gröbner system of $\langle F\rangle:\langle 1\rangle$ on $\operatorname{\mathbb{V}}(t_{1})$ w.r.t. $\succ$ is

[TABLE]

where $Q^{\prime}_{0}=\{3x_{1}^{2}x_{2}+2x_{1}x_{2}^{4},x_{1}^{4}x_{2},x_{1}^{5},2x_{1}^{4}+53x_{1}^{3}x_{2}^{3},x_{1}^{3}+4x_{1}^{2}x_{2}^{3}+10x_{2}^{9}\}$. This means that if $t_{1}=0$, then $\langle Q_{0}^{\prime}\rangle$ is the primary ideal component such that $\operatorname{\mathbb{V}}(Q_{0}^{\prime})=\{O\}$. 2. (ii)

A comprehensive Gröbner system of

$\langle F\rangle:\langle 2910897t_{1}^{3}x_{1}+16385050t_{1}x_{2}+14895500,-9801t_{1}^{2}x_{1}+52855t_{1}x_{2}^{2}+48050x_{2}\rangle$

on $\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1})$ w.r.t. $\succ$ is

$\{(\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}),Q_{0}^{\prime\prime})\}$

where $Q_{0}^{\prime\prime}=\{3x_{1}^{2}x_{2}+2x_{1}x_{2}^{4},x_{1}^{4}x_{2},x_{1}^{5},2x_{1}^{4}+53x_{1}^{3}x_{2}^{3},-1331t_{1}^{3}x_{1}^{4}+(-32065t_{1}^{2}x_{2}^{2}+29150t_{1}$ $x_{2}+5300)x_{1}^{3}+21200x_{1}^{2}x_{2}^{3}+53000x_{2}^{9}\}$. Thus, if $t_{1}\neq 0$, $\langle Q_{0}^{\prime\prime}\rangle$ is the primary ideal component such that $\operatorname{\mathbb{V}}(Q_{0}^{\prime})=\{O\}$.

4.4. Local dimensions

In this subsection, we give an algorithm for computing the local dimension at $O$.

Let $F$ be a set of polynomials in $\operatorname{\mathbb{C}}[x]$. Let $E$ be a set of families of linear polynomials

[TABLE]

in $\operatorname{\mathbb{C}}(u)[x]$ where $u=\{u_{11},\ldots,u_{1n-\ell},u_{21},\ldots,u_{2\ell-n},\ldots,u_{\ell 1},\ldots,u_{\ell n-\ell}\}$, $\ell\leq n$ and $\operatorname{\mathbb{C}}(u)$ is the fields of rational functions with $u$. Then, we have the following proposition.

Proposition 4.4**.**

Let ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle\subset\operatorname{\mathbb{C}}[x]$ and $O\in\operatorname{\mathbb{V}}(F)$. Let $G$ be a basis of the ideal $\langle F\cup E\rangle:{\mathfrak{m}}^{\infty}$ in $\operatorname{\mathbb{C}}(u)[x]$. Let $\ell$ be the minimum number that satisfies

“there exists $g\in G$ such that $g(O)\neq 0$.”

Then, $\dim_{O}(\operatorname{\mathbb{V}}(F))=\ell$.

Proof.

By Theorem 4.1, if there exists $g\in G$ such that $g(O)\neq 0$, then $\dim_{O}(\operatorname{\mathbb{V}}(F\cup E))=\dim_{O}(\operatorname{\mathbb{V}}(F)\cap\operatorname{\mathbb{V}}(E))=0$. Note that $\operatorname{\mathbb{C}}(u)$ is the fields of rational functions and $E$ is a set of $\ell$ linear polynomials with $O\in\operatorname{\mathbb{V}}(E)$ Hence, $\dim_{O}(\operatorname{\mathbb{V}}(E))=n-\ell$. Since $\ell$ is the minimum number, it follows from the classical action lemma or the concept of Chevally dimension that $\dim_{O}(\operatorname{\mathbb{V}}(F))=\ell$. ∎

By this proposition, we can construct an algorithm for computing $\dim_{O}(\operatorname{\mathbb{V}}(F))$ as follows.

**Algorithm 2-3

 ** Input: $F=\{f_{1},f_{2},\ldots,f_{s}\}\subset\operatorname{\mathbb{C}}[x]$ s.t. $O\in\operatorname{\mathbb{V}}(F)$, ${\mathfrak{m}}=\langle x_{1},\ldots,x_{n}\rangle$,

           $\succ$: a term order on $\operatorname{Term}(x)$.

Output: $\dim_{O}(\operatorname{\mathbb{V}}(F))$.

BEGIN

$\ell\leftarrow 0$;  $\text{flag}\leftarrow 0$;  $E\leftarrow\emptyset$;  $U\leftarrow\emptyset$;

while $\text{flag}\neq 1$ do

   $G\leftarrow$ Compute a basis of $\langle F\cup E\rangle:{\mathfrak{m}}^{\infty}$ w.r.t. $\succ$ in $\operatorname{\mathbb{C}}(U)[x]$;   /if $U=\emptyset$, then $\operatorname{\mathbb{C}}(U)=\operatorname{\mathbb{C}}$./

         if $\exists g\in G$ s.t. $g(O)\neq 0$ **then

**                 $\text{flag}\leftarrow 1$;

         **else

**                 $\ell\leftarrow\ell+1$;

                 $E\leftarrow\{x_{1}+u_{11}x_{\ell+1}+\cdots+u_{1n-\ell}x_{n},x_{2}+u_{21}x_{\ell+1}+\cdots+u_{2n-\ell}x_{n},$

                                                                      $\cdots,x_{\ell}+u_{\ell 1}x_{\ell+1}+\cdots+u_{\ell n-\ell}x_{n}\}$;

                 $U\leftarrow\{u_{11},\ldots,u_{1n-\ell},u_{21},\ldots,u_{2\ell-n},\ldots,u_{\ell 1},\ldots,u_{\ell n-\ell}\}$;

         **end-if

end-while

return** $\ell$;

END

We illustrate Algorithm 2-3 with the following examples.

Example 8*.*

Let $f=x_{1}^{4}+x_{2}^{6}+2x_{1}^{2}x_{2}^{3}\in\operatorname{\mathbb{C}}[x_{1},x_{2}]$ and $F=\{\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}}\}$. Let $\succ$ be the total degree reverse lexicographic term order with $x_{1}\succ x_{2}$.

The reduced Gröbner basis of $\langle F\rangle:\langle x_{1},x_{2}\rangle^{\infty}$ w.r.t. $\succ$, in $\operatorname{\mathbb{C}}[x_{1},x_{2}]$, is

[TABLE]

Thus, $\dim_{O}(\operatorname{\mathbb{V}}(F))\neq 0$.

Next, let us consider the case $\ell=1$. Let $E=\{x_{1}+u_{11}x_{2}\}$. The reduced Gröbner basis of $\langle F\cup E\rangle:\langle x_{1},x_{2}\rangle^{\infty}$ w.r.t. $\succ$ is, $\operatorname{\mathbb{C}}(u_{11})[x_{1},x_{2}]$ is

[TABLE]

Hence, as $\dim_{O}(\operatorname{\mathbb{V}}(F\cup E))=0$, we obtain $\dim_{O}(\operatorname{\mathbb{V}}(F))=1$.

Example 9*.*

Let $f=x_{1}^{3}+x_{2}x_{3}^{2}+2x_{1}^{2}x_{2}^{2}+x_{1}x_{2}^{4}\in\operatorname{\mathbb{C}}[x_{1},x_{2},x_{3}]$ and $F=\{\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}}\}$. Let $\succ$ be the total degree reverse lexicographic term order with $x_{1}\succ x_{2}\succ x_{3}$.

The reduced Gröbner basis of $\langle F\rangle:\langle x_{1},x_{2},x_{3}\rangle^{\infty}$ w.r.t. $\succ$, in $\operatorname{\mathbb{C}}[x_{1},x_{2},x_{3}]$, is

[TABLE]

Thus, $\dim_{O}(\operatorname{\mathbb{V}}(F))\neq 0$.

Next, let us consider the case $\ell=1$. Let $E=\{x_{1}+u_{11}x_{2}+u_{12}x_{3}\}$. The reduced Gröbner basis of $\langle F\cup E\rangle:\langle x_{1},x_{2},x_{3}\rangle^{\infty}$ w.r.t. $\succ$, in $\operatorname{\mathbb{C}}(u_{11},u_{12})[x_{1},x_{2},x_{3}]$ is

[TABLE]

Hence, as $\dim_{O}(\operatorname{\mathbb{V}}(F\cup E))=0$, we obtain $\dim_{O}(\operatorname{\mathbb{V}}(F))=1$.

Algorithm 2-3 can be generalized to parametric cases. The key of the generalized method is to compute comprehensive Gröbner systems in $(\operatorname{\mathbb{C}}(u)[t])[x]$. We illustrate the method with the following example.

Example 10*.*

Let $f=x_{1}^{3}+x_{1}x_{3}^{2}+t_{1}x_{1}x_{2}^{3}+x_{2}^{3}x_{3}\in\operatorname{\mathbb{C}}[t_{1}][x_{1},x_{2},x_{3}]$ and $F=\{\frac{\partial f}{\partial x_{1}},\frac{\partial f}{\partial x_{2}},\frac{\partial f}{\partial x_{3}}\}$ where $t_{1}$ is a parameter. Let $\succ$ be the total degree reverse lexicographic term order with $x_{1}\succ x_{2}\succ x_{3}$.

1. (i)

Let us consider the case $\ell=0$. A comprehensive Gröbner system of $\langle F\rangle:\langle x_{1},x_{2},x_{3}\rangle^{\infty}$ w.r.t. $\succ$ is

[TABLE]

Hence, if the parameter $t_{1}$ belongs to $\operatorname{\mathbb{C}}\backslash\operatorname{\mathbb{V}}(t_{1}^{2}+1)$, then $\dim_{O}(\operatorname{\mathbb{V}}(F))=0$. 2. (ii)

Next, let us consider the case $\ell=1$. Let $E=\{x_{1}+u_{11}x_{2}+u_{12}x_{3}\}$. A comprehensive Gröbner system of $\langle F\cup E\rangle:\langle x_{1},x_{2},x_{3}\rangle^{\infty}\subset(\operatorname{\mathbb{C}}(u_{11},u_{12})[t_{1}])[x_{1},x_{2},x_{3}]$ w.r.t. $\succ$, on the stratum $\operatorname{\mathbb{V}}(t_{1}^{2}+1)$, is

$\{(\operatorname{\mathbb{V}}(t_{1}^{2}+1),\{(u_{12}^{6}+3u_{12}^{4}+3u_{12}^{2}+1)x_{3}+(-2u_{12}^{3}+6u_{12})u_{11}^{3}t_{1}+(-6u_{12}^{2}+2)u_{11}^{3},(-u_{12}^{4}-2u_{12}^{2}-1)x_{2}+(-2u_{12}^{2}+2)u_{11}^{2}t_{1}-4u_{12}u_{11}^{2},(u_{12}^{6}+3u_{12}^{4}+3u_{12}^{2}+1)x_{1}+(-6u_{12}^{2}+2)u_{11}^{3}t_{1}+(2u_{12}^{3}-6u_{12})u_{11}^{3}\})\}.$

The first polynomial $(u_{12}^{6}+3u_{12}^{4}+3u_{12}^{2}+1)x_{3}+(-2u_{12}^{3}+6u_{12})u_{11}^{3}t_{1}+(-6u_{12}^{2}+2)u_{11}^{3}$ is not zero at the origin $O$. Hence, the local dimension of $\operatorname{\mathbb{V}}(F)$ is equal to 1 on the stratum $\operatorname{\mathbb{V}}(t_{1}^{2}+1)$.

5. Comparisons

Here we give results of the benchmark tests. All algorithms in this paper have been implemented in the computer algebra system Risa/Asir [16]. All tests presented in Table 1, have been performed on a machine [OS: Windows 10 (64bit), CPU: Intel(R) Core i9-7900 CPU @ 3.30 GHz, RAM: 128 GB] and the computer algebra system Risa/Asir version 20150126 [16]. The time is given in second (CPU time). In Table 1, “$<$0.0156” means it takes less than 0.0156 seconds, and “$>3h$” means it takes more than 3 hours.

We use the total degree reverse lexicographic term order with $x\succ y\succ z$ (or $x\succ y$) in the benchmark tests. We use the following 10 polynomials.

 

         $f_{1}=x^{3}+xz^{2}+axy^{3}+y^{3}z+xy^{4}$

         $f_{2}=x^{3}y+ay^{15}+bxy^{11}+xy^{12}$

         $f_{3}=x^{4}y+y^{8}+axy^{8}+bx^{2}y^{4}$

         $f_{4}=x^{3}y+ay^{4}+y^{3}+y^{8}x+by^{6}$

         $f_{5}=x^{4}+yz^{5}+y^{4}+ax^{4}z+y^{2}z^{7}+z^{4}$

         $f_{6}=x^{5}y^{3}+z^{8}+axz^{8}+y^{6}z+byz^{5}$

         $f_{7}=(x^{2}y+z^{4}+y^{5})^{2}+ay^{6}z^{4}+y^{4}z^{6}$

         $f_{8}=x^{10}+x^{5}y^{3}+ay^{6}+3y^{14}+bx^{10}y^{5}+xy^{14}$

         $f_{9}=x^{5}+yz^{4}+y^{3}+ax^{5}y+bx^{2}y^{7}+z^{4}$

         $f_{10}=x^{6}+yz^{7}+ax^{3}y^{4}+y^{10}+x^{2}y^{5}z^{4}$

 

where $x,y,z$ are variables and $a,b$ are parameters.

As is evident from Table 1, in Problem 2, 3, 4, 5, 9, 10, 11, Algorithm 2-2 results in better performances in contrast to Algorithm 1 and Algorithm 2-1. In Problem 6, 7, 8, Algorithm 1 results in better performances in contrast to Algorithm 2-1 and Algorithm 2-2. Hence, we cannot say that which one is the best in general. However, as Algorithm 2-2 returns all results within 230 seconds, it is better to utilize Algorithm 2-2 in general.

If $\operatorname{\mathbb{V}}(F)$ has an irreducible component $V_{0}=\{O\}$, then the ideal $\langle F\rangle$ can be written as $\langle F\rangle=Q_{0}\cap Q_{1}\cap\cdots\cap Q_{\nu}$ where $Q_{0},Q_{1},\ldots,Q_{\nu}$ are distinct primary ideals and $\operatorname{\mathbb{V}}(Q_{0})=V_{0}$. Actually, Algorithm 1 computes the ideal $Q_{0}$ and its dimension. In contrast, Algorithm 2-1 and 2-2 compute the ideal $Q_{1}\cap\cdots\cap Q_{\nu}$. Hence, if the structure of $Q_{0}$ is complicated, then we can expect that the computation cost of Algorithm 2-2 is lower than that of Algorithm 1.

In the realm of symbolic computation, the standard basis is regarded as a classical or typical tool to handle ideals in local rings. However, to the best of our knowledge, no effective algorithm for computing standard bases of parametric ideals is known. In order to treat local dimensions for parametric cases, we utilize comprehensive Gröbner systems.

To conclude this paper, we emphasize again that even though the problems considered in the present paper are local in nature, the proposed algorithms resolve the problems in polynomial rings and they are free from standard bases and Mora’s reduction (Tangent cone algorithm [4, 8]).

Appendix A Ideal quotients with parameters

Several algorithms for computing a basis of an ideal quotient in a polynomial ring are introduced in some textbooks (cf. [1, 4]). As, in general, the algorithms utilize Gröbner basis computation, the algorithms can be naturally extended to the parametric cases by utilizing comprehensive Gröbner systems (see the appendix of [11]).

Here, we briefly describe an efficient algorithm for computing ideal quotients with parameters, that utilizes a comprehensive Gröbner system of a module.

Let e$${}_{1}=(1,0) and e$${}_{2}=(0,1). Then, \{\mbox{{\boldmathe}{}{1}},\mbox{{\boldmathe}{}{2}}\} is a free basis of $(\operatorname{\mathbb{C}}[x])^{2}$. Let $\succ$ be a term order on $\operatorname{Term}(x)$ and $\succ_{\bf m}$ be a POT (position over term) module order on $(\operatorname{\mathbb{C}}[x])^{2}$ with \mbox{{\boldmathe}{}{1}}>\mbox{{\boldmathe}{}{2}} and $\succ$. The following theorems are from [4, 6].

Theorem A.1**.**

Let $f_{1},\ldots,f_{s},q$ be non-zero polynomials in $\operatorname{\mathbb{C}}[x]$. Suppose $F\subset(\operatorname{\mathbb{C}}[x])^{2}$ is a $\operatorname{\mathbb{C}}[x]$-module generated by \{f_{1}\cdot\mbox{{\boldmathe}{}{1}},f_{2}\cdot\mbox{{\boldmathe}{}{1}},\ldots,f_{s}\cdot\mbox{{\boldmathe}{}{1}},q\cdot\mbox{{\boldmathe}{}{1}}-\mbox{{\boldmathe}{}{2}}\} and $G$ is a minimal Gröbner basis of $F$ w.r.t. $\succ_{\bf m}$. Set H=\{h\in\operatorname{\mathbb{C}}[x]|h\cdot\mbox{{\boldmathe}{}{2}}\in G\}. Then, $\langle f_{1},\ldots,f_{s}\rangle:\langle q\rangle=\langle H\rangle$.

There exists algorithms and implementations for computing a comprehensive Gröbner system of a given module with parameters (cf. [6, 9]). Hence, we are able to obtain a comprehensive Gröbner system of an ideal quotient with parameters.

Theorem A.2**.**

Let $y=\{y_{2},\ldots,y_{r}\}$ be new variables such that $t\cap x=\emptyset$. Let $f_{1},\ldots,f_{s},q_{1},\ldots,q_{r}$ be non-zero polynomials in $\operatorname{\mathbb{C}}[x]$. Set $q=q_{1}+y_{3}q_{2}+\cdots+y_{r}q_{r}$ and let $G$ be a Gröbner basis of the ideal quotient $\langle f_{1},\ldots,f_{s}\rangle:\langle q\rangle$ w.r.t. a block term order such that $y\gg x$ in $\operatorname{\mathbb{C}}[y,x]$. Then, $\langle f_{1},\ldots,f_{s}\rangle:\langle q_{1},\ldots,q_{r}\rangle=\langle G\cap\operatorname{\mathbb{C}}[x]\rangle$.

As we know how to compute a comprehensive Gröbner systems of $\langle f_{1},\ldots,f_{s}\rangle:\langle q\rangle$, Theorem A.2 also can be generalized to the parametric cases, too.

An algorithm for computing a comprehensive Gröbner system of an ideal quotient is the following.

Algorithm A (ideal quotients with parameters)

Input: $f_{1},\ldots,f_{s},q_{1},\ldots,q_{r}\in\operatorname{\mathbb{C}}[t][x]$ ($\forall\bar{t}\in\operatorname{\mathbb{C}}^{m}$, $1\leq\exists i\leq s$ s.t. $\sigma_{\bar{t}}(f_{i})\neq 0$).

          $\succ$: a block term order with $y\gg x$ on $\operatorname{Term}(x\cup y)$.

          $\succ_{{\bf m}}$: a POT module order on $(\operatorname{\mathbb{C}}[x])^{2}$ with \mbox{{\boldmathe}{}{1}}>\mbox{{\boldmathe}{}{2}} and $\succ$.

Output: ${\mathcal{Q}}$: a comprehensive Gröbner systems of $\langle f_{1},\ldots,f_{s}\rangle:\langle q_{1},\ldots,q_{r}\rangle$ w.r.t. $\succ$.

BEGIN

${\mathcal{Q}}\leftarrow\emptyset$;

$q\leftarrow q_{1}+y_{2}q_{2}+\cdots+y_{r}q_{r}$;

F\leftarrow\{f_{1}\cdot\mbox{{\boldmathe}{}{1}},f_{2}\cdot\mbox{{\boldmathe}{}{1}},\ldots,f_{s}\cdot\mbox{{\boldmathe}{}{1}},q\cdot\mbox{{\boldmathe}{}{1}}-\mbox{{\boldmathe}{}_{2}}\};

${\mathcal{G}}\leftarrow$ Compute a comprehensive Gröbner system of $\langle F\rangle$ w.r.t. $\succ_{{\bf m}}$ in $(\operatorname{\mathbb{C}}[t][x])^{2}$;

while ${\mathcal{G}}\neq\emptyset$ **do

**    Select $(\operatorname{\mathbb{A}},G)$ from ${\mathcal{G}}$; ${\mathcal{G}}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}},G)\}$;

    H\leftarrow\{h\in\operatorname{\mathbb{C}}[t][y,x]|h\cdot\mbox{{\boldmathe}{}_{2}}\in G\};

    ${\mathcal{Q}}\leftarrow{\mathcal{Q}}\cup\{(\operatorname{\mathbb{A}},H\cap\operatorname{\mathbb{C}}[t][x])\}$;

**end-while

return** ${\mathcal{Q}}$;

END

Let $I,J$ be ideals in $\operatorname{\mathbb{C}}[x]$. Since \Bigl{(}I:J\Bigr{)}:J=I:J^{2}, $I:J^{\infty}$ can be obtained by utilizing the algorithm above. Our implementation for saturation with parameters is given by the following algorithm.

Algorithm B (saturation with parameters)

Input: $f_{1},\ldots,f_{s},q_{1},\ldots,q_{r}\in\operatorname{\mathbb{C}}[t][x]$ ($\forall\bar{t}\in\operatorname{\mathbb{C}}^{m}$, $1\leq\exists i\leq s$ s.t.  $\sigma_{\bar{t}}(f_{i})\neq 0$).

          $\succ$ : a block term order with $y\gg x$ on $\operatorname{Term}(x\cup y)$.

Output: ${\mathcal{Q}}$: a comprehensive Gröbner systems of $\langle f_{1},\ldots,f_{s}\rangle:\langle q_{1},\ldots,q_{r}\rangle^{\infty}$ w.r.t. $\succ$.

BEGIN

${\mathcal{Q}}\leftarrow\emptyset$;

${\mathcal{G}}\leftarrow$ Compute a comprehensive Gröbner system of $\langle f_{1},\ldots,f_{s}\rangle:\langle q_{1},\ldots,q_{r}\rangle$ w.r.t. $\succ$;

while ${\mathcal{G}}\neq\emptyset$ **do

**    Select $(\operatorname{\mathbb{A}},G)$ from ${\mathcal{G}}$; ${\mathcal{G}}\leftarrow{\mathcal{G}}\backslash\{(\operatorname{\mathbb{A}},G)\}$;

    ${\mathcal{G}}^{\prime}\leftarrow$ Compute a comprehensive Gröbner system of $\langle G\rangle:\langle q_{1},\ldots,q_{r}\rangle$ on $\operatorname{\mathbb{A}}$ w.r.t. $\succ$;

    while ${\mathcal{G}}^{\prime}\neq\emptyset$ **do

**           Select $(\operatorname{\mathbb{A}}^{\prime},G^{\prime})$ from ${\mathcal{G}}^{\prime}$; ${\mathcal{G}}^{\prime}\leftarrow{\mathcal{G}}^{\prime}\backslash\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

                if $G=G^{\prime}$ do

                       ${\mathcal{Q}}\leftarrow{\mathcal{Q}}\cup\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

                **else

**                       ${\mathcal{G}}\leftarrow{\mathcal{G}}\cup\{(\operatorname{\mathbb{A}}^{\prime},G^{\prime})\}$;

                **end-if

**    **end-while

end-while**

return ${\mathcal{Q}}$;

END

As $\operatorname{\mathbb{C}}[t][x]$ is a Noetherian ring and an algorithm for computing comprehensive Gröbner systems always terminates, Algorithm B terminates.