F0, T1, and T2 should all be matrices over some common polynomial ring.
Each element of the sequence (F,R,G,C) is a list of matrices in increasing powers of the deformation parameter specified by DefParam. Their sums satisfy the deformation equation transpose flatten ((sum F)*(sum R))+(sum C)*(sum G)==0 up to powers of the deformation parameter equal to the length of F. Furthermore, F_0=F0, R_0=gens ker F0, the columns of C_0 are multiples of those of T2 and F_1 consists of first order perturbations corresponding to T1. Thus, if T1 and T2 are tangent and obstruction spaces for some deformation functor, then F and G represent a versal family and equations for a versal base space.
Several options are available to control the termination of the calculation. The calculation will terminate at the very latest after reaching order equal to the option HighestOrder, which has default value 20. If this order is reached, a warning message is generated. If PolynomialCheck is set to true, as is the default, then the algorithm will check if the present solution lifts to infinite order and terminate if this is the case. If SanityCheck is set to true, as is the default, then the algorithm will check that the present solution really does solve the deformation equation, and terminate with an error if this is not the case.
The option Verbose may be used to control the verbosity of the output. Its value should be an integer, with higher values corresponding to more verbose output. Default value is 0.
The option SmartLift is also available, which controls whether the algorithm spends extra time trying to find liftings which introduce no new obstructions at the next highest order. By default, this option is enabled. The option CorrectionMatrix may be used to control which liftings are considered.
After each step of lifting, the solution (F,R,G,C) to the deformation equation is cached. By default, it is stored in F0.cache#VersalDeformationResults but may stored elsewhere by setting the option CacheName to something other than null.
We may use this method to compute local multigraded Hilbert schemes. Here, we consider the Borel fixed ideal for the multigraded Hilbert scheme of the diagonal in a product of three projective planes, see [CS10]:
i1 : S=QQ[x1,x2,x3,y1,y2,y3,z1,z2,z3,Degrees=> {{1,0,0},{1,0,0},{1,0,0},{0,1,0},{0,1,0},{0,1,0},{0,0,1},{0,0,1},{0,0,1}}]; |
i2 : I=ideal {y1*z2, x1*z2, y2*z1, y1*z1, x2*z1, x1*z1, x1*y2, x2*y1, x1*y1, x2*y2*z2}; o2 : Ideal of S |
i3 : (F,R,G,C)=versalDeformation(gens I,normalMatrix({0,0,0},gens I), CT^2({0,0,0},gens I),Verbose=>2); Calculating first order relations Calculating standard expressions for obstructions Starting lifting Order 2 Order 3 Order 4 Order 5 Order 6 Solution is polynomial |
Local equations for the multigraded Hilbert scheme are
i4 : T=ring first G; |
i5 : sum G o5 = | -t_2t_3t_4+t_2t_4t_7+t_1t_3t_8-t_1t_7t_8-t_1t_3t_13+t_2t_3t_13+t_1t_7t | t_1t_3t_4-t_2t_3t_4-t_1t_7t_8+t_2t_7t_8-t_1t_3t_13+t_2t_3t_13+t_1t_7t_ | -t_1t_3t_16+t_2t_7t_16+t_1t_14t_16-t_2t_14t_16+t_3t_15t_16-t_7t_15t_16 | -t_1t_3t_18+t_2t_7t_18+t_1t_14t_18-t_2t_14t_18+t_3t_15t_18-t_7t_15t_18 | -t_2t_4t_17+t_1t_8t_17-t_1t_13t_17+t_2t_13t_17+t_4t_15t_17-t_8t_15t_17 | t_2t_4t_18-t_1t_8t_18+t_1t_13t_18-t_2t_13t_18-t_4t_15t_18+t_8t_15t_18 | -t_3t_4t_17+t_7t_8t_17+t_3t_13t_17-t_7t_13t_17+t_4t_14t_17-t_8t_14t_17 | -t_3t_4t_16+t_7t_8t_16+t_3t_13t_16-t_7t_13t_16+t_4t_14t_16-t_8t_14t_16 ------------------------------------------------------------------------ _13-t_2t_7t_13+t_3t_4t_15-t_4t_7t_15-t_3t_8t_15+t_7t_8t_15 | 13-t_2t_7t_13-t_1t_4t_14+t_2t_4t_14+t_1t_8t_14-t_2t_8t_14 | | | | | | | 8 1 o5 : Matrix T <--- T |
At this point, the multigraded Hilbert scheme has 7 irreducible components:
i6 : # primaryDecomposition ideal sum G o6 = 7 |
We may use this method to compute versal deformations of modules. Here, we consider versal deformations for a torsion free rank one module of an E6 singularity, see [Si01]:
i7 : S=QQ[x,y]/ideal {x^4+y^3}; |
i8 : f= matrix {{x,-y^2},{y,x^3}}; 2 2 o8 : Matrix S <--- S |
i9 : (F,R,G,C)=versalDeformation(f,CT^1(f),extMatrix(f),Verbose=>2); Calculating first order relations Calculating standard expressions for obstructions Starting lifting Order 2 Order 3 Order 4 Order 5 Solution is polynomial |
The cokernel of f is the module M1. Its versal deformation has a singular curve as its base space:
i10 : T=ring first G; |
i11 : ideal sum G 4 3 o11 = ideal (t - t , 0) 2 1 o11 : Ideal of T |
We may also compute local Quot schemes for modules with zero-dimensional support:
i12 : S = QQ[a, b]; |
i13 : f = matrix{{a, b, 0, 0}, {0, 0, a, b}}; 2 4 o13 : Matrix S <--- S |
i14 : (F,R,G,C)=versalDeformation(f,normalMatrix(f),extMatrix(f),Verbose=>2); Calculating first order relations Calculating standard expressions for obstructions Starting lifting Order 2 Order 3 Solution is polynomial |
We are considering the local Quot scheme for the surjection of S^4 to coker f. Equations for the base space are:
i15 : T=ring first G; |
i16 : ideal sum G o16 = ideal (t t - t t , t t - t t + t t - t t , - t t + t t - t t + 4 5 2 7 2 3 1 4 4 6 2 8 3 5 1 7 6 7 ----------------------------------------------------------------------- t t , - t t + t t ) 5 8 4 5 2 7 o16 : Ideal of T |