RE: [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:address@hidden]
>Sent: Tuesday, April 30, 2002 9:03 AM
>To: address@hidden
>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 <address@hidden>
>WWW: http://www.tat.physik.uni-tuebingen.de/~rguenth/
>The GLAME Project: http://www.glame.de/