SyFi 0.3
|
#include <CrouzeixRaviart.h>
Public Member Functions | |
CrouzeixRaviart () | |
CrouzeixRaviart (Polygon &p, unsigned int order=1) | |
virtual | ~CrouzeixRaviart () |
void | compute_basis_functions () |
def | __init__ |
def | compute_basis_functions |
Public Attributes | |
this | |
Static Private Attributes | |
dictionary | __swig_setmethods__ = {} |
tuple | __setattr__ = lambdaself,name,value:_swig_setattr(self, CrouzeixRaviart, name, value) |
dictionary | __swig_getmethods__ = {} |
tuple | __getattr__ = lambdaself,name:_swig_getattr(self, CrouzeixRaviart, name) |
__repr__ = _swig_repr | |
__swig_destroy__ = _SyFi.delete_CrouzeixRaviart | |
__del__ = lambdaself:None; |
Proxy of C++ SyFi::CrouzeixRaviart class
Definition at line 12 of file CrouzeixRaviart.h.
SyFi::CrouzeixRaviart::CrouzeixRaviart | ( | ) |
Definition at line 15 of file CrouzeixRaviart.cpp.
References SyFi::StandardFE::order.
: StandardFE() { order = 1; }
SyFi::CrouzeixRaviart::CrouzeixRaviart | ( | Polygon & | p, |
unsigned int | order = 1 |
||
) |
Definition at line 20 of file CrouzeixRaviart.cpp.
References compute_basis_functions().
: StandardFE(p, order) { compute_basis_functions(); }
virtual SyFi::CrouzeixRaviart::~CrouzeixRaviart | ( | ) | [inline, virtual] |
Definition at line 17 of file CrouzeixRaviart.h.
{}
def SyFi::CrouzeixRaviart::__init__ | ( | self, | |
args | |||
) |
__init__(self) -> CrouzeixRaviart __init__(self, Polygon p, unsigned int order = 1) -> CrouzeixRaviart __init__(self, Polygon p) -> CrouzeixRaviart
Reimplemented from SyFi::StandardFE.
Definition at line 2275 of file SyFi.py.
02275 : 02276 """ 02277 __init__(self) -> CrouzeixRaviart 02278 __init__(self, Polygon p, unsigned int order = 1) -> CrouzeixRaviart 02279 __init__(self, Polygon p) -> CrouzeixRaviart 02280 """ 02281 this = _SyFi.new_CrouzeixRaviart(*args) 02282 try: self.this.append(this) 02283 except: self.this = this
void SyFi::CrouzeixRaviart::compute_basis_functions | ( | ) | [virtual] |
Reimplemented from SyFi::StandardFE.
Definition at line 25 of file CrouzeixRaviart.cpp.
References SyFi::bernstein(), SyFi::StandardFE::description, SyFi::dirac(), SyFi::StandardFE::dofs, SyFi::Line::integrate(), SyFi::Triangle::line(), SyFi::StandardFE::Ns, SyFi::StandardFE::order, SyFi::StandardFE::p, SyFi::Polygon::str(), SyFi::sub(), SyFi::t, demos::geom_test::tetrahedron, SyFi::Tetrahedron::triangle(), demos::crouzeixraviart::triangle, SyFi::Polygon::vertex(), SyFi::x, SyFi::y, and SyFi::z.
Referenced by check_CrouzeixRaviart(), SyFi::VectorCrouzeixRaviart::compute_basis_functions(), and CrouzeixRaviart().
{ // remove previously computed basis functions and dofs Ns.clear(); dofs.clear(); if ( p == NULL ) { throw(std::logic_error("You need to set a polygon before the basisfunctions can be computed")); } if (order != 1) { throw(std::logic_error("Only Crouziex-Raviart elements of order 1 is possible")); } // see e.g. Brezzi and Fortin book page 116 for the definition if ( p->str().find("ReferenceLine") != string::npos ) { cout <<"Can not define the Raviart-Thomas element on a line"<<endl; } else if ( p->str().find("Triangle") != string::npos ) { description = "CrouzeixRaviart_2D"; Triangle& triangle = (Triangle&)(*p); // create the polynomial space GiNaC::ex polynom_space = bernstein(1, triangle, "a"); GiNaC::ex polynom = polynom_space.op(0); GiNaC::lst variables = GiNaC::ex_to<GiNaC::lst>(polynom_space.op(1)); GiNaC::ex basis = polynom_space.op(2); // create the dofs GiNaC::symbol t("t"); for (int j=0; j< 3; j++) { // solve the linear system to compute // each of the basis functions GiNaC::lst equations; for (int i=0; i< 3; i++) { Line line = triangle.line(i); // GiNaC::ex dofi = line.integrate(polynom); GiNaC::lst midpoint = GiNaC::lst( (line.vertex(0).op(0) + line.vertex(1).op(0))/2, (line.vertex(0).op(1) + line.vertex(1).op(1))/2); dofs.insert(dofs.end(), midpoint); // GiNaC::ex dofi = polynom.subs( x == midpoint.op(0)).subs( y == midpoint.op(1)); GiNaC::ex dofi = line.integrate(polynom); equations.append( dofi == dirac(i,j)); if (j == 1) { // GiNaC::lst d = GiNaC::lst(line.vertex(0) , line.vertex(1)); dofs.insert(dofs.end(), midpoint); } } GiNaC::ex sub = lsolve(equations, variables); GiNaC::ex Ni = polynom.subs(sub); Ns.insert(Ns.end(),Ni); } } else if ( p->str().find("Tetrahedron") != string::npos ) { description = "CrouzeixRaviart_3D"; Tetrahedron& tetrahedron = (Tetrahedron&)(*p); GiNaC::ex polynom_space = bernstein(1, tetrahedron, "a"); GiNaC::ex polynom = polynom_space.op(0); GiNaC::lst variables = GiNaC::ex_to<GiNaC::lst>(polynom_space.op(1)); GiNaC::ex basis = polynom_space.op(2); GiNaC::ex bernstein_pol; GiNaC::symbol t("t"); // dofs related to edges for (int j=0; j< 4; j++) { GiNaC::lst equations; for (int i=0; i< 4; i++) { Triangle triangle = tetrahedron.triangle(i); GiNaC::lst midpoint = GiNaC::lst( (triangle.vertex(0).op(0) + triangle.vertex(1).op(0) + triangle.vertex(2).op(0))/3, (triangle.vertex(0).op(1) + triangle.vertex(1).op(1) + triangle.vertex(2).op(1))/3, (triangle.vertex(0).op(2) + triangle.vertex(1).op(2) + triangle.vertex(2).op(2))/3 ); GiNaC::ex dofi = polynom.subs(x == midpoint.op(0)).subs(y == midpoint.op(1)).subs(z == midpoint.op(2)); equations.append( dofi == dirac(i,j)); if ( j == 1 ) { // GiNaC::lst d = GiNaC::lst(triangle.vertex(0), triangle.vertex(1), triangle.vertex(2)); dofs.insert(dofs.end(), midpoint); } } GiNaC::ex sub = lsolve(equations, variables); GiNaC::ex Ni = polynom.subs(sub); Ns.insert(Ns.end(),Ni); } } }
def SyFi::CrouzeixRaviart::compute_basis_functions | ( | self | ) |
compute_basis_functions(self)
Reimplemented from SyFi::StandardFE.
SyFi::CrouzeixRaviart::__del__ = lambdaself:None; [static, private] |
Reimplemented from SyFi::StandardFE.
tuple SyFi::CrouzeixRaviart::__getattr__ = lambdaself,name:_swig_getattr(self, CrouzeixRaviart, name) [static, private] |
Reimplemented from SyFi::StandardFE.
SyFi::CrouzeixRaviart::__repr__ = _swig_repr [static, private] |
Reimplemented from SyFi::StandardFE.
tuple SyFi::CrouzeixRaviart::__setattr__ = lambdaself,name,value:_swig_setattr(self, CrouzeixRaviart, name, value) [static, private] |
Reimplemented from SyFi::StandardFE.
SyFi::CrouzeixRaviart::__swig_destroy__ = _SyFi.delete_CrouzeixRaviart [static, private] |
Reimplemented from SyFi::StandardFE.
dictionary SyFi::CrouzeixRaviart::__swig_getmethods__ = {} [static, private] |
Reimplemented from SyFi::StandardFE.
dictionary SyFi::CrouzeixRaviart::__swig_setmethods__ = {} [static, private] |
Reimplemented from SyFi::StandardFE.
Reimplemented from SyFi::StandardFE.