Hi everyone,
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.
code:
const Tensor<1, 3> axis({1.0, 0.0, 0.0});
const Point<3> axial_point(0.0, 0.0, 0.0);
const CylindricalManifold<3> cylinder(axis, axial_point);
GridGenerator::cylinder(triangulation, 2.5, 2.5);
for (auto cell : triangulation.active_cell_iterators())
cell->set_all_manifold_ids(2);
triangulation.set_manifold(2, cylinder);
--
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/696e28f8-775d-438d-a066-86fc5abd78aen%40googlegroups.com.
