Struct Problem

Source

struct Problem {Show 14 fields
    dims: (usize, usize, usize),
    vs: [f64; 3],
    w_mag: Vec<f64>,
    w_phs: Vec<f64>,
    t_mag: Vec<f64>,
    t_phs: Vec<f64>,
    d: Vec<f64>,
    wg: [Vec<f64>; 3],
    csf: Vec<usize>,
    drp: f64,
    drn: f64,
    cf_mhz: f64,
    lambda: f64,
    lambda_csf: f64,
}

Expand description

The arrays of one run, on the padded grid.

Fields§

§dims: (usize, usize, usize)§vs: [f64; 3]§w_mag: Vec<f64>

R2′ data weight wei_mag·m and field data weight m (M_cat’s two halves)

§w_phs: Vec<f64>§t_mag: Vec<f64>

R2′.*mask where R2′ > 0, and the local field in Hz (target_’s two halves)

§t_phs: Vec<f64>§d: Vec<f64>

Dipole kernel, MEDI’s dipole_kernel(..., 'kspace')

§wg: [Vec<f64>; 3]

MEDI edge mask of the magnitude, 1 = not an edge (wG; wG_r2p is the same)

§csf: Vec<usize>§drp: f64§drn: f64§cf_mhz: f64§lambda: f64§lambda_csf: f64

Implementations§

Source §

impl Problem

Source

fn dconv(&self, fft: &mut Fft3dWorkspace, buf: &mut [Complex64])

Dconv_p in place: buf ← (CF/1e6)·ifftn(D.*fftn(buf))

Source

fn add_csf(&self, x: &[f64], out: &mut [f64])

λ_CSF·LT(LT(|x|)) added to out, with LT(v) = Mask_CSF.*(v − mean(v(Mask_CSF)))

Source

fn add_reg( &self, ws: &mut Workspace, x: &[f64], diff: &[f64], w_single: &[Vec<f64>; 3], w_tot: &[Vec<f64>; 3], out: &mut [f64], )

The regulariser of one source: λ·∇ᵀ(W_single ∇x) + 2λ·∇ᵀ(W_tot ∇(x − y)) + CSF(x) with x − y given as diff, added to out.

Source

fn apply( &self, ws: &mut Workspace, dx: &[f64], pos: bool, w_single: &[Vec<f64>; 3], w_tot: &[Vec<f64>; 3], out: &mut [f64], )

The CG operator of one source (Ap or An): data term plus regulariser at dx with the other source’s step zero.

Source

fn gradients( &self, ws: &mut Workspace, xp: &[f64], xn: &[f64], w_pos: &[Vec<f64>; 3], w_neg: &[Vec<f64>; 3], w_tot: &[Vec<f64>; 3], bp: &mut [f64], bn: &mut [f64], )

bp = fidelity_derivative_p + reg_p and bn likewise, at (xp, xn).

Auto Trait Implementations§

§

impl UnwindSafe for Problem

Blanket Implementations§

Source §

impl<T> Any for T
where T: 'static + ?Sized,

Source §

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more

Source §

impl<T> Borrow<T> for T
where T: ?Sized,

Source §

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more

Source §

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source §

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more

§

impl<T> Downcast for T
where T: Any,

§

fn into_any(self: Box<T>) -> Box<dyn Any>

Converts Box<dyn Trait> (where Trait: Downcast) to Box<dyn Any>, which can then be downcast into Box<dyn ConcreteType> where ConcreteType implements Trait.

§

fn into_any_rc(self: Rc<T>) -> Rc<dyn Any>

Converts Rc<Trait> (where Trait: Downcast) to Rc<Any>, which can then be further downcast into Rc<ConcreteType> where ConcreteType implements Trait.

§

fn as_any(&self) -> &(dyn Any + 'static)

Converts &Trait (where Trait: Downcast) to &Any. This is needed since Rust cannot generate &Any’s vtable from &Trait’s.

§

fn as_any_mut(&mut self) -> &mut (dyn Any + 'static)

Converts &mut Trait (where Trait: Downcast) to &Any. This is needed since Rust cannot generate &mut Any’s vtable from &mut Trait’s.

§

impl<T> DowncastSend for T
where T: Any + Send,

§

fn into_any_send(self: Box<T>) -> Box<dyn Any + Send>

Converts Box<Trait> (where Trait: DowncastSend) to Box<dyn Any + Send>, which can then be downcast into Box<ConcreteType> where ConcreteType implements Trait.

§

impl<T> DowncastSync for T
where T: Any + Send + Sync,

§

fn into_any_sync(self: Box<T>) -> Box<dyn Any + Send + Sync>

Converts Box<Trait> (where Trait: DowncastSync) to Box<dyn Any + Send + Sync>, which can then be downcast into Box<ConcreteType> where ConcreteType implements Trait.

§

fn into_any_arc(self: Arc<T>) -> Arc<dyn Any + Send + Sync> ⓘ

Converts Arc<Trait> (where Trait: DowncastSync) to Arc<Any>, which can then be downcast into Arc<ConcreteType> where ConcreteType implements Trait.

Source §

impl<T> From<T> for T

Source §

fn from(t: T) -> T

Returns the argument unchanged.

Source §

impl<T, U> Into for T
where U: From<T>,

Source §

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source §

impl<T> IntoEither for T

Source §

fn into_either(self, into_left: bool) -> Either<Self, Self> ⓘ

Converts self into a Left variant of Either<Self, Self> if into_left is true. Converts self into a Right variant of Either<Self, Self> otherwise. Read more

Source §

fn into_either_with<F>(self, into_left: F) -> Either<Self, Self> ⓘ
where F: FnOnce(&Self) -> bool,

Converts self into a Left variant of Either<Self, Self> if into_left(&self) returns true. Converts self into a Right variant of Either<Self, Self> otherwise. Read more

§