On 05/04/2011 12:59 PM, mark florisson wrote:
On 4 May 2011 12:45, Dag Sverre Seljebotn<d.s.seljeb...@astro.uio.no> wrote:
On 05/04/2011 12:00 PM, mark florisson wrote:
There are two remaining issue. The first is warnings for potentially
uninitialized variables for prange(). When you do
for i in prange(start, stop, step): ...
it generates code like
nsteps = (stop - start) / step;
#pragma omp parallel for lastprivate(i)
for (temp = 0; temp< nsteps; temp++) {
i = start + temp * step;
...
}
So here it will complain about 'i' being potentially uninitialized, as
it might not be assigned to in the loop. However, simply assigning 0
to 'i' can't work either, as you expect zero iterations not to touch
it. So for now, we have a bunch of warnings, as I don't see a
__attribute__ to suppress it selectively.
Isn't this is orthogonal to OpenMP -- even if it said "range", your testcase
could get such a warning? If so, the fix is simply to initialize i in your
testcase code.
No, the problem is that 'i' needs to be lastprivate, and 'i' is
assigned to in the loop body. It's irrelevant whether 'i' is assigned
to before the loop. I think this is the case because the spec says
that lastprivate variables will get the value of the private variable
of the last sequential iteration, but it cannot at compile time know
whether there might be zero iterations, which I believe the spec
doesn't have anything to say about. So basically we could guard
against it by checking if nsteps> 0, but the compiler doesn't detect
this, so it will still issue a warning even if 'i' is initialized (the
warning is at the place of the lastprivate declaration).
Ah. But this is then more important than I initially thought it was. You
are saying that this is the case:
cdef int i = 0
with nogil:
for i in prange(n):
...
print i # garbage when n == 0?
It would be in the interest of less semantic differences w.r.t. range to
deal better with this case.
Will it silence the warning if we make "i" firstprivate as well as
lastprivate? firstprivate would only affect the case of zero iterations,
since we overwrite with NaN if the loop is entered...
Dag
_______________________________________________
cython-devel mailing list
cython-devel@python.org
http://mail.python.org/mailman/listinfo/cython-devel
