Wolfgang Bangerth

Documentation: Appropriately render other languages

1

Following #4665, we have places where we ask for mark-up of gnuplot, ksh, and prmfiles. sphinx likely doesn't understand these all yet, see the list at https://pygments.org/languages/, though gnuplot is...

Confusion about deviators

2

In the documentation, we explain how when we deal with compressible strains, we always compute ``` eps - 1/3 trace(eps)*I ``` even in for 2d models where one might think...

Limit the growth of the time step size

1

@tjhei dug up a masters thesis at https://dc.uwm.edu/cgi/viewcontent.cgi?article=1405&context=etd that is about time step choices for BDF2. Among the conclusions there is that one should not let the time step size...

We currently compute the time step based on a CFL number that defines `h` as the shortest edge length of a cell: https://github.com/geodynamics/aspect/blob/main/source/time_stepping/convection_time_step.cc#L48-L75 That seemed natural to FEM people like...

Better document the 'mr' parameter to 'metrop_select'.

1

I found myself confused why `metrop_select()` has this comparison: ``` if np.isfinite(mr) and np.log(uniform()) < mr: ``` If `mr` was just a ratio of probabilities, as one would typically see...

nvgetcommunicator returning pointer to MPI_Comm leads to complicated designs.

4

For the deal.II wrappers of the N_Vector interface, we've had several lengthy discussions about how to deal with the `ops.nvgetcommunicator` callback. The basic problem is that that interface requires one...

This is a necessary follow-up to https://github.com/dealii/dealii/pull/16543: As mentioned in https://github.com/dealii/dealii/pull/16543#discussion_r1468146758, we can't use `std::abort()` on devices. Instead, we must use `Kokkos::abort("text");`. For reference, `Assert` is currently defined in the...

[WIP] Add Tpetra tests.

22

This is definitely a WIP patch: 284 tests for the Tpetra wrappers, or which currently 41 succeed. I created these tests by cloning the existing `tests/trilinos/` directory and replacing `TrilinosWrappers::`...

MG two-level/matrix-free infrastructure and std::complex data types.

3

I tried to instantiate the matrix-free stuff with `std::complex` vectors, and utterly failed :-) Some of the minor edits necessary are captured in #16717, but the major blocker is that...

Import BOOST 1.84.

4

First draft. Probably many things still wrong with this. In reference to #16618. Fixes https://github.com/dealii/dealii/issues/16662.