PierUgit

Some code (to be polished) for set/get diagonals (work for arbitrary matrices -not necessarily square- and arbitrary diagonal number) ```fortran !---------------------------------------- function sget_diagonal(A,k) result(d) integer, parameter :: wp=kind(0.0) real(wp), intent(in)...

> I don't use openmp much, but I believe every time there is a `static` (`save`) variable somewhere, that must be declared `THREADPRIVATE`, otherwise all threads will write to it,...