i checked again in menus to generate hex-dominant mesh, it seems some option is missing in Mesh.Recombine3DAll to be activated (yes/1) instead of default (no/0). Another meshing algorithm setting also in 2D R-Tree options.
result without recombine 3D in latest Gmsh (v4.13.1)
result without recombine 3D in PrePoMax almost identical in pattern and division,
