C++ (Qt)
    while (src != srcEnd)
    {
        // присваиваем временные имена
        const float& x0 = *(src++);
        const float& y0 = *(src++);
        const float& z0 = *(src++);
 
        const float& x1 = *(src++);
        const float& y1 = *(src++);
        const float& z1 = *(src++);
 
        const float& x2 = *(src++);
        const float& y2 = *(src++);
        const float& z2 = *(src++);
 
        // вычисляем векторы
        float a[3] = {x1 - x0, y1 - y0, z1 - z0},
              b[3] = {x2 - x0, y2 - y0, z2 - z0};
 
        const int x = 0,
                  y = 1,
                  z = 2;
        // вычисление и запись перпендикуляра
        *dest = a[y] * b[z] - a[z] * b[y];
        dest[1] = a[z] * b[x] - a[x] * b[z];
        dest[2] = a[x] * b[y] - a[y] * b[x];
 
        // вычисление длины и нормализация
        float inv_length = 1.0 / sqrt(*dest * *dest +
                                      dest[y] * dest[y] +
                                      dest[z] * dest[z]);
        *dest *= inv_length;
        dest[y] *= inv_length;
        dest[z] *= inv_length;
 
        // копирование во все вершины
        memcpy(dest + 3, dest, 3 * sizeof(float));
        memcpy(dest + 6, dest, 3 * sizeof(float));
 
        dest += 9;
    }
 

C++ (Qt)
    int times = (srcEnd - src) / sizeof(float) / 9;
    #pragma omp parallel for
    for(int i = 0; i < times; i++)
    {
        // присваиваем временные имена
        int j = i * 9;
        const float& x0 = src[j];
        const float& y0 = src[j + 1];
        const float& z0 = src[j + 2];
 
        const float& x1 = src[j + 3];
        const float& y1 = src[j + 4];
        const float& z1 = src[j + 5];
 
        const float& x2 = src[j + 6];
        const float& y2 = src[j + 7];
        const float& z2 = src[j + 8];
 
        // вычисляем векторы
        float a[3] = {x1 - x0, y1 - y0, z1 - z0},
              b[3] = {x2 - x0, y2 - y0, z2 - z0};
 
        const int x = 0,
                  y = 1,
                  z = 2;
        // вычисление и запись перпендикуляра
        float& xNorm = dest[j] = a[y] * b[z] - a[z] * b[y];
        float& yNorm = dest[j + 1] = a[z] * b[x] - a[x] * b[z];
        float& zNorm = dest[j + 2] = a[x] * b[y] - a[y] * b[x];
 
        // вычисление длины и нормализация
        float inv_length = 1.0 / sqrt(xNorm * xNorm +
                                      yNorm * yNorm +
                                      zNorm * zNorm);
        xNorm *= inv_length;
        yNorm *= inv_length;
        zNorm *= inv_length;
 
        // копирование во все вершины
        memcpy(&xNorm + 3, &xNorm, 3 * sizeof(float));
        memcpy(&xNorm + 6, &xNorm, 3 * sizeof(float));
 
    }
 

C++ (Qt)
    while (src != srcEnd)
    {
        // присваиваем временные имена
        const float& x0 = *(src++);
        const float& y0 = *(src++);
        const float& z0 = *(src++);
 
        const float& x1 = *(src++);
        const float& y1 = *(src++);
        const float& z1 = *(src++);
 
        const float& x2 = *(src++);
        const float& y2 = *(src++);
        const float& z2 = *(src++);
 
        // вычисляем векторы
        float a[3] = {x1 - x0, y1 - y0, z1 - z0},
              b[3] = {x2 - x0, y2 - y0, z2 - z0};
 
        const int x = 0,
                  y = 1,
                  z = 2;
        // вычисление и запись перпендикуляра
        *dest = a[y] * b[z] - a[z] * b[y];
        dest[1] = a[z] * b[x] - a[x] * b[z];
        dest[2] = a[x] * b[y] - a[y] * b[x];
 
        // вычисление длины и нормализация
        float inv_length = 1.0 / sqrt(*dest * *dest +
                                      dest[y] * dest[y] +
                                      dest[z] * dest[z]);
        *dest *= inv_length;
        dest[y] *= inv_length;
        dest[z] *= inv_length;
 
        // копирование во все вершины
        memcpy(dest + 3, dest, 3 * sizeof(float));
        memcpy(dest + 6, dest, 3 * sizeof(float));
 
        dest += 9;
    }