bartorch.linop.Extract

bartorch.linop.Extract#

class bartorch.linop.Extract(start, oshape, ishape)#

Take the block of oshape that starts at start, BART’s linop_extract.

The restriction operator: its adjoint puts the block back where it came from and leaves the rest zero. This is what indexing an operator with A[...] builds.

Parameters:
  • start (sequence of int) – Where the block begins, one offset per axis of ishape.

  • oshape (tuple of int) – The block and the shape it comes out of, C order.

  • ishape (tuple of int) – The block and the shape it comes out of, C order.

__init__(start, oshape, ishape)#

Methods

A(x, **kwargs)

A x, under deepinv's name.

A_adjoint(y, **kwargs)

A^H y, under deepinv's name, recorded for autograd.

A_adjoint_A(x, **kwargs)

A^H A x, under deepinv's name, recorded for autograd.

A_dagger(y, **kwargs)

The pseudo-inverse, under deepinv's name.

__init__(start, oshape, ishape)

adjoint(y[, out])

A^H y, without recording for autograd.

cogram()

A A^H as an operator.

conj()

conj(A): conjugate the input, apply, conjugate the output.

forward(x[, out])

A x, without recording for autograd.

gram()

A^H A as an operator.

normal(x[, out])

A^H A x.

opnorm()

The spectral norm, by BART's power iteration on A^H A.

pinv(y[, damp])

(A^H A + damp I)^-1 A^H y, the damped least-squares solution.

to_nonlinear()

The same operator as a NonlinearOperator.

Attributes

H

A^H, from BART's own adjoint constructor.

T

A^T, the adjoint without the conjugation, as conj(A).H.

codim_rank

How many axes the codomain has.

codim_shape

The codomain, under pyxu's name for it; the same as oshape.

codim_size

How many elements the codomain holds.

device

dim_rank

How many axes the domain has.

dim_shape

The domain, under pyxu's name for it; the same as ishape.

dim_size

How many elements the domain holds.

ishape

oshape