<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<pre wrap="">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: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
To: <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov">"nek5000-users@lists.mcs.anl.gov"</a>
<a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov"><nek5000-users@lists.mcs.anl.gov></a>
Subject: Re: [Nek5000-users] fix projection for multiple passive
scalars
Message-ID:
<a class="moz-txt-link-rfc2396E" href="mailto:mailman.749.1508422411.19102.nek5000-users@lists.mcs.anl.gov"><mailman.749.1508422411.19102.nek5000-users@lists.mcs.anl.gov></a>
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 <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users-bounces@lists.mcs.anl.gov"><nek5000-users-bounces@lists.mcs.anl.gov></a> on behalf of <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov"><nek5000-users@lists.mcs.anl.gov></a>
Sent: Thursday, October 19, 2017 3:50:23 AM
To: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>
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
</pre>
</body>
</html>