Optimal Density Estimation via Shimazaki-Shinomoto Method // Implements the Shimazaki-Shinomoto method for optimizing the bin width of histograms and the bandwidth of kernel density estimators. The framework minimizes the expected Mean Integrated Squared Error (MISE) and supports both 1D and 2D distributions, fixed and locally adaptive estimators, bootstrap confidence intervals, and 'OpenMP'-accelerated 'C++' 'backends'. Ideally suited for time-dependent rate estimation and identifying intrinsic data structures. For more details see Shimazaki and Shinomoto (2007) doi:10.1162/neco.2007.19.6.1503 and Shimazaki and Shinomoto (2010) doi:10.1007/s10827-009-0180-4.