prime_numbers.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Pooyan Dadvand
//                   
//
 
#if !defined(KRATOS_PRIME_NUMBERS_H_INCLUDED )
#define  KRATOS_PRIME_NUMBERS_H_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream>
#include <array>
 
// Project includes
#include "includes/exception.h"
 
namespace Kratos
{
 
 
 
  class PrimeNumbers
    {
 
        static constexpr std::size_t mNumberOfPreCalculatedPrimes = 100000;
 
        static constexpr std::size_t mLargestPreCalculatedPrime = 1299689;
 
    public:
 
 
        PrimeNumbers() {}
 
      PrimeNumbers(PrimeNumbers const& rOther) = delete;
      
      virtual ~PrimeNumbers(){}
 
 
      
      PrimeNumbers& operator=(PrimeNumbers const& rOther) = delete;
 
      std::size_t operator[](std::size_t Index) {
          return GetPreCalculatedPrime(Index);
      }
 
 
 
 
      static std::size_t GetNumberOfPreCalculatedPrimes() {
          return mNumberOfPreCalculatedPrimes;
      }
      
      static std::size_t LargestPreCalculatedPrime() {
          return mLargestPreCalculatedPrime;
      }
 
      static std::size_t GetPreCalculatedPrime(std::size_t Index) {
          KRATOS_DEBUG_ERROR_IF(Index >= mNumberOfPreCalculatedPrimes) << "Index " << Index 
              << " is larger than Number of precalculated primes " << mNumberOfPreCalculatedPrimes << std::endl;
          return mPrecalculatedPrimes[Index];
      }
 
 
      static bool IsPrime(std::size_t TheNumber);
 
 
 
 
 
    private:
 
        static const std::array<std::size_t, mNumberOfPreCalculatedPrimes> mPrecalculatedPrimes;
 
 
        static std::size_t FindGreaterEqualPrecalculatedPrimeNumber(std::size_t TheNumber);
 
 
    }; // Class PrimeNumbers
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_PRIME_NUMBERS_H_INCLUDED  defined