본문 바로가기

Enginius/Matlab

Fast implementation of kdpp

Details about kdpp can be found in (http://enginius.tistory.com/494)


Implementation of kdpp using matrix inversion lemma to reduce the computational complexity. 

The speedup increases as the number of points increases (>100)

Starting kdpp where k is 50. 

Starting kdpp_fast where k is 50. 

k: 50 / toc_kdpp: 4.547e-02 sec / toc_kdpp_fast: 9.981e-03 sec (4.6 times faster) 


One important rule of thumb is to lower the length parameter! So that the kernel function is peaky! kdpp_fast utilizes matrix inversion lemma and computing the inverse of a kernel matrix might be numerically unstable if the kernel function is smooth. 


Result


code: main

code: kernel

code: kdpp

code: kdpp_fast


'Enginius > Matlab' 카테고리의 다른 글

MATLAB 책  (0) 2015.06.23
Figure size setting based on screen size  (0) 2015.04.09
Kinect 사용하기 2  (1) 2014.12.26
Spline Fitting  (0) 2014.12.22
다각형과 선분 사이의 접점  (0) 2014.12.17