def METHOD_NAME(myadc, Lvv, p, chnk_size):
    ''' Returns approximate antisymmetrized vvvv integrals (alpha/beta spin) used in unrestricted implementation'''
    naux = myadc.with_df.get_naoaux()
    nvir = Lvv.shape[1]
    ind_vv_g = np.tril_indices(nvir, k=-1)
    if chnk_size < nvir:
        Lvv_temp = np.ascontiguousarray(Lvv.T[p:p+chnk_size].reshape(-1,naux))
    else:
        Lvv_temp = np.ascontiguousarray(Lvv.T.reshape(-1,naux))
    Lvv = Lvv.reshape(naux,nvir*nvir)
    vvvv = np.dot(Lvv_temp, Lvv)
    vvvv = vvvv.reshape(-1, nvir, nvir, nvir)
    vvvv = vvvv.transpose(0,2,1,3).reshape(-1, nvir, nvir, nvir)
    vvvv -= vvvv.transpose(0,1,3,2)
    vvvv = np.ascontiguousarray(vvvv[:, :, ind_vv_g[0], ind_vv_g[1]])
    return vvvv