The arrangement of binary subwavelength structures is a promising alternative to the conventional multiheight level technique to generate computer generated holograms (CGHs). However, the current heuristic design approach leads to a slight mismatch between the target signal and experimental data. To evaluate this deviation, a diffractive beam splitter design is investigated rigorously using a finite-difference time-domain (FDTD) method. Since the use of a rigorous Maxwell-equation solver like FDTD requires a massive computational effort, an alternative scalar approach, a fast Fourier transform beam propagation method (FFT-BPM), is investigated with a substantial higher computing speed, showing still a good agreement with the FDTD simulation and experimental data. Therefore, an implementation of this scalar approach into the CGH design process offers the possibility to significantly increase the accuracy.