parallel/test/SPMDTest_MultiGrid.java

package appl.parallel.test;

import java.awt.Rectangle;
import java.io.File;
import java.io.Serializable;
import java.util.Vector;

import junit.framework.Test;
import junit.framework.TestCase;
import junit.framework.TestSuite;

import org.apache.log4j.LogManager;
import org.apache.log4j.Logger;
import org.apache.log4j.PropertyConfigurator;

import schmitzm.data.WritableGrid;
import schmitzm.data.WritableGridRaster;
import schmitzm.geotools.io.GeoImportUtil;
import appl.parallel.ComputingResourceContainer;
import appl.parallel.client.ClientDataServer;
import appl.parallel.client.RemoteEventHandler;
import appl.parallel.data.splittable.SplittableLLProxyGrid;
import appl.parallel.services.GlobalDiscoveryService;
import appl.parallel.spmd.SPMDClientController;
import appl.parallel.spmd.split.SplittableGrid;
import appl.util.RasterMetaData;
import appl.util.RasterUtil;
import edu.bonn.xulu.plugin.io.grid.array.WritableGridArrayFactory;

/**
 * This is a minor modification of {@link SPMDTest}. It is only used to test the
 * Multi-Splittable functionality of the parallelization. For the test of all
 * other features and more detailed information see {@link SPMDTest}.
 * 
 * @author Dominik Appl
 */
public class SPMDTest_MultiGrid extends TestCase implements Serializable {

        transient private final Logger LOG = LogManager.getLogger(this.getClass()
                        .getName());

        transient private ClientDataServer client;
        transient private Vector<ComputingResourceContainer> computingResources;
        transient private SPMDClientController clientController;

        transient SplittableGrid gridInput;
        transient SplittableGrid remoteInput;
        transient SplittableGrid localOutput;
        transient SplittableGrid remoteOutput;

        /**
         * Constructor for SPMDTest.
         * 
         * @param name
         */
        public SPMDTest_MultiGrid(String name) {
                super(name);
        }

        /**
         * Returns the JUnit test suite that implements the <b>SPMDTest</b>
         * definition.
         */
        public static Test suite() {
                TestSuite sPMDTest = new TestSuite("SPMDTest");
//              sPMDTest.setArbiter(DefaultTestArbiter.INSTANCE).setId(
//                              "E873D66A70977ED2A98AAC10B54511DB");
//
//              sPMDTest.addTest(new SPMDTest("testSPMD").setId(
//                              "E873D66A70977ED2E5E8ED20B54511DB").setTestInvocationId(
//                              "E873D66A70977ED2EF7E3C00B54511DB"));
                return sPMDTest;
        }

