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.
This is needed for both Inverse nodes form UFL (happens in Flux Vector Splitting for Maxwell DG) and for jacobians in sum factorization.