[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 07:38:12 CDT 2017


If you think this is a bug, please report it on GitHub.

On 19 Oct 2017, at 10:52, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:

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 <mailto:steffen.straub at kit.edu> 
Web: http://www.istm.kit.edu <http://www.istm.kit.edu> 

Registered office:
Kaiserstraße 12, 76131 Karlsruhe, Germany

KIT – The Research University in the Helmholtz Association

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov> 
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

<git-diff.txt>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171019/7e55a674/attachment.html>


More information about the Nek5000-users mailing list