Skip to content

Commit

Permalink
reduced the tolerance for checking whether a point lies internal or e…
Browse files Browse the repository at this point in the history
…xternal to a surface, by a factor of 10. fixes an issue of seed points being created outside the surface. previous tolerance allowed this when the vector from internal to seed point was very close to surface element boundaries
  • Loading branch information
Merryn Tawhai authored and Merryn Tawhai committed Mar 21, 2024
1 parent 1197dc9 commit 78d8306
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/lib/mesh_utilities.f90
Original file line number Diff line number Diff line change
Expand Up @@ -875,7 +875,7 @@ function point_internal_to_surface(num_vertices,triangles,point_xyz,vertex_xyz)
integer :: i,ncrossed,ntri,num_triangles
real(dp) :: area,area_triangle,cofm_surfaces(3),denominator,&
norm_v(3),point(3),P1(3),P2(3),P3(3),u
real(dp),parameter :: dist_tol = 1.0e-4_dp, user_tol = 1.0e-14_dp
real(dp),parameter :: dist_tol = 1.0e-5_dp, user_tol = 1.0e-14_dp
logical :: cross_any

num_triangles = count(triangles(:,:).ne.0)/3.0_dp
Expand Down

0 comments on commit 78d8306

Please sign in to comment.