Dear dealii users, I want to construct a matrix [image: E=1/2(F^TF-1)]. My approach is as follows.
FullMatrix<double> F(3,3) = 0.0; // deformation gradient FullMatrix<double> E(3,3) = 0.0; // Green strain FullMatrix<double> I (IdentityMatrix(3)); F.Tmmult(E, F); E.add(-1, I); E.add(2, E); I wonder if there is an more elegant or easier way to do the matrix construction as I need to do other similar operations. Best, Michael -- 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/CAEyr2zhNC9RpZVjD%2B4ggHp-PCBCRz9YCtBLSyxxqqeCDE-giSA%40mail.gmail.com.
