C1 continuity implementation in dolfinx

[dokken]

@bagla0

1. For a quad mesh of topological dimension:2 (cylinder surface) in 3D space, when I perform scalar assembly:

dx = Measure('dx')(domain=mesh)
dolfinx.fem.assemble_scalar(form(sum([1*dx])))

what quadrature rule it takes to integrate for curved surface. Can I modify the default rule?

You can set the quadrature degree in the metadata of the measure as done in

Total number of gauss point

The number of Gauss-points used depends on the variational form. DOLFINx uses ufl to estimate the quadrature degree of the variational form, see: This is in turn sent to basix which tabulates at the relevant points. You can use basix directly to investigate the default quadrature of any order, see:Creating and using a quadrature rule — Basix 0.8.0.0 documentation You can also by-pass the basix/ufl workflow by supplying your own quadrature rule:

or

Assembly timings for hexahedral elements slow

You are setting the quadrature degree in the wrong place (function space), not in the measure. If you fix this: def pre_hyperelastic_dolfinx(mesh, **kwargs): quadrature_degree = 4 ord = 1 elemu = element("P", mesh.basix_cell(), ord, shape=(3,)) V = functionspace( mesh, elemu, ) u = Function(V, name="u") v = TestFunction(V) du = TrialFunction(V) F = Identity(len(u)) + grad(u) C = F.T * F Ic = tr(C) J = det(F) dx = Measure("dx", domain=…