def METHOD_NAME(self, method, device):
    """Test save density matrix for instruction"""
    backend = self.backend(method=method, device=device)
    # Stabilizer test circuit
    circ = QuantumCircuit(3)
    circ.h(0)
    circ.sdg(0)
    circ.cx(0, 1)
    circ.cx(0, 2)
    # Target statevector
    target = qi.DensityMatrix(circ)
    # Add save to circuit
    label = "state"
    circ.save_density_matrix(label=label)
    # Run
    result = backend.run(transpile(circ, backend, optimization_level=0), shots=1).result()
    self.assertTrue(result.success)
    simdata = result.data(0)
    self.assertIn(label, simdata)
    value = simdata[label]
    self.assertEqual(value, target)