[Nek5000-users] dealias; mass matrix; interp_rstd

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

Dear Neks,
I solve
\partial A/\partial t = \nabla^2 A + u x Curl A,
where A = (s1,s2,s3)
I solve three scalar equations with respective components of u x Curl A as
forcing terms in qvol.
I would like to compute the dealiased product of u x Curl A
I build my code based on kov_stokes_test.usr and  subroutines in induct.f
I compute local curl, project curl and U onto fine mesh, perform product
and project the product back to coarse mesh using interp_rstd
However, I am not sure about the multiplication of mass matrices for the
product in these subroutines.

I attach my usr, rea and SIZE files.
Your suggestions are welcome.
Thanks,
Sandeep
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180523/3ca332ab/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conj_ht.rea
Type: application/octet-stream
Size: 304787 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180523/3ca332ab/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conj_ht.usr
Type: application/octet-stream
Size: 11669 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180523/3ca332ab/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SIZE
Type: application/octet-stream
Size: 4155 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180523/3ca332ab/attachment-0005.obj>