Dear all,
I have code that runs fine with periodic boundary condition for
a GridGenerator::parallelepiped. Also when I change the mesh to
a GridGenerator::subdivided_parallelepiped and has only 1 as
n_subdivisions, it also works. But when I change it to 2, I get the
following error:
*An error occurred in line <2686> of file
</home/joel/Downloads/dealii-8.2.1/source/grid/grid_tools.cc> in function*
* void dealii::GridTools::collect_periodic_faces(const CONTAINER&,
dealii::types::boundary_id, dealii::types::boundary_id, int,
std::vector<dealii::GridTools::PeriodicFacePair<typename
CONTAINER::cell_iterator> >&, const dealii::Tensor<1, CONTAINER::
space_dimension>&, const dealii::FullMatrix<double>&, const
std::vector<unsigned int>&) [with CONTAINER = dealii::DoFHandler<3>;
dealii::types::boundary_id = unsigned char; typename
CONTAINER::cell_iterator =
dealii::TriaIterator<dealii::DoFCellAccessor<dealii::DoFHandler<3>, false>
>]*
*The violated condition was: *
* pairs1.size() == pairs2.size()*
*The name and call sequence of the exception was:*
* ExcMessage ("Unmatched faces on periodic boundaries")*
*Additional Information: *
*Unmatched faces on periodic boundaries*
How can I fix it?
Best,
Joel
On Friday, August 12, 2016 at 9:35:23 AM UTC+2, Joel Davidsson wrote:
>
> Ok
>
> Thanks,
>
> Joel
>
> On Thursday, August 11, 2016 at 6:05:41 PM UTC+2, Bruno Turcksin wrote:
>>
>> Joel,
>>
>> On Thursday, August 11, 2016 at 11:29:35 AM UTC-4, Joel Davidsson wrote:
>>>
>>>
>>> I just had a question about the third option. If one used FE_Nothing,
>>> what is the boundary that is used when doing calculations? Is it the
>>> boundary of the big box or the green box? Is this handle automatically or
>>> does one need to mark the green box with boundary indicator for example.
>>>
>> It will be the boundary on the big box and you can't use boundary
>> indicator because you are not on the boundary. The boundary is a geometric
>> "quantity", it doesn't care about the type of finite elements you are
>> using. So if you want to use Dirichlet "boundary" condition on the green
>> box, you need to build the ConstraintMatrix yourself.
>>
>> Best,
>>
>> Bruno
>>
>
--
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.
