you can always partition regular triangulation by
GridTools::partition_triangulation (n_mpi_processes, triangulation);
and adjust your assembly loops to
if (cell->subdomain_id == this_mpi_process)
to parallelize with MPI. The downside is that all your processes will now
own the **complete mesh**.
There are also some other minor downsides like serial SolutionTransfer and
not straight-forward serialization/restart.
Regards,
Denis.
On Wednesday, February 22, 2017 at 6:51:34 AM UTC+1, zhen-lin Wang wrote:
>
> Hi all,
>
> I found that hp::DoFHandler does not support
> parallel::distributed::Triangulation currently. The error shows like:
> -----------
>
> An error occurred in line <1669> of file
> </Applications/deal.II.brew.app/Contents/Resources/brew/src/dealii/8.4.2/source/hp/dof_handler.cc>
>
> in function
>
> dealii::hp::DoFHandler<2, 2>::DoFHandler(const Triangulation<dim,
> spacedim> &) [dim = 2, spacedim = 2]
>
> The violated condition was:
>
> (dynamic_cast<const parallel::distributed::Triangulation< dim,
> spacedim >*> (&tria) == 0)
>
> The name and call sequence of the exception was:
> ExcMessage ("The given triangulation is parallel distributed but "
> "this class does not currently support this.")
> -----------
> I wondThe problem I am doing needs Fe_Nothing and I want to parallelize
> the code with BlockSparseMatrix and BlockVector. I follow the step55
> which using parallel::distributed::Triangulation. I wonder can I use
> regular Triangular and hp::DoFHandler but still follow step55 for the
> rest of the part?
>
> Thanks,
>
> Zhenlin Wang
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
