Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F78840820
test_data_distribution.cc
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
Fri, Aug 23, 06:54
Size
2 KB
Mime Type
text/x-c
Expires
Sun, Aug 25, 06:54 (2 d)
Engine
blob
Format
Raw Data
Handle
20111140
Attached To
rAKA akantu
test_data_distribution.cc
View Options
/**
* Copyright (©) 2014-2023 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This file is part of Akantu
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*/
/* -------------------------------------------------------------------------- */
#include "test_synchronizers_fixture.hh"
/* -------------------------------------------------------------------------- */
TEST_F
(
TestSynchronizerFixture
,
DataDistribution
)
{
auto
&
barycenters
=
this
->
mesh
->
getElementalData
<
Real
>
(
"barycenters"
);
auto
spatial_dimension
=
this
->
mesh
->
getSpatialDimension
();
barycenters
.
initialize
(
*
this
->
mesh
,
_spatial_dimension
=
_all_dimensions
,
_nb_component
=
spatial_dimension
);
this
->
initBarycenters
(
barycenters
,
*
this
->
mesh
);
auto
&
gids
=
this
->
mesh
->
getNodalData
<
UInt
>
(
"gid"
);
gids
.
resize
(
this
->
mesh
->
getNbNodes
());
for
(
auto
&&
data
:
enumerate
(
gids
))
{
std
::
get
<
1
>
(
data
)
=
std
::
get
<
0
>
(
data
);
}
this
->
distribute
();
for
(
auto
&&
ghost_type
:
ghost_types
)
{
for
(
const
auto
&
type
:
this
->
mesh
->
elementTypes
(
_all_dimensions
,
ghost_type
))
{
auto
&
barycenters
=
this
->
mesh
->
getData
<
Real
>
(
"barycenters"
,
type
,
ghost_type
);
for
(
auto
&&
data
:
enumerate
(
make_view
(
barycenters
,
spatial_dimension
)))
{
Element
element
{
type
,
std
::
get
<
0
>
(
data
),
ghost_type
};
Vector
<
Real
>
barycenter
(
spatial_dimension
);
this
->
mesh
->
getBarycenter
(
element
,
barycenter
);
auto
dist
=
(
std
::
get
<
1
>
(
data
)
-
barycenter
).
template
lpNorm
<
Eigen
::
Infinity
>
();
EXPECT_NEAR
(
dist
,
0
,
1e-7
);
}
}
}
if
(
psize
>
1
)
{
for
(
auto
&&
data
:
zip
(
gids
,
this
->
mesh
->
getGlobalNodesIds
()))
{
EXPECT_EQ
(
std
::
get
<
0
>
(
data
),
std
::
get
<
1
>
(
data
));
}
}
}
TEST_F
(
TestSynchronizerFixture
,
DataDistributionTags
)
{
this
->
distribute
();
for
(
const
auto
&
type
:
this
->
mesh
->
elementTypes
(
_all_dimensions
))
{
auto
&
tags
=
this
->
mesh
->
getData
<
Int
>
(
"tag_0"
,
type
);
auto
tags_it
=
tags
.
begin
(
1
);
auto
tags_end
=
tags
.
end
(
1
);
// The number of tags should match the number of elements on rank"
EXPECT_EQ
(
this
->
mesh
->
getNbElement
(
type
),
tags
.
size
());
}
}
Event Timeline
Log In to Comment