Added function which creates laplace operator matrix
Added function to work with matrices in band format
Added routines for improvment of memory usage
Refactoring of copy host-device functions
Implemented functions for factorization which use low amount of device memory