tested pragmas in transpose func, failed due to index counter being iterated mid loop by another thread. diag and h^2 matrices work fine under threading.