Kristoffer Carlsson

As a follow-up question, what is the intended way to use the wrappers in the `python` folder? It seems with `pip` you get the ones in `cffi` but those seem...

Thanks for the PR. We don't really have a good way to export a visualization of facesets but as a workaround I created dirichlet boundary conditions on all faces and...

> Regarding your first question, I also tried to figure out which PDE is solved, but I could not figure it our yet. Maybe @KristofferC or @fredrikekre could give a...

I did an MWE benchmark: ```jl u = [1,2,3,4,5,6,7,8,9,10] ue = [0, 0, 0, 0] idxs = [1,3,5,8] function assemble1!(ue, u, idxs) @inbounds for (i, j) in enumerate(idxs) ue[i] =...

> gather is basically a synonym to assemble so that seems confusing To me, `scatter` is more like assemble. You take something small and insert it into something big. This...

> 1h tuning and I have not run the benchmark further than tuning yet. Yikes, that seems a bit much. I think this is perhaps a little bit excessive to...

> Workload is now 8 minutes for tuning and 5 minutes for the actual benchmark That feel pretty ok.

Isn't the problem that there is ever only one `Iter` object in the array (they all refer to the same object) so `iter.x` will always be the same thing.

Oh, I misread the example completely.

Describing the way a global face is defined (`(global_cell_index, local_face_index)`) would also be good. The `local_face_index` is defined by e.g. https://github.com/KristofferC/JuAFEM.jl/blob/8224282ab4d67cb523ef342e4a6ceb1716764ada/src/interpolations.jl#L154 To give a description of it, say that we...