Hi there, I've been applying MatrixFree method to Step-29 where face integrals are involved, different from step-37, step-75 and so on. Firstly I think loop() instead of cell_loop() should be used. Secondly, similar to step-75, Chebyshev smoother and a coarsest level matrix are needed in my code. In step-75 we get diagonal vector by *MatrixFreeTools::compute_diagonal and MatrixFreeTools::compute_matrix*.
However, the face integral makes them different. Both two functions only take cell_operations(FEEvaluation), not face_operations(FEFaceEvaluation). My question is how to compute diagonal and matrix when my weak form has face integrals terms (only boundary integral, no interface integral). -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/c343562d-1d15-463a-8e26-cdbcbd14e463n%40googlegroups.com.