Playing SIMD

[Iakh via Digitalmars-d]

Here is my implementatation of SIMD find. Function returns index 
of ubyte in static 16 byte array with unique values.

------
immutable size_t ArraySize = 16;
int simdIndexOf(ubyte niddle, ref const ubyte[ArraySize] haystack)
{
     ubyte16 arr;
     arr.array = haystack[];
     ubyte16 niddles;
     niddles.array[] = niddle;
     ubyte16 result;
     result = __simd_sto(XMM.PCMPEQB, arr, niddles);
     alias Mask = ulong;
	static assert(2 * Mask.sizeof == result.sizeof);
	immutable BitsInByte = 8;

     if (Mask mask = *cast(Mask*)(result.array.ptr))
     {
         return bsf(mask) / BitsInByte;
     }
     else if (Mask mask = *cast(Mask*)(result.array.ptr + 
Mask.sizeof))
     {
         return bsf(mask) / BitsInByte + cast(int)Mask.sizeof;
     }
     else
     {
         return -1;
     }

}
------

Is it optimal and how do you implement this stuf?

Full exaple with comparation of algorithms (SIMD, naive, binary 
search):
http://dpaste.dzfl.pl/f3b8989841e3

Benchmark result on dpaste.dzfl.pl:
SIMD:   TickDuration(157000)
Binary: TickDuration(472000)
Naive:  TickDuration(437000)

At home with defult dub config "dub run --build=release":
SIMD:    TickDuration(241566)
Binary:  TickDuration(450515)
Naive:   TickDuration(450371)