Операция матричного умножения на замену кода ANSI/ISO C BLAS

Можно разработать библиотеку замещения кода для операций умножения с плавающей точкой для матрицы/матрицы и матрицы/вектора с функциями умножения sgemm определено в библиотеке MathWorks C BLAS. Используйте стороннюю библиотеку C BLAS для замены и измените требования к сборке в этом примере, чтобы указать на библиотеку. Для разработки библиотеки замещения кода используйте интерактивный или программный подход. Для получения дополнительной информации см. раздел «Разработка библиотеки замещения кода».

Интерактивная разработка библиотеки замещения кода

  1. Откройте Code Replacement Tool (crtool) из командной строки MATLAB с помощью следующей команды:

    >>crtool

  2. Составьте таблицу.

    1. В контекстном меню crtool выберите File > New Table.

    2. На правой панели назовите таблицу crl_table_cblas. Нажмите Apply.

  3. Создайте запись. В контекстном меню crtool выберите File > New entry > CBlas Operation (Atlas).

  4. Создайте параметры входа. В раскрывающемся списке Function выберите Multiply.

  5. Создайте концептуальное представление. Концептуальное представление описывает сигнатуру функции, которую вы хотите заменить. В Conceptual function подразделе crtool задайте возвращаемый аргумент y1и входной параметр, u1 и u2 с Data Type сингла и Argument Type матрица.

  6. Создайте представление реализации. Представление реализации описывает сигнатуру оптимизационной функции. В данном примере, чтобы указать, что аргументы реализации имеют тот же порядок и свойства, что и концептуальные аргументы, установите флажок Make conceptual and implementation argument types the same. Установите BLAS level равным 3(Vector).

    Задайте Name для функции замены под Function prototype.

  7. Укажите сведения о сборке. Перейдите на вкладку Build Information, чтобы открыть панель Requirements к сборке. Укажите файлы (источник, заголовок, объект), которые требуется генератору кода для замены кода. В данном примере вам не нужно указывать информацию о сборке.

  8. Проверьте и сохраните таблицу. На вкладке Mapping Information нажмите Validate entry. В контекстном меню crtool выберите File > Save table > Save.

  9. Зарегистрируйте библиотеку замещения кода. Регистрация создает библиотеку, состоящую из заданных таблиц. Выберите File > Generate registration file. В диалоговом окне Сгенерировать файл регистрации заполните следующие поля:

    Чтобы использовать библиотеку замены кода, обновите текущий сеанс работы с MATLAB командой:

    >>sl_refresh_customizations

  10. Проверьте библиотеку замещения кода. Из командной строки MATLAB откройте библиотеку с помощью Code Replacement Viewer и проверьте, что таблица и запись заданы правильно. Дополнительные сведения см. в разделе Проверка библиотеки замещения кода. Сконфигурируйте модель, чтобы использовать библиотеку замены кода, сгенерировать код и проверить, что замена происходит должным образом. Если происходит непредвиденное поведение, исследуйте журналы попадания и пропуска, чтобы устранить проблемы.

Программно разрабатывайте библиотеку замещения кода

  1. Откройте программный интерфейс из меню MATLAB, выбрав New > Function.

  2. Составьте таблицу.

    1. Создайте функцию, чтобы вызвать таблицу библиотеки замещения кода. Функция не должна иметь аргументов и возвращать объект таблицы.

    2. Создайте объект таблицы путем вызова RTW.TflTable.

    function hTable = crl_table_cblas
% Create a function to call the code replacement library table 

%% Create a table object
hTable = RTW.TflTable;

  3. Создайте запись. Поскольку этот пример заменяет функцию, создайте запись замены кода в вашей таблице, вызвав функцию entry RTW.TflCBlasEntryGenerator.

    function hTable = crl_table_cblas
% Create a code replacement library table 

%% Create a table object
hTable = RTW.TflTable;

%% Create an entry
%%%%%%%%%%%% Define library path Cblas library file%%%%%%%%%%%%%%%%%%

LibPath = fullfile('<Specify the path of your cblas library>');

if ispc
    libExt = 'lib';
elseif ismac
    libExt = 'dylib';
else
    libExt = 'so';
end


% Create table entry for dgemm32
hEntry = RTW.TflCBlasEntryGenerator;

  4. Создайте параметры входа. Поскольку эти примеры заменяют функцию, создайте параметры входа путем вызова функции setTflCFunctionEntryParameters.

    function hTable = crl_table_cblas
% Create a code replacement library table 

%% Create a table object
hTable = RTW.TflTable;

