Performance improvements in restrict and prolong (devectorization)

This is for non-BLAS types
Performance improvement in interp() (type-inferrence issue)
Compute linear indices without temporaries