Быстрый синус / косинус для ARMv7 + NEON: ищем тестеров… - PullRequest
Новая
       8

Быстрый синус / косинус для ARMv7 + NEON: ищем тестеров…

15 голосов
/ 06 декабря 2009

Может ли кто-нибудь иметь доступ к iPhone 3GS или Pandora, пожалуйста, протестируйте следующую процедуру сборки, которую я только что написал?

Предполагается, что синусы и косинусы действительно очень быстро вычисляются на векторном FPU NEON. Я знаю, что он прекрасно компилируется, но без адекватного оборудования я не могу его протестировать. Если бы вы могли просто вычислить несколько синусов и косинусов и сравнить результаты с результатами sinf () и cosf (), это действительно помогло бы.

Спасибо! 

#include <math.h>

/// Computes the sine and cosine of two angles
/// in: angles = Two angles, expressed in radians, in the [-PI,PI] range.
/// out: results = vector containing [sin(angles[0]),cos(angles[0]),sin(angles[1]),cos(angles[1])]
static inline void vsincos(const float angles[2], float results[4]) {
    static const float constants[]  = { 
    /* q1 */  0,                M_PI_2,           0,                M_PI_2,
    /* q2 */  M_PI,             M_PI,             M_PI,             M_PI,
    /* q3 */  4.f/M_PI,         4.f/M_PI,         4.f/M_PI,         4.f/M_PI,
    /* q4 */ -4.f/(M_PI*M_PI), -4.f/(M_PI*M_PI), -4.f/(M_PI*M_PI), -4.f/(M_PI*M_PI),
    /* q5 */  2.f,              2.f,              2.f,              2.f,
    /* q6 */  .225f,            .225f,            .225f,            .225f
    };  
    asm volatile(
        // Load q0 with [angle1,angle1,angle2,angle2]
        "vldmia %1, { d3 }\n\t"
        "vdup.f32 d0, d3[0]\n\t"
        "vdup.f32 d1, d3[1]\n\t"
        // Load q1-q6 with constants
        "vldmia %2, { q1-q6 }\n\t"
        // Cos(x) = Sin(x+PI/2), so
        // q0 = [angle1, angle1+PI/2, angle2, angle2+PI/2]
        "vadd.f32 q0,q0,q1\n\t"
        // if angle1+PI/2>PI, substract 2*PI
        // q0-=(q0>PI)?2*PI:0
        "vcge.f32 q1,q0,q2\n\t"
        "vand.f32 q1,q1,q2\n\t"
        "vmls.f32 q0,q1,q5\n\t"
        // q0=(4/PI)*q0 - q0*abs(q0)*4/(PI*PI)
        "vabs.f32 q1,q0\n\t"
        "vmul.f32 q1,q0,q1\n\t"
        "vmul.f32 q0,q0,q3\n\t"
        "vmul.f32 q1,q1,q4\n\t"
        "vadd.f32 q0,q0,q1\n\t"
        // q0+=.225*(q0*abs(q0) - q0)
        "vabs.f32 q1,q0\n\t"
        "vmul.f32 q1,q0,q1\n\t"
        "vsub.f32 q1,q0\n\t"
        "vmla.f32 q0,q1,q6\n\t"
        "vstmia %0, { q0 }\n\t"
        :: "r"(results), "r"(angles), "r"(constants)
        : "memory","cc","q0","q1","q2","q3","q4","q5","q6"
    );  
}

Ответы [ 2 ]

10 голосов
/ 06 декабря 2009

Только что протестировал его на моем beagleboard .. Как сказано в комментариях: Тот же процессор.

Ваш код примерно в 15 раз быстрее, чем Clib ... Отлично!

Я измерил 82 цикла для каждого вызова вашей реализации и 1260 для четырех вызовов c-lib. Обратите внимание, что я скомпилировал с мягким плавающим ABI, и мой OMAP3 является ранним кремнием, поэтому каждый вызов версии c-lib имеет остановку NEON по крайней мере 40 циклов.

Я скопировал результаты ..

http://torus.untergrund.net/code/sincos.zip

Скорее всего, счетчик производительности не будет работать на iphone.

Надеюсь, это то, что вы искали.

3 голосов
/ 06 декабря 2009

Ох, прежде чем я забуду это: Может быть, вы можете обезопасить себя немного работы ..

Взгляните на эти NEON оптимизированные математические функции:

http://code.google.com/p/math-neon/

...