When defining coarser->r and coarser->nr there is a problem with the size of the vectors. I think one can also reduce some redundancy