The design of intrinsically flat two-dimensional optical components, i.e., metasurfaces, generally requires an extensive parameter search to target the appropriate scattering properties of their constituting building blocks. Such design methodologies neglect important near-field interaction effects, playing an essential role in limiting the device performance. Optimization of transmission, phase-addressing and broadband performances of metasurfaces require new numerical tools. Additionally, uncertainties and systematic fabrication errors should be analysed. These estimations, of critical importance in the case of large production of metaoptics components, are useful to further project their deployment in industrial applications. Here, we report on a computational methodology to optimize metasurface designs. We complement this computational methodology by quantifying the impact of fabrication uncertainties on the experimentally characterized components. This analysis provides general perspectives on the overall metaoptics performances, giving an idea of the expected average behavior of a large number of devices.