Add LU based matrix inversion
Add DOPs