Hi, I have a question about the Dirichlet boundary conditions.
I'm using fe(FE_DGP<dim>(1),1, FE_FaceP<dim>(1),1) as the finite element space, FE_DGP is defined in interiors, FE_FaceP is defined on faces. When I try to set the Dirichlet boundary conditions to the system by *VectorTools::interpolate_boundary_values* to *constraints*, I got the error message: *An error occurred in line <1080> of file </var/folders/8z/hlb6vc015qjggytkxn84m6_c0000gn/T/heltai/spack-stage/spack-stage-dealii-9.2.0-747mf54t7nef2iff4la5ab3qxlrpd2jb/spack-src/source/fe/fe.cc> in function* * virtual Point<dim - 1> dealii::FiniteElement<3, 3>::unit_face_support_point(const unsigned int) const* *The violated condition was: * * unit_face_support_points.size() == this->dofs_per_face* *Additional information: * * You are trying to access the support points of a finite element that either has no support points at all, or for which the corresponding tables have not been implemented.* I thought that FE_FaceP has support points on boundaries. Does it mean that I can't set the Dirichlet boundary conditions in this way? Or did I misunderstand anything? Thank you. Best, Zhuoran -- 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/CAO_HR-hhT-jE1tVwTGkEUvL87PcjGb%2BCKZcQGuyTpO80tH%3DMtA%40mail.gmail.com.
