     func testDimension ()   {
         XCTAssertEqual(MultiPoint<Coordinate2D>(precision: precision, coordinateReferenceSystem: crs).dimension, geometryDimension)
     }
+    
+    func testBoundary() {
+        let geometry = MultiPoint<Coordinate2D>(elements: [Point<Coordinate2D>(coordinate: (x: 1.0, y: 1.0)),Point<Coordinate2D>(coordinate: (x: 2.0, y: 2.0))], precision: precision, coordinateReferenceSystem: crs).boundary()
+        let expected = MultiPoint<Coordinate2D>(precision: precision, coordinateReferenceSystem: crs)  // Empty Set
+        
+        XCTAssertTrue(geometry == expected, "\(geometry) is not equal to \(expected)")
+    }
+    
+    func testBoundary_Empty() {
+        let geometry = MultiPoint<Coordinate2D>(precision: precision, coordinateReferenceSystem: crs).boundary()
+        let expected = MultiPoint<Coordinate2D>(precision: precision, coordinateReferenceSystem: crs)  // Empty Set
+        
+        XCTAssertTrue(geometry == expected, "\(geometry) is not equal to \(expected)")
+    }
 }
 
 // MARK: - Coordinate2DM, FixedPrecision, Cartesian -

    func testDimension ()   {
        XCTAssertEqual(MultiPoint<Coordinate3DM>(precision: precision, coordinateReferenceSystem: crs).dimension, geometryDimension)
    }
}

    func testDimension ()   {
        XCTAssertEqual(MultiPoint<Coordinate3D>(precision: precision, coordinateReferenceSystem: crs).dimension, geometryDimension)
    }
    
    func testBoundary() {
        let geometry = MultiPoint<Coordinate3D>(elements: [Point<Coordinate3D>(coordinate: (x: 1.0, y: 1.0, z: 1.0)),Point<Coordinate3D>(coordinate: (x: 2.0, y: 2.0, z: 2.0))], precision: precision, coordinateReferenceSystem: crs).boundary()
        let expected = MultiPoint<Coordinate3D>(precision: precision, coordinateReferenceSystem: crs)  // Empty Set
        
        XCTAssertTrue(geometry == expected, "\(geometry) is not equal to \(expected)")
    }
    
    func testBoundary_Empty() {
        let geometry = MultiPoint<Coordinate3D>(precision: precision, coordinateReferenceSystem: crs).boundary()
        let expected = MultiPoint<Coordinate3D>(precision: precision, coordinateReferenceSystem: crs)  // Empty Set
        
        XCTAssertTrue(geometry == expected, "\(geometry) is not equal to \(expected)")
    }
}

// MARK: - Coordinate3DM, FixedPrecision, Cartesian -
