Я не могу выяснить, какую функцию использовать в MKL для вычитания двух одномерных матриц (например, Y = A — B).
Я просмотрел документацию, касающуюся умножения матриц (http://software.intel.com/sites/products/documentation/doclib/mkl_sa/11/tutorials/mkl_mmx_c/GUID-36BFBCE9-EB0A-43B0-ADAF-2B65275726EA.htm) но, кажется, нет никаких примеров сложения / вычитания.
Для плотных матриц вы всегда выделяете матрицы как вектор длины n*m
поэтому вы используете соответствующие подпрограммы уровня 1 BLAS, например: SAXPY(N,SA,SX,INCX,SY,INCY) <==> y = alpha*x+y
Других решений пока нет …