Tests and PETSc block-preconditioner integration for 3D epsilon Stokes operator, P1-P1-stab and P2-P1