Thank you so much.
Now it is working fine. :-)
Thanks & Regards,
Deepika
On Tue, Feb 21, 2023 at 8:42 PM Daniel Arndt wrote:
> Deepika,
>
>
> https://www.dealii.org/current/doxygen/deal.II/namespaceGridGenerator.html#a533c4778cbc9bcbed365dcab42ca4418
> says
>
> By default, the manifold_id
Deepika,
https://www.dealii.org/current/doxygen/deal.II/namespaceGridGenerator.html#a533c4778cbc9bcbed365dcab42ca4418
says
By default, the manifold_id is set to 0 on the boundary faces, 1 on the
boundary cells, and numbers::flat_manifold_id
Thank you very much sir for your explanation.
I have assigned material id to different domains. Now I am trying to set
spherical manifold id by using the following code
for (const auto : tri.cell_iterators())
{
double v = (cell->center()[0])*(cell->center()[1]) ;
if
At the moment, you are setting the manifold ID on the boundary of your *domain*
with the call to
tri.set_all_manifold_ids_on_bounary(4);
tri.set_manifold (4, manifold_neg);
And then you are telling the triangulation to use on it the
SphericalManifold<2> manifold_neg(center).
Then you re-do
You should set the manifold ids of the two regions to different numbers, and add the second manifold with the second id you choose, otherwise the second call is simply overwriting the first one, and they are both described by the last manifold. LucaIl giorno 20 feb 2023, alle ore 07:08, Deepika