[Dune] Dirichlet conditions on curved boundaries

Andreas Dedner a.s.dedner at warwick.ac.uk
Tue Jun 2 19:03:22 CEST 2015


Hi.
You should have said that your using PDElab and used the
corresponding list. I was assuming since you didn't mention any
package that you were using the dune grid interface directly.
Over to PDElab developers.
Best
Andreas

On 02/06/15 13:22, conf86 at web.de wrote:
> I found the code, but I have the impression, I can't figure out how to use it with mine.
> I am currently using a modified version of example02_bcextension.hh from the course-examples in dune-pdelab-howto. Do I need to fully rewrite bcextension.hh or is there a possibility to integrate IntersectionIterator into it?
>
>> Gesendet: Dienstag, 02. Juni 2015 um 00:08 Uhr
>> Von: "Andreas Dedner" <a.s.dedner at warwick.ac.uk>
>> An: conf86 at web.de
>> Betreff: Re: [Dune] Dirichlet conditions on curved boundaries
>>
>> Well not quite - you need to know which dofs lie on that intersection
>> (an edge of the boundary in
>> 2d or a face in 3d) and then you set those dofs.
>> In the dune-grid-howto there is a finiteelement.cc example which should
>> give you an idea.
>> Andreas
>>
>> On 01/06/15 23:03, conf86 at web.de wrote:
>>> Thank you for your answer.
>>>
>>>> Never directly compare doubles - use something like abs(...-1,0) < eps.
>>> Ok, now it makes sense to me not to do this.
>>>
>>>> But you have to use the intersection iterator to find the boundary (use
>>>> it->boundary())
>>>> and then you can set your data.
>>> So basically, I just need to replace the condition?
>>>
>>>    if (it->boundary()) y=1.0;
>>>
>>> Sorry, if this might be a stupid question, but I couldn't find any sample code for this.
>>>
>>> Best regards.
>>>
>>> _______________________________________________
>>> Dune mailing list
>>> Dune at dune-project.org
>>> http://lists.dune-project.org/mailman/listinfo/dune
>>





More information about the Dune mailing list