Boundary labels

Hi all,

I realise we can loop over triangles, and print Th(i).label so that know which nodes are at boundaries. However, if we use of P2 element, and loop over degrees of freedom such as:
for(int i=0; i<Vh.ndof; i++){ … }
How can we know which degrees of freedom sit on boundaries?


varf vGamma(u, v) = on(boundaryLabel, u = 1.0);
real[int] onGamma(Vh.ndof);
onGamma = vGamma(0, Vh, tgv = -1);
for(int i = 0; i < Vh.ndof; ++i) if(abs(onGamma[i]-1.0) < 0.1) cout << "on boundaryLabel" << endl;

Great! Thank you very much.