GCC Aliasing Checks w/Ограничивающие указатели

Ответ 1

Отключить, чтобы сообщить компилятору прекратить проверку псевдонимов:

добавьте строку:

#pragma GCC ivdep

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

https://gcc.gnu.org/onlinedocs/gcc-4.9.2/gcc/Loop-Specific-Pragmas.html

Ответ 2

Может ли это помочь?

void fn1(const float **restrict a, float *restrict b, int n)
{
    const float * restrict a0 = a[0];
    const float * restrict a1 = a[1];

    ASSUME_ALIGNED(a0); ASSUME_ALIGNED(a1); ASSUME_ALIGNED(b);

    for (int i = 0; i < n; ++i)
        b[i] = a0[i] + a1[i];
}

Изменить: вторая попытка:). С информацией из http://locklessinc.com/articles/vectorize/

gcc --fast-math ...

Ответ 3

Ну, а как насчет флага

-fno-strict-aliasing

?

Как я понял, вы правы, вы просто хотите знать, как отключить эти проверки? Если это все, этот параметр в командной строке gcc должен помочь вам.

EDIT:

В дополнение к вашему комментарию: запрещено ли использовать указатели ограничения типа const?

это из ISO/IEC 9899 (6.7.3.1 Формальное определение ограничения):

1.

Пусть D - объявление обычного идентификатора, которое предоставляет средство для обозначения объект P как указатель с ограничениями типа T.

4.

Во время каждого исполнения B пусть L - любое l-значение, которое имеет & L на основе P. Если L используется для доступ к значению объекта X, который он обозначает, и X также модифицируется (любым способом), то применяются следующие требования: T не должен быть const-квалифицированным. Все остальные значения используемый для доступа к значению X, также должен иметь свой адрес на основе P. Каждый доступ, который изменяет X, также следует рассматривать P для целей настоящего подпункта. Если P присваивается значение выражения указателя E, которое основано на другом ограниченном указателе объект P2, связанный с блоком B2, то либо выполнение B2 должно начинаться раньше выполнение B или выполнение B2 заканчивается до назначения. Если эти требования не выполняются, тогда поведение undefined.

И гораздо интереснее, так же как и в регистре:

6.

Переводчик может игнорировать любые или все последствия использования ограничений ограничения.

Итак, если вы не можете найти параметр команды, который заставляет gcc делать это, возможно, это невозможно, потому что по стандарту он не должен давать возможность сделать это.

Ответ 4

Извиняюсь заранее, потому что я не могу воспроизвести результаты с GCC 4.7 на своей машине, но есть два возможных решения.

• Используйте typedef для правильной компоновки * restrict * restrict. Это, по словам бывшего коллеги, который разработчик компилятора LLVM, единственное исключение из typedef ведет себя как препроцессор в C, и он существует, чтобы позволить вам поведение сглаживания.

  Я попытался это сделать ниже, но я не уверен, что мне это удалось. Пожалуйста, внимательно проверьте мою попытку.

• Используйте синтаксис, описанный в ответах на используя ограничитель ограничений с массивами переменной длины (VLA) C99.

  Я попытался это сделать ниже, но я не уверен, что мне это удалось. Пожалуйста, внимательно проверьте мою попытку.

Вот код, который я использовал для выполнения своих экспериментов, но я не смог определить окончательно, если бы какой-либо из моих предложений работал по желанию.

#define ALIGN_BYTES 32
#define ASSUME_ALIGNED(x) x = __builtin_assume_aligned(x, ALIGN_BYTES)

void fn0(const float *restrict a0, const float *restrict a1,
         float *restrict b, int n)
{
    ASSUME_ALIGNED(a0); ASSUME_ALIGNED(a1); ASSUME_ALIGNED(b);

    for (int i = 0; i < n; ++i)
        b[i] = a0[i] + a1[i];
}

#if defined(ARRAY_RESTRICT)
void fn1(const float *restrict a[restrict], float * restrict b, int n)
#elif defined(TYPEDEF_SOLUTION)
typedef float * restrict frp;
void fn1(const frp *restrict a, float *restrict b, int n)
#else
void fn1(const float *restrict *restrict a, float *restrict b, int n)
#endif
{
    //ASSUME_ALIGNED(a[0]); ASSUME_ALIGNED(a[1]); ASSUME_ALIGNED(b);

    for (int i = 0; i < n; ++i)
        b[i] = a[0][i] + a[1][i];
}

Опять же, я прошу прощения за наполовину испеченный характер этого ответа. Пожалуйста, не проголосуйте за меня, но не преуспевайте.