bartorch.linop.hstack

Contents

bartorch.linop.hstack#

bartorch.linop.hstack(ops, axis=0)#

Split the input between the operators and add up what they return.

[A B] @ [x; y] = A x + B y, which is what pylops calls HStack: the operators share a codomain, and their domains agree except along axis.

BART has no constructor for this, and it needs none: the adjoint of stacking the codomains is stacking the domains, so this is concatenate([op.H for op in ops], axis).H and stays one BART operator.

Parameters:
  • ops (sequence of LinearOperator) – At least one. All must have the same codomain.

  • axis (int) – Which axis of the domain grows, as an index into it.