Hi Wolfgang I tried DataPostprocessorTensor as soon as the commits merged in. It is easy to use, I really like it! Thank you! Unfortunately there is still a bug that I couldn't figure out, which is basically the same as my previous post: after successfully writing the output, the program aborts with the following message:
-------------------------------------------------------- An error occurred in line <104> of file </home/jie/dealii/source/base/subscriptor.cc> in function void dealii::Subscriptor::check_no_subscribers() const The violated condition was: counter == 0 Additional information: (none) Stacktrace: ----------- #0 /home/jie/dealii-latest/lib/libdeal_II.g.so.9.0.0-pre: dealii::Subscriptor::check_no_subscribers() const #1 /home/jie/dealii-latest/lib/libdeal_II.g.so.9.0.0-pre: dealii::Subscriptor::~Subscriptor() #2 /home/jie/dealii-latest/lib/libdeal_II.g.so.9.0.0-pre: dealii::DataPostprocessor<2>::~DataPostprocessor() #3 ./main: IFEM::LinearElasticSolver<2>::output(unsigned int) const #4 ./main: IFEM::LinearElasticSolver<2>::runStatics(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) #5 ./main: main -------------------------------------------------------- Aborted (core dumped) This looks like the early destruction bug in step 6, but I couldn't think of any similar mistake in my program. I basically copied your StrainPostprocessor class in data_out_postprocessor_tensor_01.cc, and used it like this: std::vector<DataComponentInterpretation::DataComponentInterpretation> > interpretation(dim, > DataComponentInterpretation::component_is_part_of_vector); > data.add_data_vector(this->dofHandler, this->solution, > std::vector<std::string>(dim, "displacement"), interpretation); > StrainPostprocessor<dim> strain; > data.add_data_vector(this->dofHandler, this->solution, strain); > > data.build_patches(); > data.write_vtu(out); Thanks Jie -- 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.