Hello Group, I have two questions regarding project_boundary_values for face elements.
1, I have a (dim)-vector face element, and I want to impose Dirichlet boundary conditions to the first component of the face element but not the rest. I can interpolate the boundary values using interpolate_boundary_values with a component_mask(std::vector<bool>) argument But in order to project the boundary values, I need to use component_mapping (std::vector<unsigned int>) argument in project_boundary_values. I found in the documentation of component_mapping in https://dealii.org/8.4.0/doxygen/deal.II/namespaceVectorTools.html#a0c684536027758309d42a546d8e6ca48 quite confusing, and no examples in the tutorial was given with a use of this argument. So, the question is: what's the equivalent version of VectorTools::interpolate_boundary_values(dof_handler, 0, ZeroFunction<dim>(dim), constraint_matrx, fe.component_mask(xx)) for project_boundary_values? 2, There is something wrong with FE_FaceP: Interpolate_boundary_values is not supported for FE_FaceP Project_boundary_values is not correctly implemented for FE_FaceP (project a ZeroFunction do not leads to a zero function) How shall I fix this bug? Best, Guosheng -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