        /**
         * @see junit.framework.TestCase#setUp()
         */
        protected void setUp() throws Exception {
                PropertyConfigurator.configure("xululog4j.cfg");
                computingResources = new GlobalDiscoveryService()
                                .getRemoteResources();

                for (ComputingResourceContainer computingRes : computingResources) {
                        System.out.println("DiscoveredRessource " + computingRes);
                }
                LOG.info("Discovery finished...");
                RemoteEventHandler handler = new RemoteEventHandler();
                client = new ClientDataServer(handler);
                clientController = new SPMDClientController(computingResources,null, client,handler);
                // ////////// Initialise data \\\\\\\\\\\\\
                WritableGridRaster baseGrid = GeoImportUtil
                                .readGridRasterFromArcInfoASCII(new File(
                                                "../Xulu-Data/minigrid2.arc")); // small data
//               "../Xulu-Data/large_clue_data/sc1gr0.0")); //large data
                gridInput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
                                new RasterMetaData(baseGrid));
                localOutput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
                                new RasterMetaData(baseGrid));
                remoteInput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
                                new RasterMetaData(baseGrid));
                remoteOutput = new SplittableLLProxyGrid(
                                new WritableGridArrayFactory(), new RasterMetaData(baseGrid));
                RasterUtil.copyInto(baseGrid, gridInput);
                RasterUtil.copyInto(baseGrid, remoteInput);
                //RasterUtil.copyInto(baseGrid, localOutput);
                //RasterUtil.copyInto(baseGrid, remoteOutput);
        }

        /**
         * @see junit.framework.TestCase#tearDown()
         */
        protected void tearDown() throws Exception {
                client.close();
                clientController.close();
                // GlobalDiscoveryService.getInstance().stopServices();
        }

        /**
         * TestSPMD
         * 
         * @throws Exception
         */
        public void testSPMD() throws Exception {
                int steps = 4;
                int neighborhoodRange = 3;
                System.out.println("Starting step 1: ");
                                
                /***********************************************************************
                 * START LOCAL ALGORITHM ******************** The simple Algorithm will
                 * calculate the average of the neighborhood und write the output to the
                 * coresponing cell in the output cell
                 **********************************************************************/

                Rectangle partition = new Rectangle(0, 0, gridInput.getWidth(),
                                gridInput.getHeight());
                if (gridInput.getWidth() < 20)
                        RasterUtil.printGrid(gridInput, 5, 1, "Input");
                
                for (int step = 0; step < steps; step++) {
//                       the overallsum is simply the sum over all elements
                        float overallSum = 0f;
                        long localStartTime = System.currentTimeMillis();
                        for (int y = 0; y < gridInput.getHeight(); y++)
                                for (int x = 0; x < gridInput.getWidth(); x++) {
                                        float tmp = gridInput.getRasterSampleAsFloat(x, y);
                                        if (!Float.isNaN(tmp))
                                                overallSum += tmp;
                                        // the local sum is simply the sum over all elements in the
                                        // neighborhood
                                        float localSum = 0;
                                        // number of cells over which the sum is calculated
                                        int noOfCells = 0;
                                        // for each cell: calculate the sum of all neighbors
                                        for (int y2 = y - neighborhoodRange; y2 <= y
                                                        + neighborhoodRange; y2++)
                                                for (int x2 = x - neighborhoodRange; x2 <= x
                                                                + neighborhoodRange; x2++)
                                                        // check if the coordinates are valid (inside the
                                                        // grid and not NaN)
                                                        if (partition.contains(x2, y2))
                                                                if (!Float.isNaN(gridInput
                                                                                .getRasterSampleAsFloat(x2, y2))) {
                                                                        localSum += gridInput
                                                                                        .getRasterSampleAsFloat(x2, y2);
                                                                        noOfCells++;
                                                                }

                                        localOutput.setRasterSample(localSum / (float) (noOfCells),
                                                        x, y);

                                }

                        System.out.print("Finished step " + (step+1) + " in "
                                        + ((System.currentTimeMillis() - localStartTime)) + " ms");
                        localStartTime = System.currentTimeMillis();
                        System.out.println(" with sum (local calculated): " + overallSum);
                        //copy values into results from output to input (for next step)
                        RasterUtil.copyInto(localOutput, gridInput);
                 }
                 if(gridInput.getWidth()<20)
                 RasterUtil.printGrid(localOutput, 5, 1, "Local Output");

                
                 
                 
                 /***********************************************************************
                 *                                                      START REMOTE ALGORITHM
                 *  ******************** All the same..now remotely********************/
                 
                 System.out
                                .println("******************* starting remote execution *****************");
                long remoteStartTime = System.currentTimeMillis();
                clientController.setNeighborhoodRange(neighborhoodRange);
//              initialize as MulitSplittable
                WritableGrid[] testSplittable = {remoteInput,remoteOutput};
                clientController.addToMultiDataSplitControl(testSplittable, "gridData");
                
//              clientController.addToSplitControl(remoteInput, "inputGrid");
//              clientController.addToSplitControl(remoteOutput, "outputGrid");
                System.out.println("After "
                                + ((System.currentTimeMillis() - remoteStartTime))
                                + "ms: finished with adding to splitcontroll");
                
                //********************** START STEP CONTROL *********************
                for (int step = 0; step < steps; step++) {
                        if(step > 0){
                                long updateStartTime = System.currentTimeMillis();
                                clientController.updateNeighborhood(remoteInput);
                                System.out.println("Update took:" + (System.currentTimeMillis()-updateStartTime));
                        }
                        Object[] sum;
                        try {
                                sum = clientController.runSPMDModelTask(
                                                new AverageNeighborhoodTestTask_MultiGrid(), neighborhoodRange);
                                //add results
                                float overallSum = 0;
                                for (int i = 0; i < sum.length; i++)
                                        overallSum += (Float) sum[i];
                                System.out.println("finished step "
                                                + (step+1) + " in "
                                                + ((System.currentTimeMillis() - remoteStartTime))
                                                + "ms " + " with sum "+
                                                + overallSum);
                                remoteStartTime = System.currentTimeMillis();
                                
                        } catch (Throwable e) {
                                // TODO Auto-generated catch block
                                e.printStackTrace();
                        }
                }
                //***************** END STEP CONTROL *****************************
                
                //After execution:
                // merge partitions and check if the local and remotely calculated
                // versions are equal
                remoteStartTime = System.currentTimeMillis();
                clientController.mergePartition(remoteOutput);
                System.out.println("Merging took another "
                                + ((System.currentTimeMillis() - remoteStartTime))
                                + "ms");
                // Thread.currentThread().sleep(10000);
                 if(remoteInput.getWidth()<20)
                         RasterUtil.printGrid(remoteOutput, 5, 1, "Remote Output");
                assertTrue(RasterUtil.checkEqual(remoteOutput, localOutput, true));
        }
}
1	mojays	2	package appl.parallel.test;
2
3			import java.awt.Rectangle;
4			import java.io.File;
5			import java.io.Serializable;
6			import java.util.Vector;
7
8			import junit.framework.Test;
9			import junit.framework.TestCase;
10			import junit.framework.TestSuite;
11
12			import org.apache.log4j.LogManager;
13			import org.apache.log4j.Logger;
14			import org.apache.log4j.PropertyConfigurator;
15
16			import schmitzm.data.WritableGrid;
17			import schmitzm.data.WritableGridRaster;
18			import schmitzm.geotools.io.GeoImportUtil;
19			import appl.parallel.ComputingResourceContainer;
20	alfonx	78	import appl.parallel.client.ClientDataServer;
21	mojays	2	import appl.parallel.client.RemoteEventHandler;
22			import appl.parallel.data.splittable.SplittableLLProxyGrid;
23			import appl.parallel.services.GlobalDiscoveryService;
24			import appl.parallel.spmd.SPMDClientController;
25			import appl.parallel.spmd.split.SplittableGrid;
26			import appl.util.RasterMetaData;
27			import appl.util.RasterUtil;
28	alfonx	78	import edu.bonn.xulu.plugin.io.grid.array.WritableGridArrayFactory;
29	mojays	2
30			/**
31			* This is a minor modification of {@link SPMDTest}. It is only used to test the
32			* Multi-Splittable functionality of the parallelization. For the test of all
33			* other features and more detailed information see {@link SPMDTest}.
34			*
35			* @author Dominik Appl
36			*/
37			public class SPMDTest_MultiGrid extends TestCase implements Serializable {
38
39			transient private final Logger LOG = LogManager.getLogger(this.getClass()
40			.getName());
41
42			transient private ClientDataServer client;
43			transient private Vector<ComputingResourceContainer> computingResources;
44			transient private SPMDClientController clientController;
45
46			transient SplittableGrid gridInput;
47			transient SplittableGrid remoteInput;
48			transient SplittableGrid localOutput;
49			transient SplittableGrid remoteOutput;
50
51			/**
52			* Constructor for SPMDTest.
53			*
54			* @param name
55			*/
56			public SPMDTest_MultiGrid(String name) {
57			super(name);
58			}
59
60			/**
61			* Returns the JUnit test suite that implements the <b>SPMDTest</b>
62			* definition.
63			*/
64			public static Test suite() {
65			TestSuite sPMDTest = new TestSuite("SPMDTest");
66			// sPMDTest.setArbiter(DefaultTestArbiter.INSTANCE).setId(
67			// "E873D66A70977ED2A98AAC10B54511DB");
68			//
69			// sPMDTest.addTest(new SPMDTest("testSPMD").setId(
70			// "E873D66A70977ED2E5E8ED20B54511DB").setTestInvocationId(
71			// "E873D66A70977ED2EF7E3C00B54511DB"));
72			return sPMDTest;
73			}
74
75			/**
76			* @see junit.framework.TestCase#setUp()
77			*/
78			protected void setUp() throws Exception {
79			PropertyConfigurator.configure("xululog4j.cfg");
80			computingResources = new GlobalDiscoveryService()
81			.getRemoteResources();
82
83			for (ComputingResourceContainer computingRes : computingResources) {
84			System.out.println("DiscoveredRessource " + computingRes);
85			}
86			LOG.info("Discovery finished...");
87			RemoteEventHandler handler = new RemoteEventHandler();
88			client = new ClientDataServer(handler);
89			clientController = new SPMDClientController(computingResources,null, client,handler);
90			// ////////// Initialise data \\\\\\\\\\\\\
91			WritableGridRaster baseGrid = GeoImportUtil
92			.readGridRasterFromArcInfoASCII(new File(
93			"../Xulu-Data/minigrid2.arc")); // small data
94			// "../Xulu-Data/large_clue_data/sc1gr0.0")); //large data
95			gridInput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
96			new RasterMetaData(baseGrid));
97			localOutput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
98			new RasterMetaData(baseGrid));
99			remoteInput = new SplittableLLProxyGrid(new WritableGridArrayFactory(),
100			new RasterMetaData(baseGrid));
101			remoteOutput = new SplittableLLProxyGrid(
102			new WritableGridArrayFactory(), new RasterMetaData(baseGrid));
103			RasterUtil.copyInto(baseGrid, gridInput);
104			RasterUtil.copyInto(baseGrid, remoteInput);
105			//RasterUtil.copyInto(baseGrid, localOutput);
106			//RasterUtil.copyInto(baseGrid, remoteOutput);
107			}
108
109			/**
110			* @see junit.framework.TestCase#tearDown()
111			*/
112			protected void tearDown() throws Exception {
113			client.close();
114			clientController.close();
115			// GlobalDiscoveryService.getInstance().stopServices();
116			}
117
118			/**
119			* TestSPMD
120			*
121			* @throws Exception
122			*/
123			public void testSPMD() throws Exception {
124			int steps = 4;
125			int neighborhoodRange = 3;
126			System.out.println("Starting step 1: ");
127
128			/***********************************************************************
129			* START LOCAL ALGORITHM ******************** The simple Algorithm will
130			* calculate the average of the neighborhood und write the output to the
131			* coresponing cell in the output cell
132			**********************************************************************/
133
134			Rectangle partition = new Rectangle(0, 0, gridInput.getWidth(),
135			gridInput.getHeight());
136			if (gridInput.getWidth() < 20)
137			RasterUtil.printGrid(gridInput, 5, 1, "Input");
138
139			for (int step = 0; step < steps; step++) {
140			// the overallsum is simply the sum over all elements
141			float overallSum = 0f;
142			long localStartTime = System.currentTimeMillis();
143			for (int y = 0; y < gridInput.getHeight(); y++)
144			for (int x = 0; x < gridInput.getWidth(); x++) {
145			float tmp = gridInput.getRasterSampleAsFloat(x, y);
146			if (!Float.isNaN(tmp))
147			overallSum += tmp;
148			// the local sum is simply the sum over all elements in the
149			// neighborhood
150			float localSum = 0;
151			// number of cells over which the sum is calculated
152			int noOfCells = 0;
153			// for each cell: calculate the sum of all neighbors
154			for (int y2 = y - neighborhoodRange; y2 <= y
155			+ neighborhoodRange; y2++)
156			for (int x2 = x - neighborhoodRange; x2 <= x
157			+ neighborhoodRange; x2++)
158			// check if the coordinates are valid (inside the
159			// grid and not NaN)
160			if (partition.contains(x2, y2))
161			if (!Float.isNaN(gridInput
162			.getRasterSampleAsFloat(x2, y2))) {
163			localSum += gridInput
164			.getRasterSampleAsFloat(x2, y2);
165			noOfCells++;
166			}
167
168			localOutput.setRasterSample(localSum / (float) (noOfCells),
169			x, y);
170
171			}
172
173			System.out.print("Finished step " + (step+1) + " in "
174			+ ((System.currentTimeMillis() - localStartTime)) + " ms");
175			localStartTime = System.currentTimeMillis();
176			System.out.println(" with sum (local calculated): " + overallSum);
177			//copy values into results from output to input (for next step)
178			RasterUtil.copyInto(localOutput, gridInput);
179			}
180			if(gridInput.getWidth()<20)
181			RasterUtil.printGrid(localOutput, 5, 1, "Local Output");
182
183
184
185
186			/***********************************************************************
187			* START REMOTE ALGORITHM
188			* ****************** All the same..now remotely******************/
189
190			System.out
191			.println("***************** starting remote execution ***************");
192			long remoteStartTime = System.currentTimeMillis();
193			clientController.setNeighborhoodRange(neighborhoodRange);
194			// initialize as MulitSplittable
195			WritableGrid[] testSplittable = {remoteInput,remoteOutput};
196			clientController.addToMultiDataSplitControl(testSplittable, "gridData");
197
198			// clientController.addToSplitControl(remoteInput, "inputGrid");
199			// clientController.addToSplitControl(remoteOutput, "outputGrid");
200			System.out.println("After "
201			+ ((System.currentTimeMillis() - remoteStartTime))
202			+ "ms: finished with adding to splitcontroll");
203
204			//******************** START STEP CONTROL *******************
205			for (int step = 0; step < steps; step++) {
206			if(step > 0){
207			long updateStartTime = System.currentTimeMillis();
208			clientController.updateNeighborhood(remoteInput);
209			System.out.println("Update took:" + (System.currentTimeMillis()-updateStartTime));
210			}
211			Object[] sum;
212			try {
213			sum = clientController.runSPMDModelTask(
214			new AverageNeighborhoodTestTask_MultiGrid(), neighborhoodRange);
215			//add results
216			float overallSum = 0;
217			for (int i = 0; i < sum.length; i++)
218			overallSum += (Float) sum[i];
219			System.out.println("finished step "
220			+ (step+1) + " in "
221			+ ((System.currentTimeMillis() - remoteStartTime))
222			+ "ms " + " with sum "+
223			+ overallSum);
224			remoteStartTime = System.currentTimeMillis();
225
226			} catch (Throwable e) {
227			// TODO Auto-generated catch block
228			e.printStackTrace();
229			}
230			}
231			//*************** END STEP CONTROL ***************************
232
233			//After execution:
234			// merge partitions and check if the local and remotely calculated
235			// versions are equal
236			remoteStartTime = System.currentTimeMillis();
237			clientController.mergePartition(remoteOutput);
238			System.out.println("Merging took another "
239			+ ((System.currentTimeMillis() - remoteStartTime))
240			+ "ms");
241			// Thread.currentThread().sleep(10000);
242			if(remoteInput.getWidth()<20)
243			RasterUtil.printGrid(remoteOutput, 5, 1, "Remote Output");
244			assertTrue(RasterUtil.checkEqual(remoteOutput, localOutput, true));
245			}
246			}