99 #define h(i,j) (hit_tileElemAtNoStride(H,2,i,j))
104 #define xback(i,j) (hit_tileElemAtNoStride(xTraceback,2,i,j))
105 #define yback(i,j) (hit_tileElemAtNoStride(yTraceback,2,i,j))
133 int main(
int argc,
char *argv[]) {
207 printf(
"=== Iteration: %d/%d ===\n\n",iter,iterations);
244 printf(
"=== Result ===\n");
265 printf(
"\nInputset too big to print to stdout\n");
269 printf(
"\nSmall inputset, it has been rewinded:\n"
325 for (i=0;i<=
size[0];i++){
331 for (j=0;j<=
size[1];j++){
344 for(i=1; i<
size[0]+1; i++){
345 for(j=1; j<size[1]+1; j++){
354 h_t max = MAX4(diag,down,right,0,&idx);
361 }
else if ( idx == 1 ) {
364 }
else if ( idx == 2 ) {
386 for (j=0;j<=size[1];j++){
393 for (i=0;i<=size[0];i++){
400 for (j=0;j<=size[1];j++){
401 printf(
"%6.1f",(
double)
h(i,j));
412 for (j=0;j<=size[1];j++){
419 for (i=0;i<=size[0];i++){
426 for (j=0;j<=size[1];j++){
450 int pos=size[0]+size[1]+1;
453 while (i>0 && j>0 &&
h(i,j) > 0){
455 int movex = (i>
xback(i,j));
456 int movey = (j>
yback(i,j));
461 if (movex && movey) {
475 i=
xback(tempi,tempj);
476 j=
yback(tempi,tempj);
#define hit_tileElemAt(var, ndims,...)
void printProteinMatch(HitTile_aa_t p, int begin, int end)
void readProtein(ProteinFile *fprotein, HitTile_aa_t *protein, int psize)
void openProtein(ProteinFile *fprotein, char *filename)
#define hit_clockContinue(c)
#define hit_tileDomainAlloc(newVarP, baseType, numDims,...)
void input_parameters(int argc, char *argv[])
HitTile_trace_t yTraceback
#define hit_clockReset(c)
void initPAM(char *filename)
void printProtein(HitTile_aa_t p, int offset, int psize)
void hit_comInit(int *pargc, char **pargv[])
HitTile_trace_t xTraceback
int main(int argc, char *argv[])
#define hit_tileFree(var)
#define hit_clockGetSeconds(c)
void phase_read_sequences()
#define hit_clockStart(c)
void phase_comp_hmatrix()