Struct BasisProjector 

pub struct BasisProjector { /* private fields */ }

Sparse isometric lift from a symmetry-reduced basis to its parent basis.

Implementations

impl BasisProjector

pub fn from_operator( projector: &(impl LinearOperator + ?Sized), tolerance: f64, ) -> Result<Self>

Convert any explicit rectangular isometry into a reusable basis projector.

This is the narrow waist for reduced spaces whose columns are not representable by a one-dimensional orbit character, including higher-dimensional finite-group representations.

pub fn from_embedding<Reduced, Parent>( reduced: &Reduced, parent: &Parent, ) -> Result<Self>

where Reduced: Basis, Parent: Basis<State = Reduced::State>,

One-hot embedding of a selected basis into a compatible parent basis.

pub fn between<Reduced, Parent>( reduced: &Reduced, parent: &Parent, ) -> Result<Self>

where Reduced: Basis, Parent: Basis<State = Reduced::State>,

Isometric lift from any basis into an explicitly selected parent basis.

The reduced basis owns the symmetry-reduction convention through Basis::reduction_image. Iterating the explicit parent states makes this equally useful for built-in and runtime symmetry sectors, fixed particle subspaces, and unrestricted parent spaces.

pub fn from_general<Parent>(basis: &GeneralBasis<Parent>) -> Result<Self>

where Parent: Basis, Parent::State: Hash + Ord + 'static,

pub const fn source_dimension(&self) -> usize

pub const fn reduced_dimension(&self) -> usize

pub fn project( &self, parent: &[Complex64], reduced: &mut [Complex64], ) -> Result<()>

Apply the adjoint projector to a parent-space vector.

pub fn lifted(&self, reduced: &[Complex64]) -> Result<Vec<Complex64>>

pub fn projected(&self, parent: &[Complex64]) -> Result<Vec<Complex64>>

pub fn lift_batch( &self, reduced: &[Vec<Complex64>], ) -> Result<Vec<Vec<Complex64>>>

pub fn project_batch( &self, parent: &[Vec<Complex64>], ) -> Result<Vec<Vec<Complex64>>>

pub fn lift_density(&self, reduced: &[Complex64]) -> Result<Vec<Complex64>>

Lift a row-major reduced density matrix as P ρ P†.

pub fn symmetry_leakage_norm( &self, operator: &(impl LinearOperator + ?Sized), ) -> Result<f64>

Frobenius norm of (I - P P†) A P, evaluated one reduced column at a time. Zero means the parent-space operator preserves this symmetry sector; no parent-space square projector is formed.

pub fn preserves_operator_symmetry( &self, operator: &(impl LinearOperator + ?Sized), tolerance: f64, ) -> Result<bool>

Trait Implementations

impl Clone for BasisProjector

fn clone(&self) -> BasisProjector

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for BasisProjector

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl LinearOperator for BasisProjector

fn shape(&self) -> (usize, usize)

Return (rows, columns).

fn format(&self) -> MatrixFormat

Report the public materialization format.

fn apply(&self, input: &[Complex64], output: &mut [Complex64]) -> Result<()>

Compute output = A * input.

fn stored_triplets(&self) -> Result<Option<Vec<(usize, usize, Complex64)>>>

Return canonical stored nonzeros when the representation already owns them.

fn is_real(&self) -> bool

Whether this operator preserves real vectors exactly.

fn apply_real(&self, input: &[f64], output: &mut [f64]) -> Result<()>

Apply this operator to a real vector.

fn apply_transpose( &self, input: &[Complex64], output: &mut [Complex64], ) -> Result<()>

Apply the algebraic transpose without conjugating either operand.

fn apply_adjoint( &self, input: &[Complex64], output: &mut [Complex64], ) -> Result<()>

Apply the conjugate transpose.

fn shifted_solver( &self, _shift: f64, ) -> Result<Option<Box<dyn ShiftedLinearSolver>>>

Prepare a reusable solver for (A - shift I) x = b, when supported.