r13088 | cdaley | 2011-01-19 17:36:24 -0600 (Wed, 19 Jan 2011) | 7 lines
Changed paths:
   M /trunk/source/Grid/GridParticles/GridParticlesMapToMesh/Paramesh/gr_ptSameProcMap.F90

Fix bug in a particle mapping routine.  An incorrect sign affects the
density that is deposited on the grid.

Bug found by Milos Milosavljevic and Chalence from University of
Texas.  They found that mass was not being conserved in their
simulations.



r13149 | cdaley | 2011-01-27 22:46:51 -0600 (Thu, 27 Jan 2011) | 7 lines
Changed paths:
   M /trunk/source/Grid/GridParticles/GridParticlesMapToMesh/Paramesh/gr_ptProlongSmear.F90

Fix an interface mismatch bug detected by intel compiler.  The
explicit interface for prolong_temp_patch has an array
recv(1:2,1:2,1:2), but the subroutine prolong_temp_patch has an array
recv(ia:ib,ja:jb,ka:kb).  Practically this will never cause a problem
because the subroutine is only called with the hard-coded values
ia=1,ib=2,ja=1,jb=2,ka=1,kb=2.



r13165 | cdaley | 2011-01-30 22:04:46 -0600 (Sun, 30 Jan 2011) | 12 lines
Changed paths:
   M /trunk/source/Grid/GridParticles/GridParticlesMapToMesh/Paramesh/PttoPt/gr_ptMapData.F90
   M /trunk/source/Grid/GridParticles/GridParticlesMapToMesh/Paramesh/PttoPt/gr_ptMapInit.F90
   M /trunk/source/Grid/GridParticles/GridParticlesMapToMesh/Paramesh/PttoPt/gr_ptStoreOffBlockCells.F90

Ensure a receive buffer is sized large enough in the non-default map
particles to mesh implementation (PttoPt).  The receive buffer used to
be sized as the maximum of the amount sent by any process, however, it
is possible for the accumulated receive size on a single process to
exceed the maximum send size.

This does not affect the default map particles to mesh implementation
(MoveSieve) because in this implementation a process can only receive
from 1 process at a time.

Unfortunately this is a bug that is present in FLASH3.3-release too.