parallel_partition.cpp 3.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566
  1. // ======================================================================== //
  2. // Copyright 2009-2017 Intel Corporation //
  3. // //
  4. // Licensed under the Apache License, Version 2.0 (the "License"); //
  5. // you may not use this file except in compliance with the License. //
  6. // You may obtain a copy of the License at //
  7. // //
  8. // http://www.apache.org/licenses/LICENSE-2.0 //
  9. // //
  10. // Unless required by applicable law or agreed to in writing, software //
  11. // distributed under the License is distributed on an "AS IS" BASIS, //
  12. // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. //
  13. // See the License for the specific language governing permissions and //
  14. // limitations under the License. //
  15. // ======================================================================== //
  16. #include "parallel_partition.h"
  17. #include "../sys/regression.h"
  18. namespace embree
  19. {
  20. struct parallel_partition_regression_test : public RegressionTest
  21. {
  22. parallel_partition_regression_test(const char* name) : RegressionTest(name) {
  23. registerRegressionTest(this);
  24. }
  25. bool run ()
  26. {
  27. bool passed = true;
  28. for (size_t i=0; i<100; i++)
  29. {
  30. /* create random permutation */
  31. size_t N = std::rand() % 1000000;
  32. std::vector<unsigned> array(N);
  33. for (unsigned i=0; i<N; i++) array[i] = i;
  34. for (auto& v : array) std::swap(v,array[std::rand()%array.size()]);
  35. size_t split = std::rand() % (N+1);
  36. /* perform parallel partitioning */
  37. size_t left_sum = 0, right_sum = 0;
  38. size_t mid = parallel_partitioning(array.data(),0,array.size(),0,left_sum,right_sum,
  39. [&] ( size_t i ) { return i < split; },
  40. [] ( size_t& sum, unsigned v) { sum += v; },
  41. [] ( size_t& sum, size_t v) { sum += v; },
  42. 128);
  43. /*serial_partitioning(array.data(),0,array.size(),left_sum,right_sum,
  44. [&] ( size_t i ) { return i < split; },
  45. [] ( size_t& left_sum, int v) { left_sum += v; });*/
  46. /* verify result */
  47. passed &= mid == split;
  48. passed &= left_sum == split*(split-1)/2;
  49. passed &= right_sum == N*(N-1)/2-left_sum;
  50. for (size_t i=0; i<split; i++) passed &= array[i] < split;
  51. for (size_t i=split; i<N; i++) passed &= array[i] >= split;
  52. }
  53. return passed;
  54. }
  55. };
  56. parallel_partition_regression_test parallel_partition_regression("parallel_partition_regression_test");
  57. }