A RetroSearch Logo

Home - News ( United States | United Kingdom | Italy | Germany ) - Football scores

Search Query:

Showing content from http://www.ncbi.nlm.nih.gov/IEB/ToolBox/CPP_DOC/doxyhtml/nw__spliced__aligner32_8cpp_source.html below:

NCBI C++ ToolKit: src/algo/align/nw/nw_spliced_aligner32.cpp Source File

75  x_Set

(

'A'

,

'A'

, 0x47);

76  x_Set

(

'A'

,

'G'

, 0x47);

77  x_Set

(

'A'

,

'T'

, 0x57);

78  x_Set

(

'A'

,

'C'

, 0x67);

79  x_Set

(

'G'

,

'A'

, 0x33);

80  x_Set

(

'G'

,

'G'

, 0x33);

81  x_Set

(

'G'

,

'T'

, 0x70);

82  x_Set

(

'G'

,

'C'

, 0x34);

83  x_Set

(

'T'

,

'A'

, 0x00);

84  x_Set

(

'T'

,

'G'

, 0x03);

85  x_Set

(

'T'

,

'T'

, 0x50);

86  x_Set

(

'T'

,

'C'

, 0x24);

87  x_Set

(

'C'

,

'A'

, 0x00);

88  x_Set

(

'C'

,

'G'

, 0x03);

89  x_Set

(

'C'

,

'T'

, 0x50);

90  x_Set

(

'C'

,

'C'

, 0x24);

100  m_Matrix

[(

unsigned

char)c1][(

unsigned char

)c2] =

val

;

113

m_Wd1(GetDefaultWd1()),

114

m_Wd2(GetDefaultWd2())

122  const char

* seq2,

size_t

len2)

124

m_Wd1(GetDefaultWd1()),

125

m_Wd2(GetDefaultWd2())

135

m_Wd1(GetDefaultWd1()),

136

m_Wd2(GetDefaultWd2())

146  switch

