This paper presents an algorithmic proof of the General Néron Desingularization theorem for morphisms with large smooth loci, extending previous one-dimensional results to higher dimensions.
Contribution
It provides a constructive, algorithmic approach to Néron Desingularization for a broader class of algebraic morphisms.
Findings
01
Algorithmic proof of the General Néron Desingularization theorem.
02
Extension of results from one-dimensional to higher-dimensional cases.
03
Uniform version applicable to morphisms with big smooth locus.
Abstract
An algorithmic proof of the General N\'eron Desingularization theorem and its uniform version is given for morphisms with big smooth locus. This generalizes the results for the one-dimensional case.
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Constructive Néron Desingularization of algebras with big smooth locus.
Zunaira Kosar, Gerhard Pfister and Dorin Popescu
Zunaira Kosar, Abdus Salam School of Mathematical Sciences,GC University, Lahore, Pakistan
Dorin Popescu, Simion Stoilow Institute of Mathematics of the Romanian Academy, Research unit 5,
University of Bucharest, P.O.Box 1-764, Bucharest 014700, Romania
An algorithmic proof of the General Néron Desingularization theorem and its uniform version is given for morphisms with big smooth locus. This generalizes the results for the one-dimensional case (cf. [10], [7]).
*Key words * : Smooth morphisms, regular morphisms
Motivated to generalize Artin’s Approximation Theorem (cf. [2]) to excellent Henselian
rings the third author developed a powerful tool, the General Néron Desingularization (cf.
[14]). This result was discussed and used later by many other authors (cf. [1], [20], [19]). The proof of the Desingularization was not constructive. In this paper
we want to give an algorithm to compute Néron Desingularization for an important special case. We begin recalling some standard definitions.
A ring morphism u:A→A′ of Noetherian rings has regular fibers if for all prime ideals P∈SpecA the ring A′/PA′ is a regular ring, i.e. its localizations are regular local rings.
It has geometrically regular fibers if for all prime ideals P∈SpecA and all finite field extensions K of the fraction field of A/P the ring K⊗A/PA′/PA′ is regular.
A flat morphism of Noetherian rings u is regular if its fibers are geometrically regular. If u is regular of finite type then u is called smooth. A localization of a smooth algebra is called essentially smooth. A Henselian Noetherian local ring A is excellent if the completion map A→A^ is regular.
Theorem 1**.**
(General Néron Desingularization, Popescu [14], [15], André [1], Swan [20], Spivakovsky [19]) Let u:A→A′ be a regular morphism of Noetherian rings and B an A-algebra of finite type. Then any A-morphism v:B→A′ factors through a smooth A-algebra C, that is v is a composite A-morphism B→C→A′.
Constructive General Néron Desingularization for the case when the rings A and A′ are one-dimensional local rings, is given in [12], [10] and [7], the two dimensional case is partially done in [11]. The purpose of this paper is to find a constructive proof for the case when rings A and A′ are of dimension m and the smooth locus of B→A′ is big. We proceed using induction on the dimension of rings, with the induction step given in Proposition 3. In the Section 2 we prove a uniform General Néron Desingularization for m-dimensional local Cohen-Macaulay rings and some consequences of it. We also give an algorithm to find a uniform General Néron Desingularization using Singular.
1. Constructive Néron Desingularization
Let u:A→A′ be a flat morphism of Noetherian local rings of dimension m. Suppose that the maximal ideal m of A generates the maximal ideal of A′, A′ is Henselian and u is a regular morphism.
Let B=A[Y]/I, Y=(Y1,…,Yn). If f=(f1,…,fr), r≤n is a system of polynomials from I then we can define the ideal Δf generated by all r×r-minors of the Jacobian matrix (∂fi/∂Yj). After Elkik [4] let HB/A be the radical of the ideal ∑f((f):I)ΔfB, where the sum is taken over all systems of polynomials f from I with r≤n.
HB/A defines the non smooth locus of B over A.
B is standard smooth over A if there exists f in I as above such that B=((f):I)ΔfB.
The aim of this section is to give an algorithmic proof of the following theorem.
Theorem 2**.**
Any A-morphism v:B→A′ such that v(HB/AA′) is mA′-primary factors through a standard smooth A-algebra B′.
To prove the above theorem we need the following proposition.
Proposition 3**.**
Let A and A′ be Noetherian local rings of dimension m and u:A→A′ be a regular morphism. Suppose that A′ is Henselian. Let B=A[Y]/I, Y=(Y1,…,Yn), f=(f1,…,fr), r≤n be a system of polynomials from I as above, (Mj)j∈[q] some r×r-minors 111 We use the notation [q]={1,…,q}. of the Jacobian matrix (∂fi/∂Yj′), (Nj)j∈[q]∈((f):I) and set P:=∑j=1qNjMj. Let v:B→A′ be an A-morphism. Suppose that
(1)
there exist an element d∈A such that d≡P modulo I and
2. (2)
there exist a smooth A-algebra D and an A-morphism ω:D→A′ such that Imv⊂Imω+d2e+1A′ and for Aˉ=A/(d2e+1) (defining e by (0:Ade)=(0:Ade+1)) the map vˉ=Aˉ⊗Av:Bˉ=B/d2e+1B→Aˉ′=A′/d2e+1A′ factors through Dˉ=D/d2e+1D.
Then there exist a B-algebra B′ which is standard smooth over A such that v factors through B′.
Proof.
Let δ:B⊗AD≅D[Y]/ID[Y]→A′ be the A-morphism given by b⊗λ→v(b)ω(λ).
First we show that
δ factors through a special B⊗AD-algebra E of finite type.
Let the map Bˉ→Dˉ is given by Y→y′+d2e+1D. Thus I(y′)≡0 modulo d2e+1D. Since vˉ factors through ωˉ we see that ωˉ(y′+d2e+1D)=yˉ. Set y~=ω(y′). We get
y−y~=v(Y)−y~∈d2e+1A′n, let us say y−y~=de+1ν for ν∈deA′n.
We have Mj=detHj, where Hj is the matrix (∂fi/∂Yj′)i∈[r],j′∈[n] completed with some (n−r) rows from 0,1. Since d≡P modulo I we get P(y′)≡d modulo d2e+1 in D because I(y′)≡0 modulo d2e+1D. Thus P(y′)=ds for some s∈D with s≡1 modulo d.
Let Gj′ be the adjoint matrix of Hj and Gj=NjGj′. We have
GjHj=HjGj=MjNj\mboxIdn
and so
[TABLE]
But Hj is the matrix (∂fi/∂Yj′)i∈[r],j′∈[n] completed with some (n−r) rows from 0,1. Especially we obtain
[TABLE]
Then tj:=Hj(y′)ν∈deA′n
satisfies
[TABLE]
and so
[TABLE]
Let
[TABLE]
where Tj=(T1,…,Tr,Tj,r+1,…,Tj,n) are new variables. The kernel of the map
φ:D[Y,T]→A′ given by Y→y, Tj→tj contains h. Since
[TABLE]
and
[TABLE]
modulo higher order terms in Yj′−yj′′, by Taylor’s formula we see that for p=maxidegfi we have
[TABLE]
modulo h where Q∈T2D[T]r. We have f(y′)=de+1b for some b∈deDr. Then
[TABLE]
is in the kernel of φ. Indeed, we have spfi=de+1gi\mboxmoduloh because of (3) and P(y′)=ds. Thus
de+1φ(g)=de+1g(t)∈(h(y,t),f(y))=(0) and g(t)∈deA′r and so g(t)∈(0:A′de+1)∩deA′=0 because (0:A′de)=(0:A′de+1), the map u being flat. Set E=D[Y,T]/(I,g,h) and let ψ:E→A′ be the map induced by φ. Clearly, v factors through ψ because v is the composed map B→B⊗AD≅D[Y]/I→EψA′.
Now we show that
there exist s′,s′′∈E such that Ess′s′′ is standard smooth over A and ψ factors through Ess′s′′.
Note that the r×r-minor s′ of (∂g/∂T) given by the first r-variables T is from srp+(T)⊂1+(d,T) because Q∈(T)2. Then V=(D[Y,T]/(h,g))ss′ is smooth over D. We claim that I⊂(h,g)D[Y,T]ss′s′′ for some other s′′∈1+(d,T)D[Y,T]. Indeed, we have PI⊂(h,g)D[Y,T]s and so P(y′+s−1de∑j=1qGj(y′)Tj)I⊂(h,g)D[Y,T]s. Since P(y′+s−1de∑j=1qGj(y′)Tj)∈P(y′)+de(T)D[Y,T]s we get
P(y′+s−1de∑j=1qGj(y′)Tj)=ds′′ for some s′′∈1+(T)D[Y,T]s. It follows that s′′I⊂((h,g):d)D[Y,T]ss′. Thus s′′IV⊂(0:Vd)∩deV=0 because (0:Vd)∩deV=0, and V is flat over D and so over A. This shows our claim. It follows that
I⊂(h,g)D[Y,T]ss′s′′. Thus Ess′s′′≅Vs′′ is a B-algebra which is also standard smooth over D and A.
As ω(s)≡1 modulo d and ψ(s′),ψ(s′′)≡1 modulo (d,t), d,t∈mA′ we see that ω(s),ψ(s′),ψ(s′′) are invertible because A′ is local. Thus ψ (and so v) factors through the standard smooth A-algebra B′=Ess′s′′.
We choose γ1,γ2,…,γm∈v(HB/A)A′∩A such that γk for k∈[m] is a system of parameters in A, and γk=∑i=1qv(bi)zi(k), where zi(k)∈A′, bi∈HB/A. Set B0=B[Z(1),…,Z(m)]/(f(1),…,f(m)), where f(k)=−γk+∑i=1qbiZi(k)∈B[Z(k)], Z(k)=(Z1(k),…,Zq(k)),
and let v0:B0→A′ be the map of B-algebras given by Z(k)→z(k).
Changing B by B0 we may suppose that γk∈HB/A.
(**[13, Lemma 3.4]**) Let B1 be the symmetric algebra SB(I/I2) of I/I2 over222Let M b e a finitely represented B-module and Bm(aij)Bn→M→0 a presentation then SB(M)=B[T1,…,Tn]/J with J=({i=1∑naijTi}j=1,…,m).* B. Then HB/AB1⊂HB1/A and (ΩB1/A)γ is free over (B1)γ for any γ∈HB/A.*
2. (2)
(**[20, Proposition 4.6]**) Suppose that (ΩB/A)γ is free over Bγ. Let I′=(I,Y′)⊂A[Y,Y′], Y′=(Y1′,…,Yn′). Then (I′/I′2)γ is free over Bγ.
3. (3)
(**[16, Corollary 5.10]**) Suppose that (I/I2)γ is free over Bγ. Then a power of γ is in ((g):I)Δg for some g=(g1,…gr), r≤n in I.
Using (1) of Lemma 4 we reduce our proof to the case when ΩBγk/A for all k∈[m] are free over Bγk respectively.
Let B1 be given by Lemma 1. The inclusion B⊂B1 has a retraction w which maps I/I2 to zero. For the reduction we change B,v by B1,vw.
Using (2) of Lemma 4 we may reduce to the case when (I/I2)γk is free over Bγk for all k∈[m].
Since ΩBγk/A is free over Bγk we see using Lemma 2 that changing I with (I,Y′)⊂A[Y,Y′] we may suppose that (I/I2)γk is free over Bγk.
Now using Using (3) of Lemma 4 we will reduce further to the case when a power of γk is in ((f(k)):I)Δf(k) for some f(k)=(f1(k),…frk(k)), rk≤n from I.
We reduced to the case when (I/I2)γk is free over Bγk. Then it is enough to use Lemma 3.
Replacing B1 by B we may assume that a power dk of γk for all k∈[m] has the form dk≡Pk=∑i=1qkMi(k)Li(k)\mboxmoduloI,
for some rk×rk minors Mi(k)
of (∂f(k)/∂Y) and Li(k)∈((f(k)):I).
The Jacobian matrix (∂f(k)/∂Y) can be completed with (n−rk) rows from An obtaining a square n matrix Hi(k) such that detHi(k)=Mi(k).
This is easy using just the integers 0,1.
Set d=dm, f=f(m), r=rm, q=qm, Mi=Mi(m), Ni=Ni(m), Aˉ=A/d2e+1, Bˉ=Aˉ⊗AB, Aˉ′=A′/(d2e+1A′), vˉ=Aˉ⊗Av. Then we have d≡∑jMjNj modulo I.
Now we will use the induction on m.
Case I:m=0
If m=0 then A and A′ are Artinian local rings and u:A→A′ is a regular morphism. Then we are done by Corollary 3.3 [13].
Case II:m>0
Suppose by the induction hypothesis that we have a standard smooth Aˉ-algebra Dˉ≅(Aˉ[Z]/(gˉ))hˉMˉ, for Z=(Z1,…,Zp),gˉ=(gˉ1,…,gˉq) with q≤p, hˉ∈Aˉ[Z] and Mˉ a q×q-minor of (∂Z∂gˉ), such that the map vˉ:Bˉ→A′ˉ factors through Dˉ, let us say vˉ is the composite map Bˉ→DˉωˉA′ˉ.
Now let g∈A[Z]q be a lifting of gˉ and M the q×q-minor of (∂Z∂g) corresponding Mˉ. Take h∈A[Z] such that h lifts hˉ.
Then D≅(A[Z]/(g))hM is a standard smooth A-algebra and by the Implicit Function Theorem the map ωˉ can be lifted to ω:D→A′ since A′ is Henselian. It follows that
Imv⊂Imω+d2e+1A′. Applying Proposition 3 we get a B-algebra C smooth over A such that v factors through C, B→C→A′.
2. A uniform Néron Desingularization
Let u:A→A′ be a regular morphism of Cohen-Macaulay local rings of dimension m. Suppose that the maximal ideal m of A generates the maximal ideal of A′, A′ is Henselian and A, A′ have the same completions.
Let B=A[Y]/I, Y=(Y1,…,Yn), and for i∈[m] let f(i)=(f1(i),…,fri(i)), ri≤n be a system of polynomials from I. Let Mi be an ri×ri-minor of the Jacobian matrix (∂f(i)/∂Y) and Ni∈((f(i)):I), Pi=NiMi. Let v:B→A′/m3k+cA′ be an A-morphism for some k,c∈N. Suppose that v(N1M1,…,NmMm)A′/m3k+cA′⊃mkA′/m3k+cA′. Let y′∈An be a lifting of v(Y) to A and let di=Pi(y′). Then (d1,…,dm)A′/m3k+cA′⊃mkA′/m3k+cA′. Note that mk⊂(d1,…,dm)A+m3k+c⊂(d1,…,dm)A+m3(3k+c)+c⊂…. Thus mk⊂(d1,…,dm)A and it follows that (d1,…,dm)A′⊃mkA′ . Since A is Cohen-Macaulay we get d={d1,…,dm} regular sequence in A. Note that (d1,…,dm) is the ideal corresponding to v(P1,…,Pm)A′ by the isomorphism A/m3k+c≅A′/m3k+cA′.
Theorem 5**.**
There exists a B-algebra C which is standard smooth over A with the following properties.
(1)
Every A-morphism v′:B→A′ with v′≡v\mboxmodulo(d13,…,dm3)A′ factors through C.
2. (2)
Every A-morphism v′:B→A′ with v′≡v\mboxmodulom3kA′ factors through C.
3. (3)
There exists an A-morphism w:C→A′ which makes the following diagram commutative
[TABLE]
Proof.
Let v′:B→A′ be an A-morphism with v′≡v\mboxmodulo(d13,…,dm3)A′. We apply induction on m.
Case I:m=1
If m=1 then A and A′ are Noetherian local rings of dimension 1 and u:A→A′ is a regular morphism. Then we are done by Theorem 2 [7], with e=1.
Case II:m>1
Now let Aˉ=A/(d13,…,dm−13), and consider the map v′ˉ=Aˉ⊗Av′:Bˉ=Aˉ⊗AB→Aˉ′=Aˉ⊗AA′. By the induction hypothesis there exists a standard smooth algebra Dˉ≅(Aˉ[Z]/(gˉ))hˉMˉ, for Z=(Z1,…,Zp),gˉ=(gˉ1,…,gˉq) with q≤p, hˉ∈Aˉ[Z] and Mˉ a q×q-minor of (∂Z∂gˉ), such that the map v′ˉ factors through Dˉ, say v′ˉ is the composite map Bˉ→Dˉω′ˉA′.
Now let g∈A[Z]q be a lifting of gˉ and M the q×q-minor of (∂Z∂g) corresponding Mˉ. Take h∈A[Z] such that h lifts hˉ.
Then D≅(A[Z]/(g))hM is a standard smooth A-algebra and by the Implicit Function Theorem the map ω′ˉ can be lifted to ω′:D→A′ since A′ is Henselian. It follows that
Imv′⊂Imω′+dm3A′. Applying Proposition 3 (with e=1) we get a B-algebra C standard smooth over A such that v′ factors through C. This proves (1) which obviously implies (2).
Now for (3) take the map w^:C≅(D[Y,T]/(I,g,h))ss′→A′/mcA′ given by (Y,T)→(y′,0). Then the composite map B→Cw^A′/mcA′
is lifted by v. Since C is standard smooth, we may lift w^ to an A-morphism w:C→A′ by the Implicit Function Theorem. Clearly, w makes the above diagram commutative.
Example 6**.**
(Rond) Let k be a field, A=k[[x]], x=(x1,x2,x3), B=A[Y]/(f), Y=(Y1,…,Y4), f=Y1Y2−Y3Y4. Then Δf=HB/A=(Y). Let p∈N and set y1′=x1p, y2′=x2p, y3′=x1x2−x3p. Then there exists y4′∈A such that f(y′)≡0 modulo (x)p2. It follows that d1=x1p, d2=x2p and d3=x3p2 belongs to Δf(y′) because x1px2p−x3p2=y3′(x1p−1x2p−1+x1p−2x2p−2x3p+…x3p2). For k=2p+p2−2 we have (x)k⊂(d1,d2,d3)⊂HB/A(y′). If f(y′)≡0 modulo (x)3k+p+1 then by Theorem 5 (3) we could get y∈A4 such that f(y)=0 and y≡y′ modulo (x)p+1. But this is not the case, since f(y′)≡0 modulo (x)p2 and we cannot apply the quoted theorem. Thus it is not a surprise that [18, Remark 4.7] says that there exist no y∈A4 such that f(y)=0 and y≡y′ modulo (x)p+1.
Corollary 7**.**
(Elkik) Let (A,m) be a Cohen-Macaulay Henselian local ring of dimension m and B=A[Y]/I, Y=(Y1,…,Yn) an A-algbra of finite type. Then for every k∈N there exist two integers
m0,p∈N such that if y′∈An satisfies mk⊂HB/A(y′) and I(y′)≡0 modulo mm for some m>m0 then there exists y∈An such that I(y)=0 and y≡y′ modulo mm−p.
Proof.
Suppose that A′=A. In the notation of Theorem 5 given k set m0=p=3k and suppose that y′∈An satisfies mk⊂HB/A(y′) and I(y′)≡0 modulo mm for some m>m0. Let v:B→A/mm be given by Y→y′. Set c=m−p. By Theorem 5 there exists a smooth
A-algebra C and a map w:C→A which makes the above digram commutative. Let y be the image of Y by the composite map B→CwA. Then I(y)=0 and y≡y′ modulo mc=mm−p.
Corollary 8**.**
With the assumptions and notation of the Theorem 5, let ρ:B→C be the structural algebra map. Then ρ induces bijections ρ∗ given by ρ∗(w)=w∘ρ, between
(1)
{w∈HomA(C,A′):w∘ρ≡v\mboxmodulo(d13,…,dm3)A′}* and
{v′∈HomA(B,A′):v′≡v\mboxmodulo(d13,…,dm3)A′}*
2. (2)
{w∈HomA(C,A′):w∘ρ≡v\mboxmodulom3kA′}* and
{v′∈HomA(B,A′):v′≡v\mboxmodulom3kA′}*
Proof.
We will use induction on the dimension of A.
(1)
Case I:m=1
If m=1 then A and A′ are Cohen-Macaulay local rings of dimension 1. Then we are done by Corollary 8 [7].
2. (2)
Case II:m>1
By Theorem 5, (1), ρ∗ is surjective. Let Now let Aˉ=A/(d13,…,dm−13), and the map v′ˉ=Aˉ⊗Av′:Bˉ=Aˉ⊗AB→Aˉ′=Aˉ⊗AA′. By the induction hypothesis there exists a standard smooth algebra Dˉ such that the maps wˉ and wˉ′ restricted to Dˉ coincide. This implies that w∣D and w′∣D lift the same map wˉ∣Dˉ. Thus w∣D=w′∣D by uniqueness in the Implicit Function Theorem.
By construction C=Ess′, E=D[Y,T]/(I,g,h) and
Hm(y′)(w(Y)−w′(Y))≡dm2(w(T)−w′(T))\mboxmoduloh. Thus dm2(w(T)−w′(T))=0 and so w∣E=w′∣E because dm is regular in A′ since dm is regular in A and u is flat. It follows that w=w′
3. (3)
Apply Theorem 5 (2) for the surjectivity. The injectivity follows from above.
Corollary 9**.**
With the assumptions and notation of the above Corollary, the following statements hold:
(1)
If there exists an A-morphism v~:B→A′
with v~≡v\mboxmodulo(d13,…,dm3)A′, then
there exists a unique A-morphism w~:C→A′ such that w~∘ρ=v~.
2. (2)
If there exists an A-morphism v~:B→A′
with v~≡v\mboxmodulom3kA′, then
there exists a unique A-morphism w~:C→A′ such that w~∘ρ=v~.
For the proof take w~=ρ∗−1(v~), where ρ∗ is defined in the Corollary 8.
By construction, C has the form (D[T]/(g))Mh, where M=det(∂gi/∂Tj)i,j∈[r] and h=s′∈A[T] satisfies w~(h)∈mA′.
Note that {w′∈HomA(C,A′):w′≡w~\mboxmodulo(d13,…,dm3)A′} is in bijection with the set of all t∈A′n such that g(t)=0 and t≡w~(T) modulo (d13,…,dm3)A′n.
Set V=(T1,…,Tr), Z=(Tr+1,…,Tn).
Thus g(U,w′(Z))=0 has a unique solution (namely U=w′(V)) in w~(Z)+(d13,…,dm3)A′n−r by the Implicit Function Theorem.
Consequently, w′(V) is uniquely defined by w′(Z), that is by the restriction w′∣A[Z].
Therefore, {w′∈HomA(C,A′):w′≡w~\mboxmodulo(d13,…,dm3)A′n} is in bijection with {w′′∈HomA(A[Z],A′):w′′≡w~∣A[Z]\mboxmodulo(d13,…,dm3)A′n−r}, the latter set being in bijection with w~(Z)+(d13,…,dm3)A′n−r, that is with (d13,…,dm3)A′n−r.
The proof of (2) goes similarly.
Theorem 11**.**
With the assumptions and notation of Corollary 8 there exist canonical bijections
(1)
[TABLE]
2. (2)
[TABLE]
For the proof apply Corollary 8 and the above lemma.
3. Algorithms
In this section we present the algorithms corresponding to the results of Sections 1 and 2.
We will use in our algorithm for uniform desingularization the following algorithm for the one dimensional case (cf. [7]):
Next we present the algorithm for uniform desingularization for the higher dimensional case.
Bibliography20
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] M. André, Cinq exposés sur la désingularisation , Handwritten manuscript Ecole Polytechnique Fédérale de Lausanne, (1991).
2[2] M. Artin, Algebraic approximation of structures over complete local rings , Publ. Math. IHES, 36 , (1969), 23-58.
3[3] W. Decker, G.-M. Greuel, G. Pfister, H. Schönemann: Singular 4-1-0 — A computer algebra system for polynomial computations. http://www.singular.uni-kl.de (2016).
4[4] R. Elkik, Solutions d’équations à coéfficients dans un anneaux hensélien , Ann. Sci. Ecole Normale Sup., 6 (1973), 553-604.
5[5] M. Greenberg, Rational points in henselian discrete valuation rings , Publ. Math. IHES, 31 , (1966), 59-64.
6[6] A. Khalid, Z. Kosar, An Easy proof of the General Neron Desingularization in dimension 1 , Bull. Math. Soc. Sci. Math. Roumanie, 59 (107) , 2016, 349-353, ar Xiv:AC/1605.04677.
7[7] A. Khalid, G. Pfister, D. Popescu, A uniform General Neron Desingularization in dimension one , J. of Algebra and Its Applications, 16 , (2018), ar Xiv:AC/1612.03416.
8[8] A. Khalid, A. Popescu, D. Popescu, Algorithms in the classical Néron Desingularization , to appear in Bull. Math. Soc. Sci. Math. Roumanie, ar Xiv:AC/1702.01445.