Вызовите повторно используемый код — многопоточный с персистентными данными (Windows Only)

Этот пример требует библиотек, которые характерны для операционной системы Microsoft® Windows® и, поэтому, запуски только на платформах Windows. Это - многопоточный пример, который использует персистентные данные. Два потока вызывают функцию MATLAB® matrix_exp с различными наборами входных данных.

Код MATLAB для этого примера

function [Y,numTimes] = matrix_exp(X) %#codegen
    %
    % The function matrix_exp computes matrix exponential 
    % of the input matrix using Taylor series and returns
    % the computed output. It also returns the number of 
    % times this function has been called.
    %
    persistent count;
    if isempty(count)
        count = 0;
    end
    count = count+1;
    
    E = zeros(size(X));
    F = eye(size(X));
    k = 1;
    while norm(E+F-E,1) > 0
        E = E + F;
        F = X*F/k;
        k = k+1;
    end
    Y = E ;
    
    numTimes = count;

Обеспечьте основную функцию

Чтобы вызвать повторно используемый код, который использует персистентные данные, обеспечьте main функция, что:

  • Включает заголовочный файл matrix_exp.h.

  • Для каждого потока, выделяет память для данных о стеке.

  • Выделяет память для персистентных данных, однажды на корневой процесс, если потоки осуществляют обмен данными, и однажды на поток в противном случае.

  • Вызывает matrix_exp_initialize функция обслуживания. Для получения дополнительной информации смотрите, Развертывают Сгенерированный код.

  • Вызовы matrix_exp.

  • Вызовы matrix_exp_terminate.

  • Освобождает память, используемую в стеке и персистентных данных.

В данном примере main.c содержит:

#include <stdio.h>
#include <stdlib.h>
#include <windows.h>
#include "matrix_exp.h"
#include "matrix_exp_initialize.h"
#include "matrix_exp_terminate.h"
#include "rtwtypes.h"
#define NUMELEMENTS (160*160)

typedef struct {
    real_T in[NUMELEMENTS];
    real_T out[NUMELEMENTS];
    real_T numTimes;
    matrix_expStackData* spillData;
} IODATA;

/*The thread_function calls the matrix_exp function written in MATLAB*/
DWORD WINAPI thread_function(PVOID dummyPtr) {
    IODATA *myIOData = (IODATA*)dummyPtr;
    matrix_exp_initialize(myIOData->spillData);
    matrix_exp(myIOData->spillData, myIOData->in, myIOData->out, &myIOData->numTimes);
    printf("Number of times function matrix_exp is called is %g\n",myIOData->numTimes);
    matrix_exp_terminate();
    return 0;
}

void main() {
    HANDLE thread1, thread2;
    IODATA data1;
    IODATA data2;
    int32_T i;

    /*Initializing data for passing to the 2 threads*/
    matrix_expPersistentData* pd1 = (matrix_expPersistentData*)calloc(1,sizeof(matrix_expPersistentData));
    matrix_expPersistentData* pd2 = (matrix_expPersistentData*)calloc(1,sizeof(matrix_expPersistentData));
    matrix_expStackData* sd1 = (matrix_expStackData*)calloc(1,sizeof(matrix_expStackData));
    matrix_expStackData* sd2 = (matrix_expStackData*)calloc(1,sizeof(matrix_expStackData));
     
    sd1->pd = pd1;
    sd2->pd = pd2;
    data1.spillData = sd1;
    data2.spillData = sd2;
        
    for (i=0;i<NUMELEMENTS;i++) {
        data1.in[i] = 1;
        data1.out[i] = 0;
        data2.in[i] = 1.1;
        data2.out[i] = 0;
    }
    
    data1.numTimes = 0;
    data2.numTimes = 0;
    
    /*Initializing the 2 threads and passing required data to the thread functions*/
    printf("Starting thread 1...\n");
    thread1 = CreateThread(NULL, 0, thread_function, (PVOID) &data1, 0, NULL);
    if (thread1 == NULL){
		perror( "Thread 1 creation failed.");
		exit(EXIT_FAILURE);
	}

    printf("Starting thread 2...\n");
    thread2 = CreateThread(NULL, 0, thread_function, (PVOID) &data2, 0, NULL);
    if (thread2 == NULL){
        perror( "Thread 2 creation failed.");
        exit(EXIT_FAILURE);
    }
    
    /*Wait for both the threads to finish execution*/
    if (WaitForSingleObject(thread1, INFINITE) != WAIT_OBJECT_0){
		perror( "Thread 1 join failed.");
		exit(EXIT_FAILURE);
	}
    
    if (WaitForSingleObject(thread2, INFINITE) != WAIT_OBJECT_0){
		perror( "Thread 2 join failed.");
		exit(EXIT_FAILURE);
	}  
    
    free(sd1);
    free(sd2);
    free(pd1);
    free(pd2);
    
    printf("Finished Execution!\n");
    exit(EXIT_SUCCESS);
    
}

Сгенерируйте повторно используемый код С

Запустите следующий скрипт в командной строке MATLAB, чтобы сгенерировать код.

% This example can only be run on Windows platforms
if ~ispc
 error...
   ('This example requires Windows-specific libraries and can only be run on Windows.');
end

% Setting the options for the Config object
% Create a code gen configuration object
cfg = coder.config('exe');

% Enable reentrant code generation
cfg.MultiInstanceCode = true;

% Compiling 
codegen -config cfg main.c -report matrix_exp.m -args ones(160,160)

Этот скрипт:

  • Генерирует сообщение об ошибке, если пример не работает на платформе Windows.

  • Создает объект настройки генерации кода для генерации исполняемого файла.

  • Включает MultiInstanceCode опция, чтобы сгенерировать допускающий повторное использование, повторно используемый код.

  • Вызывает codegen со следующими опциями:

    • -config передать в объекте cfg настройки генерации кода.

    • main.c включать этот файл в компиляцию.

    • -report создать отчет генерации кода.

    • -args задавать вход в качестве примера с классом, размером и сложностью.

Исследуйте сгенерированный код

codegen генерирует заголовочный файл matrix_exp_types.h, это задает:

  • matrix_expStackData глобальная структура, которая содержит локальные переменные, которые являются слишком большими, чтобы соответствовать на стеке и указателе на matrix_expPersistentData глобальная структура.

  • matrix_expPersistentData глобальная структура, которая содержит персистентные данные.

/*
 * matrix_exp_types.h
 *
 * Code generation for function 'matrix_exp'
 *
 */

#ifndef __MATRIX_EXP_TYPES_H__
#define __MATRIX_EXP_TYPES_H__

/* Include files */
#include "rtwtypes.h"

/* Type Definitions */
#ifndef typedef_matrix_expPersistentData
#define typedef_matrix_expPersistentData

typedef struct {
  double count;
} matrix_expPersistentData;

#endif                                 /*typedef_matrix_expPersistentData*/

#ifndef typedef_matrix_expStackData
#define typedef_matrix_expStackData

typedef struct {
  struct {
    double F[25600];
    double Y[25600];
    double X[25600];
  } f0;

  matrix_expPersistentData *pd;
} matrix_expStackData;

#endif                                 /*typedef_matrix_expStackData*/
#endif

/* End of code generation (matrix_exp_types.h) */

Запустите код

Вызовите код с помощью команды:

system('matrix_exp.exe')
Исполняемое завершение запусков и отчетов.

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