My code resolves a dynamic problem with coupled equation, so I study this problem in the modal space.
I have the matrix of MODAL solution (each column is the solution for time t*, with t* = 0,1,2,.. T) but I want to know the NODAL SOLUTION, so I wrote this code to execute inverse modal transformation (PHI is the modal matrix while BI2 is a translation coordinates matrix, such as g = BI2*PHI*y, where g is in the nodal space while y in the modal):
% y is the modal solution matrix % g is the nodal solution matrix g = zeros(size(y)); for i=1:length(y) g(:,i) = BI2*PHI*y(:,i); end
In fact, the modal inverse relation is true column by column [y(:,i) -> g(:,i)].
How can I improve my code? This for loop is very slow when I have a lot of dof.. I think I have to rewrite it, but I don't know how.
Thanks in advance! Pinco
No products are associated with this question.