I’m currently using a Python library called “scoringrules” to calculate energy and variogram scores for multivariate forecasts. It serves similar purposes as the scoringRules library in R.

I ran into a array shape problem even though my arrays both have shapes (1, 24). So I went back to the github page for a test code and I’m facing a similar issue with the test code provided by the original author. I haven’t changed anything in the test code. I’ve provided the test code’s image but I’m not sure I can provide mine given how long it is.

Can anyone help please ?

Images provided for better understanding.[[enter image description here](https://i.sstatic.net/nyn6cNPN.png)](https://i.sstatic.net/82DHWDlT.png)

Input two arrays of shape (1, 24) in command scoringrule.energy_score(obs,fcst) and expected energy score to be returned as a (1, 24) array.

For test code, don’t really know what to expect, since I don’t have a full mastery of this library.