¿Cómo redondear los elementos matriciales en sympy?

Tal como lo conocemos

from sympy import * x = sin(pi/4) y = sin(pi/5) A = Matrix([x, y]) print(x) print(A.evalf()) 

muestra

 sqrt(2)/2 Matrix([[0.707106781186548], [0.587785252292473]]) 

Asi que

 print(round(x.evalf(), 3)) print(round(y.evalf(), 3)) 

muestra

 0.707 0.588 

Pero, ¿cómo podemos redondear todos los elementos en una Matriz de manera concisa, para que

 print(roundMatrix(A, 3)) 

puede muestra

 Matrix([[0.707], [0.588]]) 

¿Por qué no utiliza el método evalf con evalf(3) como evalf(3) ?

 from sympy import * x = sin(pi/4) y = sin(pi/5) A = Matrix([x, y]) print(x) print(A.evalf(3)) 

Salida

 sqrt(2)/2 Matrix([[0.707], [0.588]])