********************************************************************/
/*
- * RCS $Id: hhfullalignment-C.h 243 2011-05-31 13:49:19Z fabian $
+ * RCS $Id: hhfullalignment-C.h 315 2016-12-15 17:18:30Z fabian $
*/
// hhfullalignment.C
#include "hhhit.h"
#include "hhhalfalignment.h"
#endif
+//#include "new_new.h" /* memory tracking */
* @brief Build full alignment -> qa->s[k][h] and ta->s[k][h]
*/
int
-FullAlignment::Build(HMM& q, Hit& hit)
+FullAlignment::Build(HMM& q, Hit& hit, char zcError[])
{
int step;
char prev_state=MM, state=MM;
}
}
else {
- fprintf(stderr,
+ sprintf(zcError,
"+-------------------------------+\n"
"| both sequences truncated left |\n"
"+-------------------------------+\n"
}
}
else{
- fprintf(stderr,
+ sprintf(zcError,
"+--------------------------------+\n"
"| both sequences truncated right |\n"
"+--------------------------------+\n"
int iq=hit.i1; // match state counter for query HMM (displayed in consensus line)
int jt=hit.j1; // match state counter for template HMM (displayed in consensus line)
- lq = new(short unsigned int[qa->n+2]);
- lt = new(short unsigned int[ta->n+2]);
+ lq = new short unsigned int[qa->n+2];
+ lt = new short unsigned int[ta->n+2];
for (k=0; k<qa->n; k++) lq[k]=qa->l[k][hit.i1];
for (k=0; k<ta->n; k++) lt[k]=ta->l[k][hit.j1];