Hello all, I am solving for a 3 dimensional displacement vector (*dim+1* dimensional) in a *dim=2* dimensional space (for a planer problem in elasticity). Hence I consider an FEvaluesExtractor:
const FEValuesExtractors::Vector u_fe; and initialize it as u_fe(0) Now I need to compute the gradient of the shape functions and the gradient of the solution vectors. Hence I use the command: for (unsigned int q_point = 0; q_point < n_q_points; ++q_point) for (unsigned int k = 0; k < dofs_per_cell; ++k) fe_values_ref[u_fe].gradient(k, q_point) Then I get 2x2 tensors. Could anyone suggest how to get gradient terms corresponding to all three vectors and finally get a 3x3 tensors, instead of 2x2. I am expecting to get a 3x3 tensor with non-zero entries in the first two columns and all zero entries in the third column. Thanks and regards, Anup. -- 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. For more options, visit https://groups.google.com/d/optout.
