327     Vec3d npa, npb;  in decompose()  local
363             npa[0] = S(0, 0),               npb[0] = S(0, 0);  in decompose()
364             npa[1] = S(0, 1) + rtM22,       npb[1] = S(0, 1) - rtM22;  in decompose()
365             npa[2] = S(0, 2) + e12 * rtM11, npb[2] = S(0, 2) - e12 * rtM11;  in decompose()
368             npa[0] = S(0, 1) + rtM22,       npb[0] = S(0, 1) - rtM22;  in decompose()
369             npa[1] = S(1, 1),               npb[1] = S(1, 1);  in decompose()
370             npa[2] = S(1, 2) - e02 * rtM00, npb[2] = S(1, 2) + e02 * rtM00;  in decompose()
373             npa[0] = S(0, 2) + e01 * rtM11, npb[0] = S(0, 2) - e01 * rtM11;  in decompose()
374             npa[1] = S(1, 2) + rtM00,       npb[1] = S(1, 2) - rtM00;  in decompose()
375             npa[2] = S(2, 2),               npb[2] = S(2, 2);  in decompose()
392     Vec3d nb = npb / norm(npb);  in decompose()