[Nek5000-users] fix projection for multiple passive scalars
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu Oct 19 09:13:13 CDT 2017
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 of nek5000-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
--
Karlsruhe Institute of Technology (KIT)
Institute of Fluid Mechanics
M.Sc. Steffen Straub
Doctoral Researcher
Kaiserstraße 10
Building 10.23
76131 Karlsruhe, Germany
Phone: +49 721 608-43027
E-mail: steffen.straub at kit.edu
Web: http://www.istm.kit.edu
Registered office:
Kaiserstraße 12, 76131 Karlsruhe, Germany
KIT – The Research University in the Helmholtz Association
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171019/d650be7b/attachment.html>
More information about the Nek5000-users
mailing list