Energy growth in all OCCA modes about laghos

[vladotomov]

Hi Noel,
Does the same command line run well with the master branch?

[dmed256]

I haven't ported over Laghos+OCCA to OCCA 1.0
One of the example problems did blow up after some time, it hasn't been resolved yet

[noelchalmers]

@vladotomov The commands I've used are the usual test cases

./laghos -p 1 -m data/square01_quad.mesh -rs 3 -tf 0.8 -no-vis 
./laghos -p 1 -m data/cube01_hex.mesh -rs 2 -tf 0.6 -no-vis 

I would believe they run in the master branch, but I can double check. @dmed256 is this the case that blows up?

The fact that the OCCA1.0 version compiles the kernels and runs is promising. Did any MFEM/Laghos kernel code change significantly between OCCA 0.2 and 1.0?

[camierjs]

I think that we've had this issue since the first version of the OCCA kernels.
Depending on the problem, there is either a slight divergence or such an increase.
You can use the CFL to try and control it, but we'll need to correct this.
Thanks Noel for the feedback.

[camierjs]

I'm closing this issue as all sub-directories have this issue.
The MFEM-4.0 version of Laghos that will be merged soon won't have this kind of behaviour by using directly BLAS functions from MFEM to compute the eigen values and compression directions directly from the devices. This should allow the subdirectories to sync with master and not 'diverge'.