replace use of inverse() with LU factorization in Natarajan regularization
incorporate Chen variable transform into Natarajan regularization for a more consistent interface