new operator: linear system solving through LU decomposition