Hello,

Recently I've been studying the problem of quadratic form diagonalization. Suppose that we have a form $F(x,y,z)$ with corresponding symmetric matrix $M$. This form is equivalent to another form $F'(x',y',z')$ which is in Smith normal form. The matrix that corresponds to $F'$ is $M'$ and it is diagonal. Note that both $F$ and $F'$ have integral coefficients. It is known that there exist two matrices $A$ and $B$ (which coefficients are integers as well) such that $M' = AMB$.

Now, do I understand correctly that $F$ and $F'$ represent exactly the same set of integers? If so, am I right that there must exist a change of variables of the form

$$ x = a_{11}x' + a_{12}y' + a_{13} z' $$ $$ y = a_{21}x' + a_{22}y' + a_{23} z' $$ $$ z = a_{31}x' + a_{32}y' + a_{33} z' $$

that transforms $F$ into $F'$? Finally, is it possible to deduce the matrix $(a_{ij})$ given $M'$, $M$, $A$ and $B$?

P.S. Here's an example. I have the following matrix $M$ of determinant 1:

$$ 5, 13, 1 $$

$$ 13, 34, 0 $$

$$ 1, 0, 35 $$

It corresponds to $M' = I$ where $I$ is an identity matrix. The Sage command M.smith_form() produced the following $A$:

$$ 0, 0, 1 $$ $$ 0, 1, 0 $$

$$ 1, 0, 0 $$

and $B$:

$$ -34, -455, 1190 $$ $$ 13, 174, -455 $$ $$ 1, 13, -34 $$

Now I need to find $(a_{ij})$.