当前位置: 首页>>代码示例>>C++>>正文


C++ alignment::empty方法代码示例

本文整理汇总了C++中alignment::empty方法的典型用法代码示例。如果您正苦于以下问题:C++ alignment::empty方法的具体用法?C++ alignment::empty怎么用?C++ alignment::empty使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在alignment的用法示例。


在下文中一共展示了alignment::empty方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。

示例1: n_seg

void
starling_read::
set_genome_align(const alignment& al) {

    assert(get_full_segment().genome_align().empty());
    assert(! al.empty());

    get_full_segment()._genome_align=al;
    const seg_id_t n_seg(apath_exon_count(al.path));
    if(n_seg<=1) return;

    // deal with segmented reads now:
    assert(! is_segmented());
    _segment_ptr.reset(new starling_segmented_read(n_seg));

    using namespace ALIGNPATH;

    seg_id_t seg_id(1);
    pos_t read_pos(0);
    pos_t ref_pos(al.pos);
    pos_t seg_start_read_pos(read_pos);
    pos_t seg_start_ref_pos(ref_pos);
    path_t seg_path;

    const unsigned as(al.path.size());
    for(unsigned i(0); i<as; ++i) {
        const path_segment& ps(al.path[i]);
        const pos_t last_read_pos(read_pos);
        if(is_segment_type_ref_length(ps.type)) ref_pos += ps.length;
        if(is_segment_type_read_length(ps.type)) read_pos += ps.length;

        if(ps.type!=SKIP) seg_path.push_back(ps);

        if(ps.type==SKIP || ((i+1)==as)) {
            const pos_t end_read_pos( (ps.type==SKIP) ?
                                      last_read_pos : read_pos );
            assert(end_read_pos>seg_start_read_pos);
            const unsigned size(end_read_pos-seg_start_read_pos);
            const read_segment rseg(size,seg_start_read_pos,this);
            _segment_ptr->set_segment(seg_id,rseg);
            alignment& new_al(get_segment(seg_id)._genome_align);
            new_al.path=seg_path;
            new_al.pos=seg_start_ref_pos;
            new_al.is_fwd_strand=al.is_fwd_strand;

            seg_id++;
            seg_start_read_pos=read_pos;
            seg_start_ref_pos=ref_pos;
            seg_path.clear();
        }
    }
}
开发者ID:DBHi-BiC,项目名称:isaac_variant_caller,代码行数:52,代码来源:starling_read.cpp


注:本文中的alignment::empty方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。