blob: 34518a5a0e3d7f6d33e13eb4770ed3b68021affc [file] [log] [blame]
// { dg-do run }
// { dg-additional-options -DMEM_SHARED { target offload_device_shared_as } }
#include <stdlib.h>
#include <time.h>
#include <map>
#include <omp.h>
// Make sure that KEY_MAX is less than N to ensure some duplicate keys.
#define N 3000
#define KEY_MAX 1000
void init (int data[], int max)
{
for (int i = 0; i < N; ++i)
data[i] = rand () % max;
}
bool validate (long long sum, int keys[], int data[])
{
long long total = 0;
for (int i = 0; i < N; ++i)
total += (long long) keys[i] * data[i];
return sum == total;
}
int main (void)
{
int keys[N], data[N];
std::multimap<int,int> _map;
srand (time (NULL));
init (keys, KEY_MAX);
init (data, RAND_MAX);
#ifndef MEM_SHARED
#pragma omp target enter data map (to: keys[ :N], data[ :N]) map (alloc: _map)
#endif
#pragma omp target
{
#ifndef MEM_SHARED
new (&_map) std::multimap<int,int> ();
#endif
for (int i = 0; i < N; ++i)
_map.insert({keys[i], data[i]});
}
long long sum = 0;
#pragma omp target teams distribute parallel for reduction (+:sum)
for (int i = 0; i < KEY_MAX; ++i)
{
auto range = _map.equal_range (i);
for (auto it = range.first; it != range.second; ++it)
sum += (long long) it->first * it->second;
}
#ifdef MEM_SHARED
/* Even with USM, memory allocated on the device (with _map.insert)
must be freed on the device. */
if (omp_get_default_device () != omp_initial_device
&& omp_get_default_device () != omp_get_num_devices ())
{
#pragma omp target
_map.clear ();
}
#endif
#ifndef MEM_SHARED
#pragma omp target
_map.~multimap ();
#endif
#ifndef MEM_SHARED
#pragma omp target exit data map (release: _map)
#endif
bool ok = validate (sum, keys, data);
return ok ? 0 : 1;
}