ExaHyPE issueshttps://gitlab.lrz.de/groups/exahype/-/issues2019-03-07T19:41:06+01:00https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/221Deeply check the public repository for CCZ42019-03-07T19:41:06+01:00Ghost UserDeeply check the public repository for CCZ4It is at https://github.com/exahype/exahype and we don't want to have the CCZ4 system in the code or any old commits. Make this sure by inspecting the code.
Cannot do it now since the repo is 70MB in size and I'm in a train with a bad w...It is at https://github.com/exahype/exahype and we don't want to have the CCZ4 system in the code or any old commits. Make this sure by inspecting the code.
Cannot do it now since the repo is 70MB in size and I'm in a train with a bad wifi.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/111Delete Logo.h from repository and include minimal PNG decoding library2018-06-15T15:13:43+02:00Ghost UserDelete Logo.h from repository and include minimal PNG decoding libraryThe ExaHyPE repository has more than 2000 commits and weights only 10MB. Now in January somebody added a **10MB C++ header file** containing an uncompressed picture. It's not only possible to read in this picture in any image reader any...The ExaHyPE repository has more than 2000 commits and weights only 10MB. Now in January somebody added a **10MB C++ header file** containing an uncompressed picture. It's not only possible to read in this picture in any image reader any more, but the repository just blew to the double size and nobody really wants to deal with such files (try to load it in a text editor).
First, the commit introducing this picture should be changed and the monster file deleted out of the repository.
Second, it should be replaced by an actual picture in GIF, JPG, PNG or similar efficient formats (will be less than **20kb**!) and a suitable decompression routine should be added to the Demonstrator. This mimics how actual offline initial data generators work. We do _not_ need to introduce a dependency on an external library but can embed an Open Source minimal decoding library. There are a lot of them, for instance one of these:
* https://github.com/nothings/stb/blob/master/stb_image.h - single file header only PNG decompressor (*public domain*, 200kB code)
* https://github.com/elanthis/upng - another micro PNG decompressor single file library (*as-is* license)
* http://lodev.org/lodepng/ - another pico PNG decomressor, 500lines (*as-is* license)
* https://github.com/hidefromkgb/gif_load/blob/master/gif_load.h - a 400 lines GIF reader (*public domain*, 20kB code)
This is really the way to go as it allows users to exchange the picture with another one _at runtime_.
I @sven volunteer to implement this.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/123Derivative plotter, part 22018-06-15T15:13:43+02:00Ghost UserDerivative plotter, part 2We wrote a plotter which access one cell ("patch") and computes the local derivatives, but it cannot plot to VTK but only reductions. We need to extend this to VTK, too.
Application context: Z4, we need to see the constraints.
This is ...We wrote a plotter which access one cell ("patch") and computes the local derivatives, but it cannot plot to VTK but only reductions. We need to extend this to VTK, too.
Application context: Z4, we need to see the constraints.
This is something I can try to do.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/70Discussion topics for a discussion ("Coding") day at Oct 23., 2016 in Munich2018-06-15T15:13:43+02:00Ghost UserDiscussion topics for a discussion ("Coding") day at Oct 23., 2016 in Munich* Kernel situation
* Trento will not be able to implement new numerics into C++ "nested loop" kernels
* We invested two weeks of finding bugs porting Fortran Z4 to ExaHyPE -- this is a dead end!
* We either need modern Fortran kern...* Kernel situation
* Trento will not be able to implement new numerics into C++ "nested loop" kernels
* We invested two weeks of finding bugs porting Fortran Z4 to ExaHyPE -- this is a dead end!
* We either need modern Fortran kernels or a good C++ linear algebra library (there are plenty of them)
* (Generic) research kernels might be in a different shape than optimized kernels!
* What is ExaHyPE?
* If it shall be a standalone *application*, then it should provide the best support for parameters, etc. as possible
* If it shall serve as a "utility", end users have to come up with individual solutions on theirselves. Both is fine, but there should be a common idea everybody is aware of.
* Workflow
* Durham cannot provide implementation support (code monkeying) for ExaHyPE applications in this intensity anymore
* Frankfurt cannot provide support for running benchmarks and (parameter) test in this intensity anymore
* Fortran situation in Applications
* MHD was converted to C++ with surprising results in stability and speed
* What's the joint aim for the team? There is more and more Fortran code (Applications get bigger), cf. the `Z4` application
* API and Code Generation
* Kernel calls should be method calls instead of templated static function calls
* Toolkit only needs to provide initial "scaffold" code. In 90% of applications, generated code is now in the repo, edited after generation and gets overwritten regularly.
* Toolkit should not distribute constants in the generated code but collect them in a unique file (`GeneratedConstants.h`). Paradigm: Generate maintainable code.
* Specification file format
* There is a lot to say here. We should aggree on a whishlist of features
* We could also agree to keep the spec files small and let the problem be solved by the users
* Treatment of external libraries and the build system
* Today is already tomorrow -- I have secretly introduced dependencies on GSL, BLAS due to the Astrophysics Initial Data codes.
* A decision should also affect the treatment of Peano.
* IO
* VTK Fileformat (alternatives)
* Does it make sense to define a tailored ExaHyPE output format ourselves?
* Efficient batch processing 3D VTK to movie plotting
* Unit tests
* Mathematically/Physically, there is little point in testing the numerical scheme with random data at every ExyHyPE `DEBUG or Asserts` run.
* Instead, simple but well understood tests with analytic solutions should be run (ie. via Jenkins). We should come up with a list of these tests (ie: Advection, ShuVortex, MHD Wave, Gauge Wave, etc.)
*Please just add points to this list as you whish*
# Assignments in Discussion at 25. Oct 2016
At Oct 25, in Garching there were Tobias, Vasco, Dominik, Alejandro and Sven. We discussed these points and even more and concluded these assignments:
* Fortran support (VV)
* Call linear kernels with temp arrays (DC)
* Non-conservative parts (SK) - done, cf. https://gitlab.lrz.de/gi26det/ExaHyPE/commit/1723512d16acbf118f43465d78d252cda55e9bc6 and similar
* Split of the repositories (SK) - tracked at gi26det/ExaHyPE#55
* Start the release page & remove binaries from GIT (VV) - as followup for gi26det/ExaHyPE#55, right?
* Find new Jour Fix time (incl. Dominic and Trento) (VV)
* Nightly convergence tests (SK+VV) - tracked at gi26det/ExaHyPE#74
* Period BCs: Organise small coding week in Durham (VV+TW)
* Replace static function + templates with std::function (SK+TW) - tracked at gi26det/ExaHyPE#75
* Extend Solver::solutionAdjustment to handle point sources (DC+VV)https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/159DMP is illegally called for ordinary ADERDG-Solver2017-11-02T18:15:15+01:00Ghost UserDMP is illegally called for ordinary ADERDG-SolverThis are actually two problems:
1. For an ordinary ADERDG-Solver (not Limiter) where the variable `dmp-observables` is not given, the default value is **not** `0` but instead just a random number (NaN or MaX or whatever for int). ...This are actually two problems:
1. For an ordinary ADERDG-Solver (not Limiter) where the variable `dmp-observables` is not given, the default value is **not** `0` but instead just a random number (NaN or MaX or whatever for int). This is very bad but solvable in the Parser for me.
2. The method `mapDiscreteMaximumPrincipleObservables` in the abstract ADERDG solver is called. This should not happen at all! Why does it happen?
![dmp](/uploads/fd5837af72fe6f92e871426712399023/dmp.png)https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/69Do memory precision tests2019-03-07T19:23:10+01:00Ghost UserDo memory precision testsTobias schreibt:
>
ich braeuchte mal Deine Hilfe (ja, natuerlich schnell, weil das in die Praesentation rein soll, aber das haste eh vermutet ;-). Die neue ExaHyPE-Version von gestern Nacht kann mit Zahldarstellungen <IEEE double Standa...Tobias schreibt:
>
ich braeuchte mal Deine Hilfe (ja, natuerlich schnell, weil das in die Praesentation rein soll, aber das haste eh vermutet ;-). Die neue ExaHyPE-Version von gestern Nacht kann mit Zahldarstellungen <IEEE double Standard arbeiten. Das ist teuer (erste Schritte suggerieren einen Faktor 5 in der Laufzeit, aber das koennte auch am Chip liegen, d.h. ich muss die Russenmaschine testen), aber es reduziert halt den Speicherbedarf einer Simulation etwas. Ich bekomm 10% schon fuer den einfachen Euler 2d mit Ordnung 3 - hoffe aber, dass das bei anderen Setups deutlich signifikanter ist. Einschalten laesst sich das Feature, indem man
>
double-compression = 0.000001
>
auf etwas zwischen 0 und 1 setzt. Das andere Flag spawn-double-compression-as-background-thread kannste vergessen, das ist meine Baustelle. Meine Frage/Bitte ist nun: Hast Du einen sinnvollen Benchmark zur Hand aus der Astrophysik und kannst Du mal
double-compression = 0.0
double-compression = 0.0001
double-compression = 0.000001
double-compression = 0.000000000001
>
damit ausprobieren und mir sagen, ob Werte grosser 0 die Loesung qualitativ verschlechtern? Ich hab keine Ahnung, was Ihr Euch typischerweise anseht (Ankunftszeiten/Amplituden/...?) deswegen brauch ich da ein qualifiziertes Auge. Wenn Du parallel noch drauf schaust, was er als memoryUsage ausgibt (macht er jetzt automatisch, wenn man nicht mit MPI uebersetzt), dann wuerde mir das sehr helfen, da eine Einschaetzung zu bekommen.
>
Nachtrag: So saehe dann eine Optimisation Sectino aus
```
optimisation
fuse-algorithmic-steps = on
fuse-algorithmic-steps-factor = 0.99
timestep-batch-factor = 0.0
skip-reduction-in-batched-time-steps = on
disable-amr-if-grid-has-been-stationary-in-previous-iteration = off
double-compression = 0.000001
spawn-double-compression-as-background-thread = off
end optimisation
```
@svenk: Machen.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/74Do nightly convergence tests2018-06-15T15:13:43+02:00Ghost UserDo nightly convergence testsAs an initiative for proving and monitoring the correctness of ExaHyPE, we should come up with a series of very simple tests which demonstrate how we debug the code and also give convergence rates.
This is something Sven and Vasco can d...As an initiative for proving and monitoring the correctness of ExaHyPE, we should come up with a series of very simple tests which demonstrate how we debug the code and also give convergence rates.
This is something Sven and Vasco can do.
A problem is that almost all tests require periodic BC, thought.
## Collection of problems
We could solve:
* Advection Equations
* With the Initial Data:
* Low order polynomials (`Q(:) = 1 + v0 * x(0) - v0 * t`). Intermediate solver steps are always analytically known
* No polynomials (`Q(:) = ICA * SIN(2*pi*(x(0) - t))` et al)
* For random matter distributions
* Implemented as conservative equation (`F(Q) = Q`)
* Implemented as nonconservative equation (`F=0, BgradQ = (1,0,0), S=0`)
* EulerFlow, SRHD, MHD: See the [List of Benchmarks](https://gitlab.lrz.de/gi26det/ExaHyPE/wikis/list%20of%20benchmarks)
See also: gi26det/ExaHyPE#73 and gi26det/ExaHyPE#64.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/100Do we see SRMHD convergence with reduced CFL Factor?2018-06-15T15:13:43+02:00Ghost UserDo we see SRMHD convergence with reduced CFL Factor?1. Diskussionpunkt: Zeitschritte bei ExaHyPE scheinen manchmal zu groß gewählt. Dies haben Dominic und ich vor Weihnachten bei SRMHD/AlfenWave gelöst mit pure-ADERDG beobachtet. Wenn wir den CFL_FACTOR = 0.9 auf 0.4 reduzieren, verschwin...1. Diskussionpunkt: Zeitschritte bei ExaHyPE scheinen manchmal zu groß gewählt. Dies haben Dominic und ich vor Weihnachten bei SRMHD/AlfenWave gelöst mit pure-ADERDG beobachtet. Wenn wir den CFL_FACTOR = 0.9 auf 0.4 reduzieren, verschwinden die lokalen Oszillationen in der Simulation.
1a) Sehen wir dann auch (endlich) Konvergenz? Muss ich endlich nachprüfen.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/112Documentation for user interfaces (generic kernels)2018-06-15T15:13:43+02:00Ghost UserDocumentation for user interfaces (generic kernels)https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/119DummyK2018-06-15T15:13:43+02:00Ghost UserDummyKQuestions:
* Do we consider point sources here?
* If so, why do we need an additional isDummyK() boolean getter here?
* If so, why is a point source not related to normal sources and needs its extra memory (see mapping Prediction?)
...Questions:
* Do we consider point sources here?
* If so, why do we need an additional isDummyK() boolean getter here?
* If so, why is a point source not related to normal sources and needs its extra memory (see mapping Prediction?)
We already allocate memory for the source terms in the temporary variables for space-time and spatial fluxes.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/147Dynamic AMR crashes if we switch to gloabl recomputation branch.2017-11-02T18:15:16+01:00Ghost UserDynamic AMR crashes if we switch to gloabl recomputation branch.[webmize](/uploads/db32921836ea432a6ce41060cfb87f97/webmize)Dynamic AMR crashes everytime the solver switches to the global
recomputation branch.
Potential reason:
* ~~We reuse the ADERDGTimeStep adapter which does update the limite...[webmize](/uploads/db32921836ea432a6ce41060cfb87f97/webmize)Dynamic AMR crashes everytime the solver switches to the global
recomputation branch.
Potential reason:
* ~~We reuse the ADERDGTimeStep adapter which does update the limiter
status again. This should not happen.~~
* ~~Erasing is triggered before the limiter status spreading has finished.
Add inertia for erasing requests.~~This is now realised by also taking the previous limiter status into account.
Problems seem to be solved:
![sod_shock_tube](/uploads/2bdedd63eadf6eb0bde5d093e34dbf88/sod_shock_tube.webm)https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/248Dynamic AMR Infinite Loop Scenario2019-09-20T15:46:24+02:00Ghost UserDynamic AMR Infinite Loop ScenarioFound a situation where Ancestors are not able to erase their children and the mesh refinement ends up
in an infinite loop.
A few affected cells
--------------------
```
47.8326 info (solverNumber:0,neighbourMergePerform...Found a situation where Ancestors are not able to erase their children and the mesh refinement ends up
in an infinite loop.
A few affected cells
--------------------
```
47.8326 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:41,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.481481],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,0,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8328 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:41,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.222222,0.481481],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,0,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.833 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:41,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.222222,0.518519],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,0,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8332 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:41,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.518519],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,0,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.857 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:36,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.222222,0.407407],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,1,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8572 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:36,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.407407],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,0,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8575 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:36,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.37037],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[1,0,1,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8577 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:36,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.222222,0.37037],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,3,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[0,0,2,0],communicationStatus:1,facewiseRefinementStatus:[-2,-2,-1,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.858 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:36,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.333333],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[2,0,1,0],communicationStatus:1,facewiseRefinementStatus:[-1,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8736 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:9,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.37037,0.222222],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,3,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[2,0,0,0],communicationStatus:1,facewiseRefinementStatus:[-1,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8739 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:9,parentCellLevel:-1,parentOffset:[-1,-1],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.333333,0.259259],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,3,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[1,0,2,0],communicationStatus:1,facewiseRefinementStatus:[-2,-2,-1,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8747 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:62,parentCellLevel:3,parentOffset:[0.222222,0.222222],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.296296,0.296296],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[4,4,4,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[1,0,1,0],communicationStatus:0,facewiseRefinementStatus:[-2,-2,-2,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8748 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:62,parentCellLevel:3,parentOffset:[0.222222,0.222222],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.259259,0.296296],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,3,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[2,0,2,1],communicationStatus:1,facewiseRefinementStatus:[-1,-2,-1,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
47.8752 info (solverNumber:0,neighbourMergePerformed:[0,0,0,0],hasCompletedTimeStep:1,parentIndex:62,parentCellLevel:3,parentOffset:[0.222222,0.222222],hasVirtualChildren:0,type:Ancestor,refinementEvent:ErasingChildrenRequested,level:4,offset:[0.296296,0.259259],size:[0.037037,0.037037],previousCorrectorTimeStamp:1.79769e+308,previousCorrectorTimeStepSize:0,correctorTimeStepSize:0,correctorTimeStamp:0.140328,predictorTimeStepSize:0,predictorTimeStamp:0.140328,solution:-1,solutionAverages:-1,solutionCompressed:-1,previousSolution:-1,previousSolutionAverages:-1,previousSolutionCompressed:-1,update:-1,updateAverages:-1,updateCompressed:-1,extrapolatedPredictor:-1,extrapolatedPredictorAverages:-1,extrapolatedPredictorCompressed:-1,fluctuation:-1,fluctuationAverages:-1,fluctuationCompressed:-1,solutionMin:-1,solutionMax:-1,facewiseAugmentationStatus:[3,4,3,4],augmentationStatus:4,previousAugmentationStatus:4,facewiseCommunicationStatus:[2,1,2,0],communicationStatus:1,facewiseRefinementStatus:[-1,-2,-1,-2],refinementFlag:0,refinementStatus:-2,previousRefinementStatus:-1,iterationsToCureTroubledCell:0,compressionState:Uncompressed,bytesPerDoFInPreviousSolution:-1,bytesPerDoFInSolution:-1,bytesPerDoFInUpdate:-1,bytesPerDoFInExtrapolatedPredictor:-1,bytesPerDoFInFluctuation:-1)
```
Spec file for reproducing the error
-----------------------------------
```
exahype-project Euler
peano-kernel-path const = ./Peano
exahype-path const = ./ExaHyPE
output-directory const = ./ApplicationExamples/Euler/Euler_ADERDG
architecture const = noarch
computational-domain
dimension const = 2
width = 1.0, 1.0, 1.0
offset = 0.0, 0.0, 0.0
end-time = 1.001
end computational-domain
shared-memory
identifier = dummy
configure = {background-tasks:4}
cores = 8
properties-file = sharedmemory.properties
end shared-memory
distributed-memory
identifier = static_load_balancing
configure = {hotspot,fair,ranks_per_node:4}
buffer-size = 64
timeout = 240
end distributed-memory
global-optimisation
fuse-algorithmic-steps = on
fuse-algorithmic-steps-factor = 0.99
spawn-predictor-as-background-thread = on
spawn-amr-background-threads = off
/* 0.0 und 0.8 sind schon mal zwei Faktoren */
disable-vertex-exchange-in-time-steps = on
time-step-batch-factor = 1.0
disable-metadata-exchange-in-batched-time-steps = on
double-compression = 0.0
spawn-double-compression-as-background-thread = off
end global-optimisation
solver ADER-DG EulerSolver_ADERDG
variables const = rho:1,j:3,E:1
order const = 3
/* 27 points: 0.05, 9 points: 0.15 */
maximum-mesh-size = 0.4
maximum-mesh-depth = 3
halo-cells = 1
time-stepping = globalfixed
type const = nonlinear
terms const = flux
optimisation const = generic,usestack
language const = C
constants = reference:entropywave
plot vtu::Cartesian::vertices::ascii ErrorPlotter
// abserrorl1percell[nvar],q[nvar],qanalytical[nvar]
variables const = 15
time = 0.0
repeat = 1e-2
output = ./errors
end plot
end solver
end exahype-project
```https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/128Dynamic AMR race conditions possible2018-06-15T15:13:43+02:00Ghost UserDynamic AMR race conditions possibleGradient based refinement criteria can introduce a race condition
in my current dynamic AMR implementation.
The user has to make sure here that the fine grid cells are not requesting
erasing while the coarse grid cell requests refi...Gradient based refinement criteria can introduce a race condition
in my current dynamic AMR implementation.
The user has to make sure here that the fine grid cells are not requesting
erasing while the coarse grid cell requests refinement.
I properly have to introduce a flag/cell type that notifies the fine grid cells that the
parent has requested refinement during this grid setup.
I was able to implement stable refinement criteria solely based on the magnitude of the target quantity.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/152Dynamic plotter registration2019-03-07T19:39:00+01:00Ghost UserDynamic plotter registrationUsers shall be able to add, comment out or remove plotters at runtime without recompiling. The plotter ordering should not be fixed.
This is not too hard to obtain, just requires changes at `KernelCalls.cpp` with having a generated plot...Users shall be able to add, comment out or remove plotters at runtime without recompiling. The plotter ordering should not be fixed.
This is not too hard to obtain, just requires changes at `KernelCalls.cpp` with having a generated plotter registration function (semi pseudocode)
```c++
Writer* kernels::getNamedWriter(std::string name, Solver& solver) {
if(name == "ConservedWriter") return new GRMHD::ConservedWriter(*static_cast<exahype::solvers::LimitingADERDGSolver*>(solver));
if(name == "usw") return new GRMHD::SomeOtherWriter(*static_cast<exahype::solvers::ADERDGSolver*>(solver));
if(name == ...
else
failure: Dont now this plotter type.
}
```
We then can replace the generated current section
```c++
exahype::plotters::RegisteredPlotters.push_back( new exahype::plotters::Plotter(0,0,parser,new GRMHD::ConservedWriter( *static_cast<exahype::solvers::LimitingADERDGSolver*>(exahype::solvers::RegisteredSolvers[0])) ));
exahype::plotters::RegisteredPlotters.push_back( new exahype::plotters::Plotter(0,1,parser,new GRMHD::ConservedWriter( *static_cast<exahype::solvers::LimitingADERDGSolver*>(exahype::solvers::RegisteredSolvers[0])) ));
exahype::plotters::RegisteredPlotters.push_back( new exahype::plotters::Plotter(0,2,parser,new GRMHD::ConservedWriter( *static_cast<exahype::solvers::LimitingADERDGSolver*>(exahype::solvers::RegisteredSolvers[0])) ));
exahype::plotters::RegisteredPlotters.push_back( new exahype::plotters::Plotter(0,3,parser,new GRMHD::ConservedWriter( *static_cast<exahype::solvers::LimitingADERDGSolver*>(exahype::solvers::RegisteredSolvers[0])) ));
exahype::plotters::RegisteredPlotters.push_back( new exahype::plotters::Plotter(0,4,parser,new GRMHD::IntegralsWriter( *static_cast<exahype::solvers::LimitingADERDGSolver*>(exahype::solvers::RegisteredSolvers[0])) ));
```
with a non-generated section (pseudocode)
```c+++
for(const int& solvernum : Parser->getSolvers()) {
int plotternum=0;
for( const string& plottername : Parser->getPlotterNamesForSolver(solvernum)) {
exahype::plotters::RegisteredPlotters.push_back(new exahype::plotters::Plotter(solvernum,plotternum++,parser,getNamedWriter(plottername, exahype::solvers::RegisteredSolvers[solvernum]));
}
}
```
This is something I can do on my own. If the Parser gives the neccessary data...https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/194Erasing of Peano grid vertices2018-03-20T16:13:04+01:00Ghost UserErasing of Peano grid verticesCurrently, I remove only my data from the heap when I erase a cell.
The Peano grid structure still persists.
I tried a variety of tricks to remove the grid structure on the fly during
the mesh refinement but always run into problems:
- ...Currently, I remove only my data from the heap when I erase a cell.
The Peano grid structure still persists.
I tried a variety of tricks to remove the grid structure on the fly during
the mesh refinement but always run into problems:
- Inconsistent adjacency indices
- Spurious erases of neighbours
- Data races with TBB
- ...
I think now it does make more sense to erase Peano vertices after the
mesh refinement. Here the refined grid is fixed. Refinement does not
interfere with erasing.
I think of either:
- Erase Peano vertices on the fly during the time stepping iterations
however I do not know the overhead and impact on the parallelisation.
- Run a few iterations of a dedicated mesh erasing adapter after the
refinement.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/303error in compilation2019-10-15T09:38:50+02:00Duo Lid.li@lmu.deerror in compilationwhen compiling master Exahype, i have a bunch of errors with GaussLobattoBasis.hsnippet:
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from ...when compiling master Exahype, i have a bunch of errors with GaussLobattoBasis.hsnippet:
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(1): error: expected a declaration
namespace kernels {
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(22): warning #12: parsing restarts here after previous syntax error
extern double nodes_1[1+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(23): error: invalid storage class for a class member
extern double nodes_2[2+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(24): error: invalid storage class for a class member
extern double nodes_3[3+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(25): error: invalid storage class for a class member
extern double nodes_4[4+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(26): error: invalid storage class for a class member
extern double nodes_5[5+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(27): error: invalid storage class for a class member
extern double nodes_6[6+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(28): error: invalid storage class for a class member
extern double nodes_7[7+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(29): error: invalid storage class for a class member
extern double nodes_8[8+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(30): error: invalid storage class for a class member
extern double nodes_9[9+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(31): error: invalid storage class for a class member
extern double nodes_10[10+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(32): error: invalid storage class for a class member
extern double nodes_11[11+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(33): error: invalid storage class for a class member
extern double nodes_12[12+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(34): error: invalid storage class for a class member
extern double nodes_13[13+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(35): error: invalid storage class for a class member
extern double nodes_14[14+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(36): error: invalid storage class for a class member
extern double nodes_15[15+1];
^
In file included from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/GaussLobattoBasis.h(20),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/aderdg/generic/Kernels.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/../Curvilinear/CurvilinearTransformation.h(6),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/include/Context/ContextCurvilinear.h(11),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.h(28),
from /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ApplicationExamples/ExaSeis/DynamicRupture/TPV29/SlipWeakeningSolver.cpp(9):
/hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./ExaHyPE/kernels/generated/GaussLobattoBasis.hsnippet(37): error: invalid storage class for a class member
extern double* nodes[15+1];
does any one have any idea?https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/104EulerFlow and SRMHD Convergence2018-06-15T15:13:43+02:00Ghost UserEulerFlow and SRMHD ConvergenceThere was a mistake in the EulerFlow and SRMHD error computation routines:
The analytical solution expressed in conserved variables was compared vs.
the numerical solution expressed in primitive variables.
I will report the errors ...There was a mistake in the EulerFlow and SRMHD error computation routines:
The analytical solution expressed in conserved variables was compared vs.
the numerical solution expressed in primitive variables.
I will report the errors for EulerFlow Moving Gauss 2D and SRMHD Alfven Wave 2D after the fix to make sure
we do not get the computation wrong again.
I further added a Gauss-Legendre quadrature based error computation
to the EulerFlow project in order to validate the Riemann quadrature values. Indeed, there is only a small
difference. Convergence is observed in both cases .
(I did not calculate the rates yet. Sven should have a script for that.)
Some debugging tips:
* **We need to make sure that always matching variables (cons.,prim.) are compared against eacht other -- also for the SRMHD convergence tests.)**
* **Even if the solver might get inaccurate in the long term, the initial data should always be better resolved on a finer mesh. We should always see
convergence in projecting the initial data!**
## EulerFlow moving Gauss 2D-- rho errors
settings: fused-algorithmic-steps=on,timestepping=globalfixed
### Gauss-Legendre quadrature:
maximum-mesh-size = 0.5
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.750382e-13 8.925639e-13 1.311617e-12 6.208367e-13 8.950533e-13
2 5.459310e-02 1.038389e-03 2.920092e-03 1.570339e-02 1.460953e-06 8.411925e-04
3 1.013872e-01 1.410669e-03 4.361094e-03 2.551770e-02 2.776317e-06 1.122370e-03
4 1.559803e-01 1.239346e-03 3.639293e-03 2.137039e-02 8.062877e-07 1.042290e-03
5 2.027744e-01 8.993649e-04 2.033884e-03 9.321796e-03 4.938160e-06 8.758699e-04
6 2.573675e-01 8.112109e-04 2.558304e-03 2.066562e-02 1.188077e-06 8.807274e-04
7 3.041615e-01 1.171406e-03 3.257146e-03 2.334808e-02 1.116484e-06 1.177089e-03
8 3.509556e-01 1.247573e-03 3.321656e-03 2.104360e-02 8.826953e-07 1.255401e-03
9 4.055487e-01 9.415592e-04 2.712697e-03 2.152044e-02 4.344998e-08 9.937017e-04
```
maximum-mesh-size = 0.15
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.772611e-13 8.949179e-13 1.507683e-12 6.206147e-13 8.976664e-13
2 5.199343e-02 1.824504e-04 1.173739e-03 1.840355e-02 1.319400e-11 1.656369e-04
3 1.013872e-01 1.463228e-04 8.135031e-04 1.585928e-02 9.383494e-12 1.586405e-04
4 1.507809e-01 1.511930e-04 7.129781e-04 7.640309e-03 6.468159e-13 1.499519e-04
5 2.001747e-01 1.799638e-04 1.039565e-03 1.571490e-02 7.177703e-12 1.593651e-04
6 2.521681e-01 2.472488e-04 1.460660e-03 2.225137e-02 1.960077e-11 2.191721e-04
7 3.015619e-01 2.123867e-04 9.078946e-04 1.156214e-02 2.549394e-11 2.037580e-04
8 3.509556e-01 1.996272e-04 1.025265e-03 1.816118e-02 3.995360e-11 2.028881e-04
9 4.003494e-01 2.200875e-04 1.407165e-03 2.126848e-02 6.046564e-10 1.847641e-04
```
maximum-mesh-size = 0.05
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.775711e-13 8.953091e-13 1.799227e-12 6.206147e-13 8.981400e-13
2 5.026031e-02 8.648108e-06 1.024604e-04 2.880816e-03 1.421085e-14 8.765676e-06
3 1.005206e-01 1.310812e-05 1.404085e-04 4.134064e-03 1.211253e-13 1.270370e-05
4 1.507809e-01 1.204385e-05 1.501490e-04 7.148238e-03 2.442491e-14 1.169658e-05
5 2.001747e-01 1.309551e-05 1.688337e-04 5.761659e-03 3.804734e-13 1.238790e-05
6 2.504350e-01 1.609669e-05 1.985765e-04 7.332591e-03 1.674216e-13 1.489681e-05
7 3.006953e-01 1.583313e-05 2.188010e-04 9.998774e-03 3.019807e-14 1.503079e-05
8 3.500891e-01 1.771167e-05 2.384444e-04 8.503132e-03 7.386314e-13 1.616133e-05
9 4.003494e-01 1.941804e-05 2.599394e-04 1.078235e-02 2.096101e-13 1.797866e-05
```
### Riemann quadrature:
maximum-mesh-size = 0.5
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.208791e-04 2.760006e-03 1.152340e-02 1.857846e-09 8.208791e-04
2 5.459310e-02 1.136775e-03 2.987134e-03 1.768960e-02 1.300473e-06 1.136775e-03
3 1.013872e-01 1.339213e-03 3.890847e-03 2.752230e-02 1.402174e-07 1.339213e-03
4 1.559803e-01 1.298609e-03 3.445100e-03 2.234650e-02 3.597158e-09 1.298609e-03
5 2.027744e-01 1.040780e-03 2.383341e-03 9.564163e-03 1.192050e-09 1.040780e-03
6 2.573675e-01 1.125659e-03 2.468552e-03 8.197815e-03 7.159261e-10 1.125659e-03
7 3.041615e-01 1.645793e-03 3.824891e-03 1.758752e-02 3.099943e-09 1.645793e-03
8 3.509556e-01 1.710023e-03 4.109781e-03 1.967037e-02 6.597448e-09 1.710023e-03
9 4.055487e-01 1.238782e-03 2.674109e-03 8.709187e-03 5.247924e-09 1.238782e-03
```
maximum-mesh-size = 0.15
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.976664e-13 9.187753e-13 1.507683e-12 6.206147e-13 8.976664e-13
2 5.199343e-02 1.656369e-04 1.006196e-03 1.840355e-02 1.319400e-11 1.656369e-04
3 1.013872e-01 1.586405e-04 8.668425e-04 1.585928e-02 9.383494e-12 1.586405e-04
4 1.507809e-01 1.499519e-04 7.055692e-04 7.640309e-03 6.468159e-13 1.499519e-04
5 2.001747e-01 1.593651e-04 8.780007e-04 1.571490e-02 7.177703e-12 1.593651e-04
6 2.521681e-01 2.191721e-04 1.259805e-03 2.225137e-02 1.960077e-11 2.191721e-04
7 3.015619e-01 2.037580e-04 8.841593e-04 1.156214e-02 2.549394e-11 2.037580e-04
8 3.509556e-01 2.028881e-04 1.060100e-03 1.816118e-02 3.995360e-11 2.028881e-04
9 4.003494e-01 1.847641e-04 1.146441e-03 2.126848e-02 6.046564e-10 1.847641e-04
```
maximum-mesh-size = 0.05
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 8.981400e-13 9.193660e-13 1.799227e-12 6.206147e-13 8.981400e-13
2 5.026031e-02 8.765676e-06 1.038935e-04 2.880816e-03 1.421085e-14 8.765676e-06
3 1.005206e-01 1.270370e-05 1.360760e-04 4.134064e-03 1.211253e-13 1.270370e-05
4 1.507809e-01 1.169658e-05 1.486293e-04 7.148238e-03 2.442491e-14 1.169658e-05
5 2.001747e-01 1.238790e-05 1.582185e-04 5.761659e-03 3.804734e-13 1.238790e-05
6 2.504350e-01 1.489681e-05 1.849682e-04 7.332591e-03 1.674216e-13 1.489681e-05
7 3.006953e-01 1.503079e-05 2.125449e-04 9.998774e-03 3.019807e-14 1.503079e-05
8 3.500891e-01 1.616133e-05 2.143398e-04 8.503132e-03 7.386314e-13 1.616133e-05
9 4.003494e-01 1.797866e-05 2.425857e-04 1.078235e-02 2.096101e-13 1.797866e-05
```
## SRMHD Alfven wave 2D -- vely errors
settings: fused-algorithmic-steps=off,timestepping=global
maximum-mesh-size = 0.5
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 1.547407e-03 2.445560e-03 4.573094e-03 2.602050e-05 1.547407e-03
2 6.015317e-02 1.111526e-03 1.657690e-03 3.769028e-03 3.275033e-06 1.111526e-03
3 1.002207e-01 1.398944e-03 1.869334e-03 4.319537e-03 2.836757e-04 1.398944e-03
```
maximum-mesh-size = 0.15
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 2.302968e-05 3.590189e-05 7.108230e-05 1.139464e-07 2.302968e-05
2 5.344862e-02 2.056201e-05 3.057057e-05 7.954177e-05 3.314441e-07 2.056201e-05
3 1.002171e-01 2.300635e-05 3.160379e-05 7.990042e-05 5.301505e-07 2.300635e-05
```
maximum-mesh-size = 0.05
```
plotindex time l1norm l2norm max min avg
1 0.000000e+00 2.913558e-07 4.523600e-07 9.036162e-07 4.721517e-10 2.913558e-07
2 5.122013e-02 3.221749e-07 4.326605e-07 1.331551e-06 5.221886e-09 3.221749e-07
3 1.002133e-01 3.809067e-07 5.525386e-07 1.662368e-06 1.168645e-09 3.809067e-07
```https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/89exa compile log should contain everything2018-06-15T15:13:44+02:00Ghost Userexa compile log should contain everything```
Make failed!
Uploading the file extended.log, please send the following link to your collaborators:
http://sprunge.us/GSQA
```
But the link only contains the output of `make`, not the toolkit invocation which is also important somet...```
Make failed!
Uploading the file extended.log, please send the following link to your collaborators:
http://sprunge.us/GSQA
```
But the link only contains the output of `make`, not the toolkit invocation which is also important sometimes. Make sure it contains that, too.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/91ExaHyPE binaries should tell about compile time flags2018-06-15T15:13:44+02:00Ghost UserExaHyPE binaries should tell about compile time flagsSometimes I have an ExaHyPE binary and don't remember with which options I compiled it. I would like to see as much information as possible, for instance:
```
$ ExaHyPE-Euler-p5 --version
This is ExaHyPE.
Compiled at Mon 19 Dec 15:43:36 ...Sometimes I have an ExaHyPE binary and don't remember with which options I compiled it. I would like to see as much information as possible, for instance:
```
$ ExaHyPE-Euler-p5 --version
This is ExaHyPE.
Compiled at Mon 19 Dec 15:43:36 CET 2016
Based on ExaHyPE repository git commit fbec392.
Created with md5sum(ExaHyPE.jar) = xxxxxxxxxx
Built with options:
COMPILER | GNU
MODE | Release
SHAREDMEM | None
DISTRIBUTEDMEM | None
Compiled with gcc v1.2.3.4
Compiled with MPI v2.3.4.5
Compiled with TBB v1.2.3.4
```
Especially the `Mode = Release|Debug|...` would be very helpful.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/254Finite Volumes Kernel Issues2019-09-20T15:43:37+02:00Ghost UserFinite Volumes Kernel IssuesA list of known issues with the FV kernels.
This list is open for extension. It may be split into separate issues
later on.
* MUSCL-Hancock implementation is second-order accurate but allows unlimited slopes
in the first ghost layer at ...A list of known issues with the FV kernels.
This list is open for extension. It may be split into separate issues
later on.
* MUSCL-Hancock implementation is second-order accurate but allows unlimited slopes
in the first ghost layer at corners and edges of the patch. This caused certain applications to crash.
* Stack allocation of multiple temporary arrays leads to Seg Faults for applications where these arrays
become huge (e.g. CCZ4). (This issue might be resolved after issue #255 has been resolved.)
* No (JM) optimisation has been performed yet. For most applications, the FV updates
are bandwith-bound indicating potential for speedups by vectorisation.