pumi-pic
pumi-pic copied to clipboard
SCOREC
support libraries for unstructured mesh particle in cell simulations on GPUs and CPUs
As a continuation of standardizing the search method (next steps of #143), this wraps the search related particle state arrays inside this class. It includes the following tasks: - [x]...
- [x] Replace `ray_intersects_triangle` with `line_segment_intersects_triangle`. - [x] Standardize* the behavior of arrays used as parameters. - [x] Move `set_new_element` function inside the functor. - [x] Add `next_element` parameter that...
`elmArea` is computed using the function `measure_elements_real`. Users may want to cache the information and avoid computing it in every search. Note: `tol` is related to `elmArea` and I have...
Minor change to `search_mesh_3d` function definition to be consistent with `search_mesh_2d` function. `search_mesh_2d` function is defined as: https://github.com/SCOREC/pumi-pic/blob/a9f7f8ec2630d66e68178c9ffe51610f2de2c15e/src/pumipic_adjacency.hpp#L1010-L1020
In the following segment of the [`trace_particle_through_mesh` function](https://github.com/SCOREC/pumi-pic/blob/a9f7f8ec2630d66e68178c9ffe51610f2de2c15e/src/pumipic_adjacency.tpp#L489-L491), the tolerance for search algorithms is computed based on the "minimum area" and it can be expensive. ```c++ const auto elmArea =...
Currently, the following procedures are used to load partitioned mesh in both `XGCm` and `GITRm` following `PUMIPic` example: - read in the full mesh - load in the partition file...
Documenting this here for future reference: - Through some tests, I was able to perform mesh adaptation for both 2D and 3D meshes in `PUMIPic` when the mesh is not...
To output the particle coordinate information in `XGCm`, particle data structure is accessed on the host by the following code. However, the code crashed with the latest PUMIPic commit [b6678b0](https://github.com/SCOREC/pumi-pic/commit/b6678b0a0b8c9ad1e143831bdc2920b944b0f5ff):...
We discovered that the `find_exit_face` function in `pumipic_adjacency.tpp` should not work when the `useBcc` flag is `true`. A test case is added here to demonstrate the failed scenario. Briefly, [the...