GEARRT: conservative flux exchanges

This MR does the same as !1411 (merged) did for Gizmo hydrodynamics: Change the RT interactions and flux exchanges to conservative ones.

Some additional attention was necessary for the RT sorts due to the sub-cycling.

This passes my tests, and I ran a couple of thousand (hundred) steps of EAGLE_12 (EAGLE_25) successfully (until something went wrong with hydro, which is a different problem.) Feel free to stress test this too if you so desire. There is a small catch though: The way to "chose" whether an RT scheme runs with symmetric MPI comms or not is a bit hidden. The MPI_SYMMETRIC_FORCE_INTERACTION macro is defined in src/part.h only when a flavor of Gizmo is selected as the hydrodynamics method. This remains unchanged: The RT will also only use symmetric MPI comms if a gizmo hydrodynamics flavor is selected.

As such, any RT scheme that runs on top of Gizmo hydrodynamics will also have symmetric flux exchanges. And vice versa: Any RT scheme running on top of SPH flavors will have SPH-like "force" loops.

@yuyttenhove , do you perhaps want to have a look at this MR too before the merge?

src/engine_maketasks.c

@@ -2777,10 +2777,12 @@ void engine_make_extra_hydroloop_tasks_mapper(void *map_data, int num_elements,
          * the interaction in the transport loop. Local cells don't have an
          * rt_sorts task. */
         if (ci->hydro.super->hydro.sorts != NULL)
-          scheduler_addunlock(sched, ci->hydro.sorts, t_rt_transport);
+          scheduler_addunlock(sched, ci->hydro.super->hydro.sorts,
+                              t_rt_transport);
         if ((ci->hydro.super != cj->hydro.super) &&
             (cj->hydro.super->hydro.sorts != NULL))
-          scheduler_addunlock(sched, cj->hydro.sorts, t_rt_transport);
+          scheduler_addunlock(sched, cj->hydro.super->hydro.sorts,
+                              t_rt_transport);
 #endif
       }
 
@@ -3603,10 +3605,12 @@ void engine_make_extra_hydroloop_tasks_mapper(void *map_data, int num_elements,
          * the interaction in the transport loop. Local cells don't have
          * an rt_sort task. */
         if (ci->hydro.super->hydro.sorts != NULL)
-          scheduler_addunlock(sched, ci->hydro.sorts, t_rt_transport);
+          scheduler_addunlock(sched, ci->hydro.super->hydro.sorts,
+                              t_rt_transport);
         if ((ci->hydro.super != cj->hydro.super) &&
             (cj->hydro.super->hydro.sorts != NULL))
-          scheduler_addunlock(sched, cj->hydro.sorts, t_rt_transport);
+          scheduler_addunlock(sched, cj->hydro.super->hydro.sorts,
+                              t_rt_transport);
 #endif
       }