Further cleanup on simd, add a simd debug case, disable std random nu…

…mber generator example

Makefile

@@ -1,12 +1,14 @@
-all: debug example
+all: debug example debug_simd
 
 debug:
-	clang++ -std=c++17 -O0 -g example.cpp -o debug
+	clang++ -std=c++17 -O0 -g -Wall example.cpp -o debug
 example:
 	clang++ -std=c++17 -O3 -Wall -march=native example.cpp -o example
 debug_simd:
-	clang++ -std=c++17 -O0 -g -march=native simd_debug.cpp -o simd_debug
+	clang++ -std=c++17 -O0 -g -Wall -march=native simd_debug.cpp -o simd_debug
+example_simd:
+	clang++ -std=c++17 -O0 -g -Wall -march=native example.cpp -o example_simd
 clean:
-	rm -rf debug example
+	rm -rf debug example debug_simd
 
 

changelist.txt

@@ -3,9 +3,7 @@ TODO:
 consider moving entirely away from std::vector of std::vector to a flattened 
 std::vector
 
-consider Nelder Mead PSO where PSO is used for last search using bounds found using Nelder Mead
-
-implement BFGS and L-BFGS-B
-
-implement Coordinate Descent and Gradient Descent
+implement L-BFGS-B, or at least BFGS-B
+implement CMAES
+implement Coordinate Descent
 

example.cpp

@@ -55,6 +55,7 @@ void print_vector(std::vector<double> &x) {
 
 // to use any C++ standard random number generator just pass in a generator
 // functor e.g. using Mersene Twister
+/*
 #include <cstdint>  // NOLINT
 #include <random>
 struct std_MT {
@@ -69,6 +70,7 @@ struct std_MT {
   }
   double operator()() { return this->distribution(this->generator); }
 };
+ */
 
 int main() {
   // define problem functor - in our case a variant of the rosenbrock function
@@ -119,19 +121,6 @@ int main() {
   de_res.print();
   print_vector(de_init);
 
-  std::cout << "Differential evolution with std::mt19937: " << std::endl;
-  // using standard library random number generators
-  std_MT std_gen;
-  // again initialize solver, this time also with the RNG
-  // if strategy is not specifieed, defaults to random
-  auto de_solver_MT = DE<Rosenbrock, std_MT, double>(prob, std_gen);
-  // reset initial state
-  de_init[0] = 2;
-  de_init[1] = 7;
-  de_res = de_solver_MT.minimize(de_init);
-  de_res.print();
-  print_vector(de_init);
-
   std::cout << "Particle Swarm Optimization with xoshiro: " << std::endl;
   // we also have a xoshiro generator
   xoshiro<double> xos_gen;