Inverse design is currently applied to obtain unusual acoustical devices based on ordered and non-ordered scatterers. Recently, flat and lenticular acoustic lenses, de-multiplexors, directional sound sources and acoustic cloaks have been designed using a variety of optimization methods like genetic algorithms, simulated annealing and shape or topology optimizations among others. For example, feasible one-directional cloaks were first designed in two- and three-dimensions using high symmetry objects like cylinders [Appl. Phys. Lett. 99, 074102 (2011)] and toroidal scatterers [Phys. Rev. Lett. 110, 124301 (2013)], respectively. Recently, an extraordinary simple one-directional acoustic cloak has been reported using a technique that combines the method of fundamental solutions with arbitrary shape scatterers [Sci. Rep. 8, 12924 (2018)]. Here, we report a further improvement of the method by combining the Boundary Element Method (BEM) with shape optimization to obtain quasi-omnidirectional cloaks in two-dimensions. The shapes of the scatterers are optimized for the cloaking of the whole setup with waves impinging on the stealth object (a cylinder) from many different directions. An important feature of the method is the possibility of including visco-thermal acoustic losses in the optimization process [J. Sound Vib. 447, 120–136 (2019)].