DNA* initSultan5(int& m, int& n)
{
SmartArrayPtr<double> *diag, *sdiag, *kap;
m = 5;
n = 6;
DNA* dna = new DNA(n, n - m);
diag = dna->getDiag();
sdiag = dna->getSDiag();
kap = dna->getKap();
diag->alloc(m);
sdiag->alloc(m);
kap->alloc(m);
for(int i = 0; i < 5; i ++) {
(*diag)[i] = 0.;
(*sdiag)[i] = 1.;
}
(*kap)[0] = 24.;
(*kap)[1] = 27.;
(*kap)[2] = 33.;
(*kap)[3] = 37.;
(*kap)[4] = 24.;
return dna;
}
DNA* initSultan3(int& m, int& n)
{
SmartArrayPtr<double> *diag, *sdiag, *kap;
m = 3;
n = 4;
DNA* dna = new DNA(n, n - m);
diag = dna->getDiag();
sdiag = dna->getSDiag();
kap = dna->getKap();
diag->alloc(3);
sdiag->alloc(3);
kap->alloc(3);
(*diag)[0] = 0.;
(*diag)[1] = 0.;
(*diag)[2] = 0.;
(*sdiag)[0] = 1.;
(*sdiag)[1] = 1.;
(*kap)[0] = 6.;
(*kap)[1] = 5.1;
(*kap)[2] = 3.9;
return dna;
}