RE: [deal.II] Matrix Multiplication

2021-07-26 Thread Michael Li
nt::push_forward(S, F);} I just undefine the variable to circumvent the Sacado error.  #define ENABLE_SACADO_FORMULATION   Thanks,Michael   From: Jean-Paul PelteretSent: Monday, July 26, 2021 12:46 AMTo: dealii@googlegroups.comSubject: Re: [deal.II] Matrix Multiplication Hi Michael, I revisited t

Re: [deal.II] Matrix Multiplication

2021-07-26 Thread Jean-Paul Pelteret
block_component); > | > > > Thanks, > Michael > > From: Jean-Paul Pelteret <mailto:jppelte...@gmail.com> > Sent: Thursday, July 22, 2021 11:03 AM > To: dealii@googlegroups.com <mailto:de

RE: [deal.II] Matrix Multiplication

2021-07-22 Thread Michael Li
ofs_per_block(dof_handler_ref, dofs_per_block,  | ~~^ 1137 |    block_component);  |      Thanks,Michael From: Jean-Paul PelteretSent: Thursday, July 22, 2021 11:03 AMTo: dealii@googlegroups.comSubj

Re: [deal.II] Matrix Multiplication

2021-07-22 Thread Jean-Paul Pelteret
; To: dealii@googlegroups.com <mailto:dealii@googlegroups.com> > Subject: Re: [deal.II] Matrix Multiplication > > Dear Michael, > > The classes that better serve your purpose are the Tensor and SymmetricTensor > classes: > https://dealii.org/current/doxygen/deal.

Re: [deal.II] Matrix Multiplication

2021-07-21 Thread Jean-Paul Pelteret
Dear Michael, The classes that better serve your purpose are the Tensor and SymmetricTensor classes: https://dealii.org/current/doxygen/deal.II/classTensor.html https://dealii.org/current/doxygen/deal.II/classSymmetricTensor.html

[deal.II] Matrix Multiplication

2021-07-21 Thread Michael Lee
Dear dealii users, I want to construct a matrix [image: E=1/2(F^TF-1)]. My approach is as follows. FullMatrix F(3,3) = 0.0; // deformation gradient FullMatrix E(3,3) = 0.0; // Green strain FullMatrix I (IdentityMatrix(3)); F.Tmmult(E, F); E.add(-1, I); E.add(2, E); I wonder if there is an more