Using STEP file and OpenCASCADE::NormalToMeshProjectionManifold gives error

[Amit Singh]

Dear all

 

I am trying to use OpenCASCADE::NormalToMeshProjectionManifold for a STEP file and an initial mesh for an ellipsoid x^2 + 16y^2 + 64z^2 = 1.

I make the step file in GMSH and I use an appropriately scaled unrefined hyper_sphere as the initial mesh. The error I am getting is

An error occurred in line <241> of file <./source/opencascade/manifold_lib.cc> in function
    dealii::Point<dim> dealii::OpenCASCADE::{anonymous}::internal_project_to_manifold(const TopoDS_Shape&, double, const dealii::ArrayView<const dealii::Point<dim> >&, const dealii::Point<
dim>&) [with int spacedim = 3]
The violated condition was:
    closest_point(sh, surrounding_points[i], tolerance) .distance(surrounding_points[i]) < std::max(tolerance * surrounding_points[i].norm(), tolerance)
Additional information:
    The point [ -0.57735 -0.144338 -0.0721688 ] is not on the manifold.

If I substitute the above point coordinates in the equation of my surface, I get

x^2 + 16y^2 + 64z^2 = 1.00000183878416

So the point should be on the manifold upto tolerance. I have tried reducing the tolerance when reading OpenCASCADE::get_shape_tolerance() but it has not helped.

Please help me debug this.

Regards,

Amit

The information contained in this electronic communication is intended solely for the individual(s) or entity to which it is addressed. It may contain proprietary, confidential and/or legally privileged information. Any review, retransmission, dissemination, printing, copying or other use of, or taking any action in reliance on the contents of this information by person(s) or entities other than the intended recipient is strictly prohibited and may be unlawful. If you have received this communication in error, please notify us by responding to this email or telephone and immediately and permanently delete all copies of this message and any attachments from your system(s). The contents of this message do not necessarily represent the views or policies of BITS Pilani.