Conditioning and gauge freedom
Poor conditioning:
use 2nd order method
use Cholesky decomposition
Gauge freedom
fix certain parameters (orientation)  or
zero out last few rows in Cholesky decomposition