Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F90739822
TestSegmentedView.hpp
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Mon, Nov 4, 08:06
Size
26 KB
Mime Type
text/x-c
Expires
Wed, Nov 6, 08:06 (2 d)
Engine
blob
Format
Raw Data
Handle
22095331
Attached To
rLAMMPS lammps
TestSegmentedView.hpp
View Options
/*
//@HEADER
// ************************************************************************
//
// Kokkos v. 2.0
// Copyright (2014) Sandia Corporation
//
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are
// met:
//
// 1. Redistributions of source code must retain the above copyright
// notice, this list of conditions and the following disclaimer.
//
// 2. Redistributions in binary form must reproduce the above copyright
// notice, this list of conditions and the following disclaimer in the
// documentation and/or other materials provided with the distribution.
//
// 3. Neither the name of the Corporation nor the names of the
// contributors may be used to endorse or promote products derived from
// this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
//
// Questions? Contact H. Carter Edwards (hcedwar@sandia.gov)
//
// ************************************************************************
//@HEADER
*/
#ifndef KOKKOS_TEST_SEGMENTEDVIEW_HPP
#define KOKKOS_TEST_SEGMENTEDVIEW_HPP
#include <gtest/gtest.h>
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <Kokkos_Core.hpp>
#if ! defined( KOKKOS_USING_EXPERIMENTAL_VIEW )
#include <Kokkos_SegmentedView.hpp>
#include <impl/Kokkos_Timer.hpp>
namespace
Test
{
namespace
Impl
{
template
<
class
ViewType
,
class
ExecutionSpace
,
int
Rank
=
ViewType
::
Rank
>
struct
GrowTest
;
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
1
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
value
+=
team_idx
+
team_member
.
team_rank
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
a
(
team_idx
+
team_member
.
team_rank
())
=
team_idx
+
team_member
.
team_rank
();
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
2
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
{
a
(
team_idx
+
team_member
.
team_rank
(),
k
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
;
}
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
3
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
4
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
2
;
m
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
5
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
2
;
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
3
;
n
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
6
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
2
;
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
3
;
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
2
;
o
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
7
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
2
;
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
3
;
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
2
;
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
4
;
p
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
+
15
*
p
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
a
.
dimension_6
();
p
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
,
p
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
+
15
*
p
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
GrowTest
<
ViewType
,
ExecutionSpace
,
8
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
GrowTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
a
.
grow
(
team_member
,
team_idx
+
team_member
.
team_size
());
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
7
;
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
3
;
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
2
;
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
3
;
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
2
;
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
4
;
p
++
)
for
(
typename
ExecutionSpace
::
size_type
q
=
0
;
q
<
3
;
q
++
)
value
+=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
+
15
*
p
+
17
*
q
;
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
a
.
dimension_6
();
p
++
)
for
(
typename
ExecutionSpace
::
size_type
q
=
0
;
q
<
a
.
dimension_7
();
q
++
)
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
,
p
,
q
)
=
team_idx
+
team_member
.
team_rank
()
+
13
*
k
+
3
*
l
+
7
*
m
+
5
*
n
+
2
*
o
+
15
*
p
+
17
*
q
;
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
,
int
Rank
=
ViewType
::
Rank
>
struct
VerifyTest
;
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
1
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
value
+=
a
(
team_idx
+
team_member
.
team_rank
());
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
2
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
3
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
4
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
5
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
6
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
7
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
a
.
dimension_6
();
p
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
,
p
);
}
}
};
template
<
class
ViewType
,
class
ExecutionSpace
>
struct
VerifyTest
<
ViewType
,
ExecutionSpace
,
8
>
{
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
typedef
typename
Policy
::
member_type
team_type
;
typedef
double
value_type
;
ViewType
a
;
VerifyTest
(
ViewType
in
)
:
a
(
in
)
{}
KOKKOS_INLINE_FUNCTION
void
operator
()
(
team_type
team_member
,
double
&
value
)
const
{
unsigned
int
team_idx
=
team_member
.
league_rank
()
*
team_member
.
team_size
();
if
((
a
.
dimension_0
()
>
team_idx
+
team_member
.
team_rank
())
&&
(
a
.
dimension
(
0
)
>
team_idx
+
team_member
.
team_rank
()))
{
for
(
typename
ExecutionSpace
::
size_type
k
=
0
;
k
<
a
.
dimension_1
();
k
++
)
for
(
typename
ExecutionSpace
::
size_type
l
=
0
;
l
<
a
.
dimension_2
();
l
++
)
for
(
typename
ExecutionSpace
::
size_type
m
=
0
;
m
<
a
.
dimension_3
();
m
++
)
for
(
typename
ExecutionSpace
::
size_type
n
=
0
;
n
<
a
.
dimension_4
();
n
++
)
for
(
typename
ExecutionSpace
::
size_type
o
=
0
;
o
<
a
.
dimension_5
();
o
++
)
for
(
typename
ExecutionSpace
::
size_type
p
=
0
;
p
<
a
.
dimension_6
();
p
++
)
for
(
typename
ExecutionSpace
::
size_type
q
=
0
;
q
<
a
.
dimension_7
();
q
++
)
value
+=
a
(
team_idx
+
team_member
.
team_rank
(),
k
,
l
,
m
,
n
,
o
,
p
,
q
);
}
}
};
template
<
typename
Scalar
,
class
ExecutionSpace
>
struct
test_segmented_view
{
typedef
test_segmented_view
<
Scalar
,
ExecutionSpace
>
self_type
;
typedef
Scalar
scalar_type
;
typedef
ExecutionSpace
execution_space
;
typedef
Kokkos
::
TeamPolicy
<
execution_space
>
Policy
;
double
result
;
double
reference
;
template
<
class
ViewType
>
void
run_me
(
ViewType
a
,
int
max_length
){
const
int
team_size
=
Policy
::
team_size_max
(
GrowTest
<
ViewType
,
execution_space
>
(
a
)
);
const
int
nteams
=
max_length
/
team_size
;
reference
=
0
;
result
=
0
;
Kokkos
::
parallel_reduce
(
Policy
(
nteams
,
team_size
),
GrowTest
<
ViewType
,
execution_space
>
(
a
),
reference
);
Kokkos
::
fence
();
Kokkos
::
parallel_reduce
(
Policy
(
nteams
,
team_size
),
VerifyTest
<
ViewType
,
execution_space
>
(
a
),
result
);
Kokkos
::
fence
();
}
test_segmented_view
(
unsigned
int
size
,
int
rank
)
{
reference
=
0
;
result
=
0
;
const
int
dim_1
=
7
;
const
int
dim_2
=
3
;
const
int
dim_3
=
2
;
const
int
dim_4
=
3
;
const
int
dim_5
=
2
;
const
int
dim_6
=
4
;
//const int dim_7 = 3;
if
(
rank
==
1
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*
,
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
rank1_view
;
run_me
<
rank1_view
>
(
rank1_view
(
"Rank1"
,
128
,
size
),
size
);
}
if
(
rank
==
2
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
**
,
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
rank2_view
;
run_me
<
rank2_view
>
(
rank2_view
(
"Rank2"
,
128
,
size
,
dim_1
),
size
);
}
if
(
rank
==
3
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*
[
7
][
3
][
2
],
Kokkos
::
LayoutRight
,
ExecutionSpace
>
rank3_view
;
run_me
<
rank3_view
>
(
rank3_view
(
"Rank3"
,
128
,
size
),
size
);
}
if
(
rank
==
4
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
****
,
Kokkos
::
LayoutRight
,
ExecutionSpace
>
rank4_view
;
run_me
<
rank4_view
>
(
rank4_view
(
"Rank4"
,
128
,
size
,
dim_1
,
dim_2
,
dim_3
),
size
);
}
if
(
rank
==
5
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*
[
7
][
3
][
2
][
3
],
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
rank5_view
;
run_me
<
rank5_view
>
(
rank5_view
(
"Rank5"
,
128
,
size
),
size
);
}
if
(
rank
==
6
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*****
[
2
],
Kokkos
::
LayoutRight
,
ExecutionSpace
>
rank6_view
;
run_me
<
rank6_view
>
(
rank6_view
(
"Rank6"
,
128
,
size
,
dim_1
,
dim_2
,
dim_3
,
dim_4
),
size
);
}
if
(
rank
==
7
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*******
,
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
rank7_view
;
run_me
<
rank7_view
>
(
rank7_view
(
"Rank7"
,
128
,
size
,
dim_1
,
dim_2
,
dim_3
,
dim_4
,
dim_5
,
dim_6
),
size
);
}
if
(
rank
==
8
)
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*****
[
2
][
4
][
3
],
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
rank8_view
;
run_me
<
rank8_view
>
(
rank8_view
(
"Rank8"
,
128
,
size
,
dim_1
,
dim_2
,
dim_3
,
dim_4
),
size
);
}
}
};
}
// namespace Impl
template
<
typename
Scalar
,
class
ExecutionSpace
>
void
test_segmented_view
(
unsigned
int
size
)
{
{
typedef
Kokkos
::
Experimental
::
SegmentedView
<
Scalar
*****
[
2
][
4
][
3
],
Kokkos
::
LayoutLeft
,
ExecutionSpace
>
view_type
;
view_type
a
(
"A"
,
128
,
size
,
7
,
3
,
2
,
3
);
double
reference
;
Impl
::
GrowTest
<
view_type
,
ExecutionSpace
>
f
(
a
);
const
int
team_size
=
Kokkos
::
TeamPolicy
<
ExecutionSpace
>::
team_size_max
(
f
);
const
int
nteams
=
(
size
+
team_size
-
1
)
/
team_size
;
Kokkos
::
parallel_reduce
(
Kokkos
::
TeamPolicy
<
ExecutionSpace
>
(
nteams
,
team_size
),
f
,
reference
);
size_t
real_size
=
((
size
+
127
)
/
128
)
*
128
;
ASSERT_EQ
(
real_size
,
a
.
dimension_0
());
ASSERT_EQ
(
7
,
a
.
dimension_1
());
ASSERT_EQ
(
3
,
a
.
dimension_2
());
ASSERT_EQ
(
2
,
a
.
dimension_3
());
ASSERT_EQ
(
3
,
a
.
dimension_4
());
ASSERT_EQ
(
2
,
a
.
dimension_5
());
ASSERT_EQ
(
4
,
a
.
dimension_6
());
ASSERT_EQ
(
3
,
a
.
dimension_7
());
ASSERT_EQ
(
real_size
,
a
.
dimension
(
0
));
ASSERT_EQ
(
7
,
a
.
dimension
(
1
));
ASSERT_EQ
(
3
,
a
.
dimension
(
2
));
ASSERT_EQ
(
2
,
a
.
dimension
(
3
));
ASSERT_EQ
(
3
,
a
.
dimension
(
4
));
ASSERT_EQ
(
2
,
a
.
dimension
(
5
));
ASSERT_EQ
(
4
,
a
.
dimension
(
6
));
ASSERT_EQ
(
3
,
a
.
dimension
(
7
));
ASSERT_EQ
(
8
,
a
.
Rank
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
1
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
2
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
3
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
4
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
5
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
6
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
7
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
{
Impl
::
test_segmented_view
<
Scalar
,
ExecutionSpace
>
test
(
size
,
8
);
ASSERT_EQ
(
test
.
reference
,
test
.
result
);
}
}
}
// namespace Test
#else
template
<
typename
Scalar
,
class
ExecutionSpace
>
void
test_segmented_view
(
unsigned
int
)
{}
#endif
#endif
/* #ifndef KOKKOS_TEST_SEGMENTEDVIEW_HPP */
Event Timeline
Log In to Comment