File: align.h

package info (click to toggle)
staden 2.0.0%2Bb11-5
  • links: PTS, VCS
  • area: main
  • in suites: trixie
  • size: 21,568 kB
  • sloc: ansic: 240,605; tcl: 65,360; cpp: 12,854; makefile: 11,201; sh: 2,952; fortran: 2,033; perl: 63; awk: 46
file content (91 lines) | stat: -rw-r--r-- 2,773 bytes parent folder | download | duplicates (5)
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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
#ifndef _ALIGN_H
#define _ALIGN_H

#include "os.h"
#include "misc.h"

/*
#ifdef __alpha
typedef long FastInt;
#else
typedef int FastInt;
#endif
*/
typedef int FastInt;
typedef FastInt align_int;

#define MAX_SCORE_MATRIX 30

/* Basic job number */
#define ALIGN_J_MASK    0x0F
#define ALIGN_J_SS	0
#define ALIGN_J_SV	1
#define ALIGN_J_BSV	2
#define ALIGN_J_SSH	3

/* How to expand the sequence up */
#define ALIGN_J_PADS	0x10

/* Which ends to penalise gaps for. Only applies to ALIGN_J_SSH */
#define ALIGN_GAP_S1	0x020
#define ALIGN_GAP_E1	0x040
#define ALIGN_GAP_S2	0x080
#define ALIGN_GAP_E2	0x100

/*
 * Initialise 128x128 weight matrix from our score matrix
 */
extern void init_W128(int **matrix,
		      char *order, int unknown);

extern int calign(void *seq1, void *seq2, int len1, int len2,
		  void *rseq1, void *rseq2, int *rlen1, int *rlen2,
		  int low_band, int high_band, int gap_open, int gap_extend,
		  int job, int is_protein, align_int *res);

extern void cdisplay(void *seq1, void *seq2, int len1, int len2, int job,
		     align_int *results, int pos1, int pos2);

extern void cexpand(void *seq1, void *seq2, int len1, int len2,
		    void *rseq1, void *rseq2, int *rlen1, int *rlen2,
		    int job, align_int *results);

extern int_f align_(void *seq1, int_f *len1, void *seq2, int_f *len2,
		    void *rseq1, int_f *rlen1, void *rseq2, int_f *rlen2,
		    int_f *low_band, int_f *high_band,
		    int_f *gap_open, int_f *gap_extend,
		    int_f *job, int_f *is_protein, int_fl seq1_l,
		    int_fl seq2_l, int_fl rseq1_l, int_fl rseq2_l);

extern FastInt align_ss(char *A, char *B, FastInt M, FastInt N,
			FastInt low, FastInt up, FastInt W[][128],
			FastInt G, FastInt H, FastInt *S,
			FastInt s1, FastInt s2, FastInt e1, FastInt e2);

extern FastInt align_ss2(char *A, char *B, FastInt M, FastInt N,
			FastInt low, FastInt up, FastInt W[][128],
			FastInt G, FastInt H, FastInt *S,
			FastInt s1, FastInt s2, FastInt e1, FastInt e2);

extern FastInt align_sv(char *A, FastInt (*B)[6], FastInt M,
			FastInt N, FastInt low, FastInt up,
			FastInt W[][128], FastInt G, FastInt H, FastInt *S,
			FastInt s1, FastInt s2, FastInt e1, FastInt e2);

extern FastInt balign_sv(char *A, FastInt (*B)[6], FastInt M,
			FastInt N, FastInt low, FastInt up,
			FastInt W[][128], FastInt G, FastInt H, FastInt *S,
			FastInt s1, FastInt s2, FastInt e1, FastInt e2);

extern void display_ss(char *A, char *B, FastInt M, FastInt N,
		       FastInt *S, FastInt AP, FastInt BP);

extern void display_ss2(char *A, char *B, FastInt M, FastInt N,
		       FastInt *S, FastInt AP, FastInt BP);

extern void display_sv(char *A, FastInt (*B)[6], FastInt M, FastInt N,
		       FastInt *S, FastInt AP, FastInt BP);

#include "align_lib_old.h"

#endif /* _ALIGN_H */