def square_matrix_scipy(matrix):
    """
    barebone scipy version of matmul

    :param matrix: scipy sparse coo matrix
    :return: matrix squared, also in scipy sparse coo format
    """
    matrix_csr = matrix.tocsr()
    return (matrix_csr * matrix_csr).tocoo()