psi = jacobianDual phi
psi = jacobianDual(phi, X, T)
Let I be an ideal of R and let phi be the presentation matrix of I as a module. The symmetric algebra of I has the form
$Sym_R(I) = R[T_0..T_m]/ideal(T*phi)$
where the T_i correspond to the generators of I. If $X = matrix\{\{x_1..x_n\}\}$, with x_i \in{} R, and ideal X contains the entries of the matrix phi, then there is a matrix psi defined over R[T_0..T_m], called the Jacobian Dual of phi with respect to X, such that $T*phi = X*psi$. (the matrix psi is generally not unique; Macaulay2 computes it using Groebner division with remainder.)
In the form psi = jacobianDual phi, a new ring ST := S[T_0..T_m] is created, and the vector X is set to the variables of R. The result is returned as a matrix over ST. To do the computation in a ring previously defined computed, use the form psi = jacobianDual(phi, X,T); in this case, the matrices phi, X, T should all be defined over the same ring ST, the matrix T should be a row of variables of ST, and the matrix phi should have entries in a subring not involving the entries of T.
If I is an ideal of grade >=1 and ideal X contains a nonzerodivisor of R (which will be automatic if I has finite projective dimension) then ideal X has grade >= 1 on the Rees algebra. Since ideal(T*phi) is contained in the defining ideal of the Rees algebra, the vector X is annihilated by the matrix psi when regarded over the Rees algebra. If also the number of relations of I is >= the number of generators of I, this implies that the maximal minors of psi annihilate the x_i as elements of the Rees algebra, and thus that the maximal minors of psi are inside the ideal of the Rees algebra. In very favorable circumstances, one may even have the equality reesIdeal I = ideal(T*phi)+ideal minors(psi): For example:
Theorem (S. Morey and B. Ulrich, Rees Algebras of Ideals with Low Codimension, Proc. Am. Math. Soc. 124 (1996) 3653–3661): Let R be a local Gorenstein ring with infinite residue field, let I be a perfect ideal of grade 2 with n generators, and let phi be the presentation matrix of I. Let ell = ell(I) be the analytic spread. Suppose that I satisfies the condition G_{ell} or, equivalently, that the np sized minors of phi have codimension >p for 1<= p < ell. The following conditions are equivalent:
1) reesAlgebra I is CohenMacaulay and I_(nell)(phi) = I_1(phi)^{nell} 2) reductionNumber I < ell and I_(n+1ell)(phi) = I_1(phi)^{n+1ell} 3) reesIdeal I = symmetricAlgebraIdeal I + minors(n, jacobianDual phi)
We start with the presentation matrix phi of an (n+1)generator perfect ideal Such that the first row consists of the n variables of the ring, and the rest of whose rows are reasonably general (in this case random quadrics):









We first compute the analytic spread ell and the reduction number r


Now we can check the condition G_{ell}, first probabilistically

and now deterministically

We now check the three equivalent conditions of the MoreyUlrich Theorem. Note that since ell = n1 in this case, the second part of conditions 1,2 is vacuously satisfied, and since r<ell, the conditions must all be satisfied. We first check that reesAlgebra I is CohenMacaulay:



Finally, we wish to see that reesIdeal I is generated by the ideal of the symmetric algebra together with the jacobian dual:

We now compute the ideal J of the symmetric algebra; the call symmetricAlgebra I would return the ideal over a different ring, so we do it by hand:







The name Jacobian Dual comes from the case where phi is a matrix of linear forms the x_i are the variables of R, and the generators of I are forms, all of the same degree D; in this case Euler's formula sum(df_i/dx_j*xj) = Df can be used to express the entries of psi in terms of the derivatives of the entries of phi, at least when the degrees of the columns of phi are nonzero in the coefficient field.
Explicitly, let x_1,...,x_n be the variables of R, and let phi be a presentation matrix for I. Since all the f_i have the same degree, if follows that, for each j, the entries phi_(i,j) will all have the same degree, say D_j = deg phi_(i,j). Let ST be the polynomial ring R[T_0..T_m], where the T_i correspond to f_i, and let X=matrix{{x_1,...,x_n}}, and T=matrix{{T_0,...,T_m}} be row matrices over ST. In this case, by Euler's formula, we may take
psi_{k,j}=(1/D_j)*sum_i(d phi_{i,j}/d x_k*T_i),
The division with remainder step is usually fast, but if this ever becomes a bottleneck it would be possible to test for the degree condition and use Euler's formula in the case where it applies.
The object jacobianDual is a method function with options.