Hallo zusammen,
ich habe mit der aktuellen THELI version 2.10.3 folgendes Problem (auch wie oben beschrieben):
- Code: Alles auswählen
+ /home/thomas/THELI/theli//bin/Linux_64//anisotropy -i /home/thomas/.theli/tmp//2016.08.27-02:29:45_science_iso_800_exp_120_B_1OFCC_ksb.cat2 -c rg 1.17236937946 2.12669561021 MAG_AUTO 11.452 16.8479 -o /home/thomas/.theli/tmp//2016.08.27-02:29:45_science_iso_800_exp_120_B_1OFCC_ksb.cat3 -j 0.5 -e 2.0
anisotropy V 1.3.0 (Sep 13 2016; 17:12:54)
ERROR(anisotropy): not enough stars after preselection
Fügt man im Code anisotropy.c einige VPRINTF(DEBUG,...) hinzu um zu sehen
welche Wert(e) nstars und e1,e2, Psm11 und Psm22, maxellipsq haben
- Code: Alles auswählen
...
VPRINTF(DEBUG, "i = %d, goodobject = %d\n", i, goodobject);
if (goodobject == 1) {
cl[i] = 1;
VPRINTF(DEBUG, "e1[%d] = %f, e2[%d] = %f\n", i, e1[i], i, e2[i]);
VPRINTF(DEBUG, "Psm11[%d] = %f, Psm22[%d] = %f\n", i, Psm11[i], i, Psm22[i]);
VPRINTF(DEBUG, "maxellipsq = %f\n", maxellipsq);
if (((e1[i] * e1[i] + e2[i] * e2[i]) <= maxellipsq) &&
Psm11[i] > 1.0e-08 && Psm22[i] > 1.0e-08) {
if (nstars == 0) {
istart = i;
istartcopy = i;
} else {
inext[ilast] = i;
inextcopy[ilast] = i;
}
ilast = i;
nstars++;
}
}
}
VPRINTF(DEBUG, "stars preselected\n");
VPRINTF(DEBUG, "nstars: %d, coeffs: %d\n", nstars, coeffs);
if (nstars < coeffs) {
syserr(FATAL, "anisotropy", "not enough stars after preselection");
}
...
Lässt man den anisotropy laufen:
- Code: Alles auswählen
number of conditions: 2
read input catalog
rg: 1.172369 - 2.126696
MAG_AUTO: 11.452000 - 16.847900
input keys read
input objects: 6154
i = 0, goodobject = 0
i = 1, goodobject = 1
e1[1] = 0.000000, e2[1] = 0.000000
Psm11[1] = 0.000000, Psm22[1] = 0.000000
maxellipsq = 1.000000
i = 2, goodobject = 1
e1[2] = 0.000000, e2[2] = 0.000000
Psm11[2] = 0.000000, Psm22[2] = 0.000000
maxellipsq = 1.000000
i = 3, goodobject = 1
e1[3] = 0.000000, e2[3] = 0.000000
Psm11[3] = 0.000000, Psm22[3] = 0.000000
maxellipsq = 1.000000
i = 4, goodobject = 1
e1[4] = 0.000000, e2[4] = 0.000000
Psm11[4] = 0.000000, Psm22[4] = 0.000000
maxellipsq = 1.000000
i = 5, goodobject = 1
e1[5] = 0.000000, e2[5] = 0.000000
Psm11[5] = 0.000000, Psm22[5] = 0.000000
maxellipsq = 1.000000
...
...
...
i = 6153, goodobject = 1
e1[6153] = 0.000000, e2[6153] = 0.000000
Psm11[6153] = 0.000000, Psm22[6153] = 0.000000
maxellipsq = 1.000000
stars preselected
nstars: 0, coeffs: 6
ERROR(anisotropy): not enough stars after preselection
Wenn ich mir die Tabelle in 2016.08.27-02:29:45_science_iso_800_exp_120_B_1OFCC_ksb.cat2 anschaue, dann kann man sehen
das alle Werte für e1,e2, Psm11 und Psm22 in jeder Zeile den Wert 0 haben, somit kann
nstars niemals hochgezählt wird. Kann mir jemand ein Tipp geben
was es mit den Werten e1, e2, Psm11, Psm22 auf sich hat und wo ich möglicherweise "unterwegs" etwas falsch gemacht bzw. vergessen habe!!
Danke & Viele Grüße
Thomas