[operators.cg] add support for Neumann boundaries in L2ProductFunctionalP1, L2ProductFunctionalQ1
[domaindescriptions] add more assertions