-
Notifications
You must be signed in to change notification settings - Fork 16
Standardize trace_particle_through_mesh #143
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: master
Are you sure you want to change the base?
Conversation
- passing all tests except smoke_test_particle It is giving a warning for MPI oversubscribing and "Structure not initalized at Particle"
- passing all tests
instead of checking if tol was passed using negative values, it is made optional. In the future, tol will be calculated outside according to the standardization plan and for this reason, we don't need it to be passed as reference or pointer
it is passed as a parameter rather than computing it every time the search is called.
it's an initial step to make the function standard with proper docuemntation. [important] there are some breaking changes and the function is under heavy development
the functor decides what happens to the particles at the boundary
next_element parameter holds the next element when a particle intersects a face. This is added to both the trace_particle_through_mesh and the functor. The functor may or may not use the value, but there's a high chance to be used. [Important] The test is failing. This is because although the particles doesn't intersect face 0 (it is along the way of the ray going forward) the closeness algorithm finds it as the intersected face.
|
@jacobmerson I will keep it draft until I bring the test case here. But please start reviewing the changes. I will commit the new test case for |
- works for particles going to a new element - or remains in the same - [to be fixed] particles leaking out
next_element becomes -1 when particle leaks
fix search_mesh in accordance with changes in trace_particle_through_mesh [future fixes] - interPoints should give the face intersection location when leaking out instead of the destination
jacobmerson
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I put some relevant comments here: #144
Mostly just a few naming things.
ray_intersects_trianglewithline_segment_intersects_triangle.set_new_elementfunction inside the functor.next_elementparameter that holds the next element and passes to the functor.find_exit_facesearch_mesh*After adding the callback functor, it is necessary for the user to know a clearly defined meaning of the particle state arrays.