def cross_validate(model, X, y, k_folds=5, metric="auto", shuffle=True):
    """Cross Validation

    Evaluates the given model using the given data
    repetitively fitting and predicting on different
    chunks (folds) from the data.
    
    Parameters:
    -----------
    model : dojo-model, the model to be evaluated
    X : matrix, shape (n_samples, n_features), the data used for evaluation
    y : vector, shape (n_samples, ), the desired labels
    k_folds : integer, optional, the number of iterations/folds
    metric : the single value error/accuracy metric, optional
    shuffle : boolean, whether to shuffle the data before
    splitting it or not
    
    Returns:
    --------
    dict_scores : dictionary with train scores and test scores
    
    """

    train_scores = []
    test_scores = []
    folds = KFolds(X, y, k=k_folds, shuffle=shuffle)

    for X_train, X_test, y_train, y_test in folds:
        model.fit(X_train, y_train)

        if metric is None or metric == "auto":
            train_scores.append(model.evaluate(X_train, y_train))
            test_scores.append(model.evaluate(X_test, y_test))
        else:
            train_scores.append(
                metric(y_train, model.predict(X_train))
            )
            test_scores.append(
                metric(y_test, model.predict(X_test))
            )

    return {
        "train_scores": np.array(train_scores),
        "test_scores": np.array(test_scores),
    }