I know I can do it component by component and then construct a matrix out of the output. But it would be nice if I could just say matrix.trig_reduce() and get a matrix with all the components trig_reduced. Thanks in advance
You can use the apply_map method which applies a function to each of the elements of the matrix. For example,
sage: m = matrix([[sin(x), cos(x)], [sin(x), cos(x)]]); m
[sin(x) cos(x)]
[sin(x) cos(x)]
sage: o = m*m.transpose(); o
[sin(x)^2 + cos(x)^2 sin(x)^2 + cos(x)^2]
[sin(x)^2 + cos(x)^2 sin(x)^2 + cos(x)^2]
sage: o.apply_map(lambda x: x.trig_reduce())
[1 1]
[1 1]
You could equivalently use attrcall
sage: o.apply_map(attrcall('trig_reduce'))
[1 1]
[1 1]
That solves my problem I have a small function written which does, but would it be very difficult to define trig_reduce(), full_simplify() etc for matrices in sage. I don't know of any sage webpage for feature requests. But this is a feature that would make things a lot easier and user friendly.
I have to do this often, so I wrote a function:
def matrix_full_simplify(mat,m,n):
matsimp=mat;
for i in range(m):
for j in range(n):
matsimp[i,j]=mat[i,j].full_simplify();
return matsimp;
Use:
amat = matrix_full_simplify(amat)
You could easily adapt this to be better or require fewer arguments and can create a similar function for vectors.
The trig one is possible, the other one is in the ticket. For posterity: As it turns out, we already have it, just not with all possible aliases! I didn't even know this.
sage: m = matrix([[sin(x), cos(x)], [sin(x), cos(x)]]); m
[sin(x) cos(x)]
[sin(x) cos(x)]
sage: o = m*m.transpose()
sage: o.simplify_trig()
[1 1]
[1 1]
Asked: 2010-11-25 17:54:20 +0100
Seen: 3,110 times
Last updated: May 16 '11
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
See also http://trac.sagemath.org/sage_trac/ticket/10552