Currently,I have the stiffness matrix A an density vector v.

If I know the constructing process of cell matrix in code,how should I calculate the derivative of stiffness matrix with density?

A is a 4131*4131 matrix,v is 1024*1 vector. After derivative calculation the result is a 4131*4131*1024 matrix and a Jacobin matrix,if I am wrong ,please provide your suggestions.

This is conceptually correct, but it is not how you would store things. That's because while you could probably store 1024 matrices of size 4131x4131, you definitely can't do that any more if you refine once or twice more.

But in practice, you also don't have to do that. The density on cell k only affects the matrix entries (i,j) for degrees of freedom i,j that are located on cell k. As a consequence, the (i,j,k) entries of the 4131x4131x1024 object is only going to be nonzero if i,j,k all live on the same cell. In other words, the resulting object is going to be *very* sparse.

Of course, this only tells you what *not* to do, but you want to know what to do. This depends on what you need this derivative for. You only tell us *that* you need the derivative, but not why. Can you outline what it is you really want to do? In other words, what do you need the derivative for?

Best
 W.

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/8ae56dea-34cf-1eb6-bce5-61e164b10dd9%40colostate.edu.

Reply via email to