Neuroinflab / kCSD-python

Kernel Current Source Density
https://kcsd-python.readthedocs.io/en/latest/
Other
15 stars 23 forks source link

Code vectorization #93

Closed abukaj closed 5 years ago

abukaj commented 5 years ago

https://github.com/Neuroinflab/kCSD-python/blob/2ce755a0b43324e16053d8a2319b053f614abcd7/kcsd/KCSD.py#L219-L220

is IMHO mathematically equivalent with:

estimation[:, t] = np.dot(estimation_table, beta)

which should be faster.

However, the obtained results may differ due to numerical issues.