suggestion

pennylane/_qubit_device.py

@@ -381,7 +381,7 @@ def execute_new(self, circuit, **kwargs):
         self.apply(circuit.operations, rotations=circuit.diagonalizing_gates, **kwargs)
 
         # generate computational basis samples
-        if self.shots is not None or circuit.is_sampled:
+        if self.shots is not None:
             self._samples = self.generate_samples()
 
         # compute the required statistics