segmentation fault while running with m=1 with your read from dataset

icthrm / cwSDTWnano

Here we proposed two novel algorithms, the Direct Subsequence Dynamic Time Warping for nanopore raw signal search (DSDTWnano) and the continuous wavelet Subsequence DTW for nanopore raw signal search (cwSDTWnano), to enable the direct subsequence inquiry and exact mapping in the nanopore raw signal datasets. The proposed algorithms are based on the idea of Subsequence-extended Dynamic Time Warping (SDTW) and directly operates on the raw signals, without any loss of information. DSDTWnano could ensure an output of highly accurate query result and cwSDTWnano is the accelerated version of DSDTWnano, with the help of seeding and multi-scale coarsening of signals that based on continuous wavelet transform (CWT).

10 stars 2 forks source link

I'm also observing the same. ./sat-query -i ../../ref.fa -p ../../cwSDTWnano_Dataset_lite/Lambda_data/Lambda_3000_signal/fa11fec2-6462-4493-98db-efe92af43a6f.signal -o out -m 0

RAw signal: cwSDTWnano_Dataset_lite/Lambda_data/Lambda_3000_signal/fa10a452-c965-406e-99b4-9f41e30d8b39.signal (from your dataset) reference genome: https://www.ncbi.nlm.nih.gov/nuccore/J02459.1?report=fasta

log:

./sat-query -i ../../ref.fa -p ../../cwSDTWnano_Dataset_lite/Lambda_data/Lambda_3000_signal/fa11fec2-6462-4493-98db-efe92af43a6f.signal -o out -m 0

Transform genomes to signal sequence...[Beginning] 2020/10/14 Wed 12:2:3 48502 genomes are readed. Transform genomes to signal sequence...[Finished] 2020/10/14 Wed 12:2:3 (time elapse: 140ms) Short reads inquiry...[Beginning] 2020/10/14 Wed 12:2:3 18867 signals are readed. Direct inquiry...[Beginning] 2020/10/14 Wed 12:2:3 Segmentation fault (core dumped) output: -1 -1 | 1.91714e-318 5.72706e-318 | 2.61057e-318 1.03816e-317 diff: 7.77098e-318

icthrm / cwSDTWnano

segmentation fault while running with m=1 with your read from dataset #5