OpenMP parallelized matrix assembly working now. I wish I could get
around the reduction slowdown but that seems to be unavoidable.