/** * Apply the pseudo-inverse operation of a linear operator: y = (A^H A + lambda I)^-1 x * Does not check that the dimensions are consistent for the linear operator * * @param op linear operator * @param lambda regularization parameter * @param dst output data * @param src input data */ void linop_norm_inv_unchecked(const struct linop_s* op, float lambda, complex float* dst, const complex float* src) { operator_p_apply_unchecked(op->norm_inv, lambda, dst, src); }
void operator_p_iter(iter_op_data* _data, float rho, float* dst, const float* src) { auto data = CAST_DOWN(iter_op_p_op, _data); operator_p_apply_unchecked(data->op, rho, (complex float*)dst, (const complex float*)src); }