Solved problem in Matrix^T * Matrix, never compute full transposed matrix.
Made On, Extend, Div, Grad and If use my own kernels for values instead of spMv operations, which was a hotspot for explicit sims'
Adding diagonalMatrix * matrix kernels to reduce calls to hotspot operation by half.