I believe your statement of the problem is correct. I'd like to
take some collection of FixedValue or a FixedFlux objects that were
used to build a set of boundary conditions and determine what faces
they are associated with (getting the values themselves is
straightforward). What I imagine eventually saving is six arrays.
All will be the same length as the array of x-face coordinates.
1. The x-coordinates of all the faces (saved anyway).
2. The y-coordinates of all the faces (saved anyway).
3. An array that is a mask of booleans showing where Dirichlet
boundary conditions are applied.
4. An array of the values of Dirichlet boundary conditions.
5. An array that is a mask of booleans showing where Neuman
boundary conditions are applied.
6. An array of the values of Neuman boundary conditions.
There's probably a better way to do this, but I want something
complete enough that anyone can pick up and solve the same problem
using whatever tools they like to check my answers.
This information would be nice for developing spatially varying
boundary conditions... Just a thought.
