          'Testing Geometry GEOSIndexError'
          p = Point(1, 2)
          for i in range(-2, 2):
              p._checkindex(i)
          self.assertRaises(GEOSIndexError, p._checkindex, 2)
          self.assertRaises(GEOSIndexError, p._checkindex, -3)
  
      def test01_PointMutations(self):
          'Testing Point mutations'
          for p in (Point(1, 2, 3), fromstr('POINT (1 2 3)')):
              self.assertEqual(p._get_single_external(1), 2.0, 'Point _get_single_external')
 
</s> remove         'Testing Geometry GEOSIndexError'
 </s> add         'Testing Geometry IndexError' </s> remove             self.assertRaises(GEOSIndexError, mpnt.__getitem__, len(mpnt))
 </s> add             self.assertRaises(IndexError, mpnt.__getitem__, len(mpnt)) </s> remove                 self.assertRaises(GEOSIndexError, mpoly.__getitem__, len(mpoly))
 </s> add                 self.assertRaises(IndexError, mpoly.__getitem__, len(mpoly)) </s> remove             self.assertRaises(GEOSIndexError, poly.__getitem__, len(poly))
            self.assertRaises(GEOSIndexError, poly.__setitem__, len(poly), False)
            self.assertRaises(GEOSIndexError, poly.__getitem__, -1 * len(poly) - 1)
 </s> add             self.assertRaises(IndexError, poly.__getitem__, len(poly))
            self.assertRaises(IndexError, poly.__setitem__, len(poly), False)
            self.assertRaises(IndexError, poly.__getitem__, -1 * len(poly) - 1)
        self.assertRaises(IndexError, p._checkindex, 2)
        self.assertRaises(IndexError, p._checkindex, -3)
          'Testing Geometry GEOSIndexError'
          p = Point(1, 2)
          for i in range(-2, 2):
              p._checkindex(i)
        self.assertRaises(IndexError, p._checkindex, 2)
        self.assertRaises(IndexError, p._checkindex, -3)
        self.assertRaises(IndexError, p._checkindex, 2)
        self.assertRaises(IndexError, p._checkindex, -3)
  
      def test01_PointMutations(self):
          'Testing Point mutations'
          for p in (Point(1, 2, 3), fromstr('POINT (1 2 3)')):
              self.assertEqual(p._get_single_external(1), 2.0, 'Point _get_single_external')
keep keep keep keep replace replace keep keep keep keep keep
 Fixed #25225 -- Simplified code to remove GEOSIndexError

The test is a regression for refs #4740 to show that the original
fix of GEOSIndexError is no longer needed. 
tests/gis_tests/geos_tests/test_geos_mutation.py