ff_energy 0.4.1

fuzzyfold's nearest neighbor free energy evaluations.
Documentation 
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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59

use crate::parameters::parameterset::StackParams;
use crate::parameters::parameterset::ExtendedStackParams;
use crate::parameters::parameterset::{P, E};

pub static STACKPARAMS_EN37: StackParams = [
    /* [cl] [i]:   AT     TA     CG     GC     GT     TG */
    /* [AT] */ [  -90,  -100,  -130,  -140,    10,    70],
    /* [TA] */ [ -100,   -60,  -150,  -130,    30,    40],
    /* [CG] */ [ -130,  -150,  -220,  -180,   -30,   -50],
    /* [GC] */ [ -140,  -130,  -180,  -220,   -60,    10],
    /* [GT] */ [   10,    30,   -30,   -60,   120,    70],
    /* [TG] */ [   70,    40,   -50,    10,    70,    50],
];

pub static STACKPARAMS_ENTH: StackParams = [
    /* [cl] [i]:   AT     TA     CG     GC     GT     TG */
    /* [AT] */ [ -530,  -720,  -580,  -780,  -170,   370],
    /* [TA] */ [ -720,  -670,  -990,  -850,   -90,   -90],
    /* [CG] */ [ -580,  -990,  -980,  -750,  -240,  -430],
    /* [GC] */ [ -780,  -850,  -750,  -790,  -240,   430],
    /* [GT] */ [ -170,   -90,  -240,  -240,   500,   800],
    /* [TG] */ [  340,   -90,  -430,   430,   800,  -670],
];

/// The parameters embedded into whatever dimension is currently 
/// used by the ViennaRNA-style stacking tables.
pub const STACK_EN37: ExtendedStackParams = {
    let mut full: ExtendedStackParams = [[None; E]; E];

    let mut i = 0;
    while i < P {
        let mut j = 0;
        while j < P {
            full[i][j] = Some(STACKPARAMS_EN37[i][j]);
            j += 1;
        }
        i += 1;
    }

    full
};

/// The parameters embedded into whatever dimension is currently 
/// used by the ViennaRNA-style stacking tables.
pub const STACK_ENTH: ExtendedStackParams = {
    let mut full: ExtendedStackParams = [[None; E]; E];

    let mut i = 0;
    while i < P {
        let mut j = 0;
        while j < P {
            full[i][j] = Some(STACKPARAMS_ENTH[i][j]);
            j += 1;
        }
        i += 1;
    }

    full
};