Performance of Matrix Division

[judober]

As discussed here, I found a way to increase the performance of A/B and B\A if the dimension of B is <= 3 and A is not a Vector. I will try to make a PR for this issue.

Additionally, I noticed that A/B does not work when B is factroized by lu. I will make annother PR to solve this.