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: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/63c0c8f2-2eab-20cc-bf0f-19b8ed32ec77%40colostate.edu.
