The increasing interest among the scientific community of reflectarray antennas in recent years leads to develop effective ad-hoc procedures for the antenna performances optimization, often based on evolutionary iterative algorithms. However, in most of engineering problems such an optimization approach could be very computational expensive, since it could require the computation of a complex fitness function thousands of times. An example of this is for instance the design of a large array or reflectarray, in which the structure to be optimized presents a lot of degrees of freedom and all concur to the performances of the whole antenna. To enhance the speed of the optimization task, in this work the optimization of a broadband reflectarray...