Main Article Content
A smooth version of Sylvester's law of inertia is presented for symmetric matrix functions of constant rank. The techniques used in the proof are constructive but the resulting numerical approaches are unstable, and therefore require stabilization. Two different stabilization techniques are suggested, one based on a descent method and one based on Newton's method. Some numerical tests are included to demonstrate the applicability of the obtained numerical methods.