void
timbre::similarity_raw(
musly_track* track,
musly_track** tracks,
int length,
float* similarities)
{
// map seed track to gaussian structure
gaussian g0;
g0.mu = &track[track_mu];
g0.covar = &track[track_covar];
g0.covar_logdet = &track[track_logdet];
// create the temporary buffer required for the Jensen-Shannon divergence
musly_track* tmp_t = track_alloc();
gaussian tmp;
tmp.mu = &tmp_t[track_mu];
tmp.covar = &tmp_t[track_covar];
tmp.covar_logdet = &tmp_t[track_logdet];
// iterate over all musly_tracks to compute the Jensen-Shannon divergence
for (int i = 0; i < length; i++) {
gaussian gi;
musly_track* track1 = tracks[i];
gi.mu = &track1[track_mu];
gi.covar = &track1[track_covar];
gi.covar_logdet = &track1[track_logdet];
similarities[i] = gs.jensenshannon(g0, gi, tmp);
}
delete[] tmp_t;
}
int tracks_in_pattern_alloc(struct xmp_module *mod, int num)
{
int i;
for (i = 0; i < mod->chn; i++) {
int t = num * mod->chn + i;
int rows = mod->xxp[num]->rows;
if (track_alloc(mod, t, rows) < 0)
return -1;
mod->xxp[num]->index[i] = t;
}
return 0;
}
int
mandelellis::similarity(
musly_track* track,
musly_trackid seed_trackid,
musly_track** tracks,
musly_trackid* trackids,
int length,
float* similarities)
{
if ((length <= 0) || !track || ! tracks || !similarities) {
return -1;
}
// map seed track to gaussian structure
gaussian g0;
g0.mu = &track[track_mu];
g0.covar = &track[track_covar];
g0.covar_inverse = &track[track_covar_inverse];
// create the temporary buffer required for the Kullback-Leibler divergence
musly_track* tmp_t = track_alloc();
gaussian tmp;
tmp.mu = &tmp_t[track_mu];
tmp.covar = &tmp_t[track_covar];
tmp.covar_inverse = &tmp_t[track_covar_inverse];
// iterate over all musly_tracks to compute the Kullback-Leibler divergence
for (int i = 0; i < length; i++) {
gaussian gi;
musly_track* track1 = tracks[i];
gi.mu = &track1[track_mu];
gi.covar = &track1[track_covar];
gi.covar_inverse = &track1[track_covar_inverse];
similarities[i] = gs.symmetric_kullbackleibler(g0, gi, tmp);
}
delete[] tmp_t;
return 0;
}
static int mgt_load(struct module_data *m, HIO_HANDLE *f, const int start)
{
struct xmp_module *mod = &m->mod;
struct xmp_event *event;
int i, j;
int ver;
int sng_ptr, seq_ptr, ins_ptr, pat_ptr, trk_ptr, smp_ptr;
int sdata[64];
LOAD_INIT();
hio_read24b(f); /* MGT */
ver = hio_read8(f);
hio_read32b(f); /* MCS */
set_type(m, "Megatracker MGT v%d.%d", MSN(ver), LSN(ver));
mod->chn = hio_read16b(f);
hio_read16b(f); /* number of songs */
mod->len = hio_read16b(f);
mod->pat = hio_read16b(f);
mod->trk = hio_read16b(f);
mod->ins = mod->smp = hio_read16b(f);
hio_read16b(f); /* reserved */
hio_read32b(f); /* reserved */
sng_ptr = hio_read32b(f);
seq_ptr = hio_read32b(f);
ins_ptr = hio_read32b(f);
pat_ptr = hio_read32b(f);
trk_ptr = hio_read32b(f);
smp_ptr = hio_read32b(f);
hio_read32b(f); /* total smp len */
hio_read32b(f); /* unpacked trk size */
hio_seek(f, start + sng_ptr, SEEK_SET);
hio_read(mod->name, 1, 32, f);
seq_ptr = hio_read32b(f);
mod->len = hio_read16b(f);
mod->rst = hio_read16b(f);
mod->bpm = hio_read8(f);
mod->spd = hio_read8(f);
hio_read16b(f); /* global volume */
hio_read8(f); /* master L */
hio_read8(f); /* master R */
for (i = 0; i < mod->chn; i++) {
hio_read16b(f); /* pan */
}
MODULE_INFO();
/* Sequence */
hio_seek(f, start + seq_ptr, SEEK_SET);
for (i = 0; i < mod->len; i++)
mod->xxo[i] = hio_read16b(f);
/* Instruments */
if (instrument_init(mod) < 0)
return -1;
hio_seek(f, start + ins_ptr, SEEK_SET);
for (i = 0; i < mod->ins; i++) {
int c2spd, flags;
if (subinstrument_alloc(mod, i , 1) < 0)
return -1;
hio_read(mod->xxi[i].name, 1, 32, f);
sdata[i] = hio_read32b(f);
mod->xxs[i].len = hio_read32b(f);
mod->xxs[i].lps = hio_read32b(f);
mod->xxs[i].lpe = mod->xxs[i].lps + hio_read32b(f);
hio_read32b(f);
hio_read32b(f);
c2spd = hio_read32b(f);
c2spd_to_note(c2spd, &mod->xxi[i].sub[0].xpo, &mod->xxi[i].sub[0].fin);
mod->xxi[i].sub[0].vol = hio_read16b(f) >> 4;
hio_read8(f); /* vol L */
hio_read8(f); /* vol R */
mod->xxi[i].sub[0].pan = 0x80;
flags = hio_read8(f);
mod->xxs[i].flg = flags & 0x03 ? XMP_SAMPLE_LOOP : 0;
mod->xxs[i].flg |= flags & 0x02 ? XMP_SAMPLE_LOOP_BIDIR : 0;
mod->xxi[i].sub[0].fin += 0 * hio_read8(f); // FIXME
hio_read8(f); /* unused */
hio_read8(f);
hio_read8(f);
hio_read8(f);
hio_read16b(f);
hio_read32b(f);
hio_read32b(f);
mod->xxi[i].nsm = !!mod->xxs[i].len;
mod->xxi[i].sub[0].sid = i;
D_(D_INFO "[%2X] %-32.32s %04x %04x %04x %c V%02x %5d\n",
i, mod->xxi[i].name,
mod->xxs[i].len, mod->xxs[i].lps, mod->xxs[i].lpe,
mod->xxs[i].flg & XMP_SAMPLE_LOOP_BIDIR ? 'B' :
mod->xxs[i].flg & XMP_SAMPLE_LOOP ? 'L' : ' ',
mod->xxi[i].sub[0].vol, c2spd);
}
/* PATTERN_INIT - alloc extra track*/
if (pattern_init(mod) < 0)
return -1;
D_(D_INFO "Stored tracks: %d", mod->trk);
/* Tracks */
for (i = 1; i < mod->trk; i++) {
int offset, rows;
uint8 b;
hio_seek(f, start + trk_ptr + i * 4, SEEK_SET);
offset = hio_read32b(f);
hio_seek(f, start + offset, SEEK_SET);
rows = hio_read16b(f);
if (track_alloc(mod, i, rows) < 0)
return -1;
//printf("\n=== Track %d ===\n\n", i);
for (j = 0; j < rows; j++) {
uint8 note, f2p;
b = hio_read8(f);
j += b & 0x03;
note = 0;
event = &mod->xxt[i]->event[j];
if (b & 0x04)
note = hio_read8(f);
if (b & 0x08)
event->ins = hio_read8(f);
if (b & 0x10)
event->vol = hio_read8(f);
if (b & 0x20)
event->fxt = hio_read8(f);
if (b & 0x40)
event->fxp = hio_read8(f);
if (b & 0x80)
f2p = hio_read8(f);
if (note == 1)
event->note = XMP_KEY_OFF;
else if (note > 11) /* adjusted to play codeine.mgt */
event->note = note + 1;
/* effects */
if (event->fxt < 0x10)
/* like amiga */ ;
else switch (event->fxt) {
case 0x13:
case 0x14:
case 0x15:
case 0x17:
case 0x1c:
case 0x1d:
case 0x1e:
event->fxt = FX_EXTENDED;
event->fxp = ((event->fxt & 0x0f) << 4) |
(event->fxp & 0x0f);
break;
default:
event->fxt = event->fxp = 0;
}
/* volume and volume column effects */
if ((event->vol >= 0x10) && (event->vol <= 0x50)) {
event->vol -= 0x0f;
continue;
}
switch (event->vol >> 4) {
case 0x06: /* Volume slide down */
event->f2t = FX_VOLSLIDE_2;
event->f2p = event->vol - 0x60;
break;
case 0x07: /* Volume slide up */
event->f2t = FX_VOLSLIDE_2;
event->f2p = (event->vol - 0x70) << 4;
break;
case 0x08: /* Fine volume slide down */
event->f2t = FX_EXTENDED;
event->f2p = (EX_F_VSLIDE_DN << 4) |
(event->vol - 0x80);
break;
case 0x09: /* Fine volume slide up */
event->f2t = FX_EXTENDED;
event->f2p = (EX_F_VSLIDE_UP << 4) |
(event->vol - 0x90);
break;
case 0x0a: /* Set vibrato speed */
event->f2t = FX_VIBRATO;
event->f2p = (event->vol - 0xa0) << 4;
break;
case 0x0b: /* Vibrato */
event->f2t = FX_VIBRATO;
event->f2p = event->vol - 0xb0;
break;
case 0x0c: /* Set panning */
event->f2t = FX_SETPAN;
event->f2p = ((event->vol - 0xc0) << 4) + 8;
break;
case 0x0d: /* Pan slide left */
event->f2t = FX_PANSLIDE;
event->f2p = (event->vol - 0xd0) << 4;
break;
case 0x0e: /* Pan slide right */
event->f2t = FX_PANSLIDE;
event->f2p = event->vol - 0xe0;
break;
case 0x0f: /* Tone portamento */
event->f2t = FX_TONEPORTA;
event->f2p = (event->vol - 0xf0) << 4;
break;
}
event->vol = 0;
/*printf("%02x %02x %02x %02x %02x %02x\n",
j, event->note, event->ins, event->vol,
event->fxt, event->fxp);*/
}
}
/* Extra track */
mod->xxt[0] = calloc(sizeof(struct xmp_track) +
sizeof(struct xmp_event) * 64 - 1, 1);
mod->xxt[0]->rows = 64;
/* Read and convert patterns */
D_(D_INFO "Stored patterns: %d", mod->pat);
hio_seek(f, start + pat_ptr, SEEK_SET);
for (i = 0; i < mod->pat; i++) {
if (pattern_alloc(mod, i) < 0)
return -1;
mod->xxp[i]->rows = hio_read16b(f);
for (j = 0; j < mod->chn; j++) {
mod->xxp[i]->index[j] = hio_read16b(f) - 1;
}
}
/* Read samples */
D_(D_INFO "Stored samples: %d", mod->smp);
for (i = 0; i < mod->ins; i++) {
if (mod->xxi[i].nsm == 0)
continue;
hio_seek(f, start + sdata[i], SEEK_SET);
if (load_sample(m, f, 0, &mod->xxs[i], NULL) < 0)
return -1;
}
return 0;
}
static int amf_load(struct module_data *m, HIO_HANDLE *f, const int start)
{
struct xmp_module *mod = &m->mod;
int i, j;
struct xmp_event *event;
uint8 buf[1024];
int *trkmap, newtrk;
int ver;
LOAD_INIT();
hio_read(buf, 1, 3, f);
ver = hio_read8(f);
hio_read(buf, 1, 32, f);
strncpy(mod->name, (char *)buf, 32);
set_type(m, "DSMI %d.%d AMF", ver / 10, ver % 10);
mod->ins = hio_read8(f);
mod->len = hio_read8(f);
mod->trk = hio_read16l(f);
mod->chn = hio_read8(f);
mod->smp = mod->ins;
mod->pat = mod->len;
if (ver == 0x0a)
hio_read(buf, 1, 16, f); /* channel remap table */
if (ver >= 0x0d) {
hio_read(buf, 1, 32, f); /* panning table */
for (i = 0; i < 32; i++) {
mod->xxc->pan = 0x80 + 2 * (int8)buf[i];
}
mod->bpm = hio_read8(f);
mod->spd = hio_read8(f);
} else if (ver >= 0x0b) {
hio_read(buf, 1, 16, f);
}
MODULE_INFO();
/* Orders */
/*
* Andre Timmermans <*****@*****.**> says:
*
* Order table: track numbers in this table are not explained,
* but as you noticed you have to perform -1 to obtain the index
* in the track table. For value 0, found in some files, I think
* it means an empty track.
*/
for (i = 0; i < mod->len; i++)
mod->xxo[i] = i;
D_(D_INFO "Stored patterns: %d", mod->pat);
mod->xxp = calloc(sizeof(struct xmp_pattern *), mod->pat + 1);
if (mod->xxp == NULL)
return -1;
for (i = 0; i < mod->pat; i++) {
if (pattern_alloc(mod, i) < 0)
return -1;
mod->xxp[i]->rows = ver >= 0x0e ? hio_read16l(f) : 64;
for (j = 0; j < mod->chn; j++) {
uint16 t = hio_read16l(f);
mod->xxp[i]->index[j] = t;
}
}
/* Instruments */
if (instrument_init(mod) < 0)
return -1;
/* Probe for 2-byte loop start 1.0 format
* in facing_n.amf and sweetdrm.amf have only the sample
* loop start specified in 2 bytes
*/
if (ver <= 0x0a) {
uint8 b;
uint32 len, start, end;
long pos = hio_tell(f);
for (i = 0; i < mod->ins; i++) {
b = hio_read8(f);
if (b != 0 && b != 1) {
ver = 0x09;
break;
}
hio_seek(f, 32 + 13, SEEK_CUR);
if (hio_read32l(f) > 0x100000) { /* check index */
ver = 0x09;
break;
}
len = hio_read32l(f);
if (len > 0x100000) { /* check len */
ver = 0x09;
break;
}
if (hio_read16l(f) == 0x0000) { /* check c2spd */
ver = 0x09;
break;
}
if (hio_read8(f) > 0x40) { /* check volume */
ver = 0x09;
break;
}
start = hio_read32l(f);
if (start > len) { /* check loop start */
ver = 0x09;
break;
}
end = hio_read32l(f);
if (end > len) { /* check loop end */
ver = 0x09;
break;
}
}
hio_seek(f, pos, SEEK_SET);
}
for (i = 0; i < mod->ins; i++) {
/*uint8 b;*/
int c2spd;
if (subinstrument_alloc(mod, i, 1) < 0)
return -1;
/*b =*/ hio_read8(f);
hio_read(buf, 1, 32, f);
instrument_name(mod, i, buf, 32);
hio_read(buf, 1, 13, f); /* sample name */
hio_read32l(f); /* sample index */
mod->xxi[i].nsm = 1;
mod->xxi[i].sub[0].sid = i;
mod->xxi[i].sub[0].pan = 0x80;
mod->xxs[i].len = hio_read32l(f);
c2spd = hio_read16l(f);
c2spd_to_note(c2spd, &mod->xxi[i].sub[0].xpo, &mod->xxi[i].sub[0].fin);
mod->xxi[i].sub[0].vol = hio_read8(f);
/*
* Andre Timmermans <*****@*****.**> says:
*
* [Miodrag Vallat's] doc tells that in version 1.0 only
* sample loop start is present (2 bytes) but the files I
* have tells both start and end are present (2*4 bytes).
* Maybe it should be read as version < 1.0.
*
* CM: confirmed with Maelcum's "The tribal zone"
*/
if (ver < 0x0a) {
mod->xxs[i].lps = hio_read16l(f);
mod->xxs[i].lpe = mod->xxs[i].len;
} else {
mod->xxs[i].lps = hio_read32l(f);
mod->xxs[i].lpe = hio_read32l(f);
}
if (ver < 0x0a) {
mod->xxs[i].flg = mod->xxs[i].lps > 0 ? XMP_SAMPLE_LOOP : 0;
} else {
mod->xxs[i].flg = mod->xxs[i].lpe > mod->xxs[i].lps ?
XMP_SAMPLE_LOOP : 0;
}
D_(D_INFO "[%2X] %-32.32s %05x %05x %05x %c V%02x %5d",
i, mod->xxi[i].name, mod->xxs[i].len, mod->xxs[i].lps,
mod->xxs[i].lpe, mod->xxs[i].flg & XMP_SAMPLE_LOOP ?
'L' : ' ', mod->xxi[i].sub[0].vol, c2spd);
}
/* Tracks */
trkmap = calloc(sizeof(int), mod->trk);
if (trkmap == NULL)
return -1;
newtrk = 0;
for (i = 0; i < mod->trk; i++) { /* read track table */
uint16 t;
t = hio_read16l(f);
trkmap[i] = t;
if (t > newtrk) newtrk = t;
/*printf("%d -> %d\n", i, t);*/
}
for (i = 0; i < mod->pat; i++) { /* read track table */
for (j = 0; j < mod->chn; j++) {
int k = mod->xxp[i]->index[j] - 1;
/* Use empty track if an invalid track is requested
* (such as in Lasse Makkonen "faster and louder")
*/
if (k < 0 || k >= mod->trk)
k = 0;
mod->xxp[i]->index[j] = trkmap[k];
/*printf("mod->xxp[%d]->info[%d].index = %d (k = %d)\n", i, j, trkmap[k], k);*/
}
}
mod->trk = newtrk; /* + empty track */
free(trkmap);
D_(D_INFO "Stored tracks: %d", mod->trk);
mod->trk++;
mod->xxt = calloc (sizeof (struct xmp_track *), mod->trk);
if (mod->xxt == NULL)
return -1;
/* Alloc track 0 as empty track */
if (track_alloc(mod, 0, 64) < 0)
return -1;
/* Alloc rest of the tracks */
for (i = 1; i < mod->trk; i++) {
uint8 t1, t2, t3;
int size;
if (track_alloc(mod, i, 64) < 0)
return -1;
size = hio_read24l(f);
/*printf("TRACK %d SIZE %d\n", i, size);*/
for (j = 0; j < size; j++) {
t1 = hio_read8(f); /* row */
t2 = hio_read8(f); /* type */
t3 = hio_read8(f); /* parameter */
/*printf("track %d row %d: %02x %02x %02x\n", i, t1, t1, t2, t3);*/
if (t1 == 0xff && t2 == 0xff && t3 == 0xff)
break;
event = &mod->xxt[i]->event[t1];
if (t2 < 0x7f) { /* note */
if (t2 > 0)
event->note = t2 + 1;
event->vol = t3;
} else if (t2 == 0x7f) { /* copy previous */
memcpy(event, &mod->xxt[i]->event[t1 - 1],
sizeof(struct xmp_event));
} else if (t2 == 0x80) { /* instrument */
event->ins = t3 + 1;
} else { /* effects */
uint8 fxp, fxt;
fxp = fxt = 0;
switch (t2) {
case 0x81:
fxt = FX_SPEED;
fxp = t3;
break;
case 0x82:
if ((int8)t3 > 0) {
fxt = FX_VOLSLIDE;
fxp = t3 << 4;
} else {
fxt = FX_VOLSLIDE;
fxp = -(int8)t3 & 0x0f;
}
break;
case 0x83:
event->vol = t3;
break;
case 0x84:
/* AT: Not explained for 0x84, pitch
* slide, value 0x00 corresponds to
* S3M E00 and 0x80 stands for S3M F00
* (I checked with M2AMF)
*/
if ((int8)t3 >= 0) {
fxt = FX_PORTA_DN;
fxp = t3;
} else if (t3 == 0x80) {
fxt = FX_PORTA_UP;
fxp = 0;
} else {
fxt = FX_PORTA_UP;
fxp = -(int8)t3;
}
break;
case 0x85:
/* porta abs -- unknown */
break;
case 0x86:
fxt = FX_TONEPORTA;
fxp = t3;
break;
/* AT: M2AMF maps both tremolo and tremor to
* 0x87. Since tremor is only found in certain
* formats, maybe it would be better to
* consider it is a tremolo.
*/
case 0x87:
fxt = FX_TREMOLO;
fxp = t3;
break;
case 0x88:
fxt = FX_ARPEGGIO;
fxp = t3;
break;
case 0x89:
fxt = FX_VIBRATO;
fxp = t3;
break;
case 0x8a:
if ((int8)t3 > 0) {
fxt = FX_TONE_VSLIDE;
fxp = t3 << 4;
} else {
fxt = FX_TONE_VSLIDE;
fxp = -(int8)t3 & 0x0f;
}
break;
case 0x8b:
if ((int8)t3 > 0) {
fxt = FX_VIBRA_VSLIDE;
fxp = t3 << 4;
} else {
fxt = FX_VIBRA_VSLIDE;
fxp = -(int8)t3 & 0x0f;
}
break;
case 0x8c:
fxt = FX_BREAK;
fxp = t3;
break;
case 0x8d:
fxt = FX_JUMP;
fxp = t3;
break;
case 0x8e:
/* sync -- unknown */
break;
case 0x8f:
fxt = FX_EXTENDED;
fxp = (EX_RETRIG << 4) | (t3 & 0x0f);
break;
case 0x90:
fxt = FX_OFFSET;
fxp = t3;
break;
case 0x91:
if ((int8)t3 > 0) {
fxt = FX_EXTENDED;
fxp = (EX_F_VSLIDE_UP << 4) |
(t3 & 0x0f);
} else {
fxt = FX_EXTENDED;
fxp = (EX_F_VSLIDE_DN << 4) |
(t3 & 0x0f);
}
break;
case 0x92:
if ((int8)t3 > 0) {
fxt = FX_PORTA_DN;
fxp = 0xf0 | (fxp & 0x0f);
} else {
fxt = FX_PORTA_UP;
fxp = 0xf0 | (fxp & 0x0f);
}
break;
case 0x93:
fxt = FX_EXTENDED;
fxp = (EX_DELAY << 4) | (t3 & 0x0f);
break;
case 0x94:
fxt = FX_EXTENDED;
fxp = (EX_CUT << 4) | (t3 & 0x0f);
break;
case 0x95:
fxt = FX_SPEED;
if (t3 < 0x21)
t3 = 0x21;
fxp = t3;
break;
case 0x96:
if ((int8)t3 > 0) {
fxt = FX_PORTA_DN;
fxp = 0xe0 | (fxp & 0x0f);
} else {
fxt = FX_PORTA_UP;
fxp = 0xe0 | (fxp & 0x0f);
}
break;
case 0x97:
fxt = FX_SETPAN;
fxp = 0x80 + 2 * (int8)t3;
break;
}
event->fxt = fxt;
event->fxp = fxp;
}
}
}
/* Samples */
D_(D_INFO "Stored samples: %d", mod->smp);
for (i = 0; i < mod->ins; i++) {
if (load_sample(m, f, SAMPLE_FLAG_UNS, &mod->xxs[i], NULL) < 0)
return -1;
}
m->quirk |= QUIRK_FINEFX;
return 0;
}
static int chip_load(struct module_data *m, HIO_HANDLE *f, const int start)
{
struct xmp_module *mod = &m->mod;
struct mod_header mh;
uint8 *tidx;
int i, j, tnum;
LOAD_INIT();
if ((tidx = calloc(1, 1024)) == NULL) {
goto err;
}
hio_read(&mh.name, 20, 1, f);
hio_read16b(f);
for (i = 0; i < 31; i++) {
hio_read(&mh.ins[i].name, 22, 1, f);
mh.ins[i].size = hio_read16b(f);
mh.ins[i].finetune = hio_read8(f);
mh.ins[i].volume = hio_read8(f);
mh.ins[i].loop_start = hio_read16b(f);
mh.ins[i].loop_size = hio_read16b(f);
}
hio_read(&mh.magic, 4, 1, f);
mh.len = hio_read8(f);
mh.restart = hio_read8(f);
hio_read(tidx, 1024, 1, f);
hio_read16b(f);
mod->chn = 4;
mod->ins = 31;
mod->smp = mod->ins;
mod->len = mh.len;
mod->pat = mh.len;
mod->rst = mh.restart;
tnum = 0;
for (i = 0; i < mod->len; i++) {
mod->xxo[i] = i;
for (j = 0; j < 4; j++) {
int t = tidx[2 * (4 * i + j)];
if (t > tnum)
tnum = t;
}
}
mod->trk = tnum + 1;
strncpy(mod->name, (char *)mh.name, 20);
set_type(m, "Chiptracker");
MODULE_INFO();
if (instrument_init(mod) < 0)
goto err2;
for (i = 0; i < mod->ins; i++) {
struct xmp_instrument *xxi = &mod->xxi[i];
struct xmp_sample *xxs = &mod->xxs[i];
struct xmp_subinstrument *sub;
if (subinstrument_alloc(mod, i, 1) < 0)
goto err2;
sub = &xxi->sub[0];
xxs->len = 2 * mh.ins[i].size;
xxs->lps = mh.ins[i].loop_start;
xxs->lpe = xxs->lps + 2 * mh.ins[i].loop_size;
xxs->flg = mh.ins[i].loop_size > 1 ? XMP_SAMPLE_LOOP : 0;
sub->fin = (int8) (mh.ins[i].finetune << 4);
sub->vol = mh.ins[i].volume;
sub->pan = 0x80;
sub->sid = i;
if (xxs->len > 0)
xxi->nsm = 1;
instrument_name(mod, i, mh.ins[i].name, 22);
}
if (pattern_init(mod) < 0)
goto err2;
for (i = 0; i < mod->len; i++) {
if (pattern_alloc(mod, i) < 0)
goto err2;
mod->xxp[i]->rows = 64;
for (j = 0; j < 4; j++) {
int t = tidx[2 * (4 * i + j)];
mod->xxp[i]->index[j] = t;
}
}
/* Load and convert tracks */
D_(D_INFO "Stored tracks: %d", mod->trk);
for (i = 0; i < mod->trk; i++) {
if (track_alloc(mod, i, 64) < 0)
goto err2;
for (j = 0; j < 64; j++) {
struct xmp_event *event = &mod->xxt[i]->event[j];
uint8 e[4];
hio_read(e, 1, 4, f);
if (e[0] && e[0] != 0xa8)
event->note = 13 + e[0] / 2;
event->ins = e[1];
event->fxt = e[2] & 0x0f;
event->fxp = e[3];
}
}
m->quirk |= QUIRK_MODRNG;
/* Load samples */
D_(D_INFO "Stored samples: %d", mod->smp);
for (i = 0; i < mod->smp; i++) {
if (mod->xxs[i].len == 0)
continue;
if (load_sample(m, f, SAMPLE_FLAG_FULLREP, &mod->xxs[i], NULL) < 0)
goto err2;
}
free(tidx);
return 0;
err2:
free(tidx);
err:
return -1;
}
int
main(int argc, char *argv[])
{
const char *peak_file1 = "peak1.dat";
const char *prof_file1 = "prof1.dat";
const char *peak_file2 = "peak2.dat";
const char *prof_file2 = "prof2.dat";
satdata_mag *data = NULL;
current_data sat1, sat2;
peak_workspace *peak_workspace_p;
struct timeval tv0, tv1;
sat1.n = 0;
sat2.n = 0;
while (1)
{
int c;
int option_index = 0;
static struct option long_options[] =
{
{ "curr_file", required_argument, NULL, 'j' },
{ "curr_file2", required_argument, NULL, 'k' },
{ "swarm_file", required_argument, NULL, 's' },
{ 0, 0, 0, 0 }
};
c = getopt_long(argc, argv, "j:k:s:", long_options, &option_index);
if (c == -1)
break;
switch (c)
{
case 'j':
fprintf(stderr, "main: reading %s...", optarg);
read_lc(optarg, &sat1);
fprintf(stderr, "done (%zu profiles read)\n", sat1.n);
break;
case 'k':
fprintf(stderr, "main: reading %s...", optarg);
read_lc(optarg, &sat2);
fprintf(stderr, "done (%zu profiles read)\n", sat2.n);
break;
case 's':
fprintf(stderr, "main: reading %s...", optarg);
gettimeofday(&tv0, NULL);
data = satdata_swarm_read_idx(optarg, 0);
gettimeofday(&tv1, NULL);
if (!data)
exit(1);
fprintf(stderr, "done (%zu points read, %g seconds)\n",
data->n, time_diff(tv0, tv1));
break;
default:
print_help(argv);
exit(1);
break;
}
}
if (sat1.n == 0)
{
print_help(argv);
exit(1);
}
peak_workspace_p = peak_alloc(NCURR);
if (data)
{
track_workspace *track_p = track_alloc();
fprintf(stderr, "main: separating into tracks...");
track_init(data, NULL, track_p);
fprintf(stderr, "done\n");
/* north hemisphere peak finding */
analyze_hemisphere(-3.0, 23.0, peak_file1, prof_file1, data, track_p, &sat1);
/* south hemisphere peak finding */
analyze_hemisphere(-23.0, 0.0, peak_file2, prof_file2, data, track_p, &sat1);
satdata_mag_free(data);
track_free(track_p);
}
else if (sat2.n > 0)
{
const char *data_file1 = "data.dat.north";
const char *corr_file1 = "corr.dat.north";
const char *data_file2 = "data.dat.south";
const char *corr_file2 = "corr.dat.south";
fprintf(stderr, "main: searching for northern J peaks in satellite 1...");
find_J_peaks(5.0, 23.0, -100.0, peak_file1, prof_file1, NULL, NULL, &sat1);
fprintf(stderr, "done\n");
fprintf(stderr, "main: searching for northern J peaks in satellite 2...");
find_J_peaks(5.0, 23.0, -100.0, peak_file2, prof_file2, NULL, NULL, &sat2);
fprintf(stderr, "done\n");
correlateJ(data_file1, corr_file1, &sat1, &sat2);
fprintf(stderr, "main: data printed to %s\n", data_file1);
fprintf(stderr, "main: correlation data printed to %s\n", corr_file1);
fprintf(stderr, "main: searching for northern J peaks in satellite 2...");
find_J_peaks(-23.0, -5.0, -100.0, peak_file1, prof_file1, NULL, NULL, &sat1);
fprintf(stderr, "done\n");
fprintf(stderr, "main: searching for northern J peaks in satellite 2...");
find_J_peaks(-23.0, -5.0, -100.0, peak_file2, prof_file2, NULL, NULL, &sat2);
fprintf(stderr, "done\n");
correlateJ(data_file2, corr_file2, &sat1, &sat2);
fprintf(stderr, "main: data printed to %s\n", data_file2);
fprintf(stderr, "main: correlation data printed to %s\n", corr_file2);
}
peak_free(peak_workspace_p);
return 0;
} /* main() */
