Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F80375967
fe_engine_inline_impl.hh
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 31, 10:09
Size
6 KB
Mime Type
text/x-c++
Expires
Mon, Sep 2, 10:09 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
20355760
Attached To
rAKA akantu
fe_engine_inline_impl.hh
View Options
/**
* @file fe_engine_inline_impl.hh
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date creation: Tue Jul 20 2010
* @date last modification: Fri Dec 11 2020
*
* @brief Implementation of the inline functions of the FEEngine Class
*
*
* @section LICENSE
*
* Copyright (©) 2010-2021 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 "element_class.hh"
#include "fe_engine.hh"
#include "mesh.hh"
/* -------------------------------------------------------------------------- */
#include "element_type_conversion.hh"
/* -------------------------------------------------------------------------- */
#ifndef AKANTU_FE_ENGINE_INLINE_IMPL_HH_
#define AKANTU_FE_ENGINE_INLINE_IMPL_HH_
namespace
akantu
{
/* -------------------------------------------------------------------------- */
inline
Real
FEEngine
::
getElementInradius
(
const
Matrix
<
Real
>
&
coord
,
ElementType
type
)
{
Real
inradius
=
0
;
#define GET_INRADIUS(type) inradius = ElementClass<type>::getInradius(coord);
AKANTU_BOOST_ALL_ELEMENT_SWITCH
(
GET_INRADIUS
);
#undef GET_INRADIUS
return
inradius
;
}
/* -------------------------------------------------------------------------- */
inline
Real
FEEngine
::
getElementInradius
(
const
Element
&
element
)
const
{
auto
spatial_dimension
=
mesh
.
getSpatialDimension
();
auto
positions
=
make_view
(
mesh
.
getNodes
(),
spatial_dimension
).
begin
();
auto
connectivity
=
mesh
.
getConnectivities
().
get
(
element
);
Matrix
<
Real
>
coords
(
spatial_dimension
,
connectivity
.
size
());
for
(
auto
&&
data
:
zip
(
connectivity
,
coords
))
{
Vector
<
Real
>
(
std
::
get
<
1
>
(
data
))
=
Vector
<
Real
>
(
positions
[
std
::
get
<
0
>
(
data
)]);
}
return
getElementInradius
(
coords
,
element
.
type
);
}
/* -------------------------------------------------------------------------- */
inline
InterpolationType
FEEngine
::
getInterpolationType
(
ElementType
type
)
{
return
convertType
<
ElementType
,
InterpolationType
>
(
type
);
}
/* -------------------------------------------------------------------------- */
/// @todo rewrite this function in order to get the cohesive element
/// type directly from the facet
#if defined(AKANTU_COHESIVE_ELEMENT)
inline
ElementType
FEEngine
::
getCohesiveElementType
(
ElementType
type
)
{
ElementType
ctype
;
#define GET_COHESIVE_TYPE(type) \
ctype = CohesiveFacetProperty<type>::cohesive_type;
AKANTU_BOOST_ALL_ELEMENT_SWITCH
(
GET_COHESIVE_TYPE
);
#undef GET_COHESIVE_TYPE
return
ctype
;
}
#else
inline
ElementType
FEEngine
::
getCohesiveElementType
(
__attribute__
((
unused
))
ElementType
type_facet
)
{
return
_not_defined
;
}
#endif
/* -------------------------------------------------------------------------- */
#if defined(AKANTU_IGFEM)
}
// akantu
#include "igfem_helper.hh"
namespace
akantu
{
inline
Vector
<
ElementType
>
FEEngine
::
getIGFEMElementTypes
(
ElementType
type
)
{
#define GET_IGFEM_ELEMENT_TYPES(type) \
return IGFEMHelper::getIGFEMElementTypes<type>();
AKANTU_BOOST_REGULAR_ELEMENT_SWITCH
(
GET_IGFEM_ELEMENT_TYPES
);
#undef GET_IGFEM_ELEMENT_TYPES
}
#endif
/* -------------------------------------------------------------------------- */
template
<
typename
T
>
void
FEEngine
::
extractNodalToElementField
(
const
Mesh
&
mesh
,
const
Array
<
T
>
&
nodal_f
,
Array
<
T
>
&
elemental_f
,
ElementType
type
,
GhostType
ghost_type
,
const
Array
<
UInt
>
&
filter_elements
)
{
AKANTU_DEBUG_IN
();
UInt
nb_nodes_per_element
=
Mesh
::
getNbNodesPerElement
(
type
);
UInt
nb_degree_of_freedom
=
nodal_f
.
getNbComponent
();
UInt
nb_element
=
mesh
.
getNbElement
(
type
,
ghost_type
);
UInt
*
conn_val
=
mesh
.
getConnectivity
(
type
,
ghost_type
).
storage
();
if
(
filter_elements
!=
empty_filter
)
{
nb_element
=
filter_elements
.
size
();
}
elemental_f
.
resize
(
nb_element
);
T
*
nodal_f_val
=
nodal_f
.
storage
();
T
*
f_val
=
elemental_f
.
storage
();
UInt
*
el_conn
;
for
(
UInt
el
=
0
;
el
<
nb_element
;
++
el
)
{
if
(
filter_elements
!=
empty_filter
)
{
el_conn
=
conn_val
+
filter_elements
(
el
)
*
nb_nodes_per_element
;
}
else
{
el_conn
=
conn_val
+
el
*
nb_nodes_per_element
;
}
for
(
UInt
n
=
0
;
n
<
nb_nodes_per_element
;
++
n
)
{
UInt
node
=
*
(
el_conn
+
n
);
std
::
copy
(
nodal_f_val
+
node
*
nb_degree_of_freedom
,
nodal_f_val
+
(
node
+
1
)
*
nb_degree_of_freedom
,
f_val
);
f_val
+=
nb_degree_of_freedom
;
}
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
template
<
typename
T
>
void
FEEngine
::
filterElementalData
(
const
Mesh
&
mesh
,
const
Array
<
T
>
&
elem_f
,
Array
<
T
>
&
filtered_f
,
ElementType
type
,
GhostType
ghost_type
,
const
Array
<
UInt
>
&
filter_elements
)
{
AKANTU_DEBUG_IN
();
UInt
nb_element
=
mesh
.
getNbElement
(
type
,
ghost_type
);
if
(
nb_element
==
0
)
{
filtered_f
.
resize
(
0
);
return
;
}
UInt
nb_degree_of_freedom
=
elem_f
.
getNbComponent
();
UInt
nb_data_per_element
=
elem_f
.
size
()
/
nb_element
;
if
(
filter_elements
!=
empty_filter
)
{
nb_element
=
filter_elements
.
size
();
}
filtered_f
.
resize
(
nb_element
*
nb_data_per_element
);
T
*
elem_f_val
=
elem_f
.
storage
();
T
*
f_val
=
filtered_f
.
storage
();
UInt
el_offset
;
for
(
UInt
el
=
0
;
el
<
nb_element
;
++
el
)
{
if
(
filter_elements
!=
empty_filter
)
{
el_offset
=
filter_elements
(
el
);
}
else
{
el_offset
=
el
;
}
std
::
copy
(
elem_f_val
+
el_offset
*
nb_data_per_element
*
nb_degree_of_freedom
,
elem_f_val
+
(
el_offset
+
1
)
*
nb_data_per_element
*
nb_degree_of_freedom
,
f_val
);
f_val
+=
nb_degree_of_freedom
*
nb_data_per_element
;
}
AKANTU_DEBUG_OUT
();
}
}
// namespace akantu
#endif
/* AKANTU_FE_ENGINE_INLINE_IMPL_HH_ */
Event Timeline
Log In to Comment