1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
// -----------------------------------------------------------------------------------------------------
// Copyright (c) 2006-2020, Knut Reinert & Freie Universität Berlin
// Copyright (c) 2016-2020, Knut Reinert & MPI für molekulare Genetik
// This file may be used, modified and/or redistributed under the terms of the 3-clause BSD-License
// shipped with this file and also available at: https://github.com/seqan/seqan3/blob/master/LICENSE.md
// -----------------------------------------------------------------------------------------------------
#include <gtest/gtest.h>
#include <seqan3/alignment/configuration/align_config_edit.hpp>
#include <seqan3/alignment/configuration/align_config_method.hpp>
TEST(align_cfg_edit, is_hamming)
{
auto scheme = seqan3::get<seqan3::align_cfg::scoring_scheme>(seqan3::align_cfg::edit_scheme).scheme;
for (unsigned i = 0; i < decltype(scheme)::matrix_size; ++i)
{
for (unsigned j = 0; j < decltype(scheme)::matrix_size; ++j)
{
if (i == j)
EXPECT_EQ((scheme.score(seqan3::assign_rank_to(i, seqan3::dna15{}),
seqan3::assign_rank_to(j, seqan3::dna15{}))), 0);
else
EXPECT_EQ((scheme.score(seqan3::assign_rank_to(i, seqan3::dna15{}),
seqan3::assign_rank_to(j, seqan3::dna15{}))), -1);
}
}
}
TEST(align_cfg_edit, is_simple_gap)
{
using seqan3::get;
auto gap_cost_scheme = get<seqan3::align_cfg::gap_cost_affine>(seqan3::align_cfg::edit_scheme);
EXPECT_EQ(gap_cost_scheme.extension_score, -1);
EXPECT_EQ(gap_cost_scheme.open_score, 0);
}
|