%% Create an entry
%%%%%%%%%%%% Define library path Cblas library file%%%%%%%%%%%%%%%%%%

LibPath = fullfile('<Specify the path of your cblas library>');

if ispc
    libExt = 'lib';
elseif ismac
    libExt = 'dylib';
else
    libExt = 'so';
end


% Create table entry for dgemm32
hEntry = RTW.TflCBlasEntryGenerator;

%% Create entry parameters
hEntry.setTflCOperationEntryParameters( ...
    'Key', 'RTW_OP_MUL', ...
    'Priority', 100, ...
    'ImplementationName', 'cblas_sgemm', ...
    'ImplementationSourceFile', 'ccblas.c', ...
    'ImplementationHeaderFile', 'ccblas.h', ...
    'AdditionalIncludePaths',   {LibPath}, ...
    'AdditionalLinkObjs',       {['libcblas' libExt]}, ...
    'AdditionalLinkFlags',        {'-lgfortran'},...
    'SideEffects', true);

  5. Создайте концептуальное представление. Концептуальное представление описывает сигнатуру функции, которую вы хотите заменить. Чтобы явным образом задать свойства аргумента, вызовите функцию createAndAddConceptualArg. Чтобы задать матричный аргумент в вызове функции, используйте класс аргумента RTW.TflArgMatrix и укажите базовый тип и размерности, для которых допустим аргумент. Этот тип записи таблицы поддерживает область значений измерений, заданных в формате [Dim1Min Dim2Min ... DimNMin; Dim1Max Dim2Max ... DimNMax]. Для примера [2 2; inf inf] означает двумерную матрицу размера 2x2 или больше.

    function hTable = crl_table_cblas
% Create a code replacement library table 

%% Create a table object
hTable = RTW.TflTable;

%% Create an entry
%%%%%%%%%%%% Define library path Cblas library file%%%%%%%%%%%%%%%%%%

LibPath = fullfile('<Specify the path of your cblas library>');

if ispc
    libExt = 'lib';
elseif ismac
    libExt = 'dylib';
else
    libExt = 'so';
end


% Create table entry for dgemm32
hEntry = RTW.TflCBlasEntryGenerator;

%% Create entry parameters
hEntry.setTflCOperationEntryParameters( ...
    'Key', 'RTW_OP_MUL', ...
    'Priority', 100, ...
    'ImplementationName', 'cblas_sgemm', ...
    'ImplementationSourceFile', 'ccblas.c', ...
    'ImplementationHeaderFile', 'ccblas.h', ...
    'AdditionalIncludePaths',   {LibPath}, ...
    'AdditionalLinkObjs',       {['libcblas' libExt]}, ...
    'AdditionalLinkFlags',      {'-lgfortran'},...
    'SideEffects', true);

%% Create the conceptual representation
arg = RTW.TflArgMatrix('y1', 'RTW_IO_OUTPUT',  'single');
arg.DimRange = [2 2; Inf Inf];
hEnt.addConceptualArg(arg);

arg = RTW.TflArgMatrix('u1', 'RTW_IO_INPUT',  'single');
arg.DimRange = [2 1; Inf Inf];
hEnt.addConceptualArg(arg);

arg = RTW.TflArgMatrix('u2', 'RTW_IO_INPUT',  'single');
arg.DimRange = [1 2; Inf Inf];
hEnt.addConceptualArg(arg);

  6. Создайте представление реализации. Представление реализации описывает сигнатуру оптимизационной функции. Чтобы указать, что аргументы реализации имеют тот же порядок и свойства, что и концептуальные аргументы, вызовите функцию getTflArgFromString. Пример кода конфигурирует специальные аргументы реализации, которые требуются для dgemm и dgemv замена функций. Методы удобства setReturn и addArgument задает, является ли аргумент возврата значением или аргументом, и добавляет аргумент к массиву аргументов реализации элемента. Добавьте полную запись в таблицу путем вызова функции addEntry.

    function hTable = crl_table_cblas
% Create a code replacement library table 

%% Create a table object
hTable = RTW.TflTable;

%% Create an entry
%%%%%%%%%%%% Define library path Cblas library file%%%%%%%%%%%%%%%%%%

LibPath = fullfile('<Specify the path of your cblas library>');

if ispc
    libExt = 'lib';
elseif ismac
    libExt = 'dylib';
else
    libExt = 'so';
end


% Create table entry for dgemm32
hEntry = RTW.TflCBlasEntryGenerator;

