Efficient matrix inversion for Nd arrays, similar to mtimesx operations

I have a square array of function values that I need to invert. Currently I have to invert every slice of the function in a loop. Is there a way to do this operation more efficiently in a manner similar to that used in mtimesx ? Thanks.