Function assemble_vector1

cpp/demo/custom_kernel/main.cpp:164–176

/// @return l2 norm squared of the vector.
template <std::floating_point T>
double assemble_vector1(const mesh::Geometry<T>& g, const fem::DofMap& dofmap,
                        auto kernel, const std::vector<std::int32_t>& cells)
{
  la::Vector<T> b(dofmap.index_map, 1);
  md::mdspan<const T, md::extents<std::size_t, md::dynamic_extent, 3>> x(
      g.x().data(), g.x().size() / 3, 3);
  common::Timer timer("Assembler1 lambda (vector)");
  fem::impl::assemble_cells<1>([](auto, auto, auto, auto) {}, b.array(),
                               g.dofmaps().front(), x, cells,
                               {dofmap.map(), 1, cells}, kernel, {}, {}, {});
  b.scatter_rev(std::plus<T>());
  return la::squared_norm(b);
}

/// @brief Assemble P1 mass matrix and a RHS vector using element kernel
/// approaches.