import std;
import ldc.simd;
import ldc.gccbuiltins_x86;

//alias vec4 = float4;
alias vec4 = __vector(float[4]);

enum sse = false;

struct mat4
{
   vec4[4] data;
   alias data this;

   this(float[16] vals...)
   {
       static foreach (i; 0 .. 4)
           static foreach (j; 0 .. 4)
               this[i][j] = vals[i*4+j];
   }

   mat4 T() const @property
   {
       mat4 r;
       static foreach (i; 0 .. 4)
           static foreach (j; 0 .. 4)
               r[i][j] = this[j][i];
       return r;
   }

   mat4 opBinary(string op)(auto ref const mat4 o) const
       if (op == "*")
   {
       auto ot = o.T;
       mat4 r;
       static foreach (i; 0 .. 4)
           static foreach (j; 0 .. 4)
           {
               static if (sse)
                   r[i][j] = __builtin_ia32_dpps(this[i], ot[j], cast(ubyte)0xff)[0];
               else
               {
                   r[i][j] = 0;
                   static foreach (k; 0 .. 4)
                       r[i][j] += this[i][k]…