Sorry, I fogot the subjet I found that performing programing the matrix product myself I get the desired result (Sum_k A(i,k)*B(k,i)). But is there a more elegant way? Javier Ros On Tue, 2006-03-07 at 19:31 +0100, Javier Ros Ganuza wrote:
Hello! :)
When performing matrix multiplications expressions get expanded like in the following Example.
Is there a way to avoid expansion, that is to leave the parenthesis
Thanks
Javier Ros
Example: ----------------- matrix aux1
[[-d331*sin(q3)*dq3-d332*cos(q3)*dq3],[d331*cos(q3)*dq3-d332*sin(q3)*dq3],[0]]
matrix R23
[[cos(q3),-sin(q3),0],[sin(q3),cos(q3),0],[0,0,1]]
matrix R12.mul(aux1)
[[d332*sin(q2)*sin(q3)*dq3-d331*cos(q2)*sin(q3)*dq3-d332*cos(q3)*cos(q2)*dq3-d331*sin(q2)*cos(q3)*dq3],[-d332*cos(q2)*sin(q3)*dq3-d331*sin(q2)*sin(q3)*dq3+d331*cos(q3)*cos(q2)*dq3-d332*sin(q2)*cos(q3)*dq3],[0]] -----------------
What I want is to obtain something like
cos(q3)*(-d331*sin(q3)*dq3-d332*cos(q3)*dq3)+.... and so on for the first (0,0) matrix element.
_______________________________________________ GiNaC-list mailing list GiNaC-list@ginac.de https://www.cebix.net/mailman/listinfo/ginac-list