If someone is interested, I also wrote some code for rank 4 tensors since the pcout version only outputs everything to a single line:
template<int dim> void SolidMechanics<dim>::print_tensor (const SymmetricTensor<4, dim>& tensor, char* name) { pcout << std::endl << name << std::endl; for (unsigned int i = 0; i < dim; ++i) { for (unsigned int j = 0; j < dim; ++j) { for (unsigned int k = 0; k < dim; ++k) { for (unsigned int l = 0; l < dim; ++l) { pcout << tensor[i][j][k][l] << " "; } } pcout << std::endl; } } pcout << std::endl; } BR, Seyed Ali Mohseni -- 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.