changed it to refresh LU factorization if Newton iteration >  6