BUG: Problem with Actual Implementation of Scipy (1.10.1) 

There is a known bug in Scipy that creates the following problem when using matrix arrays which affects this package, notably when using the "minimize" method: 

ValueError: `f0` passed has more than 1 dimension.

I am not certain what could be the solutions to this problem that avoids going back to a prior version of Scipy (1.4.1)

