Vladimir Tomov

Have you looked at the Memory manager section in the [website docs](https://mfem.org/gpu-support/)?

> Since rhs is the result of ForcePA->Mult(one, rhs), rhs is allocated in GPU's memory, isn't it? If `rhs` is allocated on the device, and `ForcePA` supports GPU execution, then...

Can you show the code that prints `B`? Did you do a `HostRead()` before the print? There should be no difference between the treatment of `rhs` and `B`, they are...

The differences seem too big, I'd guess it's a bug.

Strange indeed. Are you running this in your branch? Can I see the code somehow?

The problem was that the corresponding mass kernels for the velocity and energy matrices were not in MFEM. Try [this PR](https://github.com/mfem/mfem/pull/4210) for the MFEM build, for Q5Q4 2D, and let...

@sungho91, sure, more orders can be added. - The first parameter is the number of dofs in each direction. For Q5Q4 we have 6 (for the Q5 H1 space) or...

@tzanio can we remove the `in-review` label?

Hi @sungho91 the methods `-lo 3` and `-lo 4` are not implemented for NC meshes and won't work. We haven't tested Remhos with NC meshes, but I expect some combinations...

I'm on travel this week, but I'll take a look when I come back.