[Nek5000-users] Possible bug

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

Hi Frank,

You'll find many such instances in nek, since we are only
concerned with whether the space is allocated, and not
so much about its type, which matters only at point of use
(subject to caveats of sufficient space and proper byte
alignment).   In this particular case, wk is a work array
used inside mfi_gets - only pm1 is returned as useful
output.

Paul


On Mon, 12 Apr 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hello all,

In "ic.f", line 2030 there is "real*4 wk(lwk)" in the declarations for
subroutine "mfi_gets".  However, at line 2489 there is "call
mfi_gets(pm1,wk,lwk,.false.)" where "wk" is double precison.

Cheers,
Frank
-- 
Frank Herbert Muldoon, Ph.D. Mechanical Engineering
Technische Universität Wien (Technical University of Vienna)
Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
Mechanics and Heat Transfer)
Resselgasse 3
1040 Wien
Tel: +4315880132232
Fax: +4315880132299 
Cell:+436765203470
fmuldoo (skype)
http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html

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