Some part of Sage are using multiprocessing, like tensor calculus on manifolds, which is turned on by Parallelism().set(nproc=n) where n is the number of (effective) cores in your CPU (this works with all versions of Sage, including the distribution binaries). But I am not sure if matrix operations use it.