[pooma-dev] Field semantics

[James Crotinger]

I have to say that I'm surprised that this wasn't caught by some assertion.

You don't give all of the details, but it looks like the problem was due to
a domain mismatch. When you say 

   T * rh

This operation is, I believe, only performed on the "physical" domain - i.e.
it is not extended into the guards. But your p(It) explicitly says to assign
into the guards (since It is the "totalDomain", which includes the guards).
Now, given that the domain of the RHS is not the same as the domain of the
LHS, this ought to raise an exception somewhere. 

	Jim



>-----Original Message-----
>From: Richard Guenther [mailto:rguenth at tat.physik.uni-tuebingen.de]
>Sent: Tuesday, April 30, 2002 9:03 AM
>To: pooma-dev at pooma.codesourcery.com
>Subject: [pooma-dev] Field semantics
>
>Hi!
>
>Just spent the last two days debugging my HD simulation to figure out
>that (p,T,rh, etc. being Fields with 1 guard cell)
>
>Interval<1> It = rh.totalDomain();
>p(It) = T * rh;
>
>is not the same as
>
>p(It) = T(It) * rh(It);
>
>why is this so? Do I generally need to specify the domain on
>both sides of the equations?
>
>Richard.
>
>--
>Richard Guenther <richard.guenther at uni-tuebingen.de>
>WWW: http://www.tat.physik.uni-tuebingen.de/~rguenth/
>The GLAME Project: http://www.glame.de/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sourcerytools.com/pipermail/pooma-dev/attachments/20020430/d7c2af0f/attachment.html>