Page MenuHomec4science
Files F61315672

ElementQuad4.hpp
No OneTemporary
Actions

File Metadata

Created
Sun, May 5, 21:46

ElementQuad4.hpp
View Options

/* =================================================================================================
(c - GPLv3) T.W.J. de Geus (Tom) | tom@geus.me | www.geus.me | github.com/tdegeus/GooseFEM
================================================================================================= */
#include <Eigen/Eigen>
#include <pybind11/pybind11.h>
#include <pybind11/eigen.h>
#include <pyxtensor/pyxtensor.hpp>
#include "../include/GooseFEM/GooseFEM.h"
// =================================================================================================
namespace py = pybind11;
// =================================================================================================
void init_ElementQuad4(py::module &m)
{
py::class_<GooseFEM::Element::Quad4::Quadrature>(m, "Quadrature")
.def(py::init<const xt::xtensor<double,3> &>(), "Quadrature", py::arg("x"))
.def(py::init<const xt::xtensor<double,3> &, const xt::xtensor<double,2> &, const xt::xtensor<double,1> &>(), "Quadrature", py::arg("x"), py::arg("xi"), py::arg("w"))
.def("nelem", &GooseFEM::Element::Quad4::Quadrature::nelem, "Return number of elements" )
.def("nne" , &GooseFEM::Element::Quad4::Quadrature::nne , "Return number of nodes per element" )
.def("ndim" , &GooseFEM::Element::Quad4::Quadrature::ndim , "Return number of dimensions" )
.def("nip" , &GooseFEM::Element::Quad4::Quadrature::nip , "Return number of integration points")
.def("DV", py::overload_cast<size_t>(&GooseFEM::Element::Quad4::Quadrature::DV, py::const_), "Integration point volume (qtensor)")
.def("DV", py::overload_cast<> (&GooseFEM::Element::Quad4::Quadrature::DV, py::const_), "Integration point volume (qscalar)")
.def("GradN_vector" , py::overload_cast<const xt::xtensor<double,3> &>(&GooseFEM::Element::Quad4::Quadrature::GradN_vector , py::const_), "Dyadic product, returns 'qtensor'", py::arg("elemvec"))
.def("GradN_vector_T" , py::overload_cast<const xt::xtensor<double,3> &>(&GooseFEM::Element::Quad4::Quadrature::GradN_vector_T , py::const_), "Dyadic product, returns 'qtensor'", py::arg("elemvec"))
.def("SymGradN_vector", py::overload_cast<const xt::xtensor<double,3> &>(&GooseFEM::Element::Quad4::Quadrature::SymGradN_vector, py::const_), "Dyadic product, returns 'qtensor'", py::arg("elemvec"))
.def("Int_N_scalar_NT_dV" , py::overload_cast<const xt::xtensor<double,2> &>(&GooseFEM::Element::Quad4::Quadrature::Int_N_scalar_NT_dV , py::const_), "Integration, returns 'elemmat'", py::arg("qscalar"))
.def("Int_gradN_dot_tensor2_dV", py::overload_cast<const xt::xtensor<double,4> &>(&GooseFEM::Element::Quad4::Quadrature::Int_gradN_dot_tensor2_dV, py::const_), "Integration, returns 'elemvec'", py::arg("qtensor"))
.def("__repr__", [](const GooseFEM::Element::Quad4::Quadrature &){ return "<GooseFEM.Element.Quad4.Quadrature>"; });
}
// -------------------------------------------------------------------------------------------------
void init_ElementQuad4Gauss(py::module &m)
{
m.def("nip", &GooseFEM::Element::Quad4::Gauss::nip, "Return number of integration point" );
m.def("xi" , &GooseFEM::Element::Quad4::Gauss::xi , "Return integration point coordinates");
m.def("w" , &GooseFEM::Element::Quad4::Gauss::w , "Return integration point weights" );
}
// -------------------------------------------------------------------------------------------------
void init_ElementQuad4Nodal(py::module &m)
{
m.def("nip", &GooseFEM::Element::Quad4::Nodal::nip, "Return number of integration point" );
m.def("xi" , &GooseFEM::Element::Quad4::Nodal::xi , "Return integration point coordinates");
m.def("w" , &GooseFEM::Element::Quad4::Nodal::w , "Return integration point weights" );
}
// =================================================================================================

Event Timeline

c4science · Help