miscellaneous.hpp

Go to the documentation of this file.

// This file is part of GQCG-GQCP.
//
// Copyright (C) 2017-2020  the GQCG developers
//
// GQCG-GQCP is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// GQCG-GQCP is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with GQCG-GQCP.  If not, see <http://www.gnu.org/licenses/>.
 
#pragma once
 
 
#include <algorithm>
#include <cstdlib>
#include <fstream>
#include <functional>
#include <iterator>
#include <string>
#include <vector>
 
 
namespace GQCP {
 
 
template <typename T>
size_t findElementIndex(const std::vector<T>& vector, const T& value) {
 
    const auto& it = std::find(vector.begin(), vector.end(), value);  // 'it' for iterator
 
    // Check if the value was found
    if (it == vector.end()) {
        throw std::out_of_range("findElementIndex(const std::vector<T>&, const T&): the given value was not found in the given vector");
    }
 
    return std::distance(vector.begin(), it);  // the 'difference' between two iterators is an index
}
 
 
std::vector<std::vector<size_t>> generatePartitionsOf(const size_t n, const size_t k);
 
std::vector<std::vector<size_t>> generateUniquePartitionsOf(const size_t n, const size_t k);
 
size_t grayCodeOf(const size_t S);
 
size_t matrixIndexMajor(const size_t v, const size_t cols, const size_t skipped = 0);
 
size_t matrixIndexMinor(const size_t v, const size_t cols, const size_t skipped = 0);
 
void printExecutionTime(const std::string& method_name, const std::function<void()>& function);
 
size_t strictTriangularRootOf(const size_t x);
 
size_t triangularRootOf(const size_t x);
 
std::ifstream validateAndOpen(const std::string& filename, const std::string& extension);
 
size_t vectorIndex(const size_t i, const size_t j, const size_t cols, size_t skipped = 0);
 
 
}  // namespace GQCP