I realised today that sympy's matrix symbols in a vector (for multiple partial element-wise derivatives over a covariance matrix) was incompatible with scipy's cdist during the the optimize.minimize process, as it assumes actual use of numbers when making the function call (fair enough).I realised today that sympy's matrix symbols in