On Monday, November 07, 2011 10:02:40 AM Margarita Satraki wrote:
> Hi,
>
> Can you please advise me on how I can change the following code so that the
> volume_mesh (like the triangulation) also follows the curved boundary?
>
> Thanks,
>
> Margarita
>
> Triangulation<spacedim> volume_mesh;
> GridGenerator::half_hyper_ball(volume_mesh);
> static HyperBallBoundary<dim,spacedim> surface_description;
> triangulation.set_boundary (0, surface_description);
> { std::set<unsigned char> boundary_ids;
> boundary_ids.insert (0);
> GridTools::extract_boundary_mesh (volume_mesh,
> triangulation,boundary_ids);
> }
> triangulation.refine_global(3);
> volume_mesh.refine_global(3);
Obvioulsy you need to set a boundary description also for volume_mesh. step-1
would be the first place where this would be described...
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
