157 std::span<const std::int32_t> entities,
158 std::span<scalar_type> values,
159 std::array<std::size_t, 2> vshape)
const
162 if (mesh.topology()->dim() == _x_ref.second[1])
164 else if (mesh.topology()->dim() == _x_ref.second[1] + 1)
167 throw std::runtime_error(
"Invalid dimension of evaluation points.");
175 auto x_dofmap = mesh.geometry().dofmap();
178 auto& cmap = mesh.geometry().cmap();
180 std::size_t num_dofs_g = cmap.dim();
181 auto x_g = mesh.geometry().x();
184 std::vector<geometry_type> coord_dofs(3 * num_dofs_g);
186 int num_argument_dofs = 1;
187 std::span<const std::uint32_t> cell_info;
188 std::function<void(std::span<scalar_type>, std::span<const std::uint32_t>,
191 = [](std::span<scalar_type>, std::span<const std::uint32_t>,
197 if (_argument_function_space)
200 = _argument_function_space->dofmap()->element_dof_layout().num_dofs();
201 auto element = _argument_function_space->element();
202 num_argument_dofs *= _argument_function_space->dofmap()->bs();
204 if (element->needs_dof_transformations())
206 mesh.topology_mutable()->create_entity_permutations();
207 cell_info = std::span(mesh.topology()->get_cell_permutation_info());
209 = element->template dof_transformation_right_fn<scalar_type>(
215 std::function<
const std::int32_t*(std::span<const std::int32_t>,
218 = []([[maybe_unused]] std::span<const std::int32_t> entities,
219 [[maybe_unused]] std::size_t idx) {
return nullptr; };
223 = [](std::span<const std::int32_t> entities, std::size_t idx)
224 {
return entities.data() + 2 * idx + 1; };
229 std::vector<scalar_type> values_local(size0 * num_argument_dofs, 0);
230 for (std::size_t e = 0; e < entities.size() / estride; ++e)
232 std::int32_t entity = entities[e * estride];
233 auto x_dofs = MDSPAN_IMPL_STANDARD_NAMESPACE::submdspan(
234 x_dofmap, entity, MDSPAN_IMPL_STANDARD_NAMESPACE::full_extent);
235 for (std::size_t i = 0; i < x_dofs.size(); ++i)
237 std::copy_n(std::next(x_g.begin(), 3 * x_dofs[i]), 3,
238 std::next(coord_dofs.begin(), 3 * i));
241 const scalar_type* coeff_cell = coeffs.data() + e * cstride;
242 const int* entity_index = get_entity_index(entities, e);
244 std::ranges::fill(values_local, 0);
245 _fn(values_local.data(), coeff_cell, constant_data.data(),
246 coord_dofs.data(), entity_index,
nullptr);
247 post_dof_transform(values_local, cell_info, e, size0);
248 for (std::size_t j = 0; j < values_local.size(); ++j)
249 values[e * vshape[1] + j] = values_local[j];