I am using the ublas matrices with mkl. The documentation of mkl suggests/requires that allocated arrays are aligned on 16 byte boundaries.
The matrix structures within ublas delegate the memory management in the std allocator which at the very end reduces to new.
I know that one way to address this is write my own allocator for the unbounded_array that lies on the basis of the storage structures.
However, this is an error-prone job and it is suggested that this is avoided.
Any other solutions? – hopefully more straightforward ?
Thank you in advance for your help,
ps: using visual studio 2010 on winn7/visa 64 bit and mkl 10.3.