(splice_type) {

152

eInvalidSpliceTypeIndex,

165  const Uint1

where (len_dif < 0? 0: 1);

166  const size_t

shift (

abs

(len_dif) / 2);

167  const size_t

band (

abs

(len_dif) + 2*(

max

(

data

->m_len1,

data

->m_len2)/20 + 1));

175  const TScore

cds_penalty_extra = -2e-6;

177  const size_t

N1 =

data

->m_len1 + 1;

178  const size_t

N2 =

data

->m_len2 + 1;

180

vector<TScore> stl_rowV (N2), stl_rowF (N2);

181  TScore

* rowV = &stl_rowV[0];

182  TScore

* rowF = &stl_rowF[0];

190  const char

* seq1 =

m_Seq1

+

data

->m_offset1 - 1;

191  const char

* seq2 =

m_Seq2

+

data

->m_offset2 - 1;

195  bool

bFreeGapLeft1 =

data

->m_esf_L1 &&

data

->m_offset1 == 0;

196  bool

bFreeGapRight1 =

data

->m_esf_R1 &&

199  bool

bFreeGapLeft2 =

data

->m_esf_L2 &&

data

->m_offset1 == 0;

200  bool

bFreeGapRight2 =

data

->m_esf_R2 &&

205  TScore

wg1 = wgleft1, ws1 = wsleft1;

221

jAllDonors[

st

] = &stl_jAllDonors[

st

*N2];

222

vAllDonors[

st

] = &stl_vAllDonors[

st

*N2];

230

vector<size_t> stl_del_start(N2);

231  size_t

* del_start = &stl_del_start[0];

234  const Uint1

* dnr_acc_matrix = g_dnr_acc_matrix.GetMatrix();

239  for

(k = 0; k < N2; k++) {

240

rowV[k] = rowF[k] = kInfMinus;

246  if

(cds_start < cds_stop) {

247

cds_start -=

data

->m_offset1;

248

cds_stop -=

data

->m_offset1;

251  size_t i

, j = 0, k0;

253  for

(

i

= 0;

i

< N1; ++

i

, j = 0) {

255

V =

i

> 0? (V0 += wsleft2) : 0;

259

ci =

i

> 0? seq1[

i

]:

'N'

;

262

jTail[

st

] = jHead[

st

] = 0;

263

vBestDonor[

st

] = kInfMinus;

266  if

(

i

== N1 - 1 && bFreeGapRight1) {

274  unsigned char

d1 = seq2[1], d2 = seq2[2];

275  Uint1

dnr_type = 0xF0 & dnr_acc_matrix[(size_t(d1)<<8)|d2];

278

jAllDonors[

st

][jTail[

st

]] = j;

279  if

(dnr_type & (0x10 <<

st

)) {

280

vAllDonors[

st

][jTail[

st

]] =

285

vAllDonors[

st

][jTail[

st

]] = V +

m_Wd2

;

291  if

(cds_start <=

i

&&

i

< cds_stop) {

293  if

(

i

!= 0 || ! bFreeGapLeft1) {

294

ws1 += cds_penalty_extra;

296  if

(j != 0 || ! bFreeGapLeft2) {

297

ws2 += cds_penalty_extra;

301  for

(j = 1; j < N2; ++j, ++k) {

303  G

= pV[j] + sm[ci][(

unsigned

char)seq2[j]];

315  if

(j == N2 - 1 && bFreeGapRight2) {

351  if

(jTail[

st

] > jHead[

st

]) {

353  if

(vAllDonors[

st

][jHead[

st

]] > vBestDonor[

st

]) {

354

vBestDonor[

st

] = vAllDonors[

st

][jHead[

st

]];

355

jBestDonor[

st

] = jAllDonors[

st

][jHead[

st

]];

364  unsigned char

c1 = seq2[j-1], c2 = seq2[j];

365  Uint1

acc_mask = 0x0F & dnr_acc_matrix[(size_t(c1)<<8)|c2];

367  if

(acc_mask & (0x01 <<

st

)) {

376

dnr_pos = k0 + jBestDonor[

st

];

383

dnr_pos = k0 + jBestDonor[

st

];

392

backtrace_matrix[k] =

static_cast<unsigned int>

(

type

);

396  unsigned char

d1 = seq2[j+1], d2 = seq2[j+2];

397  Uint1

dnr_mask = 0xF0 & dnr_acc_matrix[(size_t(d1)<<8)|d2];

399  if

( dnr_mask & (0x10 <<

st

) ) {

403  if

(V > vBestDonor[

st

]) {

404

jAllDonors[

st

][jTail[

st

]] = j;

405

vAllDonors[

st

][jTail[

st

]] = V;

410  if

(v > vBestDonor[

st

]) {

411

jAllDonors[

st

][jTail[

st

]] = j;

412

vAllDonors[

st

][jTail[

st

]] = v;

419  if

(v > vBestDonor[

st

]) {

420

jAllDonors[

st

][jTail[

st

]] = j;

421

vAllDonors[

st

][jTail[

st

]] = v;

455  const size_t

N1 =

data

->m_len1 + 1;

456  const size_t

N2 =

data

->m_len2 + 1;

458  data

->m_transcript.clear();

459  data

->m_transcript.reserve(N1 + N2);

461  size_t

k = N1*N2 - 1;

462  size_t

i1 =

data

->m_offset1 +

data

->m_len1 - 1;

463  size_t

i2 =

data

->m_offset2 +

data

->m_len2 - 1;

465  const Uint4

mask_jump = 0x3FFFFFFF;

466  const Uint4

mask_type = ~mask_jump;

470  Uint4

Key = backtrace_matrix[k];

479  Uint4

k2 = (Key & mask_jump);

492  Uint8

kdel = k / N2, k2del = k2 / N2;

493  Uint8

kmod = k % N2, k2mod = k2 % N2;

498  else if

(kmod == k2mod) {

509  for

(; k > k2; k -= decr) {

510  data

->m_transcript.push_back(ts);

530  size_t

start1,

size_t

start2)

const 544  const size_t

dim = transcript.size();

548  const char

* p1 =

m_Seq1

+ start1;

549  const char

* p2 =

m_Seq2

+ start2;

556  switch

( transcript[0] ) {

564

state1 = 1; state2 = 0; score +=

m_Wg

;

568

state1 = 0; state2 = 0;

572

state1 = 0; state2 = 1; score +=

m_Wg

;

581  for

(

size_t i

= 0;

i

< dim; ++

i

) {

592  unsigned char

c1 = *p1;

593  unsigned char

c2 = *p2;

602  if

(state1 != 1) score +=

m_Wg

;

603

state1 = 1; state2 = 0;

610  if

(state2 != 1) score +=

m_Wg

;

611

state1 = 0; state2 = 1;

634

state1 = 2; state2 = 0;

648  for

(

Int8 i

= dim - 1;

i

>= 0; --

i

) {

658  for

(

Int8 i

= dim - 1;

i

>= 0; --

i

) {

668  for

(

size_t i

= 0;

i

< dim; ++

i

) {

678  for

(

size_t i

= 0;

i

< dim; ++

i

) {

const Uint1 * GetMatrix() const

void x_Set(char c1, char c2, Uint1 val)

void SetBand(size_t band)

virtual ETranscriptSymbol x_GetDiagTS(size_t i1, size_t i2) const

SNCBIFullScoreMatrix m_ScoreMatrix

void SetShift(Uint1 where, size_t offset)

virtual TScore x_Align(CNWAligner::SAlignInOut *data)

vector< ETranscriptSymbol > TTranscript

virtual TScore x_Align(CNWAligner::SAlignInOut *data)

TScore m_Wi[splice_type_count_32]

void x_DoBackTrace(const Uint4 *backtrace_matrix, CNWAligner::SAlignInOut *data)

virtual TScore ScoreFromTranscript(const TTranscript &transcript, size_t start1=kMax_UInt, size_t start2=kMax_UInt) const

const size_t splice_type_count_32

static TScore GetDefaultWi(unsigned char splice_type)

#define NCBI_THROW(exception_class, err_code, message)

Generic macro to throw an exception, given the exception class, error code and message string.

uint8_t Uint1

1-byte (8-bit) unsigned integer

uint32_t Uint4

4-byte (32-bit) unsigned integer

int64_t Int8

8-byte (64-bit) signed integer

uint64_t Uint8

8-byte (64-bit) unsigned integer

#define END_NCBI_SCOPE

End previously defined NCBI scope.

#define BEGIN_NCBI_SCOPE

Define ncbi namespace.

const unsigned char g_nwspl32_acceptor[splice_type_count_32][2]

const unsigned char g_nwspl32_donor[splice_type_count_32][2]

#define NCBI_FSM_DIM

Recommended approach: unpack and index directly.

int TNCBIScore

data types

static SLJIT_INLINE sljit_ins st(sljit_gpr r, sljit_s32 d, sljit_gpr x, sljit_gpr b)

const char g_msg_InconsistentArguments[]

const char g_msg_InvalidTranscriptSymbol[]

const char g_msg_InvalidSpliceTypeIndex[]

const char g_msg_InvalidBacktraceData[]

int g(Seg_Gsm *spe, Seq_Mtf *psm, Thd_Gsm *tdg)


RetroSearch is an open source project built by @garambo | Open a GitHub Issue

Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo

HTML: 3.2 | Encoding: UTF-8 | Version: 0.7.4