Three-dimensional vector beam-propagation method for second harmonic generation analysis