Hi Zhenlin,

I concur with Denis - I regularly use the method that he's outlined in 
conjunction with hp::DoFHandler.

Regards,
Jean-Paul

On Wednesday, February 22, 2017 at 9:53:26 AM UTC+1, Denis Davydov wrote:
>
> 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.

Reply via email to