1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
use crate::infer::*;
use crate::internal::*;
#[derive(Debug, Clone, new, Hash)]
pub struct RmDims {
pub axes: Vec<usize>,
}
tract_linalg::impl_dyn_hash!(RmDims);
impl RmDims {
fn compute_shape<D: DimLike>(&self, input: &[D]) -> TVec<D> {
input
.iter()
.enumerate()
.filter(|(ix, _d)| !self.axes.contains(ix))
.map(|(_ix, d)| d.clone())
.collect()
}
fn eval_t<T: Datum>(&self, input: Arc<Tensor>) -> TractResult<TVec<Arc<Tensor>>> {
let shape = self.compute_shape(input.shape());
Ok(tvec![input.into_tensor().into_array::<T>()?.into_shape(&*shape)?.into_arc_tensor()])
}
}
impl Op for RmDims {
fn name(&self) -> Cow<str> {
"RmDims".into()
}
fn info(&self) -> TractResult<Vec<String>> {
Ok(vec![format!("axes: {:?}", self.axes)])
}
op_hir!();
not_a_typed_op!();
}
impl StatelessOp for RmDims {
fn eval(&self, mut inputs: TVec<Arc<Tensor>>) -> TractResult<TVec<Arc<Tensor>>> {
let input = args_1!(inputs);
dispatch_datum!(Self::eval_t(input.datum_type())(self, input))
}
}
impl InferenceRulesOp for RmDims {
fn rules<'r, 'p: 'r, 's: 'r>(
&'s self,
s: &mut Solver<'r>,
inputs: &'p [TensorProxy],
outputs: &'p [TensorProxy],
) -> InferenceResult {
check_output_arity(&outputs, 1)?;
s.equals(&outputs[0].datum_type, &inputs[0].datum_type)?;
s.equals(&outputs[0].rank, (&inputs[0].rank).bex() - self.axes.len() as i32)?;
for axis in &self.axes {
s.equals(&inputs[0].shape[*axis], 1.to_dim())?;
}
s.given(&inputs[0].shape, move |s, shape| {
let output_shape = self.compute_shape(&shape);
s.equals(&outputs[0].shape, output_shape)
})
}
fn to_typed(
&self,
_source: &InferenceModel,
node: &InferenceNode,
target: &mut TypedModel,
mapping: &HashMap<OutletId, OutletId>,
) -> TractResult<TVec<OutletId>> {
let mut wire = mapping[&node.inputs[0]];
let mut axes = self.axes.clone();
axes.sort();
for axis in axes.into_iter().rev() {
wire = target.wire_node(
format!("{}-axis-{}", node.name, axis),
AxisOp::Rm(axis),
&[wire],
)?[0];
}
Ok(tvec!(wire))
}
as_op!();
}