I use AWR (Axiem, Anaylst), CST, and Sonnet. I'm satisfied with accuracy; need more speed. It's at a point where material and manufacturing tolerances are more of an issue, so I need to run parameter sweeps, yield analysis, etc. So multi-threading, GPU acceleration, and distributed computing help that.
Working on a thick patch antenna today. Axiem (2.5D infinite substrate) said 1.57 GHz. Analyst says 1.51 GHz. It actually measured 1.53 GHz. That's a 3% error on an 8% 12 dB RL bandwidth antenna. Those simulated results are after mesh convergence, so finer meshing won't help.
So I use 2-3 solvers on a single problem. It's important that the geometry translation between them be quick and easy. That's where AWR comes into play; everything plugs into it. 3D solver (modeling) is horrible for pushing polygons during the design process. 2.5D solvers are awesome for this, but then you need to sometimes push the geometry to 3D solver.