r.mapcalc: tests for seed option and -s flag (testing https://trac.osgeo.org/grass/changeset/61353)

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@61355 15284696-431f-4ddb-bdfa-cd5b030d7da7

raster/r.mapcalc/testsuite/test_r_mapcalc.py

@@ -0,0 +1,156 @@
+import grass.gunittest
+from grass.gunittest.gmodules import SimpleModule
+
+cell_seed_500 = """\
+north: 20
+south: 10
+east: 25
+west: 15
+rows: 10
+cols: 10
+2 12 183 135 117 176 138 117 182 40
+157 70 195 1 30 125 122 74 108 104
+163 146 82 164 67 63 60 3 147 193
+151 173 144 173 128 77 141 6 144 193
+180 56 14 121 124 107 46 107 90 60
+177 57 12 104 178 13 167 104 163 187
+55 133 27 114 195 60 78 35 49 11
+55 138 25 105 32 58 47 18 139 32
+24 74 36 71 131 124 87 161 160 154
+136 45 48 146 9 182 69 12 35 19
+"""
+
+dcell_seed_600 = """\
+north: 20
+south: 10
+east: 25
+west: 15
+rows: 10
+cols: 10
+130.790433856418332 101.3319248101491041 33.5781271447787759 37.4064724824657944 98.2794723130458152 73.9118866262841863 185.9530433718733775 74.5210037729812882 166.1178416001017695 90.9915650902159996
+109.2478664232956334 25.6499350759712215 150.9024447059825036 125.7544119036241312 66.7235333366722614 167.9375729129454271 123.1009291055983965 12.0922254083554606 59.389026967287819 113.2843489528100349
+40.0044184023145277 135.8273774212801186 71.6737798852435049 191.6223505280372876 4.1546013811569615 143.3082794522489962 177.043829294835632 115.0300571162354402 141.8985452774071234 127.8949967123638061
+93.2842559637482793 9.7471880052423856 118.1216452002055632 158.1474162140586088 67.2957262519499437 3.6524546350146849 147.0965842667525862 37.060628529871579 47.3408278816968959 66.2219633495724054
+175.5638637866295539 67.1399023507611901 162.2058392782793703 198.1586789345953719 36.474049475167746 49.2589028048889617 112.1169663235969836 22.0227597984432535 95.9169228571662131 86.7470895014531322
+93.5401613888204935 193.7821104138942587 193.8286564351004699 3.2623643889134684 94.6955247357847725 25.7099391122614307 155.592251526775442 25.3392337002970294 48.3979699868663005 99.6836079482556272
+104.16296861365457 190.7865884377180805 6.2841805474238619 49.3731395705159528 100.1903962703459285 116.927654961282343 19.8626348109264264 40.9693022766258466 81.6500759554420057 169.2220572316770131
+118.8112518721558217 55.8955021401724039 112.9150308215961331 62.6399760484719081 85.400498505854145 191.0144187084912062 124.2128169358724534 167.9341741649760706 170.6149695243870781 158.3034517206661462
+130.0453795775294736 64.1996403829061535 62.9317494959142465 175.1909990236256931 122.9624852869890361 79.9546265736285733 9.6594013716963367 114.0611338072915544 11.9371167643030809 186.9121199748369122
+3.2891990250261536 30.9245408751958379 46.4021422454598991 104.2378950097200203 47.424093232347019 73.4801303522840499 22.4778583078695213 132.870185207462697 48.1666164169167388 100.5504714442693057
+"""
+
+fcell_seed_700 = """\
+north: 20
+south: 10
+east: 25
+west: 15
+rows: 10
+cols: 10
+146.756378 192.682159 2.644822 147.270462 62.178818 192.668198 94.320778 107.710426 98.319664 114.444504
+12.995321 18.026272 151.590958 5.249451 197.266708 103.663635 115.424088 28.01062 78.555168 62.912098
+164.053619 154.652039 98.536011 44.601639 85.322289 168.383957 44.93845 128.62262 89.910591 107.242188
+111.182487 63.080284 177.791473 47.439354 42.451859 72.396568 170.597778 170.622742 141.88858 105.126854
+120.76828 148.581085 42.124866 56.432236 164.652176 98.094009 60.741329 66.286987 187.847427 160.120056
+50.530689 179.090652 138.114014 138.629211 193.147903 172.861481 133.72728 108.720459 103.508438 28.81559
+39.653179 101.948265 35.744762 25.570076 78.767021 154.600616 144.907684 82.370148 116.378654 18.218494
+35.587288 66.534409 65.744408 186.476959 137.081116 151.379272 48.261463 8.323328 130.432739 53.346546
+152.67189 15.512391 146.049072 185.276245 34.417141 127.522453 124.54998 52.08218 167.141342 87.771118
+69.0522 43.57811 63.15279 68.677063 74.202805 97.429077 167.123199 19.892767 120.593437 190.960815
+"""
+
+
+class TestRandFunction(grass.gunittest.TestCase):
+
+    # TODO: replace by unified handing of maps
+    to_remove = []
+
+    @classmethod
+    def setUpClass(cls):
+        cls.use_temp_region()
+        cls.runModule('g.region', n=20, s=10, e=25, w=15, res=1)
+
+    @classmethod
+    def tearDownClass(cls):
+        cls.del_temp_region()
+        cls.runModule('g.remove', rast=','.join(cls.to_remove))
+
+    def rinfo_contains_number(self, raster, number):
+        """Test that r.info stdandard output for raster contains a given number
+        """
+        rinfo = SimpleModule('r.info', map=raster)
+        self.runModule(rinfo)
+        self.assertIn(str(number), rinfo.outputs.stdout)
+
+    def test_seed_not_required(self):
+        """Test that seed is not required when rand() is not used"""
+        self.assertModule('r.mapcalc', expression='nonrand_cell = 200')
+        self.to_remove.append('nonrand_cell')
+
+    def test_seed_required(self):
+        """Test that seed is required when rand() is used"""
+        self.assertModuleFail('r.mapcalc', expression='rand_x = rand(1, 200)')
+        # TODO: assert map not exists but it would be handy here
+
+    def test_seed_cell(self):
+        """Test given seed with CELL against reference map"""
+        seed = 500
+        self.runModule('r.in.ascii', input='-', stdin=cell_seed_500,
+                       output='rand_cell_ref')
+        self.to_remove.append('rand_cell_ref')
+        self.assertModule('r.mapcalc', seed=seed,
+                          expression='rand_cell = rand(1, 200)')
+        self.to_remove.append('rand_cell')
+        # this assert is using r.mapcalc but we are testing different
+        # functionality than used by assert
+        self.assertRastersNoDifference(actual='rand_cell',
+                                       reference='rand_cell_ref',
+                                       precision=0)
+        self.rinfo_contains_number('rand_cell', seed)
+
+    def test_seed_dcell(self):
+        """Test given seed with DCELL against reference map"""
+        seed = 600
+        self.runModule('r.in.ascii', input='-', stdin=dcell_seed_600,
+                       output='rand_dcell_ref')
+        self.to_remove.append('rand_dcell_ref')
+        self.assertModule('r.mapcalc', seed=seed,
+                          expression='rand_dcell = rand(1.0, 200.0)')
+        self.to_remove.append('rand_dcell')
+        # this assert is using r.mapcalc but we are testing different
+        # functionality than used by assert
+        self.assertRastersNoDifference(actual='rand_dcell',
+                                       reference='rand_dcell_ref',
+                                       precision=0.00000000000001)
+        self.rinfo_contains_number('rand_dcell', seed)
+
+    def test_seed_fcell(self):
+        """Test given seed with FCELL against reference map"""
+        seed = 700
+        self.runModule('r.in.ascii', input='-', stdin=fcell_seed_700,
+                       output='rand_fcell_ref')
+        self.to_remove.append('rand_fcell_ref')
+        self.assertModule('r.mapcalc', seed=seed,
+                          expression='rand_fcell = rand(float(1), 200)')
+        self.to_remove.append('rand_fcell')
+        # this assert is using r.mapcalc but we are testing different
+        # functionality than used by assert
+        self.assertRastersNoDifference(actual='rand_fcell',
+                                       reference='rand_fcell_ref',
+                                       precision=0.000001)
+        self.rinfo_contains_number('rand_fcell', seed)
+
+    def test_auto_seed(self):
+        """Test that two runs with -s does not give same maps"""
+        self.assertModule('r.mapcalc', flags='s',
+                          expression='rand_auto_1 = rand(1., 2)')
+        self.to_remove.append('rand_auto_1')
+        self.assertModule('r.mapcalc', flags='s',
+                          expression='rand_auto_2 = rand(1., 2)')
+        self.to_remove.append('rand_auto_2')
+        self.assertRastersDifference(
+            'rand_auto_1', 'rand_auto_2',
+            statistics=dict(min=-1, max=1, mean=0),
+            precision=0.5)  # low precision, we have few cells
+
+if __name__ == '__main__':
+    grass.gunittest.test()