[Nek5000-users] fix projection for multiple passive scalars

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

Hi Neks,


thanks for dealing with the projection of passive scalars. I saw that there were 2 commits enabling projection for passive scalars
(586e675926e14ca89982adf7dfc3a4c682c6d8c7) and disabling them again
(71573b1904e27a9fa15bff13752292545d3d2943) because of too much memory consumption.

If I understand the code correctly, there is no simple flag to enable projection for passive scalars in the rea file.
Can you tell what needs to be changed in the code in order to enable them?


Best,
Steffen Straub



Message: 5
Date: Thu, 19 Oct 2017 14:13:13 +0000
From:nek5000-users at lists.mcs.anl.gov
To:"nek5000-users at lists.mcs.anl.gov"
	<nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] fix projection for multiple passive
	scalars
Message-ID:
	<mailman.749.1508422411.19102.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="windows-1252"


Yes, this is a bug.


Sorry - Because of another issue (now resolved), I hadn't been

attending to this part of projection yet (since it was essentially

broken before).


I'll try to get this fixed asap.


Paul


________________________________
From: Nek5000-users<nek5000-users-bounces at lists.mcs.anl.gov>  on behalf ofnek5000-users at lists.mcs.anl.gov  <nek5000-users at lists.mcs.anl.gov>
Sent: Thursday, October 19, 2017 3:50:23 AM
To:nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] fix projection for multiple passive scalars

Hi Neks,

in order to use projection (p94!=0) for one temperature field and at
least one additional passive scalar, I believe that the variables
napproxt and approxt need to be extended by the size of ldimt. See the
attached git-diff.txt for my changes to the code.

When I adjust the code like that, projection for veloctity, temperature
and all passive scalars works, otherwise it fails with:

...
22  Error Hmholtz TEMP     100          NaN          NaN 1.0000E-08
22proj_ortho:    1   2 TEMP   Detect rank deficiency: NaN         NaN
22  Project PS 1                  NaN          NaN          NaN 1   0
22  Error Hmholtz PS 1     100          NaN          NaN 1.0000E-08
22proj_ortho:    1   2 PS 1   Detect rank deficiency: NaN         NaN
22  Project PS 2                  NaN          NaN          NaN 1   0
22  Error Hmholtz PS 2     100          NaN          NaN 1.0000E-08
22proj_ortho:    1   2 PS 2   Detect rank deficiency: NaN         NaN
...

Can one of the developers please confirm and if confirmed include it in
the master branch?


Best,
Steffen Straub

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171114/ff79dd8f/attachment.html>