%% Create entry parameters
hEntry.setTflCOperationEntryParameters( ...
    'Key', 'RTW_OP_MUL', ...
    'Priority', 100, ...
    'ImplementationName', 'cblas_sgemm', ...
    'ImplementationSourceFile', 'ccblas.c', ...
    'ImplementationHeaderFile', 'ccblas.h', ...
    'AdditionalIncludePaths',   {LibPath}, ...
    'AdditionalLinkObjs',       {['libcblas' libExt]}, ...
    'AdditionalLinkFlags',      {'-lgfortran'},...
    'SideEffects', true);

%% Create the conceptual representation

arg = RTW.TflArgMatrix('y1', 'RTW_IO_OUTPUT',  'single');
arg.DimRange = [2 2; Inf Inf];
hEnt.addConceptualArg(arg);

arg = RTW.TflArgMatrix('u1', 'RTW_IO_INPUT',  'single');
arg.DimRange = [2 1; Inf Inf];
hEnt.addConceptualArg(arg);

arg = RTW.TflArgMatrix('u2', 'RTW_IO_INPUT',  'single');
arg.DimRange = [1 2; Inf Inf];
hEnt.addConceptualArg(arg);

%% Create the Implementation Representation
% Using RTW.TflCBlasEntryGenerator for sgemm requires the following
% implementation signature:
%
% void no_name( const integer ORDER, const integer TRANSA, const integer TRANSB, 
% const integer M, const integer N, const integer K, const double ALPHA, 
% const double* u1, const integer LDA, const double* u2, const integer LDB, 
% const double BETA, double* y1, const integer LDC );
%
% When a match occurs, the code generator computes the
% values for M, N, K, LDA, LDB, and LDC and inserts them into the
% generated code. TRANSA and TRANSB are set to 'N'.

% Specify replacement function signature
arg = hEntry.getTflArgFromString('y2','void');
arg.IOType = 'RTW_IO_OUTPUT';
hEntry.Implementation.setReturn(arg); 

arg = hEntry.getTflArgFromString('ORDER','integer',102.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('TRANSA','integer',111.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('TRANSB','integer',111.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('M','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('N','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('K','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('ALPHA','double',1.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('u1','single*');
arg.Type.BaseType.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('LDA','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('u2','single*');
arg.Type.BaseType.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('LDB','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('BETA','double',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('y1','single*');
arg.IOType = 'RTW_IO_OUTPUT';
hEntry.Implementation.addArgument(arg);

arg = hEntry.getTflArgFromString('LDC','integer',0.000000000000000);
arg.Type.ReadOnly = true;
hEntry.Implementation.addArgument(arg);


%% Add the entry to the table
addEntry(hTable, hEntry);

  7. Укажите сведения о сборке. В параметрах входа укажите файлы (заголовок, источник, объект), которые нужны генератору кода для замены кода. В данном примере информация о сборке не требуется.

  8. Проверьте и сохраните файл индивидуальной настройки. В меню MATLAB сохраните этот файл индивидуальной настройки, выбрав File > Save. В командной строке проверьте таблицу библиотеки замещения кода, вызвав ее:

    >> hTable = crl_table_cblas

  9. Зарегистрируйте библиотеку замещения кода. Регистрация создает библиотеку замещения кода путем определения имени библиотеки, таблиц замещения кода и другой информации. Создайте файл регистрации с этими спецификациями:

    function rtwTargetInfo(cm)
 
cm.registerTargetInfo(@loc_register_crl);
end
 
function this = loc_register_crl 
 
this(1) = RTW.TflRegistry; 
this(1).Name = 'CRL for matrix multiplication for Mathworks CBlas code’;
this(1).TableList = {'crl_table_cblas.m'}; % table created in this example
this(1).TargetHWDeviceType = {'*'};
this(1).Description = 'Example code replacement library';

end

    Чтобы использовать библиотеку замены кода, обновите текущий сеанс работы с MATLAB командой:

    >>sl_refresh_customizations

  10. Проверьте библиотеку замещения кода. Из командной строки MATLAB откройте библиотеку с помощью Code Replacement Viewer и проверьте, что таблица и запись заданы правильно. Дополнительные сведения см. в разделе Проверка библиотеки замещения кода. Сконфигурируйте модель, чтобы использовать библиотеку замены кода, сгенерировать код и проверить, что замена происходит должным образом. Если происходит непредвиденное поведение, исследуйте журналы попадания и пропуска, чтобы устранить проблемы.

Похожие темы

Embedded Coder

Поддержка

Для просмотра документации необходимо авторизоваться на сайте