/*******************************************************
 * Copyright (c) 2014, ArrayFire
 * All rights reserved.
 *
 * This file is distributed under 3-clause BSD license.
 * The complete license agreement can be obtained at:
 * http://arrayfire.com/licenses/BSD-3-Clause
 ********************************************************/

#include <af/array.h>
#include <af/statistics.h>
#include "error.hpp"

namespace af {

array cov(const array& X, const array& Y, const bool isbiased) {
    const af_var_bias bias =
        (isbiased ? AF_VARIANCE_SAMPLE : AF_VARIANCE_POPULATION);
    return cov(X, Y, bias);
}

array cov(const array& X, const array& Y, const af_var_bias bias) {
    af_array temp = 0;
    AF_THROW(af_cov_v2(&temp, X.get(), Y.get(), bias));
    return array(temp);
}

}  // namespace af