def initialization_VBGMM(n_components, points_data, points_test=None,
    covariance_type='full'):
    """
    This method returns an array of k means and an array of k covariances (dim,dim)
    which will be used in order to initialize an EM algorithm
    
    Parameters
    ----------
    points : an array (n_points,dim)
    
    k : int
        the number of clusters
        
    Returns
    -------
    means : an array (n_components,dim)
        The initial means computed
    
    cov : an array (n_components,dim,dim)
        The initial covariances computed
        
    log_weights : an array (n_components,)
        The initial weights (log) computed
        
    log_assignements : an array (n_points,n_components)
        The log of the soft assignements according to VBGMM
        
    """
    from .VBGMM import VariationalGaussianMixture
    GM = VariationalGaussianMixture(n_components)
    GM.fit(points_data, points_test, patience=0)
    log_assignements = GM.predict_log_resp(points_data)
    return GM.means, GM.cov, GM.log_weights, log_assignements
