[Nek5000-users] how to implement integration like (u-x, p_x) and (p_x, p_x) in Pn-Pn-2 setting

[nek5000-users at lists.mcs.anl.gov]

Sorry -

for u*p, you need to do something different...

I suggest map pr to pm1 (pressure on mesh 1), then:

       n = nx1*ny1*nz1*nelv
       val = glsc3(vx,bm1,pm1,n)

would give you what you want.

In the code, I find that you can map pr to pm1 as:

c    Map pressure onto mesh 1   (vxx and vyy are used as work arrays)
       call mappr(pm1,pr,vxx,vyy)


--- Paul