I think the error is like Marc said. And if I call refine_globle() before set_manifold(), the error will go away.
Best, LU 在2023年3月23日星期四 UTC+8 11:48:50<Wolfgang Bangerth> 写道： > On 3/22/23 01:47, RUIJIAN LU wrote: > > > > I am a beginner to dealii, now I want to solve an equation on a > cylindrical > > coordinate system. I create a triangulation by cylinder() and use > > CylindricalManifold for it, but it has some error. How can I do with it? > > > > error: > > Number dealii::Physics::VectorRelations::signed_angle(const > dealii::Tensor<1, > > dim, Number>&, const dealii::Tensor<1, dim, Number>&, const > dealii::Tensor<1, > > dim, Number>&) [with int spacedim = 3; Number = double] > > The violated condition was: > > std::abs(axis * a) < 1.e-12 * a.norm() && std::abs(axis * b) < > 1.e-12 * > > b.norm() > > Additional information: > > The vectors are not perpendicular to the axial vector. > > Can you show what the backtrace is that gets you to this place? It would > be > important to understand what these vectors a, b, and axis are. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > > -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/73e89a64-2545-463a-8cd6-258a70d757f2n%40googlegroups.com.
