avoiding using push_back in constructing FE interpolation matrix
added support for blocked storage when available