subject to change — @дневники: асоциальная сеть

среда, 28 декабря 2011

19:20

Доступ к записи ограничен

Пау-чок

Закрытая запись, не предназначенная для публичного просмотра

URL

понедельник, 26 декабря 2011

01:41

Доступ к записи ограничен

Пау-чок

Закрытая запись, не предназначенная для публичного просмотра

URL

четверг, 08 декабря 2011

01:04

Доступ к записи ограничен

Пау-чок

Закрытая запись, не предназначенная для публичного просмотра

URL

четверг, 03 ноября 2011

15:44

Пау-чок

Решил попринимать глицин. По-студенчески - две таблеточки на ночь под язык. Чтоб то, что у меня вместо мозга, хоть немного начало работать в нужных направлениях.

Боги! Какие начали сниться сны! Это что-то!

@темы: личное

URL

пятница, 07 октября 2011

00:54

Пау-чок

В продолжение предыдущего поста. В основном - для себя.
Логарифм по основанию двойки. Погрешность результата ±0.008. Отрабатывает во флеше (естественно, не в таком виде) где-то в 4.5 раз быстрее. В сях тестил - получилось быстрее раз в 6.5.

float log2_K = 0.344285328223736808439525702852;
float log2_K_dbl_and_128 = 128.688570656447473616879051405704;

float inline log2_(float x){
  union{
    float f; //32 bit float
    uint32_t i; //32 bit unsigned int
  }y,m,Ex;
  uint32_t E;
  y.f=x;
  E=y.i&0x7F800000;
  Ex.i=((E^0x7F800000)-0x800000)&0x7F800000;
  Ex.f*=x;
  m.i=(y.i&0x7FFFFF)|0x3F800000;
  y.f=m.f+(float)(E>>23)-log2_K*Ex.f*(Ex.f-3)-log2_K_dbl_and_128;
  return y.f;
}

В принципе, можно получить и бо́льшую точность. Тем более, что запас по быстродействию в позволяет.
Но мне стало лень =)

@музыка: Rock-online: Smokie - Who The X Is Alice (living Next Door To Alice)

@темы: программерское

URL

среда, 05 октября 2011

21:43

Пау-чок

Вот то, о чём я и писал. "Проблема википедии".
Многим программистам известен довольно нетривиальный алгоритм приближённого вычисления "инверсного квадратного корня" - т.е. функции y(x)=1/sqrt(x). Ну, или по крайней мере известно о его существовании, а дальнейшее гуглится на раз.
Рассматривая возможность применения алгоритма применительно к флешу... В общем...
Немного программерской магии.
Функция "быстрого инверсного квадратного корня", как она представлена в коде Quake III, имеет вид:

float Q_rsqrt( float number )
{
        long i;
        float x2, y;
        const float threehalfs = 1.5F;
        x2 = number * 0.5F;
        y  = number;
        i  = * ( long * ) &y;                       // evil floating point bit level hacking
        i  = 0x5f3759df - ( i >> 1 );               // what the fuck?
        y  = * ( float * ) &i;
        y  = y * ( threehalfs - ( x2 * y * y ) );   // 1st iteration
//    y  = y * ( threehalfs - ( x2 * y * y ) );   // 2nd iteration, this can be removed
        return y;
}

На первый взгляд вообще не понятно, как эта функция работает и почему результат вообще получается верным.
Основой алгоритма является строка, к которой кто-то из разработчиков id Software приписал вполне себе закономерный комментарий: "what the fuck?". Да, там происходит магия =) Но магия эта, как оказалось, вполне объяснима с математической точки зрения. Правда, что должно было упасть на голову изобретателю метода, чтобы начать мыслить в таком направлении, я даже представить не могу...

Покопавшись в нескольких гримуарах, описывающих ритуалы призвания магических чисел типа выше упомянутого 0x5f3759df, и кое-что даже усвоив, не смог устоять и попробовал скастовать что-то подобное.
Получилось не ахти, конечно. Т.к. с какого-то перепуга я решил приняться за функцию y(x)=1/x. Ну, подумалось мне, во флеше деление работает не ахти как быстро. А вдруг удастся ускорить?
В общем-то получилось вот что, дающее относительную погрешность не больше ±6.5e-4%.

float inline inv(float x){
  union{
    float f; //32 bit float
    uint32_t i; //32 bit unsigned int
  } y;
  y.f=x;
  y.i=0x7EF311C3-y.i;
  y.f*=2-y.f*x;
  y.f*=2-y.f*x;
  return y.f;
}

Практической пользы от этой функции чуть меньше, чем никакой =) Она отрабатывает раза в два-три медленнее, чем простое деление =)
Однако же сам факт того, что можно получить число 1/x вообще не используя операцию деления, а лишь чуть-чуть пошаманив, меня весьма вдохновляет =)

@темы: программерское

URL

воскресенье, 02 октября 2011