quantrs2_ml/torchquantum/gates/two_qubit/special/
tqch_traits.rs1use crate::error::{MLError, Result};
14use crate::torchquantum::{
15 CType, NParamsEnum, OpHistoryEntry, TQDevice, TQModule, TQOperator, TQParameter, WiresEnum,
16};
17use scirs2_core::ndarray::{Array2, ArrayD, IxDyn};
18
19use super::types::TQCH;
20
21impl Default for TQCH {
22 fn default() -> Self {
23 Self::new()
24 }
25}
26
27impl TQModule for TQCH {
28 fn forward(&mut self, _qdev: &mut TQDevice) -> Result<()> {
29 Err(MLError::InvalidConfiguration(
30 "Use apply() instead of forward() for operators".to_string(),
31 ))
32 }
33 fn parameters(&self) -> Vec<TQParameter> {
34 Vec::new()
35 }
36 fn n_wires(&self) -> Option<usize> {
37 Some(2)
38 }
39 fn set_n_wires(&mut self, _n_wires: usize) {}
40 fn is_static_mode(&self) -> bool {
41 self.static_mode
42 }
43 fn static_on(&mut self) {
44 self.static_mode = true;
45 }
46 fn static_off(&mut self) {
47 self.static_mode = false;
48 }
49 fn name(&self) -> &str {
50 "CH"
51 }
52}
53
54impl TQOperator for TQCH {
55 fn num_wires(&self) -> WiresEnum {
56 WiresEnum::Fixed(2)
57 }
58 fn num_params(&self) -> NParamsEnum {
59 NParamsEnum::Fixed(0)
60 }
61 fn get_matrix(&self, _params: Option<&[f64]>) -> Array2<CType> {
62 let scale = 1.0 / std::f64::consts::SQRT_2;
63 Array2::from_shape_vec(
64 (4, 4),
65 vec![
66 CType::new(1.0, 0.0),
67 CType::new(0.0, 0.0),
68 CType::new(0.0, 0.0),
69 CType::new(0.0, 0.0),
70 CType::new(0.0, 0.0),
71 CType::new(1.0, 0.0),
72 CType::new(0.0, 0.0),
73 CType::new(0.0, 0.0),
74 CType::new(0.0, 0.0),
75 CType::new(0.0, 0.0),
76 CType::new(scale, 0.0),
77 CType::new(scale, 0.0),
78 CType::new(0.0, 0.0),
79 CType::new(0.0, 0.0),
80 CType::new(scale, 0.0),
81 CType::new(-scale, 0.0),
82 ],
83 )
84 .unwrap_or_else(|_| Array2::eye(4).mapv(|x| CType::new(x, 0.0)))
85 }
86 fn apply(&mut self, qdev: &mut TQDevice, wires: &[usize]) -> Result<()> {
87 self.apply_with_params(qdev, wires, None)
88 }
89 fn apply_with_params(
90 &mut self,
91 qdev: &mut TQDevice,
92 wires: &[usize],
93 _params: Option<&[f64]>,
94 ) -> Result<()> {
95 if wires.len() < 2 {
96 return Err(MLError::InvalidConfiguration(
97 "CH gate requires exactly 2 wires".to_string(),
98 ));
99 }
100 let matrix = self.get_matrix(None);
101 qdev.apply_two_qubit_gate(wires[0], wires[1], &matrix)?;
102 if qdev.record_op {
103 qdev.record_operation(OpHistoryEntry {
104 name: "ch".to_string(),
105 wires: wires.to_vec(),
106 params: None,
107 inverse: self.inverse,
108 trainable: false,
109 });
110 }
111 Ok(())
112 }
113 fn has_params(&self) -> bool {
114 false
115 }
116 fn trainable(&self) -> bool {
117 false
118 }
119 fn inverse(&self) -> bool {
120 self.inverse
121 }
122 fn set_inverse(&mut self, inverse: bool) {
123 self.inverse = inverse;
124 }
125}