Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F79069975
test_fe_engine_precomputation.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
Sat, Aug 24, 14:38
Size
4 KB
Mime Type
text/x-c++
Expires
Mon, Aug 26, 14:38 (2 d)
Engine
blob
Format
Raw Data
Handle
20011044
Attached To
rAKA akantu
test_fe_engine_precomputation.cc
View Options
/**
* @file test_fe_engine_precomputation.cc
*
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date creation: Mon Jun 14 2010
* @date last modification: Mon Feb 19 2018
*
* @brief test of the fem class
*
*
* Copyright (©) 2010-2018 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* 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 "py_aka_array.hh"
#include "test_fe_engine_fixture.hh"
/* -------------------------------------------------------------------------- */
#include <pybind11/embed.h>
#include <pybind11/numpy.h>
/* -------------------------------------------------------------------------- */
using
namespace
akantu
;
namespace
py
=
pybind11
;
using
namespace
py
::
literals
;
template
<
class
T
>
decltype
(
auto
)
make_proxy
(
Array
<
T
>
&
array
)
{
return
detail
::
ArrayProxy
<
T
>
(
array
);
}
template
<
typename
type_
>
class
TestFEMPyFixture
:
public
TestFEMFixture
<
type_
>
{
using
parent
=
TestFEMFixture
<
type_
>
;
public
:
void
SetUp
()
override
{
parent
::
SetUp
();
const
auto
&
connectivities
=
this
->
mesh
->
getConnectivity
(
this
->
type
);
const
auto
&
nodes
=
this
->
mesh
->
getNodes
().
begin
(
this
->
dim
);
coordinates
=
std
::
make_unique
<
Array
<
Real
>>
(
connectivities
.
size
(),
connectivities
.
getNbComponent
()
*
this
->
dim
);
for
(
auto
&&
tuple
:
zip
(
make_view
(
connectivities
,
connectivities
.
getNbComponent
()),
make_view
(
*
coordinates
,
this
->
dim
,
connectivities
.
getNbComponent
())))
{
const
auto
&
conn
=
std
::
get
<
0
>
(
tuple
);
const
auto
&
X
=
std
::
get
<
1
>
(
tuple
);
for
(
auto
s
:
arange
(
conn
.
size
()))
{
Vector
<
Real
>
(
X
(
s
))
=
Vector
<
Real
>
(
nodes
[
conn
(
s
)]);
}
}
}
void
TearDown
()
override
{
parent
::
TearDown
();
coordinates
.
reset
(
nullptr
);
}
protected
:
std
::
unique_ptr
<
Array
<
Real
>>
coordinates
;
};
TYPED_TEST_SUITE
(
TestFEMPyFixture
,
fe_engine_types
,
);
TYPED_TEST
(
TestFEMPyFixture
,
Precompute
)
{
SCOPED_TRACE
(
std
::
to_string
(
this
->
type
));
this
->
fem
->
initShapeFunctions
();
const
auto
&
N
=
this
->
fem
->
getShapeFunctions
().
getShapes
(
this
->
type
);
const
auto
&
B
=
this
->
fem
->
getShapeFunctions
().
getShapesDerivatives
(
this
->
type
);
const
auto
&
j
=
this
->
fem
->
getIntegrator
().
getJacobians
(
this
->
type
);
// Array<Real> ref_N(this->nb_quadrature_points_total, N.getNbComponent());
// Array<Real> ref_B(this->nb_quadrature_points_total, B.getNbComponent());
Array
<
Real
>
ref_j
(
this
->
nb_quadrature_points_total
,
j
.
getNbComponent
());
auto
ref_N
(
N
);
auto
ref_B
(
B
);
py
::
module
py_engine
=
py
::
module
::
import
(
"py_engine"
);
auto
py_shape
=
py_engine
.
attr
(
"Shapes"
)(
py
::
str
(
std
::
to_string
(
this
->
type
)));
auto
kwargs
=
py
::
dict
(
"N"
_a
=
ref_N
,
"B"
_a
=
ref_B
,
"j"
_a
=
ref_j
,
"X"
_a
=
*
this
->
coordinates
,
"Q"
_a
=
this
->
fem
->
getIntegrationPoints
(
this
->
type
));
auto
ret
=
py_shape
.
attr
(
"precompute"
)(
**
kwargs
);
auto
check
=
[
&
](
auto
&
ref_A
,
auto
&
A
,
const
auto
&
id
)
{
SCOPED_TRACE
(
std
::
to_string
(
this
->
type
)
+
" "
+
id
);
for
(
auto
&&
n
:
zip
(
make_view
(
ref_A
,
ref_A
.
getNbComponent
()),
make_view
(
A
,
A
.
getNbComponent
())))
{
auto
diff
=
(
std
::
get
<
0
>
(
n
)
-
std
::
get
<
1
>
(
n
)).
template
norm
<
L_inf
>
();
EXPECT_NEAR
(
0.
,
diff
,
1e-10
);
}
};
check
(
ref_N
,
N
,
"N"
);
check
(
ref_B
,
B
,
"B"
);
check
(
ref_j
,
j
,
"j"
);
}
Event Timeline
Log In to Comment