... Sorry, there is a mistake in one of the last code lines. It should be:
// append cell if it is cut by the interface
if (surface_vertices.size() > 0)
cut_cells.emplace_back(cell);
Magdalena schrieb am Freitag, 16. Dezember 2022 um 11:09:40 UTC+1:
> Dear Praveen,
>
> if you have an
Dear Praveen,
if you have an implicit description of your curve, let's say a level-set
like DoF-vector,
VectorType level_set_vector;
where a value "contour_value" represents your surface, a helpful tool
could be the Marching Cube Algorithm (
Hi,
I think Step 85 does this by using a (global) level set function.
Otherwise, you can use RTrees of bounding boxes
with
https://www.dealii.org/developer/doxygen/deal.II/classGridTools_1_1Cache.html#aca2782d6e93b5a0033c046b57904c67f
combined with boost::geometry::index to identify