|
SyFi
0.3
|
#include <ArnoldFalkWintherWeakSym.h>
Public Member Functions | |
| ArnoldFalkWintherWeakSymSigma () | |
| ArnoldFalkWintherWeakSymSigma (Polygon &p, int order=1) | |
| virtual | ~ArnoldFalkWintherWeakSymSigma () |
| virtual void | compute_basis_functions () |
Definition at line 26 of file ArnoldFalkWintherWeakSym.h.
Definition at line 32 of file ArnoldFalkWintherWeakSym.cpp.
References SyFi::StandardFE::description.
: StandardFE() { description = "ArnoldFalkWintherWeakSymSigma"; }
| SyFi::ArnoldFalkWintherWeakSymSigma::ArnoldFalkWintherWeakSymSigma | ( | Polygon & | p, |
| int | order = 1 |
||
| ) |
Definition at line 37 of file ArnoldFalkWintherWeakSym.cpp.
References compute_basis_functions().
: StandardFE(p, order) { compute_basis_functions(); }
| virtual SyFi::ArnoldFalkWintherWeakSymSigma::~ArnoldFalkWintherWeakSymSigma | ( | ) | [inline, virtual] |
Definition at line 31 of file ArnoldFalkWintherWeakSym.h.
{}
| void SyFi::ArnoldFalkWintherWeakSymSigma::compute_basis_functions | ( | ) | [virtual] |
Reimplemented from SyFi::StandardFE.
Definition at line 42 of file ArnoldFalkWintherWeakSym.cpp.
References SyFi::Nedelec2Hdiv::compute_basis_functions(), SyFi::StandardFE::description, SyFi::StandardFE::dof(), SyFi::StandardFE::dofs, SyFi::istr(), SyFi::StandardFE::N(), SyFi::StandardFE::nbf(), SyFi::StandardFE::Ns, SyFi::StandardFE::order, SyFi::StandardFE::p, SyFi::StandardFE::set_order(), and SyFi::StandardFE::set_polygon().
Referenced by ArnoldFalkWintherWeakSymSigma(), and main().
{
// remove previously computed basis functions and dofs
Ns.clear();
dofs.clear();
if ( order < 1 )
{
throw(std::logic_error("Arnold-Falk-Winther elements must be of order 1 or higher."));
}
if ( p == NULL )
{
throw(std::logic_error("You need to set a polygon before the basisfunctions can be computed"));
}
description = istr("ArnoldFalkWintherWeakSymSigma_", order) + "_3D";
Nedelec2Hdiv fe;
fe.set_order(order);
fe.set_polygon(*p);
fe.compute_basis_functions();
for (int d=0; d<3; d++)
{
for (unsigned int i=0; i<fe.nbf(); i++)
{
GiNaC::matrix Nmat = GiNaC::matrix(3,3);
Nmat(d,0) = fe.N(i).op(0);
Nmat(d,1) = fe.N(i).op(1);
Nmat(d,2) = fe.N(i).op(2);
Ns.insert(Ns.end(), Nmat);
dofs.insert(dofs.end(),GiNaC::lst(fe.dof(i).op(0), fe.dof(i).op(1), d));
}
}
}