Dear all, I'm trying to use a Matrix-based solution for a general 3D FE elastic problem (just like example-8). In that example the authors used a linear scaler solution for calculating the stiffness matrix of each element; but, due to some reasons I want to use a matrix-based solution so that the stiffness matrix of each element can be calculated as follows: K_cell = B^T . D . B, where B matrix follows from the diļ¬erential operator for strain calculation and D is the usual elasticity matrix for 3D isotropic problems. The size of the B and D matrices are 6X3 and 6X6 respectively. I defined these matrices as "vector" of "FullMatrix<double>" and try to do operate on them by FullMatrix operators as follows: B[i].Tmmult(tmpMat,D); tmpMat.mmult(BDB,B[j]); and it works correctly; however, it works rather slowly. My question is, is "FullMatrix" the most efficient and fast way in dealii to evaluate such like matrix calculations? Are there any alternatives?
Best regards, Masoud -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/711e55f4-df0c-4673-8c29-b9cba6cc2185n%40googlegroups.com.
