Неблокирующие алгоритмы на C++, многопоточность

atomic_int global_min; // атомарное значение, можем одновременно читать и записывать из разных потоков
int global_min_index; // не атомарное значение
 
// функция потока
void thread(int* array, int start, int stop)
{
  // цикл по заданным элементам
  for (int pos=start; pos<stop; pos++)
  {
    // сравнение с атомиком
    if (array[pos] < global_min)
    {
      // << lock ??? на запись нужно как-то заблокировать, что бы два потока не лезли одновременно
 
      // повторная проверка из-за возможности пересечения нескольких потоков
      // если какой-то поток уже внёс свои изменения, то нужно перепроверить минимум
      //  если всё ещё минимум, то обновить
      //  если минимум устарел, то ничего не делаем
      if (array[pos] < global_min)
      {
        global_min = array[pos];// запись в атомик
        global_min_index = pos;
      }
 
      // << unlock ???
    }
  }
}

// поиск минимума для части диапазона по аналогии с кодом джава
for (int pos=start; pos<stop; pos++)
{
  int prev = atomic_min_index;
  do{
    if (array[pos] >= array[prev])
      break;
  } while (!atomic_min_index.compare_exchange_strong(prev, pos));
}

#include <thread>
#include <future>
#include <random>
#include <algorithm>
 
template <size_t N>
void fill(int(&array)[N], int lower_bound, int upper_bound) noexcept {
    static std::mt19937 gen { std::random_device()() }; 
    std::uniform_int_distribution<int> uid(lower_bound, upper_bound);
 
    std::generate(array, array + N, [&uid] { return uid(gen); });
}
 
template <size_t N>
void print(int(&array)[N]) noexcept {
    for (size_t i = 0; i < N; i++) {
        printf("%d ", array[i]);
    }
}
 
int find_max(int* begin, int* end) noexcept {
    return *std::max_element(begin, end);
}
 
int main() {
    int data[100];
    fill(data, 1, 99);
    print(data);
 
    auto worker1 = std::async(find_max, data,      data + 50 );
    auto worker2 = std::async(find_max, data + 50, data + 100);
 
    int max_from_thread_1 = worker1.get();
    int max_from_thread_2 = worker2.get();
 
    printf("\nmax value = %d", std::max(max_from_thread_1, max_from_thread_2));
}