Hi all,
Instead of defining a 'first' and 'second' boundary with the help of two boundary_ids this function defines a 'left' boundary as all faces with local face index 2*dimension and boundary indicator b_id and, similarly, a 'right' boundary consisting of all face with local face index 2*dimension+1 and boundary indicator b_id. Faces with coordinates only differing in the direction component are identified.
I am confused as to what could 'dimension' mean here. Neither the template, nor the arguments have any variable called dimension. For dimension=3, this would mean faces with local ids 6 and 7 would be looked for. They don't exist right??
I was wondering if this is a typo and 'dimension' should actually be 'direction'. Any clarification would be greatly helpful.
Thanking in anticipation
Vachan