If x is a matrix, the default returns mean-zero and variance-one columns.
If x is a simple_triplet_matrix, columns are scaled by sdev(x) but left unshifted; the function returns a simple_triplet_matrix. If shift and scale are specified AND x is a matrix, columns are shifted by -m and divided by s.