Implement new Space time predictor without (variable number of) picard loops
The task is relatively easy: Diff Dumbsers fortran prototype (in repository) and check out the changes in the generic kernels.
From Dumbser, 11. März 2018 um 11:20:
Wäre natürlich super, wenn Du meinen Fortran Code in C übersetzen könntest. Wenn es Probleme gibt, machen wir wieder eine Skype session. Das Format der Schleifen und der nötigen Berechnungen ist quasi dasselbe wie für alle anderen Rechnungen im space-time predictor, d.h. da kann man sehr viel übernehmen. Nur, dass man den Zeitindex nicht mehr mitschleppen muss, sondern nur im Raum arbeiten kann. Der einzige Kernel der geändert werden muss ist der space-time predictor (in 2D und 3D).
Ich würde nur den second und third order initial guess implementieren, siehe den Code in SpaceTimePredictor.f90 unter
#ifdef SECOND_ORDER_INITIAL_GUESS
#ifdef THIRD_ORDER_INITIAL_GUESS
Ich habe mich diese Woche auf Scaling und den Vergleich Runge-Kutta DG / ADER-DG konzentriert, d.h. an dem 2D FO-fCCZ4 habe ich nicht mehr weitergearbeitet. Ich will erst das GRMHD Paper vom Tisch haben.