def hadamard_complex(x_re, x_im, y_re, y_im):
    """Hadamard product for complex vectors"""
    result_re = x_re * y_re - x_im * y_im
    result_im = x_re * y_im + x_im * y_re
    return result_re, result_im