fix substantive error: need |det(J)| in all integrals, and note det(J)^2 = |det(J)|^2
record *all* boundary segments into .poly file