BNN optimization.

See merge request !12

Includes input energy parameter in the model and adds non-linearities

See merge request !11

Handle multi-pulse data

See merge request !10

Speed up prediction and outlier detection.

See merge request !9

Using threading for parallel processing speeds up prediction by a factor of 100: most of the time consumed was on passing data to the processes. Switched from the EllipticEnvelope to an IQR-based outlier detection method -- faster and does not assume Gaussianity. Since this happens after the PCA, the data is already decorrelated.