Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F86222558
Vector.h
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, Oct 5, 01:20
Size
7 KB
Mime Type
text/x-c++
Expires
Mon, Oct 7, 01:20 (2 d)
Engine
blob
Format
Raw Data
Handle
21374160
Attached To
rGOOSEFEM GooseFEM
Vector.h
View Options
/* =================================================================================================
(c - GPLv3) T.W.J. de Geus (Tom) | tom@geus.me | www.geus.me | github.com/tdegeus/GooseFEM
================================================================================================= */
#ifndef GOOSEFEM_VECTOR_H
#define GOOSEFEM_VECTOR_H
// -------------------------------------------------------------------------------------------------
#include "GooseFEM.h"
// =========================================== GooseFEM ============================================
namespace
GooseFEM
{
// -------------------------------------------------------------------------------------------------
class
Vector
{
private:
// information
xt
::
xtensor
<
size_t
,
2
>
m_conn
;
// connectivity [nelem, nne ]
xt
::
xtensor
<
size_t
,
2
>
m_dofs
;
// DOF-numbers per node [nnode, ndim]
xt
::
xtensor
<
size_t
,
2
>
m_part
;
// DOF-numbers per node, after partitioning [nnode, ndim]
xt
::
xtensor
<
size_t
,
1
>
m_iiu
;
// DOF-numbers that are unknown [nnu]
xt
::
xtensor
<
size_t
,
1
>
m_iip
;
// DOF-numbers that are prescribed [nnp]
// dimensions
size_t
m_nelem
;
// number of elements
size_t
m_nne
;
// number of nodes per element
size_t
m_nnode
;
// number of nodes
size_t
m_ndim
;
// number of dimensions
size_t
m_ndof
;
// number of DOFs
size_t
m_nnu
;
// number of unknown DOFs
size_t
m_nnp
;
// number of prescribed DOFs
public:
// notation:
// "nodevec" - nodal vectors - [nnode, ndim]
// "elemvec" - nodal vectors stored per element - [nelem, nne, ndim]
// "dofval" - DOF values - [ndof]
// "dofval_u" - DOF values (Unknown) - [nnu]
// "dofval_p" - DOF values (Prescribed) - [nnp]
// constructor
Vector
(){};
Vector
(
const
xt
::
xtensor
<
size_t
,
2
>
&
conn
,
const
xt
::
xtensor
<
size_t
,
2
>
&
dofs
);
Vector
(
const
xt
::
xtensor
<
size_t
,
2
>
&
conn
,
const
xt
::
xtensor
<
size_t
,
2
>
&
dofs
,
const
xt
::
xtensor
<
size_t
,
1
>
&
iip
);
// dimensions
size_t
nelem
()
const
;
// number of elements
size_t
nne
()
const
;
// number of nodes per element
size_t
nnode
()
const
;
// number of nodes
size_t
ndim
()
const
;
// number of dimensions
size_t
ndof
()
const
;
// number of DOFs
size_t
nnu
()
const
;
// number of unknown DOFs
size_t
nnp
()
const
;
// number of prescribed DOFs
// DOF lists
xt
::
xtensor
<
size_t
,
2
>
dofs
()
const
;
// DOFs
xt
::
xtensor
<
size_t
,
1
>
iiu
()
const
;
// unknown DOFs
xt
::
xtensor
<
size_t
,
1
>
iip
()
const
;
// prescribed DOFs
// convert vectors (overwrite entries that occur more that once) -- no allocation
void
asDofs
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs_u
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs_u
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs_p
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asDofs_p
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
asNode
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval
,
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
void
asNode
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
,
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
void
asNode
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
void
asElement
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval
,
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
void
asElement
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
,
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
void
asElement
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
// assemble vectors (adds entries that occur more that once) -- no allocation
void
assembleDofs
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleDofs
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleDofs_u
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleDofs_u
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleDofs_p
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleDofs_p
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
void
assembleNode
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
,
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
// convert vectors (overwrite entries that occur more that once)
xt
::
xtensor
<
double
,
1
>
asDofs
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs_u
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs_u
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs_p
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
asDofs_p
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
2
>
asNode
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
xt
::
xtensor
<
double
,
2
>
asNode
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
)
const
;
xt
::
xtensor
<
double
,
2
>
asNode
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
3
>
asElement
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval
)
const
;
xt
::
xtensor
<
double
,
3
>
asElement
(
const
xt
::
xtensor
<
double
,
1
>
&
dofval_u
,
const
xt
::
xtensor
<
double
,
1
>
&
dofval_p
)
const
;
xt
::
xtensor
<
double
,
3
>
asElement
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
// assemble vectors (adds entries that occur more that once)
xt
::
xtensor
<
double
,
1
>
assembleDofs
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
assembleDofs
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
1
>
assembleDofs_u
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
assembleDofs_u
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
1
>
assembleDofs_p
(
const
xt
::
xtensor
<
double
,
2
>
&
nodevec
)
const
;
xt
::
xtensor
<
double
,
1
>
assembleDofs_p
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
xt
::
xtensor
<
double
,
2
>
assembleNode
(
const
xt
::
xtensor
<
double
,
3
>
&
elemvec
)
const
;
};
// -------------------------------------------------------------------------------------------------
}
// namespace ...
// =================================================================================================
#endif
Event Timeline
Log In to Comment