def __init__(self, cell, degree):
assert degree == 2, "Only defined for degree 2"
Ps = ONSymTensorPolynomialSet(cell, degree)
Ls = ArnoldWintherNCDual(cell, degree)
mapping = "double contravariant piola"
super(ArnoldWintherNC, self).__init__(Ps, Ls, degree, mapping=mapping)
def __init__(self, cell, degree):
assert degree >= 0, "Regge start at degree 0!"
# shape functions
Ps = ONSymTensorPolynomialSet(cell, degree)
# degrees of freedom
Ls = ReggeDual(cell, degree)
# mapping under affine transformation
mapping = "double covariant piola"
super(Regge, self).__init__(Ps, Ls, degree, mapping=mapping)
def __init__(self, cell, degree):
assert degree >= 0, "Hellan-Herrmann-Johnson starts at degree 0!"
# shape functions
Ps = ONSymTensorPolynomialSet(cell, degree)
# degrees of freedom
Ls = HellanHerrmannJohnsonDual(cell, degree)
# mapping under affine transformation
mapping = "double contravariant piola"
super(HellanHerrmannJohnson, self).__init__(Ps, Ls, degree,
mapping=mapping)