Comments (7)
Finding alignment is not so hard! I know for each block its vertical deltas for every cell and horizontal delta of last cell in block. This means that if I store vertical and horizontal deltas for all blocks, I can reconstruct path of alignment. For each cell on path, I will have to compute part of block it belongs to.
I implemented this, it slows whole algorithm about 3, 4 times and takes a lot of memory.
from edlib.
Interesting: if query is larger then target, then optimal path will not go through so many blocks and finding alignment will be faster! In HW and SHW this does not affect us because we can not swap query and target, but in NW we can so this should be considered. On the other side, in NW query and target will usually be of similar length, so this does not really help
from edlib.
I should speed up alignment, and also make it consume much less memory.
Options I have for that: Hirschbergs algorithm, or something similar to what they did in SSW.
SSW idea could work pretty well for HW and SHW, since query is small there.
For NW I do not believe it will be very useful, so this is where I could consider Hirschberg.
from edlib.
I could also just try to use less space, by doing something smart hm. But big problem is that I do not want to allocate memory too often. I could use something like vector, to allocate memory dinamically?
from edlib.
I implemented SSW idea of finding alignment. SHW and HW now use NW alignment to find alignment. NW alignment is still implemented like before. Next step is to make NW use less memory, or implement it with Hischberg. I did some tests and it seems like Hirschberg method can be much slower in some cases, so I will not try it so soon.
from edlib.
I should change finding of alignment for NW so it uses only space needed for band, not the whole matrix.
from edlib.
I will create new issue for this!
from edlib.
Related Issues (20)
- Is there a suffix mode? HOT 3
- align two arrays of strings HOT 7
- anchored alignment HOT 2
- Error when query or target is the empty string and task='path'
- memory leak identified in myersCalcEditDistanceSemiGlobal
- Mention Perl binding in Readme HOT 2
- finding all alignments within the string HOT 5
- Cigar results interpretation
- Missing table data in paper
- Possible to customize costs used for gap opening, gap extension, indel and substitutions? HOT 3
- Wheels for python 3.10 HOT 2
- [Documentation] (1) mention ruby-bindings and (2) mention "why use edlib" in the main README HOT 1
- results different with usearch and vsearch global alignment
- Wheels for macos python 3.11 HOT 1
- Choosing longest or shortest alignment HOT 1
- How to report all possibale positions not only the best one ? HOT 1
- Cannot install edlib for Python 3.11 on Windows HOT 1
- Perl Wrapper around Edlib
- Multiple locations but one CIGAR?
- Documentation outdated (v1.1.2 vs v1.2.7)
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from edlib.