default to substitution in matrix computation