Intel MKL provides an additional GEMM variant, GEMMT, that updates only the upper or lower triangular part of the result matrix. This would be a great addition to OpenBLAS.

https://software.intel.com/en-us/node/590047