Matrix Aware divison

MatrixDivison.py

import numpy as np

# matrix-aware division
def mdiv(x, y):

	# everything is float
	x = 1.0*x
	y = 1.0*y

	# if y is a matrix, do matrix division
	if isinstance(y, np.ndarray):

		# turn x into an array of the right size if it isnt already
		if not isinstance(x, np.ndarray):
			# x time sidentity
			x = x * np.identity(y.shape[0])
		
		# invert y
		y = np.linalg.inv(y)

		#dot with x
		return y.dot(x)

	# else do normal division
	else:
		return x / y