ExaHyPE-Engine issueshttps://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues2018-06-15T15:13:43+02:00https://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/267Domain sizes of TOV star2019-03-07T17:26:15+01:00Ghost UserDomain sizes of TOV starInterior: `[-3,+3]`, `maximum-mesh-size = 0.23`, DG order 3
Plotting repeat: `integrals=0.1, volumetric=0.5`Interior: `[-3,+3]`, `maximum-mesh-size = 0.23`, DG order 3
Plotting repeat: `integrals=0.1, volumetric=0.5`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/210Efficent output format for large simulations2018-03-02T13:23:47+01:00Leonhard RannabauerEfficent output format for large simulationsAs soon as output reaches a certain size plotting is our major bottleneck. (a mesh of ~40 mio dofs with 14 unknowns for pvtu on supermuc takes 1h per plot)
We should look into different approaches like the ASYNC lib https://github.com/TU...As soon as output reaches a certain size plotting is our major bottleneck. (a mesh of ~40 mio dofs with 14 unknowns for pvtu on supermuc takes 1h per plot)
We should look into different approaches like the ASYNC lib https://github.com/TUM-I5/ASYNC which sacrifices a singe thread per rank on output.
For supermuc: We should also start to look into parallel file systems like LUSTRE https://en.wikipedia.org/wiki/Lustre_(file_system).Leonhard RannabauerLeonhard Rannabauerhttps://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/302error when use dynamic AMR in LimitedADERDG2020-11-20T15:10:48+01:00Duo Lid.li@lmu.deerror when use dynamic AMR in LimitedADERDGparameter tarch::parallel::Node::getInstance().getRank(): 17
parameter "if the neighbour data buffer is empty, you have perhaps forgotten to call releaseMessages() on the heap in the traversal before": if the neighbour data buffer is emp...parameter tarch::parallel::Node::getInstance().getRank(): 17
parameter "if the neighbour data buffer is empty, you have perhaps forgotten to call releaseMessages() on the heap in the traversal before": if the neighbour data buffer is empty, you have perhaps forgotten to call releaseMessages() on the heap in the traversal before
ExaHyPE-Diffuse: /hppfs/work/pr63qo/di52lak2/ExaHyPE_deadlock/ExaHyPE-Engine/./Peano/peano/../peano/heap/BoundaryDataExchanger.cpph:282: void peano::heap::BoundaryDataExchanger<Data, SendReceiveTaskType, VectorContainer>::freeAllMessagesInDeployBuffer() [with Data = char, SendReceiveTaskType = peano::heap::SendReceiveTask<char>, VectorContainer = std::vector<char, std::allocator<char>>]: Assertion `false' failed.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/232EulerFV headers2018-07-25T12:41:20+02:00Ghost UserEulerFV headersHi,
in EulerFV there seems to be a mismatch/inconsistency between the MyEulerSolver.cpp and MyEulerSolver.h
on the repository this is the function used:
void EulerFV::MyEulerSolver::eigenvalues(const double* const Q, const int normalNo...Hi,
in EulerFV there seems to be a mismatch/inconsistency between the MyEulerSolver.cpp and MyEulerSolver.h
on the repository this is the function used:
void EulerFV::MyEulerSolver::eigenvalues(const double* const Q, const int normalNonZeroIndex, double* lambda)
but there is no such signature in the header generated by the toolkit. I assume const int normalNonZeroIndex, const int d and const int dIndex refer to the same variable.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/199Extend batching to global time stepping2017-11-21T14:41:53+01:00Ghost UserExtend batching to global time steppingI am currently working on extending the batching from "globalfixed" to "global" time stepping.
The idea is here to only reduce the time step size and other data in the last iteration of the batch.
If something went wrong, we then roll ba...I am currently working on extending the batching from "globalfixed" to "global" time stepping.
The idea is here to only reduce the time step size and other data in the last iteration of the batch.
If something went wrong, we then roll back to the state at the beginning of the batch.
This WiP issue is for tracking my progress and issues I ran into.
Reduction of time step size
---------------------------
While we will not update the currently used CFL time step size in every single iteration
in the planned batched global time stepping scheme,
we will search the minimum CFL time step size over all batch iterations.
This allows us to check if the used time step size did violate the CFL constraint
during the execution of the batch.
The consequence might then be to use the newly determined minimum batch time step size as
time step size for a rerun of the whole batch.
Basic modifications to the fused time stepping algorithm
--------------------------------------------------------
- We will store the previous solution only at the beginning of a batch.
The same happens with the previous time stamp and time step size.
- TBC...
Rollbacks:
----------
TBC...
https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/164Extrapolate space-time polynomials and wrap boundary condition imposition in ...2017-07-12T20:24:40+02:00Ghost UserExtrapolate space-time polynomials and wrap boundary condition imposition in time-integralThis will further be a first step for local time steppingThis will further be a first step for local time steppinghttps://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/214Faster mesh refinement2018-03-02T11:23:39+01:00Ghost UserFaster mesh refinementProblem 1 (Costly operation is performed while no parallelism available)
------------------------------------------------------------------------
Peano's shared-memory parallelism is based on identifying regular subgrids in the
tree. If...Problem 1 (Costly operation is performed while no parallelism available)
------------------------------------------------------------------------
Peano's shared-memory parallelism is based on identifying regular subgrids in the
tree. If a new cells is introduced to the tree it might not be yet be identified
as part of a regular subgrid and the operations performed on the cell
are thus parallelised.
Solution: If imposing initial conditions or evaluating a refinement criterion is too
expensive we could let the user choose to perform it as background task.
This might lead to more mesh setup iteration but potentially to a better
exploitation of the available cores. It should also benefit the hiding of MPI communication
during the mesh setup.
Problem 2 (Overall concurrency)
-------------------------------
ExaHyPE's regular shared-memory parallelism during the mesh setup is currently further limited as
multiple cells might write to the same vertex in order to set refinement events.
Solution: We should be able to solve this by inverting the control. The
vertex checks in touchVertexLastTime if any cell has set a refinement event,
and refines if that is the case. This would increase the concurrency of the
enterCell operations.
Problem 3 (Memory)
-------------------------------
ExaHyPE's initial mesh setup is performed at the beginning by a single rank.
Gradually more and more ranks are added. In order to prevent that
any of the ranks runs out of memory during the initial mesh setup,
it might make sense to only temporarily allocate memory, impose initial
conditions, evaluate the refinement criterion, and then free the memory again (better:
recylce it).
After the initial mesh setup, we would then allocate memory on all ranks
and impose initial conditions.
Problem 4 (Load Balancing)
-------------------------------
The load-balancing does currently only count the number of cells.
It does not take the different cell types in ExaHyPE's grid into account.
The helper cell types Descendant and Ancestor have way less work to do than
the compute cells of type Cell.https://gitlab.lrz.de/exahype/ExaHyPE-Engine/-/issues/272Find out exact properties of the Cactus schemes2019-03-18T14:04:53+01:00Ghost UserFind out exact properties of the Cactus schemes* WhiskyTHC: Find details; probably in papers; maybe find comparisons with other codes
* FIL/IllinoisGRMHD: Run a 2nd order FV Cowling star, in order to compare different codes and see how they compare* WhiskyTHC: Find details; probably in papers; maybe find comparisons with other codes
* FIL/IllinoisGRMHD: Run a 2nd order FV Cowling star, in order to compare different codes and see how they compare