sim.field.grid.DoubleGrid2D Maven / Gradle / Ivy
Show all versions of mason Show documentation
/*
Copyright 2006 by Sean Luke and George Mason University
Licensed under the Academic Free License version 3.0
See the file "LICENSE" for more information
*/
package sim.field.grid;
import sim.util.*;
/**
A wrapper for 2D arrays of doubles.
This object expects that the 2D arrays are rectangular. You are encouraged to access the array
directly. The object
implements all of the Grid2D interface. See Grid2D for rules on how to properly implement toroidal
or hexagonal grids.
The width and height of the object are provided to avoid having to say field[x].length, etc.
*/
public /*strictfp*/ class DoubleGrid2D extends AbstractGrid2D
{
private static final long serialVersionUID = 1;
public double[/**x*/][/**y*/] field;
public double[][] getField() { return field; }
public DoubleGrid2D (int width, int height)
{
this.width = width;
this.height = height;
field = new double[width][height];
}
public DoubleGrid2D (int width, int height, double initialValue)
{
this(width,height);
setTo(initialValue);
}
public DoubleGrid2D (DoubleGrid2D values)
{
setTo(values);
}
public DoubleGrid2D(double[][] values)
{
setTo(values);
}
/** Sets location (x,y) to val */
public final void set(final int x, final int y, final double val)
{
field[x][y] = val;
}
/** Returns the element at location (x,y) */
public final double get(final int x, final int y)
{
return field[x][y];
}
/** Sets all the locations in the grid the provided element */
public final DoubleGrid2D setTo(final double thisMuch)
{
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;x fieldx[y]) min = fieldx[y];
}
return min;
}
/** Returns the mean value stored in the grid */
public final double mean()
{
long count = 0;
double mean = 0;
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;xtoNoMoreThanThisMuch are changed to toNoMoreThanThisMuch.
Returns the modified grid.
*/
public final DoubleGrid2D upperBound(final double toNoMoreThanThisMuch)
{
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;x toNoMoreThanThisMuch)
fieldx[y] = toNoMoreThanThisMuch;
}
return this;
}
/** Thresholds the grid so that values smaller than toNoLowerThanThisMuch are changed to toNoLowerThanThisMuch
Returns the modified grid.
*/
public final DoubleGrid2D lowerBound(final double toNoLowerThanThisMuch)
{
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;xwithThisMuch
Returns the modified grid.
*/
public final DoubleGrid2D add(final double withThisMuch)
{
final int width = this.width;
final int height = this.height;
if (withThisMuch==0.0) return this;
double[] fieldx = null;
for(int x=0;xbyThisMuch
Returns the modified grid.
*/
public final DoubleGrid2D multiply(final double byThisMuch)
{
if (byThisMuch==1.0) return this;
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;x 0.0)
// fieldx[y] = /*Strict*/Math.floor(fieldx[y]);
//else
// fieldx[y] = /*Strict*/Math.ceil(fieldx[y]);
}
return this;
}
/** Sets each value in the grid to rint(value). That is, each value
is rounded to the closest integer value. If two integers are the same
distance, the value is rounded to the even integer.
Returns the modified grid.
*/
public final DoubleGrid2D rint()
{
double[] fieldx = null;
final int width = this.width;
final int height = this.height;
for(int x=0;x0 && y>0 && xThen places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsMaxDistance(...)
* Returns the result DoubleBag.
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may only run in two modes: toroidal or bounded. Unbounded lookup is not permitted, and so
* this function is deprecated: instead you should use the other version of this function which has more functionality.
* If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
The origin -- that is, the (x,y) point at the center of the neighborhood -- is always included in the results.
*
*
This function is equivalent to: getNeighborsMaxDistance(x,y,dist,toroidal ? Grid2D.TOROIDAL : Grid2D.BOUNDED, true, result, xPos, yPos);
*
* @deprecated
*/
public void getNeighborsMaxDistance( final int x, final int y, final int dist, final boolean toroidal, DoubleBag result, IntBag xPos, IntBag yPos )
{
getMooreNeighbors(x, y, dist, toroidal ? TOROIDAL : BOUNDED, true, result, xPos, yPos);
}
/**
* Gets all neighbors of a location that satisfy max( abs(x-X) , abs(y-Y) ) <= dist, This region forms a
* square 2*dist+1 cells across, centered at (X,Y). If dist==1, this
* is equivalent to the so-called "Moore Neighborhood" (the eight neighbors surrounding (X,Y)), plus (X,Y) itself.
* Places each x and y value of these locations in the provided IntBags xPos and yPos, clearing the bags first.
*
*
Then places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsMaxDistance(...)
* Returns the result DoubleBag.
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may be run in one of three modes: Grid2D.BOUNDED, Grid2D.UNBOUNDED, and Grid2D.TOROIDAL. If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. If "unbounded", then the neighbors are not so restricted. Note that unbounded
* neighborhood lookup only makes sense if your grid allows locations to actually be outside this box. For example,
* SparseGrid2D permits this but ObjectGrid2D and DoubleGrid2D and IntGrid2D and DenseGrid2D do not. Finally if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
You can also opt to include the origin -- that is, the (x,y) point at the center of the neighborhood -- in the neighborhood results.
*/
public DoubleBag getMooreNeighbors( final int x, final int y, final int dist, int mode, boolean includeOrigin, DoubleBag result, IntBag xPos, IntBag yPos )
{
if( xPos == null )
xPos = new IntBag();
if( yPos == null )
yPos = new IntBag();
getMooreLocations( x, y, dist, mode, includeOrigin, xPos, yPos );
return getObjectsAtLocations(xPos,yPos,result);
}
/**
* Gets all neighbors of a location that satisfy abs(x-X) + abs(y-Y) <= dist. This region forms a diamond
* 2*dist+1 cells from point to opposite point inclusive, centered at (X,Y). If dist==1 this is
* equivalent to the so-called "Von-Neumann Neighborhood" (the four neighbors above, below, left, and right of (X,Y)),
* plus (X,Y) itself.
*
*
Places each x and y value of these locations in the provided IntBags xPos and yPos, clearing the bags first.
* Then places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsHamiltonianDistance(...)
* Returns the result DoubleBag (constructing one if null had been passed in).
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may only run in two modes: toroidal or bounded. Unbounded lookup is not permitted, and so
* this function is deprecated: instead you should use the other version of this function which has more functionality.
* If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
The origin -- that is, the (x,y) point at the center of the neighborhood -- is always included in the results.
*
*
This function is equivalent to: getNeighborsHamiltonianDistance(x,y,dist,toroidal ? Grid2D.TOROIDAL : Grid2D.BOUNDED, true, result, xPos, yPos);
*
* @deprecated
*/
public void getNeighborsHamiltonianDistance( final int x, final int y, final int dist, final boolean toroidal, DoubleBag result, IntBag xPos, IntBag yPos )
{
getVonNeumannNeighbors(x, y, dist, toroidal ? TOROIDAL : BOUNDED, true,result, xPos, yPos);
}
/**
* Gets all neighbors of a location that satisfy abs(x-X) + abs(y-Y) <= dist. This region forms a diamond
* 2*dist+1 cells from point to opposite point inclusive, centered at (X,Y). If dist==1 this is
* equivalent to the so-called "Von-Neumann Neighborhood" (the four neighbors above, below, left, and right of (X,Y)),
* plus (X,Y) itself.
*
*
Places each x and y value of these locations in the provided IntBags xPos and yPos, clearing the bags first.
* Then places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsHamiltonianDistance(...)
* Returns the result DoubleBag (constructing one if null had been passed in).
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may be run in one of three modes: Grid2D.BOUNDED, Grid2D.UNBOUNDED, and Grid2D.TOROIDAL. If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. If "unbounded", then the neighbors are not so restricted. Note that unbounded
* neighborhood lookup only makes sense if your grid allows locations to actually be outside this box. For example,
* SparseGrid2D permits this but ObjectGrid2D and DoubleGrid2D and IntGrid2D and DenseGrid2D do not. Finally if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
You can also opt to include the origin -- that is, the (x,y) point at the center of the neighborhood -- in the neighborhood results.
*/
public DoubleBag getVonNeumannNeighbors( final int x, final int y, final int dist, int mode, boolean includeOrigin, DoubleBag result, IntBag xPos, IntBag yPos )
{
if( xPos == null )
xPos = new IntBag();
if( yPos == null )
yPos = new IntBag();
getVonNeumannLocations( x, y, dist, mode, includeOrigin, xPos, yPos );
return getObjectsAtLocations(xPos,yPos,result);
}
/**
* Gets all neighbors located within the hexagon centered at (X,Y) and 2*dist+1 cells from point to opposite point
* inclusive.
* If dist==1, this is equivalent to the six neighbors immediately surrounding (X,Y),
* plus (X,Y) itself.
*
*
Places each x and y value of these locations in the provided IntBags xPos and yPos, clearing the bags first.
* Then places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsHamiltonianDistance(...)
* Returns the result DoubleBag (constructing one if null had been passed in).
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may only run in two modes: toroidal or bounded. Unbounded lookup is not permitted, and so
* this function is deprecated: instead you should use the other version of this function which has more functionality.
* If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
The origin -- that is, the (x,y) point at the center of the neighborhood -- is always included in the results.
*
*
This function is equivalent to: getNeighborsHexagonalDistance(x,y,dist,toroidal ? Grid2D.TOROIDAL : Grid2D.BOUNDED, true, result, xPos, yPos);
*
* @deprecated
*/
public void getNeighborsHexagonalDistance( final int x, final int y, final int dist, final boolean toroidal, DoubleBag result, IntBag xPos, IntBag yPos )
{
getHexagonalNeighbors(x, y, dist, toroidal ? TOROIDAL : BOUNDED, true, result, xPos, yPos);
}
/**
* Gets all neighbors located within the hexagon centered at (X,Y) and 2*dist+1 cells from point to opposite point
* inclusive.
* If dist==1, this is equivalent to the six neighbors immediately surrounding (X,Y),
* plus (X,Y) itself.
*
*
Places each x and y value of these locations in the provided IntBags xPos and yPos, clearing the bags first.
* Then places into the result DoubleBag any Objects which fall on one of these locations, clearning it first.
* Note that the order and size of the result DoubleBag may not correspond to the X and Y bags. If you want
* all three bags to correspond (x, y, object) then use getNeighborsAndCorrespondingPositionsHamiltonianDistance(...)
* Returns the result DoubleBag (constructing one if null had been passed in).
* null may be passed in for the various bags, though it is more efficient to pass in a 'scratch bag' for
* each one.
*
* This function may be run in one of three modes: Grid2D.BOUNDED, Grid2D.UNBOUNDED, and Grid2D.TOROIDAL. If "bounded",
* then the neighbors are restricted to be only those which lie within the box ranging from (0,0) to (width, height),
* that is, the width and height of the grid. If "unbounded", then the neighbors are not so restricted. Note that unbounded
* neighborhood lookup only makes sense if your grid allows locations to actually be outside this box. For example,
* SparseGrid2D permits this but ObjectGrid2D and DoubleGrid2D and IntGrid2D and DenseGrid2D do not. Finally if "toroidal",
* then the environment is assumed to be toroidal, that is, wrap-around, and neighbors are computed in this fashion. Toroidal
* locations will not appear multiple times: specifically, if the neighborhood distance is so large that it wraps completely around
* the width or height of the box, neighbors will not be counted multiple times. Note that to ensure this, subclasses may need to
* resort to expensive duplicate removal, so it's not suggested you use so unreasonably large distances.
*
*
You can also opt to include the origin -- that is, the (x,y) point at the center of the neighborhood -- in the neighborhood results.
*/
public DoubleBag getHexagonalNeighbors( final int x, final int y, final int dist, int mode, boolean includeOrigin, DoubleBag result, IntBag xPos, IntBag yPos )
{
if( xPos == null )
xPos = new IntBag();
if( yPos == null )
yPos = new IntBag();
getHexagonalLocations( x, y, dist, mode, includeOrigin, xPos, yPos );
return getObjectsAtLocations(xPos,yPos,result);
}
public DoubleBag getRadialNeighbors( final int x, final int y, final int dist, int mode, boolean includeOrigin,DoubleBag result, IntBag xPos, IntBag yPos )
{
return getRadialNeighbors(x, y, dist, mode, includeOrigin, Grid2D.ANY, true, result, xPos, yPos);
}
public DoubleBag getRadialNeighbors( final int x, final int y, final int dist, int mode, boolean includeOrigin, int measurementRule, boolean closed, DoubleBag result, IntBag xPos, IntBag yPos )
{
if( xPos == null )
xPos = new IntBag();
if( yPos == null )
yPos = new IntBag();
getRadialLocations( x, y, dist, mode, includeOrigin, measurementRule, closed, xPos, yPos );
return getObjectsAtLocations(xPos,yPos,result);
}
// For each location, puts all such objects into the result DoubleBag. Modifies
// the xPos and yPos bags so that each position corresponds to the equivalent result in
// in the result DoubleBag.
void reduceObjectsAtLocations(final IntBag xPos, final IntBag yPos, DoubleBag result)
{
if (result==null) result = new DoubleBag();
else result.clear();
for( int i = 0 ; i < xPos.numObjs ; i++ )
{
assert sim.util.LocationLog.it(this, new Int2D(xPos.objs[i],yPos.objs[i]));
double val = field[xPos.objs[i]][yPos.objs[i]] ;
result.add( val );
}
}
/* For each location, puts all such objects into the result DoubleBag. Returns the result DoubleBag.
If the provided result DoubleBag is null, one will be created and returned. */
DoubleBag getObjectsAtLocations(final IntBag xPos, final IntBag yPos, DoubleBag result)
{
if (result==null) result = new DoubleBag();
else result.clear();
for( int i = 0 ; i < xPos.numObjs ; i++ )
{
assert sim.util.LocationLog.it(this, new Int2D(xPos.objs[i],yPos.objs[i]));
double val = field[xPos.objs[i]][yPos.objs[i]] ;
result.add( val );
}
return result;
}
}