Implement vectorized matrix inversion

This is needed for both Inverse nodes form UFL (happens in Flux Vector Splitting for Maxwell DG) and for jacobians in sum factorization.