Tedious object allocation and deallocation

Using the old C++ (i.e. C++03) does not allow the use of move constructor and move assignment, which is very efficient in evaluating mathematical expressions with many temporary objects. The problem is how to add these new features into `matrix` and be able to compile it with CUDA 7.0
