Hej, I found the issue ... I forgot the line
constraints.distribute ( sol );
after the solution routine ... so it was just mainly a visualization issue.
Thanks for the quick response!
El jueves, 22 de septiembre de 2016, 19:26:28 (UTC+2), Daniel Arndt
escribió:
>
> Juan,
>
> [...]
>> 1) My first question is about how the "direction" works, is it 1
>> correspond to "x" and 2 corresponds to "y"?
>>
> No, direction 0/1/2 means that the support_points of the dofs to be
> identified just differ in the x-/y-/z-component.
> What you are doing should be perfectly fine.
>
>
>> 2) When I do this I obtain fairly reasonable solutions, however I am
>> getting a weird Dirichlet-type of boundary conditions on boundaries 1 and 3
>> (see the attached figures), that distort the fields.
>> I seem unable to find what is wrong, so I wonder if anyone has
>> experienced this kind of behavior.
>>
> No, I haven't observed this problem. Can you check in your
> ConstraintMatrix that you don't use any other constraints? For doing so use
> ConstraintMatrix.print(std::cout)
> and check if there is a DoF that is not constrained to another one.
>
> Best,
> Daniel
>
--
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.
