Учитывая плотную матрицу A формы nxn, хранящуюся в мажорной колонке, константу c и единичную матрицу I (также nxn). Я хотел бы рассчитать A + c* I в рамках BLAS. Это можно сделать, используя SAXPY с соответствующим приращением.
Тем не менее, я должен сделать это несколько раз для пакета матриц A. Каждый отдельный A в пакете имеет nxn, а c одинаков для всех A. Пакет матриц A хранится в основном формате столбцов в «сцепленном» массиве, как это типично, например, для пакетного GEMM в BLAS. A для каждого экземпляра пакета находится с фиксированным смещением адреса (с шагом). Каков наиболее эффективный способ вычисления A + c* I в BLAS (последовательным) пакетным способом?