Struct DiscreteFld 

pub struct DiscreteFld { /* private fields */ }

Order-independent accumulator for deriving a fragment-length distribution (and library format) from uniquely-mapped proper pairs.

Unlike FragmentLengthDistribution — which accumulates in log space and is required for the online phase (it folds in per-fragment forgetting mass and is read per fragment) — this stores plain integer counts per length, bucketed by orientation. Integer increments are commutative, so the tallies are independent of thread / chunk order; the FLD is then built once, deterministically (fixed length order), in finish. Keeping it a separate type means the online FLD’s float/log-space accumulation is never disturbed, and no runtime dispatch is needed.

Implementations

impl DiscreteFld

pub fn new(fld_max: usize) -> Self

Create an accumulator covering raw fragment lengths [0, fld_max].

pub fn add(&self, len: usize, is_fw: bool, mate_fw: bool)

Record one uniquely-mapped proper pair from its mate strands: fragment length, orientation bucket, and observed format (derived from is_fw / mate_fw, mirroring the RAD reader’s rad_frag_format so the mapping pass and the RAD-derive path agree). Works in sketch mode, where no precomputed LibraryFormat is available. Thread-safe and order-independent.

pub fn count(&self) -> u64

Total unique proper pairs seen across both orientation buckets.

pub fn finish( &self, fld_mean: f64, fld_sd: f64, ) -> (FragmentLengthDistribution, Option<LibraryFormat>)

Build the FLD from the majority-orientation bucket (deterministically, in fixed length order) and infer the library format from the format tally. fld_mean/fld_sd seed the prior. Returns the cached FLD and the detected format (None when no proper pairs were seen).

Trait Implementations

impl Debug for DiscreteFld

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

Formats the value using the given formatter.