Getting cppsim to compile from sources alongside the MindQuantum sources

build.bat

@@ -372,7 +372,6 @@ call %SCRIPTDIR%\dos\build_cmake_option.bat BUILD_TESTING !enable_tests!
 call %SCRIPTDIR%\dos\build_cmake_option.bat CLEAN_3RDPARTY_INSTALL_DIR !do_clean_3rdparty!
 call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_ANALYZER !enable_analyzer!
 call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_CMAKE_DEBUG !cmake_debug_mode!
-call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_CPPSIM !enable_cppsim!
 call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_CUDA !enable_gpu!
 call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_CXX_EXPERIMENTAL !enable_cxx!
 call %SCRIPTDIR%\dos\build_cmake_option.bat ENABLE_GITEE !enable_gitee!
@@ -504,7 +503,6 @@ rem ============================================================================
   echo   /CleanBuildDir      Delete build directory before building
   echo   /CleanCache         Re-run CMake with a clean CMake cache
   echo   /CleanVenv          Delete Python virtualenv before building
-  echo   /CppSim             (experimental) Enable the use of cppsim to generate simulation kernels
   echo   /Cxx                (experimental) Enable MindQuantum C++ support
   echo   /Debug              Build in debug mode
   echo   /Delocate           Delocate the binary wheels after build is finished

build.ps1

@@ -24,7 +24,6 @@ Param(
     [switch]$CleanCache,
     [switch]$CleanVenv,
     [ValidateNotNullOrEmpty()][string]$Config,
-    [switch]$CppSim,
     [ValidateNotNullOrEmpty()][string]$CudaArch,
     [switch]$Cxx,
     [switch]$Debug,
@@ -488,9 +487,6 @@ Do not use the CMake registry to find packages
 .PARAMETER Config
 Path to INI configuration file with default values for the parameters
 
-.PARAMETER CppSim
-(experimental) Enable the use of cppsim to generate simulation kernels
-
 .PARAMETER Cxx
 (experimental) Enable MindQuantum C++ support
 

build_locally.bat

@@ -132,11 +132,6 @@ rem ============================================================================
     shift & shift & goto :initial
   )
 
