INNER CODE UNIT · Rust
map_params_from_sketch
bluenote-1577/skani · src/chain.rs:88
pub fn map_params_from_sketch <'a>(
ref_sketch: &Sketch,
amino_acid: bool,
command_params: &CommandParams,
model_opt: &'a Option<GBDT>
) -> MapParams<'a> {
let max_gap_length = if amino_acid{D_MAX_GAP_LENGTH_AAI} else {D_MAX_GAP_LENGTH};
let anchor_score = if amino_acid{D_ANCHOR_SCORE_AAI} else {D_ANCHOR_SCORE_ANI};
let min_anchors = if amino_acid{D_MIN_ANCHORS_AAI} else {D_MIN_ANCHORS_ANI};
let min_length_cover = if amino_acid{MIN_LENGTH_COVER_AAI} else {MIN_LENGTH_COVER};
let fragment_length = fragment_length_formula(ref_sketch.total_sequence_length, amino_acid);
let length_cutoff = fragment_length;
let mut frac_cover_cutoff = command_params.min_aligned_frac;
if frac_cover_cutoff < 0.{
if amino_acid {
frac_cover_cutoff = D_FRAC_COVER_CUTOFF_AA.parse::<f64>().unwrap()/100.;
} else {
frac_cover_cutoff = D_FRAC_COVER_CUTOFF.parse::<f64>().unwrap()/100.;