Source code for sfepy.discrete.fem.fields_positive

from sfepy.discrete.fem.fields_base import FEField, H1Mixin
from sfepy.discrete.fem.fields_nodal import GlobalNodalLikeBasis

[docs]class H1BernsteinVolumeField(H1Mixin, GlobalNodalLikeBasis, FEField): """ Bernstein basis approximation with positive-only basis function values. """ family_name = 'volume_H1_bernstein'
[docs] def create_basis_context(self): """ Create the context required for evaluating the field basis. """ # Hack for tests to pass - the reference coordinates are determined # from vertices only - we can use the Lagrange basis context for the # moment. The true context for Field.evaluate_at() is not implemented. gps = self.gel.poly_space mesh = self.create_mesh(extra_nodes=False) ctx = geo_ctx = gps.create_context(self.cmesh, 0, 1e-15, 100, 1e-8) ctx.geo_ctx = geo_ctx return ctx
[docs]class H1BernsteinSurfaceField(H1Mixin, GlobalNodalLikeBasis, FEField): family_name = 'surface_H1_bernstein'