void Genome::copy(Genome *dest) const
{
copyDimensions(dest);
copySequence(dest);
copyTopSegments(dest);
copyBottomSegments(dest);
copyMetadata(dest);
}
value powell_3D_align (value sa, value sb, value sc, value ra, value rb, \
value rc, value mm, value go, value ge)
{
CAMLparam5(sa, sb, sc, mm, go);
CAMLxparam4(ge, ra, rb, rc);
struct seq *csa, *csb, *csc;
struct seq *cra, *crb, *crc;
misCost = Int_val(mm);
startInsert = Int_val(go);
continueInsert = Int_val(ge);
startDelete = startInsert;
continueDelete = continueInsert;
Seq_custom_val(csa,sa);
Seq_custom_val(csb,sb);
Seq_custom_val(csc,sc);
Seq_custom_val(cra,ra);
Seq_custom_val(crb,rb);
Seq_custom_val(crc,rc);
assert(misCost != 0 && startInsert>=0 && continueInsert>0);
copySequence (csa, Astr);
copySequence (csb, Bstr);
copySequence (csc, Cstr);
Alen=strlen(Astr);
Blen=strlen(Bstr);
Clen=strlen(Cstr);
setup();
CAMLreturn(Val_int(doUkk(cra, crb, crc)));
}
