Reference data and output of simulator function are square matrices
We want to use trace norm of the difference between the matrices as a distance function
Currently, since the keep_all summary statistic reshapes everything into 1D array, we'd have to reshape everything back into a matrix in the distance function, before taking the trace norm. The limitation of all summary statistics having to produce a 1D array has to be removed.
Consider the following use case:
Currently, since the
keep_all
summary statistic reshapes everything into 1D array, we'd have to reshape everything back into a matrix in the distance function, before taking the trace norm. The limitation of all summary statistics having to produce a 1D array has to be removed.