hagenw
commented
8 years ago OK, I was now able to use only one Dirac impulse. I changed findnearestneighbour() in a way that it will return now only as many neighbors as available.
Would be nice if you could test it, too.
In order to improve memory consumption and speeding up calculation I reduced the number of stored dirac pulses returned by
dummy_irs()to only 4. It principle it should also work by using only 1, but this would need some changes in the other functions. We could still think about if this would be worth the effort. I tested the new version at the example from the README and the difference between the results were <10^-18.