Home Online Manual
Back: adjointIdeal
Forward: paraPlaneCurve
Up: paraplanecurves_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

D.5.4.2 invertBirMap

Procedure from library paraplanecurves.lib (see paraplanecurves_lib).

invertBirMap(phi, I); phi ideal, I ideal

The ideal phi in the basering R represents a birational map of the variety given by the ideal I in R to its image in projective space P = PP^(size(phi)-1).

The procedure might fail or give a wrong output if phi does not define a birational map.

ring, the coordinate ring of P, with an ideal named J and an ideal named psi.
The ideal J defines the image of phi.
The ideal psi gives the inverse of phi.
Note that the entries of psi should be considered as representatives of classes in the quotient ring R/J.

We compute the ideal I(G) in R**S of the graph G of phi.
The ideal J is given by the intersection of I(G) with S.
The map psi is given by a relation mod J of those relations in I(G) which are linear in the variables of R.

LIB "paraplanecurves.lib";
ring R = 0,(x,y,z),dp;
poly f = y^8-x^3*(z+x)^5;
ideal adj = adjointIdeal(f);
def Rn = invertBirMap(adj,ideal(f));
==> // 'invertBirMap' created a ring together with two ideals J and psi.
==> // Supposing you typed, say,  def RPn = invertBirMap(phi,I);
==> // you may access the ideals by typing
==> //      setring RPn; J; psi;
==> J[1]=y(5)*y(6)-y(4)*y(7)
==> J[2]=y(4)*y(6)-y(3)*y(7)
==> J[3]=y(2)*y(6)-y(1)*y(7)
==> J[4]=y(4)*y(5)-y(2)*y(7)
==> J[5]=y(3)*y(5)-y(1)*y(7)
==> J[6]=y(1)*y(5)-y(7)^2
==> J[7]=y(4)^2-y(1)*y(7)
==> J[8]=y(3)*y(4)-y(1)*y(6)
==> J[9]=y(2)*y(4)-y(7)^2
==> J[10]=y(1)*y(4)-y(6)*y(7)
==> J[11]=y(2)*y(3)-y(6)*y(7)
==> J[12]=y(1)*y(3)-y(6)^2
==> J[13]=y(2)^2-y(5)*y(7)
==> J[14]=y(1)*y(2)-y(4)*y(7)
==> J[15]=y(1)^2-y(3)*y(7)
==> J[16]=y(1)*y(6)^2-y(3)^2*y(7)
==> J[17]=y(6)^4-y(3)^3*y(7)
==> psi[1]=-y(6)^2
==> psi[2]=-y(4)*y(7)
==> psi[3]=y(6)^2-y(5)*y(7)