1 |
mojays |
2 |
package appl.parallel.test; |
2 |
|
|
|
3 |
|
|
import java.awt.Rectangle; |
4 |
|
|
|
5 |
|
|
import schmitzm.data.WritableGrid; |
6 |
|
|
import appl.parallel.spmd.AbstractSPMDTask; |
7 |
|
|
import appl.parallel.spmd.SPMDServerInterface; |
8 |
|
|
import appl.util.RasterUtil; |
9 |
|
|
|
10 |
|
|
/** |
11 |
|
|
* @author Dominik Appl |
12 |
|
|
*/ |
13 |
|
|
public class AverageNeighborhoodTestTask extends AbstractSPMDTask { |
14 |
|
|
|
15 |
|
|
/** |
16 |
|
|
* |
17 |
|
|
*/ |
18 |
|
|
private static final long serialVersionUID = 13L; |
19 |
|
|
|
20 |
|
|
/* (non-Javadoc) |
21 |
|
|
* @see appl.parallel.spmd.SPMDTask#run() |
22 |
|
|
*/ |
23 |
|
|
public Object run(Object...parameters) { |
24 |
|
|
|
25 |
|
|
SPMDServerInterface controller = getSPMDServerController(); |
26 |
|
|
WritableGrid gridInput = (WritableGrid) controller |
27 |
|
|
.getPartition("inputGrid"); |
28 |
|
|
WritableGrid gridOutput = (WritableGrid) controller |
29 |
|
|
.getPartition("outputGrid"); |
30 |
|
|
int neighborhoodRange = (Integer) parameters[0]; |
31 |
|
|
/********************** START REMOTE ALGORITHM ********************* |
32 |
|
|
* start the calculation (notice that only the outer loop variables have changed) |
33 |
|
|
* *****************************************************************/ |
34 |
|
|
float overallSum = 0; |
35 |
|
|
Rectangle partition = controller.getLocalBounds(); |
36 |
|
|
//System.out.println("290,64: " + gridInput.getRasterSampleAsFloat(controller.getLocalCalcMinX()+289,64)); |
37 |
|
|
for (int y = controller.getLocalCalcMinY(); y < controller |
38 |
|
|
.getLocalCalcMaxY() + 1; y++) |
39 |
|
|
for (int x = controller.getLocalCalcMinX(); x < controller |
40 |
|
|
.getLocalCalcMaxX() + 1; x++) { |
41 |
|
|
float tmp = gridInput.getRasterSampleAsFloat(x,y); |
42 |
|
|
if(!Float.isNaN(tmp)) |
43 |
|
|
overallSum += tmp; |
44 |
|
|
//the local sum is simply the sum over all elements in the neighborhood |
45 |
|
|
float localSum = 0; |
46 |
|
|
//number of cells over which the sum is calculated |
47 |
|
|
int noOfCells = 0; |
48 |
|
|
//for each cell: calculate the sum of all neighbors |
49 |
|
|
for (int y2 = y - neighborhoodRange; y2 <= y |
50 |
|
|
+ neighborhoodRange; y2++) |
51 |
|
|
for (int x2 = x - neighborhoodRange; x2 <= x |
52 |
|
|
+ neighborhoodRange; x2++) |
53 |
|
|
//check if the coordinates are valid (inside the grid and not NaN) |
54 |
|
|
if (partition.contains(x2, y2)) |
55 |
|
|
if (!Float.isNaN(gridInput |
56 |
|
|
.getRasterSampleAsFloat(x2, y2))) { |
57 |
|
|
localSum += gridInput.getRasterSampleAsFloat( |
58 |
|
|
x2, y2); |
59 |
|
|
noOfCells++; |
60 |
|
|
} |
61 |
|
|
gridOutput.setRasterSample(localSum / (float) (noOfCells), x,y); |
62 |
|
|
|
63 |
|
|
} |
64 |
|
|
//copy values into results from output to input (for next step) |
65 |
|
|
RasterUtil.copyInto(gridOutput, gridInput); |
66 |
|
|
return overallSum; |
67 |
|
|
} |
68 |
|
|
|
69 |
|
|
/* (non-Javadoc) |
70 |
|
|
* @see appl.parallel.spmd.AbstractSPMDTask#init() |
71 |
|
|
*/ |
72 |
|
|
@Override |
73 |
|
|
public void init() { |
74 |
|
|
// TODO Auto-generated method stub |
75 |
|
|
|
76 |
|
|
} |
77 |
|
|
|
78 |
|
|
|
79 |
|
|
|
80 |
|
|
} |