Inverts an array of matrices so that: inva[,,i] = inv(a[,,i])
mult__inv(a, whichoutput, nargout)
"nargout == 1" returns inva -> matrix array: a^(-1) (N x N x T)
"nargout == 2" returns inva -> matrix array: a^(-1) (N x N x T) and logabsdet -> vector of log(abs(det(a))) (1 x T)
matrix array (N x N x T)
if = 'lobabsdet' computes the log(abs(det(a))) only (text).
"nargout == 1" or "nargout == 2"(scalar)
This function is modified version of the "mult__inv" function by Le and Singleton (2018).
"A Small Package of Matlab Routines for the Estimation of Some Term Structure Models."
(Euro Area Business Cycle Network Training School - Term Structure Modelling).
Available at: https://cepr.org/40029