I'm using Comsol script to solve coefficient form PDE. I assign boundary conditions using e.g. fem.bnd.h={1 0 0 0} and fem.bnd.r={1 0 0 0}. However, depending on relative positions of various structures, the Comsol arbitrarily changes the numbering order of the boundaries in my geometry, which...