-  if /I "%1" == "/CppSim" (
-    set enable_cppsim=1
-    shift & goto :initial
-  )
-
   if /I "%1" == "/Cxx" (
     set enable_cxx=1
     shift & goto :initial
@@ -348,7 +343,6 @@ call %SCRIPTDIR%\dos\build_locally_cmake_option.bat BUILD_TESTING !enable_tests!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat CLEAN_3RDPARTY_INSTALL_DIR !do_clean_3rdparty!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_ANALYZER !enable_analyzer!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_CMAKE_DEBUG !cmake_debug_mode!
-call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_CPPSIM !enable_cppsim!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_CUDA !enable_gpu!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_CXX_EXPERIMENTAL !enable_cxx!
 call %SCRIPTDIR%\dos\build_locally_cmake_option.bat ENABLE_DOCUMENTATION !do_docs!
@@ -546,7 +540,6 @@ exit /B 0
   echo   /CleanCache         Re-run CMake with a clean CMake cache
   echo   /CleanVenv          Delete Python virtualenv before building
   echo   /ConfigureOnly      Stop after the CMake configure and generation steps (ie. before building MindQuantum)
-  echo   /CppSim             (experimental) Enable the use of cppsim to generate simulation kernels
   echo   /Cxx                (experimental) Enable MindQuantum C++ support
   echo   /Debug              Build in debug mode
   echo   /DebugCMake         Enable debugging mode for CMake configuration step

build_locally.ps1

@@ -27,7 +27,6 @@ Param(
     [ValidateNotNullOrEmpty()][string]$Config,
     [Alias("C")][switch]$Configure,
     [switch]$ConfigureOnly,
-    [switch]$CppSim,
     [ValidateNotNullOrEmpty()][string]$CudaArch,
     [switch]$Cxx,
     [switch]$Debug,
@@ -244,7 +243,6 @@ $cmake_args = @('-DIN_PLACE_BUILD:BOOL=ON'
                 "-DENABLE_ANALYZER:BOOL={0}" -f $CMAKE_BOOL[$enable_analyzer]
                 "-DENABLE_PROJECTQ:BOOL={0}" -f $CMAKE_BOOL[$enable_projectq]
                 "-DENABLE_CMAKE_DEBUG:BOOL={0}" -f $CMAKE_BOOL[$cmake_debug_mode]
-                "-DENABLE_CPPSIM:BOOL={0}" -f $CMAKE_BOOL[$enable_cppsim]
                 "-DENABLE_CUDA:BOOL={0}" -f $CMAKE_BOOL[$enable_gpu]
                 "-DENABLE_CXX_EXPERIMENTAL:BOOL={0}" -f $CMAKE_BOOL[$enable_cxx]
                 "-DENABLE_DOCUMENTATION:BOOL={0}" -f $CMAKE_BOOL[$do_docs]
@@ -439,9 +437,6 @@ Path to INI configuration file with default values for the parameters
 .PARAMETER ConfigureOnly
 Stop after the CMake configure and generation steps (ie. before building MindQuantum)
 
-.PARAMETER CppSim
-(experimental) Enable the use of cppsim to generate simulation kernels
-
 .PARAMETER Cxx
 (experimental) Enable MindQuantum C++ support
 

build_locally.sh

@@ -169,7 +169,6 @@ cmake_args=(-DIN_PLACE_BUILD:BOOL=ON
             -DENABLE_PROJECTQ:BOOL="${CMAKE_BOOL[$enable_projectq]}"
             -DENABLE_CMAKE_DEBUG:BOOL="${CMAKE_BOOL[$cmake_debug_mode]}"
             -DENABLE_CUDA:BOOL="${CMAKE_BOOL[$enable_gpu]}"
-            -DENABLE_CPPSIM:BOOL="${CMAKE_BOOL[$enable_cppsim]}"
             -DENABLE_CXX_EXPERIMENTAL:BOOL="${CMAKE_BOOL[$enable_cxx]}"
             -DENABLE_DOCUMENTATION:BOOL="${CMAKE_BOOL[$do_docs]}"
             -DENABLE_GITEE:BOOL="${CMAKE_BOOL[$enable_gitee]}"

ccsrc/include/cppsim/combinations.h

@@ -0,0 +1,185 @@
+#ifndef COMBINATIONS_H
+#define COMBINATIONS_H
+
+#include "gpu_support.h"
+
+#include <algorithm>
+#include <array>
+#include <cstdio>
+#include <stdint.h>
+#include <type_traits>
+
+// Iterate through the combinations using currying approach: https://stackoverflow.com/a/54508163/4063520
+class Combinations
+{
+// *****************************************************************************
+// Simple case: no sum constraint
+// *****************************************************************************
+private :
+
+	template<
+		uint32_t n0, uint32_t n1, uint32_t ...n, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void _iterate(Callable&& c)
+	{
+		for (uint32_t i = 0; i < n0; i += 2 * n1)
+		{
+			auto bind_an_argument = [i, &c](auto... args)
+			{
+				c(i, args...);
+			};
+
+			_iterate<n1, n...>(bind_an_argument);
+		}
+	}
+
+	template<
+		uint32_t n0, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void _iterate(Callable&& c)
+	{
+		for (uint32_t i = 0; i < n0; i++)
+		{
+			c(i);
+		}
+	}
+
+
+public :
+
+	template<uint32_t ...n>
+	struct Combination
+	{
+		using type = std::array<uint32_t, sizeof...(n)>;
+	};
+
+	// Iterate through all combinations.
+	// For each combination, call a user-provided function.
+	template<
+		uint32_t ...n, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void iterate(Callable&& c)
+	{
+		_iterate<n...>(c);
+	}
+
+	// Tell the length of a combination supplied by the iterator
+	// configured with the given set of template parameters.
+	template<uint32_t ...n>
+	static constexpr uint32_t length()
+	{
+		return sizeof...(n);
+	}
+
+	// Tell the number of combinations supplied by the iterator
+	// configured with the given set of template parameters.
+	template<uint32_t n0>
+	static uint32_t _popcount()
+	{
+		return n0;
+	}
+
+	// Tell the number of combinations supplied by the iterator
+	// configured with the given set of template parameters.
+	template<uint32_t n0, uint32_t n1, uint32_t ...n>
+	static uint32_t _popcount()
+	{
+		return (n0 / (2 * n1)) * _popcount<n1, n...>();
+	}
+
+	// Tell the number of combinations supplied by the iterator
+	// configured with the given set of template parameters.
+	template<uint32_t ...n>
+	static uint32_t popcount()
+	{
+		if (sizeof...(n) == 0) return 0;
+		return _popcount<n...>();
+	}
+
+	// Reverse the order of elements in a combination
+	// configured with the given set of template parameters.	
+	template<uint32_t ...n>
+	static void reverse(typename Combination<n...>::type& c)
+	{
+		std::reverse(c.begin(), c.end());
+	}
+
+// *****************************************************************************
+// Simple case: no sum constraint with a user-defined range
+// *****************************************************************************
+public :
+
+	template<
+		uint32_t n0, uint32_t n1, uint32_t ...n, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void _iterate(uint32_t* start, uint32_t& limit, Callable&& c)
+	{
+		for (uint32_t i = *start; (i < n0) && limit; i += 2 * n1)
+		{
+			// Flush starting point to zero, in order for all subsequent iterations
+			// to start from zero as usual.
+			*start = 0;
+
+			auto bind_an_argument = [i, &c](auto... args)
+			{
+				c(i, args...);
+			};
+
+			_iterate<n1, n...>(start++, limit, bind_an_argument);
+		}
+	}
+
+	template<
+		uint32_t n0, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void _iterate(uint32_t* start, uint32_t& limit, Callable&& c)
+	{
+		for (uint32_t i = *start; (i < n0) && limit; i++)
+		{
+			// Flush starting point to zero, in order for all subsequent iterations
+			// to start from zero as usual.
+			*start = 0;
+
+			c(i);
+			limit--;
+		}
+	}
+
+public :
+
+	// Iterate through combinations with specific starting point and duration.
+	// For each combination, call a user-provided function.
+	template<
+		uint32_t ...n, // max allowed sequence element value
+		class Callable
+	>
+	GPU_SUPPORT
+	static constexpr void iterate(const typename Combination<n...>::type& start_, const uint32_t limit_, Callable&& c)
+	{
+		auto start = start_;
+		uint32_t limit = limit_;
+		_iterate<n...>(start.data(), limit, c);
+	}
+
+	// Tell the number of combinations supplied by the iterator
+	// configured with the given set of template parameters.
+	template<uint32_t ...n>
+	static uint32_t popcount(const uint32_t limit)
+	{
+		// XXX Actually could be less than limit, if start is closer to the end.
+		return limit;
+	}
+};
+
+#endif // COMBINATIONS_H
+

ccsrc/include/cppsim/compiler.h

@@ -0,0 +1,19 @@
+#ifndef COMPILER_H
+#define COMPILER_H
+
+#include <string>
+
+class Compiler
+{
+public :
+
+	void* codegen(int nqubits, const std::string& source, std::string& errmsg);
+
+	Compiler();
+
+};
+
+Compiler& get_compiler();
+
+#endif // COMPILER_H
+

ccsrc/include/cppsim/cppsim_omp.hpp

@@ -0,0 +1,18 @@
+#ifndef CPPSIM_OMP_HPP
+#define CPPSIM_OMP_HPP
+
+#include <cstdint>
+
+#if defined(_OPENMP)
+#  include <omp.h>
+#endif
+
+namespace omp {
+#ifdef _MSC_VER
+using idx_t = int64_t;
+#else
+using idx_t = uint64_t;
+#endif  // _MSC_VER
+}  // namespace omp
+
+#endif /* CPPSIM_OMP_HPP */