Sign in to follow this  
Dirk Gregorius

Another matrix identity

Recommended Posts

can't you just do the multiplication yourself with a general matrix A?

The diagonal matrix scales the elements of one of the matrices and then it's just a normal multiplication of two matrices.

Share this post

Link to post
Share on other sites
Thanks for the help guys. Of course the resulting matrix is not diagonal, at least not in general. The background is solving is system DAEs like you have to to in constrained dynamics. In one step the ODE and OpenTissue just calculated the reciprocal elements of A, where A is the product as described:

A = J * W * JT // J = Jacobian, W = Inverse mass matrix

The reason that only the reciprocal is needed is that the Gauss-Seidel solver divdes by A[i][i].

So what actually is done to prepare for solving the DAE is that the equation is transformed such that the transpose of J is never needed.

For those who are interessted:

A = J * W * JT

1) W is block diagonal => W = WT
2) Matrix Multiplication is associativ

<=> A = J * ( W * JT )
<=> A = J * ( WT * JT )
<=> A = J * Transposed( W * J )

This way I save creating the transpose of the Jacobian...


Share this post

Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

Sign in to follow this