This article presents a systematic procedure to shape the body of an electrically large object in order to achieve a desired radiation/scattering pattern. In particular, the object surface is described by non-uniform rational B-splines (NURBS) and physical optics (PO) is used for electromagnetic field calculations. The NURBS control points and weights are considered as degrees of freedom for an evolutionary optimization algorithm. They are iteratively optimized according to the specific constraints of radiation/scattering problem under consideration. To speedup the optimization process, the electromagnetic codes are implemented on graphics processing unit (GPU). The accuracy of the simulation results are verified by CST commercial electromagnetic software. Furthermore, the results show simulation speedup on the order of several times to a few tens of times. This makes an iterative optimization of a surface rather practical on a single personal computer.