Apply symmetric Gaussian elimination without row exchanges. At step , let be the leading diagonal entry of the remaining symmetric Schur complement. If , stop and report that is not positive definite. If , use it to eliminate the rest of its row and column. If all steps succeed, this constructs an LDL decomposition
where is unit lower triangular and has positive diagonal.
The test is correct from first principles. If all , then for every nonzero ,
because is invertible. Conversely, if is positive definite, its first pivot is , and completing the square gives
Choosing shows that the Schur complement is positive definite. Induction forces every pivot to be positive. This is also Sylvester's criterion.
At step , updating the remaining matrix costs arithmetic operations. Hence the total is
which proves the existence of the required algorithm.
Introduce a Lagrange multiplier for the normalization and vary
For a smooth variation with on , integration by parts gives
The fundamental lemma of the calculus of variations therefore yields the Euler-Lagrange equation
Multiplying by and integrating, while using
the divergence theorem and on the boundary give
The normalization is one, so the multiplier equals the stationary value:
The integrand has an order-two pole at and a simple pole at . Write it near zero as , where
Its residue at zero is
At the residue is
The residue theorem now gives
The condition says that is a regular value of . The regular level set theorem therefore makes a smooth submanifold of codimension one in , hence a smooth surface.
For
the gradient is
On , , so and cannot both vanish. Thus there, and the given set is a smooth surface.
Not every smooth surface in is a global zero set. Every set is closed because is continuous, whereas the open unit disk
is a smooth surface but is not closed in .
The map
is a nonzero linear functional: otherwise would lie in the radical of a bilinear form, contradicting nondegeneracy. Its kernel is , so the rank-nullity theorem gives . Antisymmetry gives , hence .
Suppose is orthogonal to every vector of . Since , it is also orthogonal to , and therefore to
For a nondegenerate bilinear form, : both sides have dimension one and the latter is contained in the former. Thus , proving that the restriction to is nondegenerate.
The space has dimension and again carries a nondegenerate antisymmetric form. Induction, starting from the zero-dimensional space, shows that is even. Therefore is even. Equivalently, every finite-dimensional symplectic vector space has even dimension.

There are unlisted articles, also show them or only show them.