Page MenuHomec4science
Files F62072643

Mesh.h
No OneTemporary
Actions

File Metadata

Created
Fri, May 10, 18:13

Mesh.h
View Options

/*
(c - GPLv3) T.W.J. de Geus (Tom) | tom@geus.me | www.geus.me | github.com/tdegeus/GooseFEM
*/
#ifndef GOOSEFEM_MESH_H
#define GOOSEFEM_MESH_H
#include "config.h"
namespace GooseFEM {
namespace Mesh {
enum class ElementType {
Quad4,
Hex8,
Tri3 };
// Renumber to lowest possible index. For example [0,3,4,2] -> [0,2,3,1]
class Renumber {
public:
// constructors
Renumber() = default;
Renumber(const xt::xarray<size_t>& dofs);
// get renumbered DOFs (same as "Renumber::apply(dofs)")
xt::xtensor<size_t, 2> get(const xt::xtensor<size_t, 2>& dofs) const;
// apply renumbering to other set
template <class T>
T apply(const T& list) const;
// get the list needed to renumber, e.g.:
// dofs_renumbered(i,j) = index(dofs(i,j))
xt::xtensor<size_t, 1> index() const;
private:
xt::xtensor<size_t, 1> m_renum;
};
// Reorder to lowest possible index, in specific order.
//
// For example for "Reorder({iiu,iip})" after reordering:
//
// iiu = xt::range<size_t>(nnu);
// iip = xt::range<size_t>(nnp) + nnu;
class Reorder {
public:
// constructors
Reorder() = default;
Reorder(const std::initializer_list<xt::xtensor<size_t, 1>> args);
// get reordered DOFs (same as "Reorder::apply(dofs)")
xt::xtensor<size_t, 2> get(const xt::xtensor<size_t, 2>& dofs) const;
// apply renumbering to other set
template <class T>
T apply(const T& list) const;
// get the list needed to reorder, e.g.:
// dofs_reordered(i,j) = index(dofs(i,j))
xt::xtensor<size_t, 1> index() const;
private:
xt::xtensor<size_t, 1> m_renum;
};
// list with DOF-numbers in sequential order
inline xt::xtensor<size_t, 2> dofs(size_t nnode, size_t ndim);
// renumber to lowest possible index (see "GooseFEM::Mesh::Renumber")
inline xt::xtensor<size_t, 2> renumber(const xt::xtensor<size_t, 2>& dofs);
// number of elements connected to each node
inline xt::xtensor<size_t, 1> coordination(const xt::xtensor<size_t, 2>& conn);
// elements connected to each node
inline std::vector<std::vector<size_t>> elem2node(const xt::xtensor<size_t, 2>& conn);
// return size of each element edge
inline xt::xtensor<double, 2> edgesize(
const xt::xtensor<double, 2>& coor, const xt::xtensor<size_t, 2>& conn, ElementType type);
// return size of each element edge: extract element-type based on shape of "conn"
inline xt::xtensor<double, 2> edgesize(
const xt::xtensor<double, 2>& coor, const xt::xtensor<size_t, 2>& conn);
} // namespace Mesh
} // namespace GooseFEM
#include "Mesh.hpp"
#endif

Event Timeline

c4science · Help