dpsoftrast.c

   1
   2 #include <stdio.h>
   3 #include <string.h>
   4 #include <math.h>
   5 #include "quakedef.h"
   6 #include "dpsoftrast.h"
   7
   8 #ifndef __cplusplus
   9 typedef qboolean bool;
  10 #endif
  11
  12 #if defined(__GNUC__)
  13 #define ALIGN(var) var __attribute__((__aligned__(16)))
  14 #elif defined(_MSC_VER)
  15 #define ALIGN(var) __declspec(align(16)) var
  16 #else
  17 #define ALIGN(var) var
  18 #endif
  19
  20 #ifdef SSE2_PRESENT
  21 #include <emmintrin.h>
  22
  23 #define MM_MALLOC(size) _mm_malloc(size, 16)
  24
  25 static void *MM_CALLOC(size_t nmemb, size_t size)
  26 {
  27         void *ptr = _mm_malloc(nmemb*size, 16);
  28         if(ptr != NULL) memset(ptr, 0, nmemb*size);
  29         return ptr;
  30 }
  31
  32 #define MM_FREE _mm_free
  33 #else
  34 #define MM_MALLOC(size) malloc(size)
  35 #define MM_CALLOC(nmemb, size) calloc(nmemb, size)
  36 #define MM_FREE free
  37 #endif
  38
  39 typedef enum DPSOFTRAST_ARRAY_e
  40 {
  41         DPSOFTRAST_ARRAY_POSITION,
  42         DPSOFTRAST_ARRAY_COLOR,
  43         DPSOFTRAST_ARRAY_TEXCOORD0,
  44         DPSOFTRAST_ARRAY_TEXCOORD1,
  45         DPSOFTRAST_ARRAY_TEXCOORD2,
  46         DPSOFTRAST_ARRAY_TEXCOORD3,
  47         DPSOFTRAST_ARRAY_TEXCOORD4,
  48         DPSOFTRAST_ARRAY_TEXCOORD5,
  49         DPSOFTRAST_ARRAY_TEXCOORD6,
  50         DPSOFTRAST_ARRAY_TEXCOORD7,
  51         DPSOFTRAST_ARRAY_TOTAL
  52 }
  53 DPSOFTRAST_ARRAY;
  54
  55 typedef struct DPSOFTRAST_Texture_s
  56 {
  57         int flags;
  58         int width;
  59         int height;
  60         int depth;
  61         int sides;
  62         DPSOFTRAST_TEXTURE_FILTER filter;
  63         int mipmaps;
  64         int size;
  65         unsigned char *bytes;
  66         int mipmap[DPSOFTRAST_MAXMIPMAPS][5];
  67 }
  68 DPSOFTRAST_Texture;
  69
  70 typedef struct DPSOFTRAST_State_User_s
  71 {
  72         int colormask[4];
  73         int blendfunc[2];
  74         int blendsubtract;
  75         int depthmask;
  76         int depthtest;
  77         int depthfunc;
  78         int scissortest;
  79         int cullface;
  80         int alphatest;
  81         int alphafunc;
  82         float alphavalue;
  83         int scissor[4];
  84         int viewport[4];
  85         float depthrange[2];
  86         float polygonoffset[2];
  87         float color[4];
  88 }
  89 DPSOFTRAST_State_User;
  90
  91 #define DPSOFTRAST_MAXSUBSPAN 16
  92
  93 typedef ALIGN(struct DPSOFTRAST_State_Draw_Span_s
  94 {
  95         int start; // pixel index
  96         int length; // pixel count
  97         int startx; // usable range (according to pixelmask)
  98         int endx; // usable range (according to pixelmask)
  99         unsigned char mip[DPSOFTRAST_MAXTEXTUREUNITS]; // texcoord to screen space density values (for picking mipmap of textures)
 100         unsigned char *pixelmask; // true for pixels that passed depth test, false for others
 101         // [0][n][] is start interpolant values (projected)
 102         // [1][n][] is end interpolant values (projected)
 103         // [0][DPSOFTRAST_ARRAY_TOTAL][] is start screencoord4f
 104         // [1][DPSOFTRAST_ARRAY_TOTAL][] is end screencoord4f
 105         // NOTE: screencoord4f[3] is W (basically 1/Z), useful for depthbuffer
 106         ALIGN(float data[2][DPSOFTRAST_ARRAY_TOTAL+1][4]);
 107 }
 108 DPSOFTRAST_State_Draw_Span);
 109
 110 #define DPSOFTRAST_DRAW_MAXSPANQUEUE 1024
 111
 112 typedef struct DPSOFTRAST_State_Draw_s
 113 {
 114         int numvertices;
 115         int maxvertices;
 116         float *in_array4f[DPSOFTRAST_ARRAY_TOTAL];
 117         float *post_array4f[DPSOFTRAST_ARRAY_TOTAL];
 118         float *screencoord4f;
 119
 120         // spans are queued in this structure for dispatch to the pixel shader,
 121         // partly to improve cache locality, partly for batching purposes, spans
 122         // are flushed before DrawTriangles returns to caller
 123         int numspans;
 124         DPSOFTRAST_State_Draw_Span spanqueue[DPSOFTRAST_DRAW_MAXSPANQUEUE];
 125 }
 126 DPSOFTRAST_State_Draw;
 127
 128 #define DPSOFTRAST_VALIDATE_FB 1
 129 #define DPSOFTRAST_VALIDATE_DEPTHFUNC 2
 130 #define DPSOFTRAST_VALIDATE_BLENDFUNC 4
 131 #define DPSOFTRAST_VALIDATE_DRAW (DPSOFTRAST_VALIDATE_FB | DPSOFTRAST_VALIDATE_DEPTHFUNC | DPSOFTRAST_VALIDATE_BLENDFUNC)
 132
 133 typedef enum DPSOFTRAST_BLENDMODE_e
 134 {
 135         DPSOFTRAST_BLENDMODE_OPAQUE,
 136         DPSOFTRAST_BLENDMODE_ALPHA,
 137         DPSOFTRAST_BLENDMODE_ADDALPHA,
 138         DPSOFTRAST_BLENDMODE_ADD,
 139         DPSOFTRAST_BLENDMODE_INVMOD,
 140         DPSOFTRAST_BLENDMODE_MUL,
 141         DPSOFTRAST_BLENDMODE_MUL2,
 142         DPSOFTRAST_BLENDMODE_SUBALPHA,
 143         DPSOFTRAST_BLENDMODE_PSEUDOALPHA,
 144         DPSOFTRAST_BLENDMODE_TOTAL
 145 }
 146 DPSOFTRAST_BLENDMODE;
 147
 148 typedef ALIGN(struct DPSOFTRAST_State_s
 149 {
 150         // DPSOFTRAST_VALIDATE_ flags
 151         int validate;
 152
 153         int fb_colormask;
 154         int fb_width;
 155         int fb_height;
 156         unsigned int *fb_depthpixels;
 157         unsigned int *fb_colorpixels[4];
 158
 159         const float *pointer_vertex3f;
 160         const float *pointer_color4f;
 161         const unsigned char *pointer_color4ub;
 162         const float *pointer_texcoordf[DPSOFTRAST_MAXTEXCOORDARRAYS];
 163         int stride_vertex;
 164         int stride_color;
 165         int stride_texcoord[DPSOFTRAST_MAXTEXCOORDARRAYS];
 166         int components_texcoord[DPSOFTRAST_MAXTEXCOORDARRAYS];
 167         DPSOFTRAST_Texture *texbound[DPSOFTRAST_MAXTEXTUREUNITS];
 168
 169         int shader_mode;
 170         int shader_permutation;
 171         ALIGN(float uniform4f[DPSOFTRAST_UNIFORM_TOTAL*4]);
 172         int uniform1i[DPSOFTRAST_UNIFORM_TOTAL];
 173
 174         // derived values (DPSOFTRAST_VALIDATE_FB)
 175         int fb_clearscissor[4];
 176         int fb_viewport[4];
 177         int fb_viewportscissor[4];
 178         ALIGN(float fb_viewportcenter[4]);
 179         ALIGN(float fb_viewportscale[4]);
 180
 181         // derived values (DPSOFTRAST_VALIDATE_DEPTHFUNC)
 182         int fb_depthfunc;
 183
 184         // derived values (DPSOFTRAST_VALIDATE_BLENDFUNC)
 185         int fb_blendmode;
 186
 187         int texture_max;
 188         int texture_end;
 189         int texture_firstfree;
 190         DPSOFTRAST_Texture *texture;
 191
 192         int bigendian;
 193
 194         // error reporting
 195         const char *errorstring;
 196
 197         DPSOFTRAST_State_User user;
 198
 199         DPSOFTRAST_State_Draw draw;
 200 }
 201 DPSOFTRAST_State);
 202
 203 DPSOFTRAST_State dpsoftrast;
 204
 205 extern int dpsoftrast_test;
 206
 207 #define DPSOFTRAST_DEPTHSCALE (1024.0f*1048576.0f)
 208 #define DPSOFTRAST_DEPTHOFFSET (128.0f)
 209 #define DPSOFTRAST_BGRA8_FROM_RGBA32F(r,g,b,a) (((int)(r * 255.0f + 0.5f) << 16) | ((int)(g * 255.0f + 0.5f) << 8) | (int)(b * 255.0f + 0.5f) | ((int)(a * 255.0f + 0.5f) << 24))
 210 #define DPSOFTRAST_DEPTH32_FROM_DEPTH32F(d) ((int)(DPSOFTRAST_DEPTHSCALE * (1-d)))
 211 #define DPSOFTRAST_DRAW_MAXSPANLENGTH 256
 212
 213 void DPSOFTRAST_RecalcFB(void)
 214 {
 215         // calculate framebuffer scissor, viewport, viewport clipped by scissor,
 216         // and viewport projection values
 217         int x1, x2, x3, x4, x5, x6;
 218         int y1, y2, y3, y4, y5, y6;
 219         x1 = dpsoftrast.user.scissor[0];
 220         x2 = dpsoftrast.user.scissor[0] + dpsoftrast.user.scissor[2];
 221         x3 = dpsoftrast.user.viewport[0];
 222         x4 = dpsoftrast.user.viewport[0] + dpsoftrast.user.viewport[2];
 223         y1 = dpsoftrast.fb_height - dpsoftrast.user.scissor[1] - dpsoftrast.user.scissor[3];
 224         y2 = dpsoftrast.fb_height - dpsoftrast.user.scissor[1];
 225         y3 = dpsoftrast.fb_height - dpsoftrast.user.viewport[1] - dpsoftrast.user.viewport[3];
 226         y4 = dpsoftrast.fb_height - dpsoftrast.user.viewport[1];
 227         if (!dpsoftrast.user.scissortest) {x1 = 0;y1 = 0;x2 = dpsoftrast.fb_width;y2 = dpsoftrast.fb_height;}
 228         if (x1 < 0) x1 = 0;
 229         if (x2 > dpsoftrast.fb_width) x2 = dpsoftrast.fb_width;
 230         if (x3 < 0) x1 = 0;
 231         if (x4 > dpsoftrast.fb_width) x4 = dpsoftrast.fb_width;
 232         if (y1 < 0) y1 = 0;
 233         if (y2 > dpsoftrast.fb_height) y2 = dpsoftrast.fb_height;
 234         if (y3 < 0) y1 = 0;
 235         if (y4 > dpsoftrast.fb_height) y4 = dpsoftrast.fb_height;
 236         x5 = x1;if (x5 < x3) x5 = x3;
 237         x6 = x2;if (x6 > x4) x4 = x4;
 238         y5 = y1;if (y5 < y3) y5 = y3;
 239         y6 = y2;if (y6 > y4) y6 = y4;
 240         dpsoftrast.fb_clearscissor[0] = x1;
 241         dpsoftrast.fb_clearscissor[1] = y1;
 242         dpsoftrast.fb_clearscissor[2] = x2 - x1;
 243         dpsoftrast.fb_clearscissor[3] = y2 - y1;
 244         dpsoftrast.fb_viewport[0] = x3;
 245         dpsoftrast.fb_viewport[1] = y3;
 246         dpsoftrast.fb_viewport[2] = x4 - x3;
 247         dpsoftrast.fb_viewport[3] = y4 - y3;
 248         dpsoftrast.fb_viewportscissor[0] = x5;
 249         dpsoftrast.fb_viewportscissor[1] = y5;
 250         dpsoftrast.fb_viewportscissor[2] = x6 - x5;
 251         dpsoftrast.fb_viewportscissor[3] = y6 - y5;
 252         dpsoftrast.fb_viewportcenter[1] = dpsoftrast.user.viewport[0] + 0.5f * dpsoftrast.user.viewport[2] - 0.5f;
 253         dpsoftrast.fb_viewportcenter[2] = dpsoftrast.fb_height - dpsoftrast.user.viewport[1] - 0.5f * dpsoftrast.user.viewport[3] - 0.5f;
 254         dpsoftrast.fb_viewportcenter[3] = 0.5f;
 255         dpsoftrast.fb_viewportcenter[0] = 0.0f;
 256         dpsoftrast.fb_viewportscale[1] = 0.5f * dpsoftrast.user.viewport[2];
 257         dpsoftrast.fb_viewportscale[2] = -0.5f * dpsoftrast.user.viewport[3];
 258         dpsoftrast.fb_viewportscale[3] = 0.5f;
 259         dpsoftrast.fb_viewportscale[0] = 0.0f;
 260 }
 261
 262 void DPSOFTRAST_RecalcDepthFunc(void)
 263 {
 264         dpsoftrast.fb_depthfunc = dpsoftrast.user.depthtest ? dpsoftrast.user.depthfunc : GL_ALWAYS;
 265 }
 266
 267 int blendmodetable[][4] =
 268 {
 269         {DPSOFTRAST_BLENDMODE_OPAQUE, GL_ONE, GL_ZERO, false},
 270         {DPSOFTRAST_BLENDMODE_ALPHA, GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA, false},
 271         {DPSOFTRAST_BLENDMODE_ADDALPHA, GL_SRC_ALPHA, GL_ONE, false},
 272         {DPSOFTRAST_BLENDMODE_ADD, GL_ONE, GL_ONE, false},
 273         {DPSOFTRAST_BLENDMODE_INVMOD, GL_ZERO, GL_ONE_MINUS_SRC_COLOR, false},
 274         {DPSOFTRAST_BLENDMODE_MUL, GL_ZERO, GL_SRC_COLOR, false},
 275         {DPSOFTRAST_BLENDMODE_MUL, GL_DST_COLOR, GL_ZERO, false},
 276         {DPSOFTRAST_BLENDMODE_MUL2, GL_DST_COLOR, GL_SRC_COLOR, false},
 277         {DPSOFTRAST_BLENDMODE_PSEUDOALPHA, GL_ONE, GL_ONE_MINUS_SRC_ALPHA, false},
 278         {DPSOFTRAST_BLENDMODE_SUBALPHA, GL_SRC_COLOR, GL_ONE, true}
 279 };
 280
 281 void DPSOFTRAST_RecalcBlendFunc(void)
 282 {
 283         int i;
 284         dpsoftrast.fb_blendmode = DPSOFTRAST_BLENDMODE_OPAQUE;
 285         for (i = 0;i < (int)(sizeof(blendmodetable) / sizeof(blendmodetable[0]));i++)
 286         {
 287                 if (dpsoftrast.user.blendfunc[0] == blendmodetable[i][1] && dpsoftrast.user.blendfunc[1] == blendmodetable[i][2] && dpsoftrast.user.blendsubtract == blendmodetable[i][3])
 288                 {
 289                         dpsoftrast.fb_blendmode = blendmodetable[i][0];
 290                         break;
 291                 }
 292         }
 293 }
 294
 295 #define DPSOFTRAST_ValidateQuick(f) ((dpsoftrast.validate & (f)) ? (DPSOFTRAST_Validate(f), 0) : 0)
 296
 297 void DPSOFTRAST_Validate(int mask)
 298 {
 299         mask &= dpsoftrast.validate;
 300         if (!mask)
 301                 return;
 302         if (mask & DPSOFTRAST_VALIDATE_FB)
 303         {
 304                 dpsoftrast.validate &= ~DPSOFTRAST_VALIDATE_FB;
 305                 DPSOFTRAST_RecalcFB();
 306         }
 307         if (mask & DPSOFTRAST_VALIDATE_DEPTHFUNC)
 308         {
 309                 dpsoftrast.validate &= ~DPSOFTRAST_VALIDATE_DEPTHFUNC;
 310                 DPSOFTRAST_RecalcDepthFunc();
 311         }
 312         if (mask & DPSOFTRAST_VALIDATE_BLENDFUNC)
 313         {
 314                 dpsoftrast.validate &= ~DPSOFTRAST_VALIDATE_BLENDFUNC;
 315                 DPSOFTRAST_RecalcBlendFunc();
 316         }
 317 }
 318
 319 DPSOFTRAST_Texture *DPSOFTRAST_Texture_GetByIndex(int index)
 320 {
 321         if (index >= 1 && index < dpsoftrast.texture_end && dpsoftrast.texture[index].bytes)
 322                 return &dpsoftrast.texture[index];
 323         return NULL;
 324 }
 325
 326 int DPSOFTRAST_Texture_New(int flags, int width, int height, int depth)
 327 {
 328         int w;
 329         int h;
 330         int d;
 331         int size;
 332         int s;
 333         int texnum;
 334         int mipmaps;
 335         int sides = (flags & DPSOFTRAST_TEXTURE_FLAG_CUBEMAP) ? 6 : 1;
 336         int texformat = flags & DPSOFTRAST_TEXTURE_FORMAT_COMPAREMASK;
 337         DPSOFTRAST_Texture *texture;
 338         if (width*height*depth < 1)
 339         {
 340                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: width, height or depth is less than 1";
 341                 return 0;
 342         }
 343         if (width > DPSOFTRAST_TEXTURE_MAXSIZE || height > DPSOFTRAST_TEXTURE_MAXSIZE || depth > DPSOFTRAST_TEXTURE_MAXSIZE)
 344         {
 345                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: texture size is too large";
 346                 return 0;
 347         }
 348         switch(texformat)
 349         {
 350         case DPSOFTRAST_TEXTURE_FORMAT_BGRA8:
 351         case DPSOFTRAST_TEXTURE_FORMAT_RGBA8:
 352         case DPSOFTRAST_TEXTURE_FORMAT_ALPHA8:
 353                 break;
 354         case DPSOFTRAST_TEXTURE_FORMAT_DEPTH:
 355                 if (flags & DPSOFTRAST_TEXTURE_FLAG_CUBEMAP)
 356                 {
 357                         dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FORMAT_DEPTH only permitted on 2D textures";
 358                         return 0;
 359                 }
 360                 if (depth != 1)
 361                 {
 362                         dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FORMAT_DEPTH only permitted on 2D textures";
 363                         return 0;
 364                 }
 365                 if ((flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP) && (texformat == DPSOFTRAST_TEXTURE_FORMAT_DEPTH))
 366                 {
 367                         dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FORMAT_DEPTH does not permit mipmaps";
 368                         return 0;
 369                 }
 370                 break;
 371         }
 372         if (depth != 1 && (flags & DPSOFTRAST_TEXTURE_FLAG_CUBEMAP))
 373         {
 374                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FLAG_CUBEMAP can not be used on 3D textures";
 375                 return 0;
 376         }
 377         if (depth != 1 && (flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP))
 378         {
 379                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FLAG_MIPMAP can not be used on 3D textures";
 380                 return 0;
 381         }
 382         if (depth != 1 && (flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP))
 383         {
 384                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FLAG_MIPMAP can not be used on 3D textures";
 385                 return 0;
 386         }
 387         if ((flags & DPSOFTRAST_TEXTURE_FLAG_CUBEMAP) && (flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP))
 388         {
 389                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: DPSOFTRAST_TEXTURE_FLAG_MIPMAP can not be used on cubemap textures";
 390                 return 0;
 391         }
 392         if ((width & (width-1)) || (height & (height-1)) || (depth & (depth-1)))
 393         {
 394                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_New: dimensions are not power of two";
 395                 return 0;
 396         }
 397         // find first empty slot in texture array
 398         for (texnum = dpsoftrast.texture_firstfree;texnum < dpsoftrast.texture_end;texnum++)
 399                 if (!dpsoftrast.texture[texnum].bytes)
 400                         break;
 401         dpsoftrast.texture_firstfree = texnum + 1;
 402         if (dpsoftrast.texture_max <= texnum)
 403         {
 404                 // expand texture array as needed
 405                 if (dpsoftrast.texture_max < 1024)
 406                         dpsoftrast.texture_max = 1024;
 407                 else
 408                         dpsoftrast.texture_max *= 2;
 409                 dpsoftrast.texture = (DPSOFTRAST_Texture *)realloc(dpsoftrast.texture, dpsoftrast.texture_max * sizeof(DPSOFTRAST_Texture));
 410         }
 411         if (dpsoftrast.texture_end <= texnum)
 412                 dpsoftrast.texture_end = texnum + 1;
 413         texture = &dpsoftrast.texture[texnum];
 414         memset(texture, 0, sizeof(*texture));
 415         texture->flags = flags;
 416         texture->width = width;
 417         texture->height = height;
 418         texture->depth = depth;
 419         texture->sides = sides;
 420         w = width;
 421         h = height;
 422         d = depth;
 423         size = 0;
 424         mipmaps = 0;
 425         w = width;
 426         h = height;
 427         d = depth;
 428         for (;;)
 429         {
 430                 s = w * h * d * sides * 4;
 431                 texture->mipmap[mipmaps][0] = size;
 432                 texture->mipmap[mipmaps][1] = s;
 433                 texture->mipmap[mipmaps][2] = w;
 434                 texture->mipmap[mipmaps][3] = h;
 435                 texture->mipmap[mipmaps][4] = d;
 436                 size += s;
 437                 mipmaps++;
 438                 if (w * h * d == 1 || !(flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP))
 439                         break;
 440                 if (w > 1) w >>= 1;
 441                 if (h > 1) h >>= 1;
 442                 if (d > 1) d >>= 1;
 443         }
 444         texture->mipmaps = mipmaps;
 445         texture->size = size;
 446
 447         // allocate the pixels now
 448         texture->bytes = (unsigned char *)MM_CALLOC(1, size);
 449
 450         return texnum;
 451 }
 452 void DPSOFTRAST_Texture_Free(int index)
 453 {
 454         DPSOFTRAST_Texture *texture;
 455         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 456         if (texture->bytes)
 457                 MM_FREE(texture->bytes);
 458         texture->bytes = NULL;
 459         memset(texture, 0, sizeof(*texture));
 460         // adjust the free range and used range
 461         if (dpsoftrast.texture_firstfree > index)
 462                 dpsoftrast.texture_firstfree = index;
 463         while (dpsoftrast.texture_end > 0 && dpsoftrast.texture[dpsoftrast.texture_end-1].bytes == NULL)
 464                 dpsoftrast.texture_end--;
 465 }
 466 void DPSOFTRAST_Texture_CalculateMipmaps(int index)
 467 {
 468         int i, x, y, z, w, layer0, layer1, row0, row1;
 469         unsigned char *o, *i0, *i1, *i2, *i3;
 470         DPSOFTRAST_Texture *texture;
 471         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 472         if (texture->mipmaps <= 1)
 473                 return;
 474         for (i = 1;i < texture->mipmaps;i++)
 475         {
 476                 for (z = 0;z < texture->mipmap[i][4];z++)
 477                 {
 478                         layer0 = z*2;
 479                         layer1 = z*2+1;
 480                         if (layer1 >= texture->mipmap[i-1][4])
 481                                 layer1 = texture->mipmap[i-1][4]-1;
 482                         for (y = 0;y < texture->mipmap[i][3];y++)
 483                         {
 484                                 row0 = y*2;
 485                                 row1 = y*2+1;
 486                                 if (row1 >= texture->mipmap[i-1][3])
 487                                         row1 = texture->mipmap[i-1][3]-1;
 488                                 o =  texture->bytes + texture->mipmap[i  ][0] + 4*((texture->mipmap[i  ][3] * z      + y   ) * texture->mipmap[i  ][2]);
 489                                 i0 = texture->bytes + texture->mipmap[i-1][0] + 4*((texture->mipmap[i-1][3] * layer0 + row0) * texture->mipmap[i-1][2]);
 490                                 i1 = texture->bytes + texture->mipmap[i-1][0] + 4*((texture->mipmap[i-1][3] * layer0 + row1) * texture->mipmap[i-1][2]);
 491                                 i2 = texture->bytes + texture->mipmap[i-1][0] + 4*((texture->mipmap[i-1][3] * layer1 + row0) * texture->mipmap[i-1][2]);
 492                                 i3 = texture->bytes + texture->mipmap[i-1][0] + 4*((texture->mipmap[i-1][3] * layer1 + row1) * texture->mipmap[i-1][2]);
 493                                 w = texture->mipmap[i][2];
 494                                 if (layer1 > layer0)
 495                                 {
 496                                         if (texture->mipmap[i-1][2] > 1)
 497                                         {
 498                                                 // average 3D texture
 499                                                 for (x = 0;x < w;x++, o += 4, i0 += 8, i1 += 8, i2 += 8, i3 += 8)
 500                                                 {
 501                                                         o[0] = (i0[0] + i0[4] + i1[0] + i1[4] + i2[0] + i2[4] + i3[0] + i3[4] + 4) >> 3;
 502                                                         o[1] = (i0[1] + i0[5] + i1[1] + i1[5] + i2[1] + i2[5] + i3[1] + i3[5] + 4) >> 3;
 503                                                         o[2] = (i0[2] + i0[6] + i1[2] + i1[6] + i2[2] + i2[6] + i3[2] + i3[6] + 4) >> 3;
 504                                                         o[3] = (i0[3] + i0[7] + i1[3] + i1[7] + i2[3] + i2[7] + i3[3] + i3[7] + 4) >> 3;
 505                                                 }
 506                                         }
 507                                         else
 508                                         {
 509                                                 // average 3D mipmap with parent width == 1
 510                                                 for (x = 0;x < w;x++, o += 4, i0 += 8, i1 += 8)
 511                                                 {
 512                                                         o[0] = (i0[0] + i1[0] + i2[0] + i3[0] + 2) >> 2;
 513                                                         o[1] = (i0[1] + i1[1] + i2[1] + i3[1] + 2) >> 2;
 514                                                         o[2] = (i0[2] + i1[2] + i2[2] + i3[2] + 2) >> 2;
 515                                                         o[3] = (i0[3] + i1[3] + i2[3] + i3[3] + 2) >> 2;
 516                                                 }
 517                                         }
 518                                 }
 519                                 else
 520                                 {
 521                                         if (texture->mipmap[i-1][2] > 1)
 522                                         {
 523                                                 // average 2D texture (common case)
 524                                                 for (x = 0;x < w;x++, o += 4, i0 += 8, i1 += 8)
 525                                                 {
 526                                                         o[0] = (i0[0] + i0[4] + i1[0] + i1[4] + 2) >> 2;
 527                                                         o[1] = (i0[1] + i0[5] + i1[1] + i1[5] + 2) >> 2;
 528                                                         o[2] = (i0[2] + i0[6] + i1[2] + i1[6] + 2) >> 2;
 529                                                         o[3] = (i0[3] + i0[7] + i1[3] + i1[7] + 2) >> 2;
 530                                                 }
 531                                         }
 532                                         else
 533                                         {
 534                                                 // 2D texture with parent width == 1
 535                                                 o[0] = (i0[0] + i1[0] + 1) >> 1;
 536                                                 o[1] = (i0[1] + i1[1] + 1) >> 1;
 537                                                 o[2] = (i0[2] + i1[2] + 1) >> 1;
 538                                                 o[3] = (i0[3] + i1[3] + 1) >> 1;
 539                                         }
 540                                 }
 541                         }
 542                 }
 543         }
 544 }
 545 void DPSOFTRAST_Texture_UpdatePartial(int index, int mip, const unsigned char *pixels, int blockx, int blocky, int blockwidth, int blockheight)
 546 {
 547         DPSOFTRAST_Texture *texture;
 548         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 549
 550         // FIXME IMPLEMENT
 551
 552         dpsoftrast.errorstring = "DPSOFTRAST_Texture_UpdatePartial: Not implemented.";
 553 }
 554 void DPSOFTRAST_Texture_UpdateFull(int index, const unsigned char *pixels)
 555 {
 556         DPSOFTRAST_Texture *texture;
 557         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 558
 559         memcpy(texture->bytes, pixels, texture->mipmap[0][1]);
 560         DPSOFTRAST_Texture_CalculateMipmaps(index);
 561 }
 562 int DPSOFTRAST_Texture_GetWidth(int index, int mip)
 563 {
 564         DPSOFTRAST_Texture *texture;
 565         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return 0;
 566         return texture->mipmap[mip][2];
 567 }
 568 int DPSOFTRAST_Texture_GetHeight(int index, int mip)
 569 {
 570         DPSOFTRAST_Texture *texture;
 571         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return 0;
 572         return texture->mipmap[mip][3];
 573 }
 574 int DPSOFTRAST_Texture_GetDepth(int index, int mip)
 575 {
 576         DPSOFTRAST_Texture *texture;
 577         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return 0;
 578         return texture->mipmap[mip][4];
 579 }
 580 unsigned char *DPSOFTRAST_Texture_GetPixelPointer(int index, int mip)
 581 {
 582         DPSOFTRAST_Texture *texture;
 583         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return 0;
 584         return texture->bytes + texture->mipmap[mip][0];
 585 }
 586 void DPSOFTRAST_Texture_Filter(int index, DPSOFTRAST_TEXTURE_FILTER filter)
 587 {
 588         DPSOFTRAST_Texture *texture;
 589         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 590         if (!(texture->flags & DPSOFTRAST_TEXTURE_FLAG_MIPMAP) && filter > DPSOFTRAST_TEXTURE_FILTER_LINEAR)
 591         {
 592                 dpsoftrast.errorstring = "DPSOFTRAST_Texture_Filter: requested filter mode requires mipmaps";
 593                 return;
 594         }
 595         texture->filter = filter;
 596 }
 597
 598 void DPSOFTRAST_SetRenderTargets(int width, int height, unsigned int *depthpixels, unsigned int *colorpixels0, unsigned int *colorpixels1, unsigned int *colorpixels2, unsigned int *colorpixels3)
 599 {
 600         dpsoftrast.fb_width = width;
 601         dpsoftrast.fb_height = height;
 602         dpsoftrast.fb_depthpixels = depthpixels;
 603         dpsoftrast.fb_colorpixels[0] = colorpixels0;
 604         dpsoftrast.fb_colorpixels[1] = colorpixels1;
 605         dpsoftrast.fb_colorpixels[2] = colorpixels2;
 606         dpsoftrast.fb_colorpixels[3] = colorpixels3;
 607 }
 608 void DPSOFTRAST_Viewport(int x, int y, int width, int height)
 609 {
 610         dpsoftrast.user.viewport[0] = x;
 611         dpsoftrast.user.viewport[1] = y;
 612         dpsoftrast.user.viewport[2] = width;
 613         dpsoftrast.user.viewport[3] = height;
 614         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_FB;
 615 }
 616 void DPSOFTRAST_ClearColor(float r, float g, float b, float a)
 617 {
 618         int i, x1, y1, x2, y2, w, h, x, y;
 619         unsigned int *p;
 620         unsigned int c;
 621         DPSOFTRAST_Validate(DPSOFTRAST_VALIDATE_FB);
 622         x1 = dpsoftrast.fb_clearscissor[0];
 623         y1 = dpsoftrast.fb_clearscissor[1];
 624         x2 = dpsoftrast.fb_clearscissor[2];
 625         y2 = dpsoftrast.fb_clearscissor[1] + dpsoftrast.fb_clearscissor[3];
 626         w = x2 - x1;
 627         h = y2 - y1;
 628         if (w < 1 || h < 1)
 629                 return;
 630         // FIXME: honor dpsoftrast.fb_colormask?
 631         c = DPSOFTRAST_BGRA8_FROM_RGBA32F(r,g,b,a);
 632         for (i = 0;i < 4;i++)
 633         {
 634                 if (!dpsoftrast.fb_colorpixels[i])
 635                         continue;
 636                 for (y = y1;y < y2;y++)
 637                 {
 638                         p = dpsoftrast.fb_colorpixels[i] + y * dpsoftrast.fb_width;
 639                         for (x = x1;x < x2;x++)
 640                                 p[x] = c;
 641                 }
 642         }
 643 }
 644 void DPSOFTRAST_ClearDepth(float d)
 645 {
 646         int x1, y1, x2, y2, w, h, x, y;
 647         unsigned int *p;
 648         unsigned int c;
 649         DPSOFTRAST_Validate(DPSOFTRAST_VALIDATE_FB);
 650         x1 = dpsoftrast.fb_clearscissor[0];
 651         y1 = dpsoftrast.fb_clearscissor[1];
 652         x2 = dpsoftrast.fb_clearscissor[2];
 653         y2 = dpsoftrast.fb_clearscissor[1] + dpsoftrast.fb_clearscissor[3];
 654         w = x2 - x1;
 655         h = y2 - y1;
 656         if (w < 1 || h < 1)
 657                 return;
 658         c = DPSOFTRAST_DEPTH32_FROM_DEPTH32F(d);
 659         for (y = y1;y < y2;y++)
 660         {
 661                 p = dpsoftrast.fb_depthpixels + y * dpsoftrast.fb_width;
 662                 for (x = x1;x < x2;x++)
 663                         p[x] = c;
 664         }
 665 }
 666 void DPSOFTRAST_ColorMask(int r, int g, int b, int a)
 667 {
 668         dpsoftrast.user.colormask[0] = r != 0;
 669         dpsoftrast.user.colormask[1] = g != 0;
 670         dpsoftrast.user.colormask[2] = b != 0;
 671         dpsoftrast.user.colormask[3] = a != 0;
 672         dpsoftrast.fb_colormask = ((-dpsoftrast.user.colormask[0]) & 0x00FF0000) | ((-dpsoftrast.user.colormask[1]) & 0x0000FF00) | ((-dpsoftrast.user.colormask[2]) & 0x000000FF) | ((-dpsoftrast.user.colormask[3]) & 0xFF000000);
 673 }
 674 void DPSOFTRAST_DepthTest(int enable)
 675 {
 676         dpsoftrast.user.depthtest = enable;
 677         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_DEPTHFUNC;
 678 }
 679 void DPSOFTRAST_ScissorTest(int enable)
 680 {
 681         dpsoftrast.user.scissortest = enable;
 682         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_FB;
 683 }
 684 void DPSOFTRAST_Scissor(float x, float y, float width, float height)
 685 {
 686         dpsoftrast.user.scissor[0] = x;
 687         dpsoftrast.user.scissor[1] = y;
 688         dpsoftrast.user.scissor[2] = width;
 689         dpsoftrast.user.scissor[3] = height;
 690         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_FB;
 691 }
 692
 693 void DPSOFTRAST_BlendFunc(int smodulate, int dmodulate)
 694 {
 695         // FIXME: validate
 696         dpsoftrast.user.blendfunc[0] = smodulate;
 697         dpsoftrast.user.blendfunc[1] = dmodulate;
 698         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_BLENDFUNC;
 699 }
 700 void DPSOFTRAST_BlendSubtract(int enable)
 701 {
 702         dpsoftrast.user.blendsubtract = enable != 0;
 703         dpsoftrast.validate |= DPSOFTRAST_VALIDATE_BLENDFUNC;
 704 }
 705 void DPSOFTRAST_DepthMask(int enable)
 706 {
 707         dpsoftrast.user.depthmask = enable;
 708 }
 709 void DPSOFTRAST_DepthFunc(int comparemode)
 710 {
 711         // FIXME: validate
 712         dpsoftrast.user.depthfunc = comparemode;
 713 }
 714 void DPSOFTRAST_DepthRange(float range0, float range1)
 715 {
 716         dpsoftrast.user.depthrange[0] = range0;
 717         dpsoftrast.user.depthrange[1] = range1;
 718 }
 719 void DPSOFTRAST_PolygonOffset(float alongnormal, float intoview)
 720 {
 721         dpsoftrast.user.polygonoffset[0] = alongnormal;
 722         dpsoftrast.user.polygonoffset[1] = intoview;
 723 }
 724 void DPSOFTRAST_CullFace(int mode)
 725 {
 726         // FIXME: validate
 727         dpsoftrast.user.cullface = mode;
 728 }
 729 void DPSOFTRAST_AlphaTest(float enable)
 730 {
 731         dpsoftrast.user.alphatest = enable;
 732 }
 733 void DPSOFTRAST_AlphaFunc(int alphafunc, float alphavalue)
 734 {
 735         // FIXME: validate
 736         dpsoftrast.user.alphafunc = alphafunc;
 737         dpsoftrast.user.alphavalue = alphavalue;
 738 }
 739 void DPSOFTRAST_Color4f(float r, float g, float b, float a)
 740 {
 741         dpsoftrast.user.color[0] = r;
 742         dpsoftrast.user.color[1] = g;
 743         dpsoftrast.user.color[2] = b;
 744         dpsoftrast.user.color[3] = a;
 745 }
 746 void DPSOFTRAST_GetPixelsBGRA(int blockx, int blocky, int blockwidth, int blockheight, unsigned char *outpixels)
 747 {
 748         int outstride = blockwidth * 4;
 749         int instride = dpsoftrast.fb_width * 4;
 750         int bx1 = blockx;
 751         int by1 = blocky;
 752         int bx2 = blockx + blockwidth;
 753         int by2 = blocky + blockheight;
 754         int bw;
 755         int bh;
 756         int x;
 757         int y;
 758         unsigned char *inpixels;
 759         unsigned char *b;
 760         unsigned char *o;
 761         if (bx1 < 0) bx1 = 0;
 762         if (by1 < 0) by1 = 0;
 763         if (bx2 > dpsoftrast.fb_width) bx2 = dpsoftrast.fb_width;
 764         if (by2 > dpsoftrast.fb_height) by2 = dpsoftrast.fb_height;
 765         bw = bx2 - bx1;
 766         bh = by2 - by1;
 767         inpixels = (unsigned char *)dpsoftrast.fb_colorpixels[0];
 768         if (dpsoftrast.bigendian)
 769         {
 770                 for (y = by1;y < by2;y++)
 771                 {
 772                         b = (unsigned char *)inpixels + (dpsoftrast.fb_height - 1 - y) * instride + 4 * bx1;
 773                         o = (unsigned char *)outpixels + (y - by1) * outstride;
 774                         for (x = bx1;x < bx2;x++)
 775                         {
 776                                 o[0] = b[3];
 777                                 o[1] = b[2];
 778                                 o[2] = b[1];
 779                                 o[3] = b[0];
 780                                 o += 4;
 781                                 b += 4;
 782                         }
 783                 }
 784         }
 785         else
 786         {
 787                 for (y = by1;y < by2;y++)
 788                 {
 789                         b = (unsigned char *)inpixels + (dpsoftrast.fb_height - 1 - y) * instride + 4 * bx1;
 790                         o = (unsigned char *)outpixels + (y - by1) * outstride;
 791                         memcpy(o, b, bw*4);
 792                 }
 793         }
 794
 795 }
 796 void DPSOFTRAST_CopyRectangleToTexture(int index, int mip, int tx, int ty, int sx, int sy, int width, int height)
 797 {
 798         int tx1 = tx;
 799         int ty1 = ty;
 800         int tx2 = tx + width;
 801         int ty2 = ty + height;
 802         int sx1 = sx;
 803         int sy1 = sy;
 804         int sx2 = sx + width;
 805         int sy2 = sy + height;
 806         int swidth;
 807         int sheight;
 808         int twidth;
 809         int theight;
 810         int sw;
 811         int sh;
 812         int tw;
 813         int th;
 814         int y;
 815         unsigned int *spixels;
 816         unsigned int *tpixels;
 817         DPSOFTRAST_Texture *texture;
 818         texture = DPSOFTRAST_Texture_GetByIndex(index);if (!texture) return;
 819         if (mip < 0 || mip >= texture->mipmaps) return;
 820         spixels = dpsoftrast.fb_colorpixels[0];
 821         swidth = dpsoftrast.fb_width;
 822         sheight = dpsoftrast.fb_height;
 823         tpixels = (unsigned int *)(texture->bytes + texture->mipmap[mip][0]);
 824         twidth = texture->mipmap[mip][2];
 825         theight = texture->mipmap[mip][3];
 826         if (tx1 < 0) tx1 = 0;
 827         if (ty1 < 0) ty1 = 0;
 828         if (tx2 > twidth) tx2 = twidth;
 829         if (ty2 > theight) ty2 = theight;
 830         if (sx1 < 0) sx1 = 0;
 831         if (sy1 < 0) sy1 = 0;
 832         if (sx2 > swidth) sx2 = swidth;
 833         if (sy2 > sheight) sy2 = sheight;
 834         tw = tx2 - tx1;
 835         th = ty2 - ty1;
 836         sw = sx2 - sx1;
 837         sh = sy2 - sy1;
 838         if (tw > sw) tw = sw;
 839         if (th > sh) th = sh;
 840         if (tw < 1 || th < 1)
 841                 return;
 842         for (y = 0;y < th;y++)
 843                 memcpy(tpixels + ((ty1 + y) * twidth + tx1), spixels + ((sy1 + y) * swidth + sx1), tw*4);
 844         if (texture->mipmaps > 1)
 845                 DPSOFTRAST_Texture_CalculateMipmaps(index);
 846 }
 847 void DPSOFTRAST_SetTexture(int unitnum, int index)
 848 {
 849         DPSOFTRAST_Texture *texture;
 850         if (unitnum < 0 || unitnum >= DPSOFTRAST_MAXTEXTUREUNITS)
 851         {
 852                 dpsoftrast.errorstring = "DPSOFTRAST_SetTexture: invalid unit number";
 853                 return;
 854         }
 855         texture = DPSOFTRAST_Texture_GetByIndex(index);
 856         if (index && !texture)
 857         {
 858                 dpsoftrast.errorstring = "DPSOFTRAST_SetTexture: invalid texture handle";
 859                 return;
 860         }
 861         dpsoftrast.texbound[unitnum] = texture;
 862 }
 863
 864 void DPSOFTRAST_SetVertexPointer(const float *vertex3f, size_t stride)
 865 {
 866         dpsoftrast.pointer_vertex3f = vertex3f;
 867         dpsoftrast.stride_vertex = stride;
 868 }
 869 void DPSOFTRAST_SetColorPointer(const float *color4f, size_t stride)
 870 {
 871         dpsoftrast.pointer_color4f = color4f;
 872         dpsoftrast.pointer_color4ub = NULL;
 873         dpsoftrast.stride_color = stride;
 874 }
 875 void DPSOFTRAST_SetColorPointer4ub(const unsigned char *color4ub, size_t stride)
 876 {
 877         dpsoftrast.pointer_color4f = NULL;
 878         dpsoftrast.pointer_color4ub = color4ub;
 879         dpsoftrast.stride_color = stride;
 880 }
 881 void DPSOFTRAST_SetTexCoordPointer(int unitnum, int numcomponents, size_t stride, const float *texcoordf)
 882 {
 883         dpsoftrast.pointer_texcoordf[unitnum] = texcoordf;
 884         dpsoftrast.components_texcoord[unitnum] = numcomponents;
 885         dpsoftrast.stride_texcoord[unitnum] = stride;
 886 }
 887
 888 void DPSOFTRAST_SetShader(unsigned int mode, unsigned int permutation)
 889 {
 890         dpsoftrast.shader_mode = mode;
 891         dpsoftrast.shader_permutation = permutation;
 892 }
 893 void DPSOFTRAST_Uniform4fARB(DPSOFTRAST_UNIFORM index, float v0, float v1, float v2, float v3)
 894 {
 895         dpsoftrast.uniform4f[index*4+0] = v0;
 896         dpsoftrast.uniform4f[index*4+1] = v1;
 897         dpsoftrast.uniform4f[index*4+2] = v2;
 898         dpsoftrast.uniform4f[index*4+3] = v3;
 899 }
 900 void DPSOFTRAST_Uniform4fvARB(DPSOFTRAST_UNIFORM index, const float *v)
 901 {
 902         dpsoftrast.uniform4f[index*4+0] = v[0];
 903         dpsoftrast.uniform4f[index*4+1] = v[1];
 904         dpsoftrast.uniform4f[index*4+2] = v[2];
 905         dpsoftrast.uniform4f[index*4+3] = v[3];
 906 }
 907 void DPSOFTRAST_UniformMatrix4fvARB(DPSOFTRAST_UNIFORM index, int arraysize, int transpose, const float *v)
 908 {
 909         int i;
 910         for (i = 0;i < arraysize;i++, index += 4, v += 16)
 911         {
 912                 if (transpose)
 913                 {
 914                         dpsoftrast.uniform4f[index*4+0] = v[0];
 915                         dpsoftrast.uniform4f[index*4+1] = v[4];
 916                         dpsoftrast.uniform4f[index*4+2] = v[8];
 917                         dpsoftrast.uniform4f[index*4+3] = v[12];
 918                         dpsoftrast.uniform4f[index*4+4] = v[1];
 919                         dpsoftrast.uniform4f[index*4+5] = v[5];
 920                         dpsoftrast.uniform4f[index*4+6] = v[9];
 921                         dpsoftrast.uniform4f[index*4+7] = v[13];
 922                         dpsoftrast.uniform4f[index*4+8] = v[2];
 923                         dpsoftrast.uniform4f[index*4+9] = v[6];
 924                         dpsoftrast.uniform4f[index*4+10] = v[10];
 925                         dpsoftrast.uniform4f[index*4+11] = v[14];
 926                         dpsoftrast.uniform4f[index*4+12] = v[3];
 927                         dpsoftrast.uniform4f[index*4+13] = v[7];
 928                         dpsoftrast.uniform4f[index*4+14] = v[11];
 929                         dpsoftrast.uniform4f[index*4+15] = v[15];
 930                 }
 931                 else
 932                 {
 933                         dpsoftrast.uniform4f[index*4+0] = v[0];
 934                         dpsoftrast.uniform4f[index*4+1] = v[1];
 935                         dpsoftrast.uniform4f[index*4+2] = v[2];
 936                         dpsoftrast.uniform4f[index*4+3] = v[3];
 937                         dpsoftrast.uniform4f[index*4+4] = v[4];
 938                         dpsoftrast.uniform4f[index*4+5] = v[5];
 939                         dpsoftrast.uniform4f[index*4+6] = v[6];
 940                         dpsoftrast.uniform4f[index*4+7] = v[7];
 941                         dpsoftrast.uniform4f[index*4+8] = v[8];
 942                         dpsoftrast.uniform4f[index*4+9] = v[9];
 943                         dpsoftrast.uniform4f[index*4+10] = v[10];
 944                         dpsoftrast.uniform4f[index*4+11] = v[11];
 945                         dpsoftrast.uniform4f[index*4+12] = v[12];
 946                         dpsoftrast.uniform4f[index*4+13] = v[13];
 947                         dpsoftrast.uniform4f[index*4+14] = v[14];
 948                         dpsoftrast.uniform4f[index*4+15] = v[15];
 949                 }
 950         }
 951 }
 952 void DPSOFTRAST_Uniform1iARB(DPSOFTRAST_UNIFORM index, int i0)
 953 {
 954         dpsoftrast.uniform1i[index] = i0;
 955 }
 956
 957 void DPSOFTRAST_Draw_LoadVertices(int firstvertex, int numvertices, bool needcolors)
 958 {
 959         int i;
 960         int j;
 961         int stride;
 962         const float *v;
 963         float *p;
 964         float *data;
 965         const unsigned char *b;
 966         dpsoftrast.draw.numvertices = numvertices;
 967         if (dpsoftrast.draw.maxvertices < dpsoftrast.draw.numvertices)
 968         {
 969                 if (dpsoftrast.draw.maxvertices < 4096)
 970                         dpsoftrast.draw.maxvertices = 4096;
 971                 while (dpsoftrast.draw.maxvertices < dpsoftrast.draw.numvertices)
 972                         dpsoftrast.draw.maxvertices *= 2;
 973                 if (dpsoftrast.draw.in_array4f[0])
 974                         MM_FREE(dpsoftrast.draw.in_array4f[0]);
 975                 data = (float *)MM_CALLOC(1, dpsoftrast.draw.maxvertices * sizeof(float[4])*(DPSOFTRAST_ARRAY_TOTAL*2 + 1));
 976                 for (i = 0;i < DPSOFTRAST_ARRAY_TOTAL;i++, data += dpsoftrast.draw.maxvertices * 4)
 977                         dpsoftrast.draw.in_array4f[i] = data;
 978                 for (i = 0;i < DPSOFTRAST_ARRAY_TOTAL;i++, data += dpsoftrast.draw.maxvertices * 4)
 979                         dpsoftrast.draw.post_array4f[i] = data;
 980                 dpsoftrast.draw.screencoord4f = data;
 981                 data += dpsoftrast.draw.maxvertices * 4;
 982         }
 983         stride = dpsoftrast.stride_vertex;
 984         v = (const float *)((unsigned char *)dpsoftrast.pointer_vertex3f + firstvertex * stride);
 985         p = dpsoftrast.draw.in_array4f[0];
 986         for (i = 0;i < numvertices;i++)
 987         {
 988                 p[0] = v[0];
 989                 p[1] = v[1];
 990                 p[2] = v[2];
 991                 p[3] = 1.0f;
 992                 p += 4;
 993                 v = (const float *)((const unsigned char *)v + stride);
 994         }
 995         if (needcolors)
 996         {
 997                 if (dpsoftrast.pointer_color4f)
 998                 {
 999                         stride = dpsoftrast.stride_color;
1000                         v = (const float *)((const unsigned char *)dpsoftrast.pointer_color4f + firstvertex * stride);
1001                         p = dpsoftrast.draw.in_array4f[1];
1002                         for (i = 0;i < numvertices;i++)
1003                         {
1004                                 p[0] = v[0];
1005                                 p[1] = v[1];
1006                                 p[2] = v[2];
1007                                 p[3] = v[3];
1008                                 p += 4;
1009                                 v = (const float *)((const unsigned char *)v + stride);
1010                         }
1011                 }
1012                 else if (dpsoftrast.pointer_color4ub)
1013                 {
1014                         stride = dpsoftrast.stride_color;
1015                         b = (const unsigned char *)((const unsigned char *)dpsoftrast.pointer_color4ub + firstvertex * stride);
1016                         p = dpsoftrast.draw.in_array4f[1];
1017                         for (i = 0;i < numvertices;i++)
1018                         {
1019                                 p[0] = b[0] * (1.0f / 255.0f);
1020                                 p[1] = b[1] * (1.0f / 255.0f);
1021                                 p[2] = b[2] * (1.0f / 255.0f);
1022                                 p[3] = b[3] * (1.0f / 255.0f);
1023                                 p += 4;
1024                                 b = (const unsigned char *)((const unsigned char *)b + stride);
1025                         }
1026                 }
1027                 else
1028                 {
1029                         v = dpsoftrast.user.color;
1030                         p = dpsoftrast.draw.in_array4f[1];
1031                         for (i = 0;i < numvertices;i++)
1032                         {
1033                                 p[0] = v[0];
1034                                 p[1] = v[1];
1035                                 p[2] = v[2];
1036                                 p[3] = v[3];
1037                                 p += 4;
1038                         }
1039                 }
1040         }
1041         for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL-2;j++)
1042         {
1043                 if (dpsoftrast.pointer_texcoordf[j])
1044                 {
1045                         stride = dpsoftrast.stride_texcoord[j];
1046                         v = (const float *)((const unsigned char *)dpsoftrast.pointer_texcoordf[j] + firstvertex * stride);
1047                         p = dpsoftrast.draw.in_array4f[j+2];
1048                         switch(dpsoftrast.components_texcoord[j])
1049                         {
1050                         case 2:
1051                                 for (i = 0;i < numvertices;i++)
1052                                 {
1053                                         p[0] = v[0];
1054                                         p[1] = v[1];
1055                                         p[2] = 0.0f;
1056                                         p[3] = 1.0f;
1057                                         p += 4;
1058                                         v = (const float *)((const unsigned char *)v + stride);
1059                                 }
1060                                 break;
1061                         case 3:
1062                                 for (i = 0;i < numvertices;i++)
1063                                 {
1064                                         p[0] = v[0];
1065                                         p[1] = v[1];
1066                                         p[2] = v[2];
1067                                         p[3] = 1.0f;
1068                                         p += 4;
1069                                         v = (const float *)((const unsigned char *)v + stride);
1070                                 }
1071                                 break;
1072                         case 4:
1073                                 for (i = 0;i < numvertices;i++)
1074                                 {
1075                                         p[0] = v[0];
1076                                         p[1] = v[1];
1077                                         p[2] = v[2];
1078                                         p[3] = v[3];
1079                                         p += 4;
1080                                         v = (const float *)((const unsigned char *)v + stride);
1081                                 }
1082                                 break;
1083                         }
1084                 }
1085         }
1086 }
1087
1088 void DPSOFTRAST_Array_Transform(float *out4f, const float *in4f, int numitems, const float *inmatrix16f)
1089 {
1090         static const float identitymatrix[4][4] = {{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}};
1091         // TODO: SIMD
1092         float matrix[4][4];
1093         int i;
1094         memcpy(matrix, inmatrix16f, sizeof(float[16]));
1095         if (!memcmp(identitymatrix, matrix, sizeof(float[16])))
1096         {
1097                 // fast case for identity matrix
1098                 memcpy(out4f, in4f, numitems * sizeof(float[4]));
1099                 return;
1100         }
1101         for (i = 0;i < numitems;i++, out4f += 4, in4f += 4)
1102         {
1103                 out4f[0] = in4f[0] * matrix[0][0] + in4f[1] * matrix[1][0] + in4f[2] * matrix[2][0] + in4f[3] * matrix[3][0];
1104                 out4f[1] = in4f[0] * matrix[0][1] + in4f[1] * matrix[1][1] + in4f[2] * matrix[2][1] + in4f[3] * matrix[3][1];
1105                 out4f[2] = in4f[0] * matrix[0][2] + in4f[1] * matrix[1][2] + in4f[2] * matrix[2][2] + in4f[3] * matrix[3][2];
1106                 out4f[3] = in4f[0] * matrix[0][3] + in4f[1] * matrix[1][3] + in4f[2] * matrix[2][3] + in4f[3] * matrix[3][3];
1107         }
1108 }
1109
1110 void DPSOFTRAST_Array_Copy(float *out4f, const float *in4f, int numitems)
1111 {
1112         memcpy(out4f, in4f, numitems * sizeof(float[4]));
1113 }
1114
1115 #ifdef SSE2_PRESENT
1116 static __m128 DPSOFTRAST_Draw_ProjectVertex(__m128 v)
1117 {
1118         __m128 viewportcenter = _mm_load_ps(dpsoftrast.fb_viewportcenter), viewportscale = _mm_load_ps(dpsoftrast.fb_viewportscale);
1119         __m128 w = _mm_rcp_ss(_mm_shuffle_ps(v, v, _MM_SHUFFLE(2, 1, 0, 3)));
1120         v = _mm_add_ps(viewportcenter, _mm_mul_ps(_mm_mul_ps(viewportscale, w), _mm_shuffle_ps(w, w, _MM_SHUFFLE(0, 0, 0, 0))));
1121         v = _mm_move_ss(v, w);
1122         v = _mm_shuffle_ps(v, v, _MM_SHUFFLE(0, 3, 2, 1));
1123         return v;
1124 }
1125 #endif
1126
1127 void DPSOFTRAST_Draw_ProjectVertices(float *out4f, const float *in4f, int numitems)
1128 {
1129 #ifdef SSE2_PRESENT
1130         // NOTE: this is used both as a whole mesh transform function and a
1131         // per-triangle transform function (for clipped triangles), accordingly
1132         // it should not crash on divide by 0 but the result of divide by 0 is
1133         // unimportant...
1134         // TODO: SIMD
1135         int i;
1136         __m128 viewportcenter = _mm_load_ps(dpsoftrast.fb_viewportcenter), viewportscale = _mm_load_ps(dpsoftrast.fb_viewportscale);
1137         for (i = 0;i < numitems;i++)
1138         {
1139                 __m128 v = _mm_load_ps(in4f), w = _mm_rcp_ss(_mm_shuffle_ps(v, v, _MM_SHUFFLE(2, 1, 0, 3)));
1140                 v = _mm_add_ps(viewportcenter, _mm_mul_ps(_mm_mul_ps(viewportscale, w), _mm_shuffle_ps(w, w, _MM_SHUFFLE(0, 0, 0, 0))));
1141                 v = _mm_move_ss(v, w);
1142                 _mm_store_ps(out4f, _mm_shuffle_ps(v, v, _MM_SHUFFLE(0, 3, 2, 1)));
1143                 in4f += 4;
1144                 out4f += 4;
1145         }
1146 #endif
1147 }
1148
1149 void DPSOFTRAST_Draw_DebugEdgePoints(const float *screen0, const float *screen1)
1150 {
1151         int i;
1152         int x;
1153         int y;
1154         int w = dpsoftrast.fb_width;
1155         int bounds[4];
1156         float v0[2], v1[2];
1157         unsigned int *pixels = dpsoftrast.fb_colorpixels[0];
1158         //const float *c4f;
1159         bounds[0] = dpsoftrast.fb_viewportscissor[0];
1160         bounds[1] = dpsoftrast.fb_viewportscissor[1];
1161         bounds[2] = dpsoftrast.fb_viewportscissor[0] + dpsoftrast.fb_viewportscissor[2];
1162         bounds[3] = dpsoftrast.fb_viewportscissor[1] + dpsoftrast.fb_viewportscissor[3];
1163         v0[0] = screen0[0];
1164         v0[1] = screen0[1];
1165         v1[0] = screen1[0];
1166         v1[1] = screen1[1];
1167         for (i = 0;i <= 128;i++)
1168         {
1169                 // check nearclip
1170                 //if (dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+3] != 1.0f)
1171                 //      continue;
1172                 x = (int)(v0[0] + (v1[0] - v0[0]) * (i/128.0f));
1173                 y = (int)(v0[1] + (v1[1] - v0[1]) * (i/128.0f));
1174                 if (x < bounds[0] || y < bounds[1] || x >= bounds[2] || y >= bounds[3])
1175                         continue;
1176                 //c4f = dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_COLOR] + element0*4;
1177                 //pixels[y*w+x] = DPSOFTRAST_BGRA8_FROM_RGBA32F(c4f[0], c4f[1], c4f[2], c4f[3]);
1178                 pixels[y*w+x] = 0xFFFFFFFF;
1179         }
1180 }
1181
1182 void DPSOFTRAST_Draw_Span_Begin(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *zf)
1183 {
1184         int x;
1185         int startx = span->startx;
1186         int endx = span->endx;
1187         float w = span->data[0][DPSOFTRAST_ARRAY_TOTAL][3];
1188         float wslope = span->data[1][DPSOFTRAST_ARRAY_TOTAL][3];
1189         float endz = 1.0f / (w + wslope * startx);
1190         for (x = startx;x < endx;)
1191         {
1192                 int nextsub = x + DPSOFTRAST_MAXSUBSPAN, endsub = nextsub - 1;
1193                 float z = endz, dz;
1194                 if(nextsub >= endx) nextsub = endsub = endx-1;
1195                 endz = 1.0f / (w + wslope * nextsub);
1196                 dz = x < nextsub ? (endz - z) / (nextsub - x) : 0.0f;
1197                 for (; x <= endsub; x++, z += dz)
1198                         zf[x] = z;
1199         }
1200 }
1201
1202 void DPSOFTRAST_Draw_Span_Finish(const DPSOFTRAST_State_Draw_Span * RESTRICT span, const float * RESTRICT in4f)
1203 {
1204         int x;
1205         int startx = span->startx;
1206         int endx = span->endx;
1207         int d[4];
1208         float a, b;
1209         unsigned char * RESTRICT pixelmask = span->pixelmask;
1210         unsigned char * RESTRICT pixel = (unsigned char *)dpsoftrast.fb_colorpixels[0];
1211         if (!pixel)
1212                 return;
1213         pixel += span->start * 4;
1214         // handle alphatest now (this affects depth writes too)
1215         if (dpsoftrast.user.alphatest)
1216                 for (x = startx;x < endx;x++)
1217                         if (in4f[x*4+3] < 0.5f)
1218                                 pixelmask[x] = false;
1219         // FIXME: this does not handle bigendian
1220         switch(dpsoftrast.fb_blendmode)
1221         {
1222         case DPSOFTRAST_BLENDMODE_OPAQUE:
1223                 for (x = startx;x < endx;x++)
1224                 {
1225                         if (!pixelmask[x])
1226                                 continue;
1227                         d[0] = (int)(in4f[x*4+2]*255.0f);if (d[0] > 255) d[0] = 255;
1228                         d[1] = (int)(in4f[x*4+1]*255.0f);if (d[1] > 255) d[1] = 255;
1229                         d[2] = (int)(in4f[x*4+0]*255.0f);if (d[2] > 255) d[2] = 255;
1230                         d[3] = (int)(in4f[x*4+3]*255.0f);if (d[3] > 255) d[3] = 255;
1231                         pixel[x*4+0] = d[0];
1232                         pixel[x*4+1] = d[1];
1233                         pixel[x*4+2] = d[2];
1234                         pixel[x*4+3] = d[3];
1235                 }
1236                 break;
1237         case DPSOFTRAST_BLENDMODE_ALPHA:
1238                 for (x = startx;x < endx;x++)
1239                 {
1240                         if (!pixelmask[x])
1241                                 continue;
1242                         a = in4f[x*4+3] * 255.0f;
1243                         b = 1.0f - in4f[x*4+3];
1244                         d[0] = (int)(in4f[x*4+2]*a+pixel[x*4+0]*b);if (d[0] > 255) d[0] = 255;
1245                         d[1] = (int)(in4f[x*4+1]*a+pixel[x*4+1]*b);if (d[1] > 255) d[1] = 255;
1246                         d[2] = (int)(in4f[x*4+0]*a+pixel[x*4+2]*b);if (d[2] > 255) d[2] = 255;
1247                         d[3] = (int)(in4f[x*4+3]*a+pixel[x*4+3]*b);if (d[3] > 255) d[3] = 255;
1248                         pixel[x*4+0] = d[0];
1249                         pixel[x*4+1] = d[1];
1250                         pixel[x*4+2] = d[2];
1251                         pixel[x*4+3] = d[3];
1252                 }
1253                 break;
1254         case DPSOFTRAST_BLENDMODE_ADDALPHA:
1255                 for (x = startx;x < endx;x++)
1256                 {
1257                         if (!pixelmask[x])
1258                                 continue;
1259                         a = in4f[x*4+3] * 255.0f;
1260                         d[0] = (int)(in4f[x*4+2]*a+pixel[x*4+0]);if (d[0] > 255) d[0] = 255;
1261                         d[1] = (int)(in4f[x*4+1]*a+pixel[x*4+1]);if (d[1] > 255) d[1] = 255;
1262                         d[2] = (int)(in4f[x*4+0]*a+pixel[x*4+2]);if (d[2] > 255) d[2] = 255;
1263                         d[3] = (int)(in4f[x*4+3]*a+pixel[x*4+3]);if (d[3] > 255) d[3] = 255;
1264                         pixel[x*4+0] = d[0];
1265                         pixel[x*4+1] = d[1];
1266                         pixel[x*4+2] = d[2];
1267                         pixel[x*4+3] = d[3];
1268                 }
1269                 break;
1270         case DPSOFTRAST_BLENDMODE_ADD:
1271                 for (x = startx;x < endx;x++)
1272                 {
1273                         if (!pixelmask[x])
1274                                 continue;
1275                         d[0] = (int)(in4f[x*4+2]*255.0f+pixel[x*4+0]);if (d[0] > 255) d[0] = 255;
1276                         d[1] = (int)(in4f[x*4+1]*255.0f+pixel[x*4+1]);if (d[1] > 255) d[1] = 255;
1277                         d[2] = (int)(in4f[x*4+0]*255.0f+pixel[x*4+2]);if (d[2] > 255) d[2] = 255;
1278                         d[3] = (int)(in4f[x*4+3]*255.0f+pixel[x*4+3]);if (d[3] > 255) d[3] = 255;
1279                         pixel[x*4+0] = d[0];
1280                         pixel[x*4+1] = d[1];
1281                         pixel[x*4+2] = d[2];
1282                         pixel[x*4+3] = d[3];
1283                 }
1284                 break;
1285         case DPSOFTRAST_BLENDMODE_INVMOD:
1286                 for (x = startx;x < endx;x++)
1287                 {
1288                         if (!pixelmask[x])
1289                                 continue;
1290                         d[0] = (int)((1.0f-in4f[x*4+2])*pixel[x*4+0]);if (d[0] > 255) d[0] = 255;
1291                         d[1] = (int)((1.0f-in4f[x*4+1])*pixel[x*4+1]);if (d[1] > 255) d[1] = 255;
1292                         d[2] = (int)((1.0f-in4f[x*4+0])*pixel[x*4+2]);if (d[2] > 255) d[2] = 255;
1293                         d[3] = (int)((1.0f-in4f[x*4+3])*pixel[x*4+3]);if (d[3] > 255) d[3] = 255;
1294                         pixel[x*4+0] = d[0];
1295                         pixel[x*4+1] = d[1];
1296                         pixel[x*4+2] = d[2];
1297                         pixel[x*4+3] = d[3];
1298                 }
1299                 break;
1300         case DPSOFTRAST_BLENDMODE_MUL:
1301                 for (x = startx;x < endx;x++)
1302                 {
1303                         if (!pixelmask[x])
1304                                 continue;
1305                         d[0] = (int)(in4f[x*4+2]*pixel[x*4+0]);if (d[0] > 255) d[0] = 255;
1306                         d[1] = (int)(in4f[x*4+1]*pixel[x*4+1]);if (d[1] > 255) d[1] = 255;
1307                         d[2] = (int)(in4f[x*4+0]*pixel[x*4+2]);if (d[2] > 255) d[2] = 255;
1308                         d[3] = (int)(in4f[x*4+3]*pixel[x*4+3]);if (d[3] > 255) d[3] = 255;
1309                         pixel[x*4+0] = d[0];
1310                         pixel[x*4+1] = d[1];
1311                         pixel[x*4+2] = d[2];
1312                         pixel[x*4+3] = d[3];
1313                 }
1314                 break;
1315         case DPSOFTRAST_BLENDMODE_MUL2:
1316                 for (x = startx;x < endx;x++)
1317                 {
1318                         if (!pixelmask[x])
1319                                 continue;
1320                         d[0] = (int)(in4f[x*4+2]*pixel[x*4+0]*2.0f);if (d[0] > 255) d[0] = 255;
1321                         d[1] = (int)(in4f[x*4+1]*pixel[x*4+1]*2.0f);if (d[1] > 255) d[1] = 255;
1322                         d[2] = (int)(in4f[x*4+0]*pixel[x*4+2]*2.0f);if (d[2] > 255) d[2] = 255;
1323                         d[3] = (int)(in4f[x*4+3]*pixel[x*4+3]*2.0f);if (d[3] > 255) d[3] = 255;
1324                         pixel[x*4+0] = d[0];
1325                         pixel[x*4+1] = d[1];
1326                         pixel[x*4+2] = d[2];
1327                         pixel[x*4+3] = d[3];
1328                 }
1329                 break;
1330         case DPSOFTRAST_BLENDMODE_SUBALPHA:
1331                 for (x = startx;x < endx;x++)
1332                 {
1333                         if (!pixelmask[x])
1334                                 continue;
1335                         a = in4f[x*4+3] * -255.0f;
1336                         d[0] = (int)(in4f[x*4+2]*a+pixel[x*4+0]);if (d[0] > 255) d[0] = 255;if (d[0] < 0) d[0] = 0;
1337                         d[1] = (int)(in4f[x*4+1]*a+pixel[x*4+1]);if (d[1] > 255) d[1] = 255;if (d[1] < 0) d[1] = 0;
1338                         d[2] = (int)(in4f[x*4+0]*a+pixel[x*4+2]);if (d[2] > 255) d[2] = 255;if (d[2] < 0) d[2] = 0;
1339                         d[3] = (int)(in4f[x*4+3]*a+pixel[x*4+3]);if (d[3] > 255) d[3] = 255;if (d[3] < 0) d[3] = 0;
1340                         pixel[x*4+0] = d[0];
1341                         pixel[x*4+1] = d[1];
1342                         pixel[x*4+2] = d[2];
1343                         pixel[x*4+3] = d[3];
1344                 }
1345                 break;
1346         case DPSOFTRAST_BLENDMODE_PSEUDOALPHA:
1347                 for (x = startx;x < endx;x++)
1348                 {
1349                         if (!pixelmask[x])
1350                                 continue;
1351                         a = 255.0f;
1352                         b = 1.0f - in4f[x*4+3];
1353                         d[0] = (int)(in4f[x*4+2]*a+pixel[x*4+0]*b);if (d[0] > 255) d[0] = 255;
1354                         d[1] = (int)(in4f[x*4+1]*a+pixel[x*4+1]*b);if (d[1] > 255) d[1] = 255;
1355                         d[2] = (int)(in4f[x*4+0]*a+pixel[x*4+2]*b);if (d[2] > 255) d[2] = 255;
1356                         d[3] = (int)(in4f[x*4+3]*a+pixel[x*4+3]*b);if (d[3] > 255) d[3] = 255;
1357                         pixel[x*4+0] = d[0];
1358                         pixel[x*4+1] = d[1];
1359                         pixel[x*4+2] = d[2];
1360                         pixel[x*4+3] = d[3];
1361                 }
1362                 break;
1363         }
1364 }
1365
1366 void DPSOFTRAST_Draw_Span_FinishBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, const unsigned char* RESTRICT in4ub)
1367 {
1368 #ifdef SSE2_PRESENT
1369         int x;
1370         int startx = span->startx;
1371         int endx = span->endx;
1372         const unsigned int * RESTRICT ini = (const unsigned int *)in4ub;
1373         unsigned char * RESTRICT pixelmask = span->pixelmask;
1374         unsigned char * RESTRICT pixel = (unsigned char *)dpsoftrast.fb_colorpixels[0];
1375         unsigned int * RESTRICT pixeli = (unsigned int *)dpsoftrast.fb_colorpixels[0];
1376         if (!pixel)
1377                 return;
1378         pixel += span->start * 4;
1379         pixeli += span->start;
1380         // handle alphatest now (this affects depth writes too)
1381         if (dpsoftrast.user.alphatest)
1382                 for (x = startx;x < endx;x++)
1383                         if (in4ub[x*4+3] < 0.5f)
1384                                 pixelmask[x] = false;
1385         // FIXME: this does not handle bigendian
1386         switch(dpsoftrast.fb_blendmode)
1387         {
1388         case DPSOFTRAST_BLENDMODE_OPAQUE:
1389                 for (x = startx;x + 4 <= endx;)
1390                 {
1391                         if (*(const unsigned int *)&pixelmask[x] == 0x01010101)
1392                         {
1393                                 _mm_storeu_si128((__m128i *)&pixeli[x], _mm_loadu_si128((const __m128i *)&ini[x]));
1394                                 x += 4;
1395                         }
1396                         else
1397                         {
1398                                 if (pixelmask[x])
1399                                         pixeli[x] = ini[x];
1400                                 x++;
1401                         }
1402                 }
1403                 for (;x < endx;x++)
1404                         if (pixelmask[x])
1405                                 pixeli[x] = ini[x];
1406                 break;
1407         case DPSOFTRAST_BLENDMODE_ALPHA:
1408         #define FINISHBLEND(blend2, blend1) \
1409                 for (x = startx;x + 2 <= endx;x += 2) \
1410                 { \
1411                         __m128i src, dst; \
1412                         switch (*(const unsigned short*)&pixelmask[x]) \
1413                         { \
1414                         case 0x0101: \
1415                                 src = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ini[x]), _mm_setzero_si128()); \
1416                                 dst = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixeli[x]), _mm_setzero_si128()); \
1417                                 blend2; \
1418                                 _mm_storel_epi64((__m128i *)&pixeli[x], _mm_packus_epi16(dst, dst)); \
1419                                 continue; \
1420                         case 0x0100: \
1421                                 src = _mm_unpacklo_epi8(_mm_cvtsi32_si128(ini[x+1]), _mm_setzero_si128()); \
1422                                 dst = _mm_unpacklo_epi8(_mm_cvtsi32_si128(pixeli[x+1]), _mm_setzero_si128()); \
1423                                 blend1; \
1424                                 pixeli[x+1] = _mm_cvtsi128_si32(_mm_packus_epi16(dst, dst));  \
1425                                 continue; \
1426                         case 0x0001: \
1427                                 src = _mm_unpacklo_epi8(_mm_cvtsi32_si128(ini[x]), _mm_setzero_si128()); \
1428                                 dst = _mm_unpacklo_epi8(_mm_cvtsi32_si128(pixeli[x]), _mm_setzero_si128()); \
1429                                 blend1; \
1430                                 pixeli[x] = _mm_cvtsi128_si32(_mm_packus_epi16(dst, dst)); \
1431                                 continue; \
1432                         } \
1433                         break; \
1434                 } \
1435                 for(;x < endx; x++) \
1436                 { \
1437                         __m128i src, dst; \
1438                         if (!pixelmask[x]) \
1439                                 continue; \
1440                         src = _mm_unpacklo_epi8(_mm_cvtsi32_si128(ini[x]), _mm_setzero_si128()); \
1441                         dst = _mm_unpacklo_epi8(_mm_cvtsi32_si128(pixeli[x]), _mm_setzero_si128()); \
1442                         blend1; \
1443                         pixeli[x] = _mm_cvtsi128_si32(_mm_packus_epi16(dst, dst)); \
1444                 }
1445
1446                 FINISHBLEND({
1447                         __m128i blend = _mm_shufflehi_epi16(_mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3));
1448                         dst = _mm_add_epi16(dst, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(src, dst), 4), _mm_slli_epi16(blend, 4)));
1449                 }, {
1450                         __m128i blend = _mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3));
1451                         dst = _mm_add_epi16(dst, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(src, dst), 4), _mm_slli_epi16(blend, 4)));
1452                 });
1453                 break;
1454         case DPSOFTRAST_BLENDMODE_ADDALPHA:
1455                 FINISHBLEND({
1456                         __m128i blend = _mm_shufflehi_epi16(_mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3));
1457                         dst = _mm_add_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(src, blend), 8));
1458                 }, {
1459                         __m128i blend = _mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3));
1460                         dst = _mm_add_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(src, blend), 8));
1461                 });
1462                 break;
1463         case DPSOFTRAST_BLENDMODE_ADD:
1464                 FINISHBLEND({ dst = _mm_add_epi16(src, dst); }, { dst = _mm_add_epi16(src, dst); });
1465                 break;
1466         case DPSOFTRAST_BLENDMODE_INVMOD:
1467                 FINISHBLEND({
1468                         dst = _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(dst, src), 8));
1469                 }, {
1470                         dst = _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(dst, src), 8));
1471                 });
1472                 break;
1473         case DPSOFTRAST_BLENDMODE_MUL:
1474                 FINISHBLEND({ dst = _mm_srli_epi16(_mm_mullo_epi16(src, dst), 8); }, { dst = _mm_srli_epi16(_mm_mullo_epi16(src, dst), 8); });
1475                 break;
1476         case DPSOFTRAST_BLENDMODE_MUL2:
1477                 FINISHBLEND({ dst = _mm_srli_epi16(_mm_mullo_epi16(src, dst), 7); }, { dst = _mm_srli_epi16(_mm_mullo_epi16(src, dst), 7); });
1478                 break;
1479         case DPSOFTRAST_BLENDMODE_SUBALPHA:
1480                 FINISHBLEND({
1481                         __m128i blend = _mm_shufflehi_epi16(_mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3));
1482                         dst = _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(src, blend), 8));
1483                 }, {
1484                         __m128i blend = _mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3));
1485                         dst = _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(src, blend), 8));
1486                 });
1487                 break;
1488         case DPSOFTRAST_BLENDMODE_PSEUDOALPHA:
1489                 FINISHBLEND({
1490                         __m128i blend = _mm_shufflehi_epi16(_mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3));
1491                         dst = _mm_add_epi16(src, _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(dst, blend), 8)));
1492                 }, {
1493                         __m128i blend = _mm_shufflelo_epi16(src, _MM_SHUFFLE(3, 3, 3, 3));
1494                         dst = _mm_add_epi16(src, _mm_sub_epi16(dst, _mm_srli_epi16(_mm_mullo_epi16(dst, blend), 8)));
1495                 });
1496                 break;
1497         }
1498 #endif
1499 }
1500
1501 void DPSOFTRAST_Draw_Span_Texture2DVarying(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float * RESTRICT out4f, int texunitindex, int arrayindex, const float * RESTRICT zf)
1502 {
1503         int x;
1504         int startx = span->startx;
1505         int endx = span->endx;
1506         int flags;
1507         float c[4];
1508         float data[4];
1509         float slope[4];
1510         float tc[2], endtc[2];
1511         float tcscale[2];
1512         unsigned int tci[2];
1513         unsigned int tci1[2];
1514         unsigned int tcimin[2];
1515         unsigned int tcimax[2];
1516         int tciwrapmask[2];
1517         int tciwidth;
1518         int filter;
1519         int mip;
1520         const unsigned char * RESTRICT pixelbase;
1521         const unsigned char * RESTRICT pixel[4];
1522         DPSOFTRAST_Texture *texture = dpsoftrast.texbound[texunitindex];
1523         // if no texture is bound, just fill it with white
1524         if (!texture)
1525         {
1526                 for (x = startx;x < endx;x++)
1527                 {
1528                         out4f[x*4+0] = 1.0f;
1529                         out4f[x*4+1] = 1.0f;
1530                         out4f[x*4+2] = 1.0f;
1531                         out4f[x*4+3] = 1.0f;
1532                 }
1533                 return;
1534         }
1535         mip = span->mip[texunitindex];
1536         pixelbase = (unsigned char *)texture->bytes + texture->mipmap[mip][0];
1537         // if this mipmap of the texture is 1 pixel, just fill it with that color
1538         if (texture->mipmap[mip][1] == 4)
1539         {
1540                 c[0] = texture->bytes[2] * (1.0f/255.0f);
1541                 c[1] = texture->bytes[1] * (1.0f/255.0f);
1542                 c[2] = texture->bytes[0] * (1.0f/255.0f);
1543                 c[3] = texture->bytes[3] * (1.0f/255.0f);
1544                 for (x = startx;x < endx;x++)
1545                 {
1546                         out4f[x*4+0] = c[0];
1547                         out4f[x*4+1] = c[1];
1548                         out4f[x*4+2] = c[2];
1549                         out4f[x*4+3] = c[3];
1550                 }
1551                 return;
1552         }
1553         filter = texture->filter & DPSOFTRAST_TEXTURE_FILTER_LINEAR;
1554         data[0] = span->data[0][arrayindex][0];
1555         data[1] = span->data[0][arrayindex][1];
1556         data[2] = span->data[0][arrayindex][2];
1557         data[3] = span->data[0][arrayindex][3];
1558         slope[0] = span->data[1][arrayindex][0];
1559         slope[1] = span->data[1][arrayindex][1];
1560         slope[2] = span->data[1][arrayindex][2];
1561         slope[3] = span->data[1][arrayindex][3];
1562         flags = texture->flags;
1563         tcscale[0] = texture->mipmap[mip][2];
1564         tcscale[1] = texture->mipmap[mip][3];
1565         tciwidth = texture->mipmap[mip][2];
1566         tcimin[0] = 0;
1567         tcimin[1] = 0;
1568         tcimax[0] = texture->mipmap[mip][2]-1;
1569         tcimax[1] = texture->mipmap[mip][3]-1;
1570         tciwrapmask[0] = texture->mipmap[mip][2]-1;
1571         tciwrapmask[1] = texture->mipmap[mip][3]-1;
1572         endtc[0] = (data[0] + slope[0]*startx) * zf[startx] * tcscale[0] - 0.5f;
1573         endtc[1] = (data[1] + slope[1]*startx) * zf[startx] * tcscale[1] - 0.5f;
1574         for (x = startx;x < endx;)
1575         {
1576                 unsigned int subtc[2];
1577                 unsigned int substep[2];
1578                 float subscale = 65536.0f/DPSOFTRAST_MAXSUBSPAN;
1579                 int nextsub = x + DPSOFTRAST_MAXSUBSPAN, endsub = nextsub - 1;
1580                 if(nextsub >= endx)
1581                 {
1582                         nextsub = endsub = endx-1;
1583                         if(x < nextsub) subscale = 65536.0f / (nextsub - x);
1584                 }
1585                 tc[0] = endtc[0];
1586                 tc[1] = endtc[1];
1587                 endtc[0] = (data[0] + slope[0]*nextsub) * zf[nextsub] * tcscale[0] - 0.5f;
1588                 endtc[1] = (data[1] + slope[1]*nextsub) * zf[nextsub] * tcscale[1] - 0.5f;
1589                 substep[0] = (endtc[0] - tc[0]) * subscale;
1590                 substep[1] = (endtc[1] - tc[1]) * subscale;
1591                 subtc[0] = tc[0] * (1<<16);
1592                 subtc[1] = tc[1] * (1<<16);
1593                 if(filter)
1594                 {
1595                         if (flags & DPSOFTRAST_TEXTURE_FLAG_CLAMPTOEDGE)
1596                         {
1597                                 for (; x <= endsub; x++, subtc[0] += substep[0], subtc[1] += substep[1])
1598                                 {
1599                                         unsigned int frac[2] = { subtc[0]&0xFFF, subtc[1]&0xFFF };
1600                                         unsigned int ifrac[2] = { 0x1000 - frac[0], 0x1000 - frac[1] };
1601                                         unsigned int lerp[4] = { ifrac[0]*ifrac[1], frac[0]*ifrac[1], ifrac[0]*frac[1], frac[0]*frac[1] };
1602                                         tci[0] = subtc[0]>>16;
1603                                         tci[1] = subtc[1]>>16;
1604                                         tci1[0] = tci[0] + 1;
1605                                         tci1[1] = tci[1] + 1;
1606                                         tci[0] = tci[0] >= tcimin[0] ? (tci[0] <= tcimax[0] ? tci[0] : tcimax[0]) : tcimin[0];
1607                                         tci[1] = tci[1] >= tcimin[1] ? (tci[1] <= tcimax[1] ? tci[1] : tcimax[1]) : tcimin[1];
1608                                         tci1[0] = tci1[0] >= tcimin[0] ? (tci1[0] <= tcimax[0] ? tci1[0] : tcimax[0]) : tcimin[0];
1609                                         tci1[1] = tci1[1] >= tcimin[1] ? (tci1[1] <= tcimax[1] ? tci1[1] : tcimax[1]) : tcimin[1];
1610                                         pixel[0] = pixelbase + 4 * (tci[1]*tciwidth+tci[0]);
1611                                         pixel[1] = pixelbase + 4 * (tci[1]*tciwidth+tci1[0]);
1612                                         pixel[2] = pixelbase + 4 * (tci1[1]*tciwidth+tci[0]);
1613                                         pixel[3] = pixelbase + 4 * (tci1[1]*tciwidth+tci1[0]);
1614                                         c[0] = (pixel[0][2]*lerp[0]+pixel[1][2]*lerp[1]+pixel[2][2]*lerp[2]+pixel[3][2]*lerp[3]) * (1.0f / 0xFF000000);
1615                                         c[1] = (pixel[0][1]*lerp[0]+pixel[1][1]*lerp[1]+pixel[2][1]*lerp[2]+pixel[3][1]*lerp[3]) * (1.0f / 0xFF000000);
1616                                         c[2] = (pixel[0][0]*lerp[0]+pixel[1][0]*lerp[1]+pixel[2][0]*lerp[2]+pixel[3][0]*lerp[3]) * (1.0f / 0xFF000000);
1617                                         c[3] = (pixel[0][3]*lerp[0]+pixel[1][3]*lerp[1]+pixel[2][3]*lerp[2]+pixel[3][3]*lerp[3]) * (1.0f / 0xFF000000);
1618                                         out4f[x*4+0] = c[0];
1619                                         out4f[x*4+1] = c[1];
1620                                         out4f[x*4+2] = c[2];
1621                                         out4f[x*4+3] = c[3];
1622                                 }
1623                         }
1624                         else
1625                         {
1626                                 for (; x <= endsub; x++, subtc[0] += substep[0], subtc[1] += substep[1])
1627                                 {
1628                                         unsigned int frac[2] = { subtc[0]&0xFFF, subtc[1]&0xFFF };
1629                                         unsigned int ifrac[2] = { 0x1000 - frac[0], 0x1000 - frac[1] };
1630                                         unsigned int lerp[4] = { ifrac[0]*ifrac[1], frac[0]*ifrac[1], ifrac[0]*frac[1], frac[0]*frac[1] };
1631                                         tci[0] = subtc[0]>>16;
1632                                         tci[1] = subtc[1]>>16;
1633                                         tci1[0] = tci[0] + 1;
1634                                         tci1[1] = tci[1] + 1;
1635                                         tci[0] &= tciwrapmask[0];
1636                                         tci[1] &= tciwrapmask[1];
1637                                         tci1[0] &= tciwrapmask[0];
1638                                         tci1[1] &= tciwrapmask[1];
1639                                         pixel[0] = pixelbase + 4 * (tci[1]*tciwidth+tci[0]);
1640                                         pixel[1] = pixelbase + 4 * (tci[1]*tciwidth+tci1[0]);
1641                                         pixel[2] = pixelbase + 4 * (tci1[1]*tciwidth+tci[0]);
1642                                         pixel[3] = pixelbase + 4 * (tci1[1]*tciwidth+tci1[0]);
1643                                         c[0] = (pixel[0][2]*lerp[0]+pixel[1][2]*lerp[1]+pixel[2][2]*lerp[2]+pixel[3][2]*lerp[3]) * (1.0f / 0xFF000000);
1644                                         c[1] = (pixel[0][1]*lerp[0]+pixel[1][1]*lerp[1]+pixel[2][1]*lerp[2]+pixel[3][1]*lerp[3]) * (1.0f / 0xFF000000);
1645                                         c[2] = (pixel[0][0]*lerp[0]+pixel[1][0]*lerp[1]+pixel[2][0]*lerp[2]+pixel[3][0]*lerp[3]) * (1.0f / 0xFF000000);
1646                                         c[3] = (pixel[0][3]*lerp[0]+pixel[1][3]*lerp[1]+pixel[2][3]*lerp[2]+pixel[3][3]*lerp[3]) * (1.0f / 0xFF000000);
1647                                         out4f[x*4+0] = c[0];
1648                                         out4f[x*4+1] = c[1];
1649                                         out4f[x*4+2] = c[2];
1650                                         out4f[x*4+3] = c[3];
1651                                 }
1652                         }
1653                 }
1654                 else if (flags & DPSOFTRAST_TEXTURE_FLAG_CLAMPTOEDGE)
1655                 {
1656                         for (; x <= endsub; x++, subtc[0] += substep[0], subtc[1] += substep[1])
1657                         {
1658                                 tci[0] = subtc[0]>>16;
1659                                 tci[1] = subtc[1]>>16;
1660                                 tci[0] = tci[0] >= tcimin[0] ? (tci[0] <= tcimax[0] ? tci[0] : tcimax[0]) : tcimin[0];
1661                                 tci[1] = tci[1] >= tcimin[1] ? (tci[1] <= tcimax[1] ? tci[1] : tcimax[1]) : tcimin[1];
1662                                 pixel[0] = pixelbase + 4 * (tci[1]*tciwidth+tci[0]);
1663                                 c[0] = pixel[0][2] * (1.0f / 255.0f);
1664                                 c[1] = pixel[0][1] * (1.0f / 255.0f);
1665                                 c[2] = pixel[0][0] * (1.0f / 255.0f);
1666                                 c[3] = pixel[0][3] * (1.0f / 255.0f);
1667                                 out4f[x*4+0] = c[0];
1668                                 out4f[x*4+1] = c[1];
1669                                 out4f[x*4+2] = c[2];
1670                                 out4f[x*4+3] = c[3];
1671                         }
1672                 }
1673                 else
1674                 {
1675                         for (; x <= endsub; x++, subtc[0] += substep[0], subtc[1] += substep[1])
1676                         {
1677                                 tci[0] = subtc[0]>>16;
1678                                 tci[1] = subtc[1]>>16;
1679                                 tci[0] &= tciwrapmask[0];
1680                                 tci[1] &= tciwrapmask[1];
1681                                 pixel[0] = pixelbase + 4 * (tci[1]*tciwidth+tci[0]);
1682                                 c[0] = pixel[0][2] * (1.0f / 255.0f);
1683                                 c[1] = pixel[0][1] * (1.0f / 255.0f);
1684                                 c[2] = pixel[0][0] * (1.0f / 255.0f);
1685                                 c[3] = pixel[0][3] * (1.0f / 255.0f);
1686                                 out4f[x*4+0] = c[0];
1687                                 out4f[x*4+1] = c[1];
1688                                 out4f[x*4+2] = c[2];
1689                                 out4f[x*4+3] = c[3];
1690                         }
1691                 }
1692         }
1693 }
1694
1695 void DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char * RESTRICT out4ub, int texunitindex, int arrayindex, const float * RESTRICT zf)
1696 {
1697 #ifdef SSE2_PRESENT
1698         int x;
1699         int startx = span->startx;
1700         int endx = span->endx;
1701         int flags;
1702         __m128 data, slope, tcscale;
1703         __m128i tcsize, tcmask, tcoffset, tcmax;
1704         __m128 tc, endtc;
1705         __m128i subtc, substep, endsubtc;
1706         int filter;
1707         int mip;
1708         unsigned int *outi = (unsigned int *)out4ub;
1709         const unsigned char * RESTRICT pixelbase;
1710         DPSOFTRAST_Texture *texture = dpsoftrast.texbound[texunitindex];
1711         // if no texture is bound, just fill it with white
1712         if (!texture)
1713         {
1714                 memset(out4ub + startx*4, 255, span->length*4);
1715                 return;
1716         }
1717         mip = span->mip[texunitindex];
1718         pixelbase = (const unsigned char *)texture->bytes + texture->mipmap[mip][0];
1719         // if this mipmap of the texture is 1 pixel, just fill it with that color
1720         if (texture->mipmap[mip][1] == 4)
1721         {
1722                 unsigned int k = *((const unsigned int *)pixelbase);
1723                 for (x = startx;x < endx;x++)
1724                         outi[x] = k;
1725                 return;
1726         }
1727         filter = texture->filter & DPSOFTRAST_TEXTURE_FILTER_LINEAR;
1728         data = _mm_load_ps(span->data[0][arrayindex]);
1729         slope = _mm_load_ps(span->data[1][arrayindex]);
1730         flags = texture->flags;
1731         tcsize = _mm_shuffle_epi32(_mm_loadu_si128((const __m128i *)&texture->mipmap[mip][0]), _MM_SHUFFLE(3, 2, 3, 2));
1732         tcmask = _mm_sub_epi32(tcsize, _mm_set1_epi32(1));
1733         tcscale = _mm_cvtepi32_ps(tcsize);
1734         data = _mm_mul_ps(_mm_shuffle_ps(data, data, _MM_SHUFFLE(1, 0, 1, 0)), tcscale);
1735         slope = _mm_mul_ps(_mm_shuffle_ps(slope, slope, _MM_SHUFFLE(1, 0, 1, 0)), tcscale);
1736         endtc = _mm_sub_ps(_mm_mul_ps(_mm_add_ps(data, _mm_mul_ps(slope, _mm_set1_ps(startx))), _mm_set1_ps(zf[startx])), _mm_set1_ps(0.5f));
1737         endsubtc = _mm_cvtps_epi32(_mm_mul_ps(endtc, _mm_set1_ps(65536.0f)));
1738         tcoffset = _mm_add_epi32(_mm_slli_epi32(_mm_shuffle_epi32(tcsize, _MM_SHUFFLE(0, 0, 0, 0)), 18), _mm_set1_epi32(4));
1739         tcmax = filter ? _mm_packs_epi32(tcmask, tcmask) : _mm_slli_epi32(tcmask, 16);
1740         for (x = startx;x < endx;)
1741         {
1742                 int nextsub = x + DPSOFTRAST_MAXSUBSPAN, endsub = nextsub - 1;
1743                 __m128 subscale = _mm_set1_ps(65536.0f/DPSOFTRAST_MAXSUBSPAN);
1744                 if(nextsub >= endx)
1745                 {
1746                         nextsub = endsub = endx-1;
1747                         if(x < nextsub) subscale = _mm_set1_ps(65536.0f / (nextsub - x));
1748                 }
1749                 tc = endtc;
1750                 subtc = endsubtc;
1751                 endtc = _mm_sub_ps(_mm_mul_ps(_mm_add_ps(data, _mm_mul_ps(slope, _mm_set1_ps(nextsub))), _mm_set1_ps(zf[nextsub])), _mm_set1_ps(0.5f));
1752                 substep = _mm_cvtps_epi32(_mm_mul_ps(_mm_sub_ps(endtc, tc), subscale));
1753                 endsubtc = _mm_cvtps_epi32(_mm_mul_ps(endtc, _mm_set1_ps(65536.0f)));
1754                 if (filter)
1755                 {
1756                         __m128i tcrange = _mm_srai_epi32(_mm_unpacklo_epi64(subtc, _mm_add_epi32(endsubtc, substep)), 16);
1757                         if (_mm_movemask_epi8(_mm_andnot_si128(_mm_cmplt_epi32(tcrange, _mm_setzero_si128()), _mm_cmplt_epi32(tcrange, tcmask))) == 0xFFFF)
1758                         {
1759                                 subtc = _mm_unpacklo_epi64(subtc, _mm_add_epi32(subtc, substep));
1760                                 substep = _mm_slli_epi32(substep, 1);
1761                                 for (; x + 1 <= endsub; x += 2, subtc = _mm_add_epi32(subtc, substep))
1762                                 {
1763                                         __m128i tci = _mm_shufflehi_epi16(_mm_shufflelo_epi16(subtc, _MM_SHUFFLE(3, 1, 3, 1)), _MM_SHUFFLE(3, 1, 3, 1)), pix1, pix2, pix3, pix4, fracm;
1764                                         tci = _mm_madd_epi16(_mm_add_epi16(tci, _mm_setr_epi32(0, 0x10000, 0, 0x10000)), tcoffset);
1765                                         pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(tci)]), _mm_setzero_si128());
1766                                         pix2 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))]), _mm_setzero_si128());
1767                                         pix3 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(2, 2, 2, 2)))]), _mm_setzero_si128());
1768                                         pix4 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(3, 3, 3, 3)))]), _mm_setzero_si128());
1769                                         fracm = _mm_srli_epi16(subtc, 1);
1770                                         pix1 = _mm_add_epi16(pix1,
1771                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1772                                                                                                                  _mm_shuffle_epi32(_mm_shufflelo_epi16(fracm, _MM_SHUFFLE(2, 2, 2, 2)), _MM_SHUFFLE(1, 0, 1, 0))));
1773                                         pix3 = _mm_add_epi16(pix3,
1774                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix4, pix3), 1),
1775                                                                                                                  _mm_shuffle_epi32(_mm_shufflehi_epi16(fracm, _MM_SHUFFLE(2, 2, 2, 2)), _MM_SHUFFLE(3, 2, 3, 2))));
1776                                         pix2 = _mm_unpacklo_epi64(pix1, pix3);
1777                                         pix4 = _mm_unpackhi_epi64(pix1, pix3);
1778                                         pix2 = _mm_add_epi16(pix2,
1779                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix4, pix2), 1),
1780                                                                                                                  _mm_shufflehi_epi16(_mm_shufflelo_epi16(fracm, _MM_SHUFFLE(0, 0, 0, 0)), _MM_SHUFFLE(0, 0, 0, 0))));
1781                                         _mm_storel_epi64((__m128i *)&outi[x], _mm_packus_epi16(pix2, _mm_shufflelo_epi16(pix2, _MM_SHUFFLE(3, 2, 3, 2))));
1782                                 }
1783                                 if (x <= endsub)
1784                                 {
1785                                         __m128i tci = _mm_shufflelo_epi16(subtc, _MM_SHUFFLE(3, 1, 3, 1)), pix1, pix2, fracm;
1786                                         tci = _mm_madd_epi16(_mm_add_epi16(tci, _mm_setr_epi32(0, 0x10000, 0, 0)), tcoffset);
1787                                         pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(tci)]), _mm_setzero_si128());
1788                                         pix2 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))]), _mm_setzero_si128());
1789                                         fracm = _mm_srli_epi16(subtc, 1);
1790                                         pix1 = _mm_add_epi16(pix1,
1791                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1792                                                                                                                  _mm_shuffle_epi32(_mm_shufflelo_epi16(fracm, _MM_SHUFFLE(2, 2, 2, 2)), _MM_SHUFFLE(1, 0, 1, 0))));
1793                                         pix2 = _mm_shuffle_epi32(pix1, _MM_SHUFFLE(3, 2, 3, 2));
1794                                         pix1 = _mm_add_epi16(pix1,
1795                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1796                                                                                                                  _mm_shufflelo_epi16(fracm, _MM_SHUFFLE(0, 0, 0, 0))));
1797                                         outi[x] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
1798                                         x++;
1799                                 }
1800                         }
1801                         else if (flags & DPSOFTRAST_TEXTURE_FLAG_CLAMPTOEDGE)
1802                         {
1803                                 for (; x <= endsub; x++, subtc = _mm_add_epi32(subtc, substep))
1804                                 {
1805                                         __m128i tci = _mm_shufflehi_epi16(_mm_shufflelo_epi16(subtc, _MM_SHUFFLE(3, 1, 3, 1)), _MM_SHUFFLE(3, 1, 3, 1)), pix1, pix2, fracm;
1806                                         tci = _mm_min_epi16(_mm_max_epi16(_mm_add_epi16(tci, _mm_setr_epi32(0, 1, 0x10000, 0x10001)), _mm_setzero_si128()), tcmax);
1807                                         tci = _mm_madd_epi16(tci, tcoffset);
1808                                         pix1 = _mm_unpacklo_epi8(_mm_unpacklo_epi32(_mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(tci)]),
1809                                                                                                                                 _mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))])),
1810                                                                                         _mm_setzero_si128());
1811                                         pix2 = _mm_unpacklo_epi8(_mm_unpacklo_epi32(_mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(2, 2, 2, 2)))]),
1812                                                                                                                                 _mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(3, 3, 3, 3)))])),
1813                                                                                         _mm_setzero_si128());
1814                                         fracm = _mm_srli_epi16(subtc, 1);
1815                                         pix1 = _mm_add_epi16(pix1,
1816                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1817                                                                                                                  _mm_shuffle_epi32(_mm_shufflelo_epi16(fracm, _MM_SHUFFLE(2, 2, 2, 2)), _MM_SHUFFLE(1, 0, 1, 0))));
1818                                         pix2 = _mm_shuffle_epi32(pix1, _MM_SHUFFLE(3, 2, 3, 2));
1819                                         pix1 = _mm_add_epi16(pix1,
1820                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1821                                                                                                                  _mm_shufflelo_epi16(fracm, _MM_SHUFFLE(0, 0, 0, 0))));
1822                                         outi[x] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
1823                                 }
1824                         }
1825                         else
1826                         {
1827                                 for (; x <= endsub; x++, subtc = _mm_add_epi32(subtc, substep))
1828                                 {
1829                                         __m128i tci = _mm_shufflehi_epi16(_mm_shufflelo_epi16(subtc, _MM_SHUFFLE(3, 1, 3, 1)), _MM_SHUFFLE(3, 1, 3, 1)), pix1, pix2, fracm;
1830                                         tci = _mm_and_si128(_mm_add_epi16(tci, _mm_setr_epi32(0, 1, 0x10000, 0x10001)), tcmax);
1831                                         tci = _mm_madd_epi16(tci, tcoffset);
1832                                         pix1 = _mm_unpacklo_epi8(_mm_unpacklo_epi32(_mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(tci)]),
1833                                                                                                                                 _mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))])),
1834                                                                                         _mm_setzero_si128());
1835                                         pix2 = _mm_unpacklo_epi8(_mm_unpacklo_epi32(_mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(2, 2, 2, 2)))]),
1836                                                                                                                                 _mm_cvtsi32_si128(*(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(3, 3, 3, 3)))])),
1837                                                                                         _mm_setzero_si128());
1838                                         fracm = _mm_srli_epi16(subtc, 1);
1839                                         pix1 = _mm_add_epi16(pix1,
1840                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1841                                                                                                                  _mm_shuffle_epi32(_mm_shufflelo_epi16(fracm, _MM_SHUFFLE(2, 2, 2, 2)), _MM_SHUFFLE(1, 0, 1, 0))));
1842                                         pix2 = _mm_shuffle_epi32(pix1, _MM_SHUFFLE(3, 2, 3, 2));
1843                                         pix1 = _mm_add_epi16(pix1,
1844                                                                                  _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 1),
1845                                                                                                                  _mm_shufflelo_epi16(fracm, _MM_SHUFFLE(0, 0, 0, 0))));
1846                                         outi[x] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
1847                                 }
1848                         }
1849                 }
1850                 else
1851                 {
1852                         subtc = _mm_unpacklo_epi64(subtc, _mm_add_epi32(subtc, substep));
1853                         substep = _mm_slli_epi32(substep, 1);
1854                         if (flags & DPSOFTRAST_TEXTURE_FLAG_CLAMPTOEDGE)
1855                         {
1856                                 for (; x + 1 <= endsub; x += 2, subtc = _mm_add_epi32(subtc, substep))
1857                                 {
1858                                         __m128i tci = _mm_min_epi16(_mm_max_epi16(subtc, _mm_setzero_si128()), tcmax);
1859                                         tci = _mm_shufflehi_epi16(_mm_shufflelo_epi16(tci, _MM_SHUFFLE(3, 1, 3, 1)), _MM_SHUFFLE(3, 1, 3, 1));
1860                                         tci = _mm_madd_epi16(tci, tcoffset);
1861                                         outi[x] = *(const int *)&pixelbase[_mm_cvtsi128_si32(tci)];
1862                                         outi[x+1] = *(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))];
1863                                 }
1864                                 if (x <= endsub)
1865                                 {
1866                                         __m128i tci = _mm_min_epi16(_mm_max_epi16(subtc, _mm_setzero_si128()), tcmax);
1867                                         tci = _mm_shufflelo_epi16(tci, _MM_SHUFFLE(3, 1, 3, 1));
1868                                         tci = _mm_madd_epi16(tci, tcoffset);
1869                                         outi[x] = *(const int *)&pixelbase[_mm_cvtsi128_si32(tci)];
1870                                         x++;
1871                                 }
1872                         }
1873                         else
1874                         {
1875                                 for (; x + 1 <= endsub; x += 2, subtc = _mm_add_epi32(subtc, substep))
1876                                 {
1877                                         __m128i tci = _mm_and_si128(subtc, tcmax);
1878                                         tci = _mm_shufflehi_epi16(_mm_shufflelo_epi16(tci, _MM_SHUFFLE(3, 1, 3, 1)), _MM_SHUFFLE(3, 1, 3, 1));
1879                                         tci = _mm_madd_epi16(tci, tcoffset);
1880                                         outi[x] = *(const int *)&pixelbase[_mm_cvtsi128_si32(tci)];
1881                                         outi[x+1] = *(const int *)&pixelbase[_mm_cvtsi128_si32(_mm_shuffle_epi32(tci, _MM_SHUFFLE(1, 1, 1, 1)))];
1882                                 }
1883                                 if (x <= endsub)
1884                                 {
1885                                         __m128i tci = _mm_and_si128(subtc, tcmax);
1886                                         tci = _mm_shufflelo_epi16(tci, _MM_SHUFFLE(3, 1, 3, 1));
1887                                         tci = _mm_madd_epi16(tci, tcoffset);
1888                                         outi[x] = *(const int *)&pixelbase[_mm_cvtsi128_si32(tci)];
1889                                         x++;
1890                                 }
1891                         }
1892                 }
1893         }
1894 #endif
1895 }
1896
1897 void DPSOFTRAST_Draw_Span_TextureCubeVaryingBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char * RESTRICT out4ub, int texunitindex, int arrayindex, const float * RESTRICT zf)
1898 {
1899         // TODO: IMPLEMENT
1900         memset(out4ub, 255, span->length*4);
1901 }
1902
1903 float DPSOFTRAST_SampleShadowmap(const float *vector)
1904 {
1905         // TODO: IMPLEMENT
1906         return 1.0f;
1907 }
1908
1909 void DPSOFTRAST_Draw_Span_MultiplyVarying(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *in4f, int arrayindex, const float *zf)
1910 {
1911         int x;
1912         int startx = span->startx;
1913         int endx = span->endx;
1914         float c[4];
1915         float data[4];
1916         float slope[4];
1917         float z;
1918         data[0] = span->data[0][arrayindex][0];
1919         data[1] = span->data[0][arrayindex][1];
1920         data[2] = span->data[0][arrayindex][2];
1921         data[3] = span->data[0][arrayindex][3];
1922         slope[0] = span->data[1][arrayindex][0];
1923         slope[1] = span->data[1][arrayindex][1];
1924         slope[2] = span->data[1][arrayindex][2];
1925         slope[3] = span->data[1][arrayindex][3];
1926         for (x = startx;x < endx;x++)
1927         {
1928                 z = zf[x];
1929                 c[0] = (data[0] + slope[0]*x) * z;
1930                 c[1] = (data[1] + slope[1]*x) * z;
1931                 c[2] = (data[2] + slope[2]*x) * z;
1932                 c[3] = (data[3] + slope[3]*x) * z;
1933                 out4f[x*4+0] = in4f[x*4+0] * c[0];
1934                 out4f[x*4+1] = in4f[x*4+1] * c[1];
1935                 out4f[x*4+2] = in4f[x*4+2] * c[2];
1936                 out4f[x*4+3] = in4f[x*4+3] * c[3];
1937         }
1938 }
1939
1940 void DPSOFTRAST_Draw_Span_Varying(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, int arrayindex, const float *zf)
1941 {
1942         int x;
1943         int startx = span->startx;
1944         int endx = span->endx;
1945         float c[4];
1946         float data[4];
1947         float slope[4];
1948         float z;
1949         data[0] = span->data[0][arrayindex][0];
1950         data[1] = span->data[0][arrayindex][1];
1951         data[2] = span->data[0][arrayindex][2];
1952         data[3] = span->data[0][arrayindex][3];
1953         slope[0] = span->data[1][arrayindex][0];
1954         slope[1] = span->data[1][arrayindex][1];
1955         slope[2] = span->data[1][arrayindex][2];
1956         slope[3] = span->data[1][arrayindex][3];
1957         for (x = startx;x < endx;x++)
1958         {
1959                 z = zf[x];
1960                 c[0] = (data[0] + slope[0]*x) * z;
1961                 c[1] = (data[1] + slope[1]*x) * z;
1962                 c[2] = (data[2] + slope[2]*x) * z;
1963                 c[3] = (data[3] + slope[3]*x) * z;
1964                 out4f[x*4+0] = c[0];
1965                 out4f[x*4+1] = c[1];
1966                 out4f[x*4+2] = c[2];
1967                 out4f[x*4+3] = c[3];
1968         }
1969 }
1970
1971 void DPSOFTRAST_Draw_Span_AddBloom(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *ina4f, const float *inb4f, const float *subcolor)
1972 {
1973         int x, startx = span->startx, endx = span->endx;
1974         float c[4], localcolor[4];
1975         localcolor[0] = subcolor[0];
1976         localcolor[1] = subcolor[1];
1977         localcolor[2] = subcolor[2];
1978         localcolor[3] = subcolor[3];
1979         for (x = startx;x < endx;x++)
1980         {
1981                 c[0] = inb4f[x*4+0] - localcolor[0];if (c[0] < 0.0f) c[0] = 0.0f;
1982                 c[1] = inb4f[x*4+1] - localcolor[1];if (c[1] < 0.0f) c[1] = 0.0f;
1983                 c[2] = inb4f[x*4+2] - localcolor[2];if (c[2] < 0.0f) c[2] = 0.0f;
1984                 c[3] = inb4f[x*4+3] - localcolor[3];if (c[3] < 0.0f) c[3] = 0.0f;
1985                 out4f[x*4+0] = ina4f[x*4+0] + c[0];
1986                 out4f[x*4+1] = ina4f[x*4+1] + c[1];
1987                 out4f[x*4+2] = ina4f[x*4+2] + c[2];
1988                 out4f[x*4+3] = ina4f[x*4+3] + c[3];
1989         }
1990 }
1991
1992 void DPSOFTRAST_Draw_Span_MultiplyBuffers(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *ina4f, const float *inb4f)
1993 {
1994         int x, startx = span->startx, endx = span->endx;
1995         for (x = startx;x < endx;x++)
1996         {
1997                 out4f[x*4+0] = ina4f[x*4+0] * inb4f[x*4+0];
1998                 out4f[x*4+1] = ina4f[x*4+1] * inb4f[x*4+1];
1999                 out4f[x*4+2] = ina4f[x*4+2] * inb4f[x*4+2];
2000                 out4f[x*4+3] = ina4f[x*4+3] * inb4f[x*4+3];
2001         }
2002 }
2003
2004 void DPSOFTRAST_Draw_Span_AddBuffers(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *ina4f, const float *inb4f)
2005 {
2006         int x, startx = span->startx, endx = span->endx;
2007         for (x = startx;x < endx;x++)
2008         {
2009                 out4f[x*4+0] = ina4f[x*4+0] + inb4f[x*4+0];
2010                 out4f[x*4+1] = ina4f[x*4+1] + inb4f[x*4+1];
2011                 out4f[x*4+2] = ina4f[x*4+2] + inb4f[x*4+2];
2012                 out4f[x*4+3] = ina4f[x*4+3] + inb4f[x*4+3];
2013         }
2014 }
2015
2016 void DPSOFTRAST_Draw_Span_MixBuffers(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *ina4f, const float *inb4f)
2017 {
2018         int x, startx = span->startx, endx = span->endx;
2019         float a, b;
2020         for (x = startx;x < endx;x++)
2021         {
2022                 a = 1.0f - inb4f[x*4+3];
2023                 b = inb4f[x*4+3];
2024                 out4f[x*4+0] = ina4f[x*4+0] * a + inb4f[x*4+0] * b;
2025                 out4f[x*4+1] = ina4f[x*4+1] * a + inb4f[x*4+1] * b;
2026                 out4f[x*4+2] = ina4f[x*4+2] * a + inb4f[x*4+2] * b;
2027                 out4f[x*4+3] = ina4f[x*4+3] * a + inb4f[x*4+3] * b;
2028         }
2029 }
2030
2031 void DPSOFTRAST_Draw_Span_MixUniformColor(const DPSOFTRAST_State_Draw_Span * RESTRICT span, float *out4f, const float *in4f, const float *color)
2032 {
2033         int x, startx = span->startx, endx = span->endx;
2034         float localcolor[4], ilerp, lerp;
2035         localcolor[0] = color[0];
2036         localcolor[1] = color[1];
2037         localcolor[2] = color[2];
2038         localcolor[3] = color[3];
2039         ilerp = 1.0f - localcolor[3];
2040         lerp = localcolor[3];
2041         for (x = startx;x < endx;x++)
2042         {
2043                 out4f[x*4+0] = in4f[x*4+0] * ilerp + localcolor[0] * lerp;
2044                 out4f[x*4+1] = in4f[x*4+1] * ilerp + localcolor[1] * lerp;
2045                 out4f[x*4+2] = in4f[x*4+2] * ilerp + localcolor[2] * lerp;
2046                 out4f[x*4+3] = in4f[x*4+3] * ilerp + localcolor[3] * lerp;
2047         }
2048 }
2049
2050
2051
2052 void DPSOFTRAST_Draw_Span_MultiplyVaryingBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *in4ub, int arrayindex, const float *zf)
2053 {
2054 #ifdef SSE2_PRESENT
2055         int x;
2056         int startx = span->startx;
2057         int endx = span->endx;
2058         __m128 data = _mm_load_ps(span->data[0][arrayindex]), slope = _mm_load_ps(span->data[1][arrayindex]);
2059         data = _mm_shuffle_ps(data, data, _MM_SHUFFLE(3, 0, 1, 2));
2060         slope = _mm_shuffle_ps(slope, slope, _MM_SHUFFLE(3, 0, 1, 2));
2061         data = _mm_add_ps(data, _mm_mul_ps(slope, _mm_set1_ps(startx)));
2062         data = _mm_mul_ps(data, _mm_set1_ps(256.0f));
2063         slope = _mm_mul_ps(slope, _mm_set1_ps(256.0f));
2064         for (x = startx;x+2 <= endx;x += 2, data = _mm_add_ps(data, slope))
2065         {
2066                 __m128i pix = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_loadl_epi64((const __m128i *)&in4ub[x*4]));
2067                 __m128i mod = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x]))), mod2;
2068                 data = _mm_add_ps(data, slope);
2069                 mod2 = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x+1])));
2070                 mod = _mm_unpacklo_epi64(_mm_packs_epi32(mod, mod), _mm_packs_epi32(mod2, mod2));
2071                 pix = _mm_mulhi_epu16(pix, mod);
2072                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix, pix));
2073         }
2074         for (;x < endx;x++, data = _mm_add_ps(data, slope))
2075         {
2076                 __m128i pix = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&in4ub[x*4]));
2077                 __m128i mod = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x])));
2078                 mod = _mm_packs_epi32(mod, mod);
2079                 pix = _mm_mulhi_epu16(pix, mod);
2080                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix, pix));
2081         }
2082 #endif
2083 }
2084
2085 void DPSOFTRAST_Draw_Span_VaryingBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, int arrayindex, const float *zf)
2086 {
2087 #ifdef SSE2_PRESENT
2088         int x;
2089         int startx = span->startx;
2090         int endx = span->endx;
2091         __m128 data = _mm_load_ps(span->data[0][arrayindex]), slope = _mm_load_ps(span->data[1][arrayindex]);
2092         data = _mm_shuffle_ps(data, data, _MM_SHUFFLE(3, 0, 1, 2));
2093         slope = _mm_shuffle_ps(slope, slope, _MM_SHUFFLE(3, 0, 1, 2));
2094         data = _mm_add_ps(data, _mm_mul_ps(slope, _mm_set1_ps(startx)));
2095         data = _mm_mul_ps(data, _mm_set1_ps(255.0f));
2096         slope = _mm_mul_ps(slope, _mm_set1_ps(255.0f));
2097         for (x = startx;x+2 <= endx;x += 2, data = _mm_add_ps(data, slope))
2098         {
2099                 __m128i pix = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x]))), pix2;
2100                 data = _mm_add_ps(data, slope);
2101                 pix2 = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x+1])));
2102                 pix = _mm_unpacklo_epi64(_mm_packs_epi32(pix, pix), _mm_packs_epi32(pix2, pix2));
2103                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix, pix));
2104         }
2105         for (;x < endx;x++, data = _mm_add_ps(data, slope))
2106         {
2107                 __m128i pix = _mm_cvtps_epi32(_mm_mul_ps(data, _mm_load1_ps(&zf[x])));
2108                 pix = _mm_packs_epi32(pix, pix);
2109                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix, pix));
2110         }
2111 #endif
2112 }
2113
2114 void DPSOFTRAST_Draw_Span_AddBloomBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *ina4ub, const unsigned char *inb4ub, const float *subcolor)
2115 {
2116 #ifdef SSE2_PRESENT
2117         int x, startx = span->startx, endx = span->endx;
2118         __m128i localcolor = _mm_shuffle_epi32(_mm_cvtps_epi32(_mm_mul_ps(_mm_loadu_ps(subcolor), _mm_set1_ps(255.0f))), _MM_SHUFFLE(3, 0, 1, 2));
2119         localcolor = _mm_shuffle_epi32(_mm_packs_epi32(localcolor, localcolor), _MM_SHUFFLE(1, 0, 1, 0));
2120         for (x = startx;x+2 <= endx;x+=2)
2121         {
2122                 __m128i pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ina4ub[x*4]), _mm_setzero_si128());
2123                 __m128i pix2 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&inb4ub[x*4]), _mm_setzero_si128());
2124                 pix1 = _mm_add_epi16(pix1, _mm_sub_epi16(pix2, localcolor));
2125                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix1, pix1));
2126         }
2127         if(x < endx)
2128         {
2129                 __m128i pix1 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&ina4ub[x*4]), _mm_setzero_si128());
2130                 __m128i pix2 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&inb4ub[x*4]), _mm_setzero_si128());
2131                 pix1 = _mm_add_epi16(pix1, _mm_sub_epi16(pix2, localcolor));
2132                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
2133         }
2134 #endif
2135 }
2136
2137 void DPSOFTRAST_Draw_Span_MultiplyBuffersBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *ina4ub, const unsigned char *inb4ub)
2138 {
2139 #ifdef SSE2_PRESENT
2140         int x, startx = span->startx, endx = span->endx;
2141         for (x = startx;x+2 <= endx;x+=2)
2142         {
2143                 __m128i pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ina4ub[x*4]), _mm_setzero_si128());
2144                 __m128i pix2 = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_loadl_epi64((const __m128i *)&inb4ub[x*4]));
2145                 pix1 = _mm_mulhi_epu16(pix1, pix2);
2146                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix1, pix1));
2147         }
2148         if(x < endx)
2149         {
2150                 __m128i pix1 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&ina4ub[x*4]), _mm_setzero_si128());
2151                 __m128i pix2 = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&inb4ub[x*4]));
2152                 pix1 = _mm_mulhi_epu16(pix1, pix2);
2153                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
2154         }
2155 #endif
2156 }
2157
2158 void DPSOFTRAST_Draw_Span_AddBuffersBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *ina4ub, const unsigned char *inb4ub)
2159 {
2160 #ifdef SSE2_PRESENT
2161         int x, startx = span->startx, endx = span->endx;
2162         for (x = startx;x+2 <= endx;x+=2)
2163         {
2164                 __m128i pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ina4ub[x*4]), _mm_setzero_si128());
2165                 __m128i pix2 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&inb4ub[x*4]), _mm_setzero_si128());
2166                 pix1 = _mm_add_epi16(pix1, pix2);
2167                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix1, pix1));
2168         }
2169         if(x < endx)
2170         {
2171                 __m128i pix1 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&ina4ub[x*4]), _mm_setzero_si128());
2172                 __m128i pix2 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&inb4ub[x*4]), _mm_setzero_si128());
2173                 pix1 = _mm_add_epi16(pix1, pix2);
2174                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
2175         }
2176 #endif
2177 }
2178
2179 void DPSOFTRAST_Draw_Span_TintedAddBuffersBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *ina4ub, const unsigned char *inb4ub, const float *inbtintbgra)
2180 {
2181 #ifdef SSE2_PRESENT
2182         int x, startx = span->startx, endx = span->endx;
2183         __m128i tint = _mm_cvtps_epi32(_mm_mul_ps(_mm_loadu_ps(inbtintbgra), _mm_set1_ps(256.0f)));
2184         tint = _mm_shuffle_epi32(_mm_packs_epi32(tint, tint), _MM_SHUFFLE(1, 0, 1, 0));
2185         for (x = startx;x+2 <= endx;x+=2)
2186         {
2187                 __m128i pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ina4ub[x*4]), _mm_setzero_si128());
2188                 __m128i pix2 = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_loadl_epi64((const __m128i *)&inb4ub[x*4]));
2189                 pix1 = _mm_add_epi16(pix1, _mm_mulhi_epu16(tint, pix2));
2190                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix1, pix1));
2191         }
2192         if(x < endx)
2193         {
2194                 __m128i pix1 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&ina4ub[x*4]), _mm_setzero_si128());
2195                 __m128i pix2 = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&inb4ub[x*4]));
2196                 pix1 = _mm_add_epi16(pix1, _mm_mulhi_epu16(tint, pix2));
2197                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
2198         }
2199 #endif
2200 }
2201
2202 void DPSOFTRAST_Draw_Span_MixBuffersBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *ina4ub, const unsigned char *inb4ub)
2203 {
2204 #ifdef SSE2_PRESENT
2205         int x, startx = span->startx, endx = span->endx;
2206         for (x = startx;x+2 <= endx;x+=2)
2207         {
2208                 __m128i pix1 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&ina4ub[x*4]), _mm_setzero_si128());
2209                 __m128i pix2 = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&inb4ub[x*4]), _mm_setzero_si128());
2210                 __m128i blend = _mm_shufflehi_epi16(_mm_shufflelo_epi16(pix2, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3));
2211                 pix1 = _mm_add_epi16(pix1, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 4), _mm_slli_epi16(blend, 4)));
2212                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix1, pix1));
2213         }
2214         if(x < endx)
2215         {
2216                 __m128i pix1 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&ina4ub[x*4]), _mm_setzero_si128());
2217                 __m128i pix2 = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&inb4ub[x*4]), _mm_setzero_si128());
2218                 __m128i blend = _mm_shufflelo_epi16(pix2, _MM_SHUFFLE(3, 3, 3, 3));
2219                 pix1 = _mm_add_epi16(pix1, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(pix2, pix1), 4), _mm_slli_epi16(blend, 4)));
2220                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix1, pix1));
2221         }
2222 #endif
2223 }
2224
2225 void DPSOFTRAST_Draw_Span_MixUniformColorBGRA8(const DPSOFTRAST_State_Draw_Span * RESTRICT span, unsigned char *out4ub, const unsigned char *in4ub, const float *color)
2226 {
2227 #ifdef SSE2_PRESENT
2228         int x, startx = span->startx, endx = span->endx;
2229         __m128i localcolor = _mm_shuffle_epi32(_mm_cvtps_epi32(_mm_mul_ps(_mm_loadu_ps(color), _mm_set1_ps(255.0f))), _MM_SHUFFLE(3, 0, 1, 2)), blend;
2230         localcolor = _mm_shuffle_epi32(_mm_packs_epi32(localcolor, localcolor), _MM_SHUFFLE(1, 0, 1, 0));
2231         blend = _mm_slli_epi16(_mm_shufflehi_epi16(_mm_shufflelo_epi16(localcolor, _MM_SHUFFLE(3, 3, 3, 3)), _MM_SHUFFLE(3, 3, 3, 3)), 4);
2232         for (x = startx;x+2 <= endx;x+=2)
2233         {
2234                 __m128i pix = _mm_unpacklo_epi8(_mm_loadl_epi64((const __m128i *)&in4ub[x*4]), _mm_setzero_si128());
2235                 pix = _mm_add_epi16(pix, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(localcolor, pix), 4), blend));
2236                 _mm_storel_epi64((__m128i *)&out4ub[x*4], _mm_packus_epi16(pix, pix));
2237         }
2238         if(x < endx)
2239         {
2240                 __m128i pix = _mm_unpacklo_epi8(_mm_cvtsi32_si128(*(const int *)&in4ub[x*4]), _mm_setzero_si128());
2241                 pix = _mm_add_epi16(pix, _mm_mulhi_epi16(_mm_slli_epi16(_mm_sub_epi16(localcolor, pix), 4), blend));
2242                 *(int *)&out4ub[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix, pix));
2243         }
2244 #endif
2245 }
2246
2247
2248
2249 void DPSOFTRAST_VertexShader_Generic(void)
2250 {
2251         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2252         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_COLOR], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_COLOR], dpsoftrast.draw.numvertices);
2253         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices);
2254         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SPECULAR)
2255                 DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1], dpsoftrast.draw.numvertices);
2256 }
2257
2258 void DPSOFTRAST_PixelShader_Generic(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2259 {
2260         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2261         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2262         unsigned char buffer_texture_lightmapbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2263         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2264         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2265         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_DIFFUSE)
2266         {
2267                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_FIRST, 2, buffer_z);
2268                 DPSOFTRAST_Draw_Span_MultiplyVaryingBGRA8(span, buffer_FragColorbgra8, buffer_texture_colorbgra8, 1, buffer_z);
2269                 if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SPECULAR)
2270                 {
2271                         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_lightmapbgra8, GL20TU_SECOND, 2, buffer_z);
2272                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
2273                         {
2274                                 // multiply
2275                                 DPSOFTRAST_Draw_Span_MultiplyBuffersBGRA8(span, buffer_FragColorbgra8, buffer_FragColorbgra8, buffer_texture_lightmapbgra8);
2276                         }
2277                         else if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
2278                         {
2279                                 // add
2280                                 DPSOFTRAST_Draw_Span_AddBuffersBGRA8(span, buffer_FragColorbgra8, buffer_FragColorbgra8, buffer_texture_lightmapbgra8);
2281                         }
2282                         else if (dpsoftrast.shader_permutation & SHADERPERMUTATION_VERTEXTEXTUREBLEND)
2283                         {
2284                                 // alphablend
2285                                 DPSOFTRAST_Draw_Span_MixBuffersBGRA8(span, buffer_FragColorbgra8, buffer_FragColorbgra8, buffer_texture_lightmapbgra8);
2286                         }
2287                 }
2288         }
2289         else
2290                 DPSOFTRAST_Draw_Span_VaryingBGRA8(span, buffer_FragColorbgra8, 1, buffer_z);
2291         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2292 }
2293
2294
2295
2296 void DPSOFTRAST_VertexShader_PostProcess(void)
2297 {
2298         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2299         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices);
2300         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1], dpsoftrast.draw.numvertices);
2301 }
2302
2303 void DPSOFTRAST_PixelShader_PostProcess(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2304 {
2305         // TODO: optimize!!  at the very least there is no reason to use texture sampling on the frame texture
2306         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2307         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2308         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2309         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2310         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_FragColorbgra8, GL20TU_FIRST, 2, buffer_z);
2311         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_BLOOM)
2312         {
2313                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_SECOND, 3, buffer_z);
2314                 DPSOFTRAST_Draw_Span_AddBloomBGRA8(span, buffer_FragColorbgra8, buffer_FragColorbgra8, buffer_texture_colorbgra8, dpsoftrast.uniform4f + DPSOFTRAST_UNIFORM_BloomColorSubtract * 4);
2315         }
2316         DPSOFTRAST_Draw_Span_MixUniformColorBGRA8(span, buffer_FragColorbgra8, buffer_FragColorbgra8, dpsoftrast.uniform4f + DPSOFTRAST_UNIFORM_ViewTintColor * 4);
2317         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SATURATION)
2318         {
2319                 // TODO: implement saturation
2320         }
2321         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GAMMARAMPS)
2322         {
2323                 // TODO: implement gammaramps
2324         }
2325         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2326 }
2327
2328
2329
2330 void DPSOFTRAST_VertexShader_Depth_Or_Shadow(void)
2331 {
2332         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2333 }
2334
2335 void DPSOFTRAST_PixelShader_Depth_Or_Shadow(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2336 {
2337         // this is never called (because colormask is off when this shader is used)
2338         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2339         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2340         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2341         memset(buffer_FragColorbgra8, 0, span->length*4);
2342         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2343 }
2344
2345
2346
2347 void DPSOFTRAST_VertexShader_FlatColor(void)
2348 {
2349         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2350         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_TexMatrixM1);
2351 }
2352
2353 void DPSOFTRAST_PixelShader_FlatColor(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2354 {
2355         int x, startx = span->startx, endx = span->endx;
2356         int Color_Ambienti[4];
2357         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2358         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2359         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2360         Color_Ambienti[2] = (int)(dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+0]*256.0f);
2361         Color_Ambienti[1] = (int)(dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+1]*256.0f);
2362         Color_Ambienti[0] = (int)(dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+2]*256.0f);
2363         Color_Ambienti[3] = (int)(dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Alpha*4+0]        *256.0f);
2364         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2365         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_COLOR, 2, buffer_z);
2366         for (x = startx;x < endx;x++)
2367         {
2368                 buffer_FragColorbgra8[x*4+0] = (buffer_texture_colorbgra8[x*4+0] * Color_Ambienti[0])>>8;
2369                 buffer_FragColorbgra8[x*4+1] = (buffer_texture_colorbgra8[x*4+1] * Color_Ambienti[1])>>8;
2370                 buffer_FragColorbgra8[x*4+2] = (buffer_texture_colorbgra8[x*4+2] * Color_Ambienti[2])>>8;
2371                 buffer_FragColorbgra8[x*4+3] = (buffer_texture_colorbgra8[x*4+3] * Color_Ambienti[3])>>8;
2372         }
2373         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2374 }
2375
2376
2377
2378 void DPSOFTRAST_VertexShader_VertexColor(void)
2379 {
2380         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2381         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_COLOR], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_COLOR], dpsoftrast.draw.numvertices);
2382         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_TexMatrixM1);
2383 }
2384
2385 void DPSOFTRAST_PixelShader_VertexColor(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2386 {
2387         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2388         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2389         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2390         int x, startx = span->startx, endx = span->endx;
2391         float Color_Ambient[4], Color_Diffuse[4];
2392         float data[4];
2393         float slope[4];
2394         float z;
2395         int arrayindex = DPSOFTRAST_ARRAY_COLOR;
2396         data[2] = span->data[0][arrayindex][0];
2397         data[1] = span->data[0][arrayindex][1];
2398         data[0] = span->data[0][arrayindex][2];
2399         data[3] = span->data[0][arrayindex][3];
2400         slope[2] = span->data[1][arrayindex][0];
2401         slope[1] = span->data[1][arrayindex][1];
2402         slope[0] = span->data[1][arrayindex][2];
2403         slope[3] = span->data[1][arrayindex][3];
2404         Color_Ambient[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+0];
2405         Color_Ambient[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+1];
2406         Color_Ambient[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+2];
2407         Color_Ambient[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Alpha*4+0];
2408         Color_Diffuse[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+0];
2409         Color_Diffuse[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+1];
2410         Color_Diffuse[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+2];
2411         Color_Diffuse[3] = 0.0f;
2412         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2413         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_COLOR, 2, buffer_z);
2414         for (x = startx;x < endx;x++)
2415         {
2416                 z = buffer_z[x];
2417                 buffer_FragColorbgra8[x*4+0] = (int)(buffer_texture_colorbgra8[x*4+0] * (Color_Ambient[0] + ((data[0] + slope[0]*x) * z) * Color_Diffuse[0]));
2418                 buffer_FragColorbgra8[x*4+1] = (int)(buffer_texture_colorbgra8[x*4+1] * (Color_Ambient[1] + ((data[1] + slope[1]*x) * z) * Color_Diffuse[1]));
2419                 buffer_FragColorbgra8[x*4+2] = (int)(buffer_texture_colorbgra8[x*4+2] * (Color_Ambient[2] + ((data[2] + slope[2]*x) * z) * Color_Diffuse[2]));
2420                 buffer_FragColorbgra8[x*4+3] = (int)(buffer_texture_colorbgra8[x*4+3] * (Color_Ambient[3] + ((data[3] + slope[3]*x) * z) * Color_Diffuse[3]));
2421         }
2422         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2423 }
2424
2425
2426
2427 void DPSOFTRAST_VertexShader_Lightmap(void)
2428 {
2429         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2430         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_TexMatrixM1);
2431         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD4], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD4], dpsoftrast.draw.numvertices);
2432 }
2433
2434 void DPSOFTRAST_PixelShader_Lightmap(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2435 {
2436 #ifdef SSE2_PRESENT
2437         unsigned char * RESTRICT pixelmask = span->pixelmask;
2438         unsigned char * RESTRICT pixel = (unsigned char *)dpsoftrast.fb_colorpixels[0] + span->start * 4;
2439         int x, startx = span->startx, endx = span->endx;
2440         __m128i Color_Ambientm, Color_Diffusem, Color_Glowm, Color_AmbientGlowm;
2441         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2442         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2443         unsigned char buffer_texture_lightmapbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2444         unsigned char buffer_texture_glowbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2445         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2446         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2447         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_COLOR, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2448         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_lightmapbgra8, GL20TU_LIGHTMAP, DPSOFTRAST_ARRAY_TEXCOORD4, buffer_z);
2449         if (dpsoftrast.user.alphatest || dpsoftrast.fb_blendmode != DPSOFTRAST_BLENDMODE_OPAQUE)
2450                 pixel = buffer_FragColorbgra8;
2451         Color_Ambientm = _mm_shuffle_epi32(_mm_cvtps_epi32(_mm_mul_ps(_mm_load_ps(&dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4]), _mm_set1_ps(256.0f))), _MM_SHUFFLE(3, 0, 1, 2));
2452         Color_Ambientm = _mm_and_si128(Color_Ambientm, _mm_setr_epi32(0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0));
2453         Color_Ambientm = _mm_or_si128(Color_Ambientm, _mm_setr_epi32(0, 0, 0, (int)(dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Alpha*4+0]*255.0f)));
2454         Color_Ambientm = _mm_packs_epi32(Color_Ambientm, Color_Ambientm);
2455         Color_Diffusem = _mm_shuffle_epi32(_mm_cvtps_epi32(_mm_mul_ps(_mm_load_ps(&dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4]), _mm_set1_ps(256.0f))), _MM_SHUFFLE(3, 0, 1, 2));
2456         Color_Diffusem = _mm_and_si128(Color_Diffusem, _mm_setr_epi32(0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0));
2457         Color_Diffusem = _mm_packs_epi32(Color_Diffusem, Color_Diffusem);
2458         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GLOW)
2459         {
2460                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_glowbgra8, GL20TU_GLOW, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2461                 Color_Glowm = _mm_shuffle_epi32(_mm_cvtps_epi32(_mm_mul_ps(_mm_load_ps(&dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4]), _mm_set1_ps(256.0f))), _MM_SHUFFLE(3, 0, 1, 2));
2462                 Color_Glowm = _mm_and_si128(Color_Glowm, _mm_setr_epi32(0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0));
2463                 Color_Glowm = _mm_packs_epi32(Color_Glowm, Color_Glowm);
2464                 Color_AmbientGlowm = _mm_unpacklo_epi64(Color_Ambientm, Color_Glowm);
2465                 for (x = startx;x < endx;x++)
2466                 {
2467                         __m128i color, lightmap, glow, pix;
2468                         if (x + 4 <= endx && *(const unsigned int *)&pixelmask[x] == 0x01010101)
2469                         {
2470                                 __m128i pix2;
2471                                 color = _mm_loadu_si128((const __m128i *)&buffer_texture_colorbgra8[x*4]);
2472                                 lightmap = _mm_loadu_si128((const __m128i *)&buffer_texture_lightmapbgra8[x*4]);
2473                                 glow = _mm_loadu_si128((const __m128i *)&buffer_texture_glowbgra8[x*4]);
2474                                 pix = _mm_add_epi16(_mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(Color_Diffusem, _mm_unpacklo_epi8(_mm_setzero_si128(), lightmap)), Color_Ambientm),
2475                                                                                                         _mm_unpacklo_epi8(_mm_setzero_si128(), color)),
2476                                                                         _mm_mulhi_epu16(Color_Glowm, _mm_unpacklo_epi8(_mm_setzero_si128(), glow)));
2477                                 pix2 = _mm_add_epi16(_mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(Color_Diffusem, _mm_unpackhi_epi8(_mm_setzero_si128(), lightmap)), Color_Ambientm),
2478                                                                                                         _mm_unpackhi_epi8(_mm_setzero_si128(), color)),
2479                                                                         _mm_mulhi_epu16(Color_Glowm, _mm_unpackhi_epi8(_mm_setzero_si128(), glow)));
2480                                 _mm_storeu_si128((__m128i *)&pixel[x*4], _mm_packus_epi16(pix, pix2));
2481                                 x += 3;
2482                                 continue;
2483                         }
2484                         if(!pixelmask[x])
2485                                 continue;
2486                         color = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&buffer_texture_colorbgra8[x*4]));
2487                         lightmap = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&buffer_texture_lightmapbgra8[x*4]));
2488                         glow = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&buffer_texture_glowbgra8[x*4]));
2489                         pix = _mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(Color_Diffusem, lightmap), Color_AmbientGlowm), _mm_unpacklo_epi64(color, glow));
2490                         pix = _mm_add_epi16(pix, _mm_shuffle_epi32(pix, _MM_SHUFFLE(3, 2, 3, 2)));
2491                         *(int *)&pixel[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix, pix));
2492                 }
2493         }
2494         else
2495         {
2496                 for (x = startx;x < endx;x++)
2497                 {
2498                         __m128i color, lightmap, pix;
2499                         if (x + 4 <= endx && *(const unsigned int *)&pixelmask[x] == 0x01010101)
2500                         {
2501                                 __m128i pix2;
2502                                 color = _mm_loadu_si128((const __m128i *)&buffer_texture_colorbgra8[x*4]);
2503                                 lightmap = _mm_loadu_si128((const __m128i *)&buffer_texture_lightmapbgra8[x*4]);
2504                                 pix = _mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(Color_Diffusem, _mm_unpacklo_epi8(_mm_setzero_si128(), lightmap)), Color_Ambientm),
2505                                                                           _mm_unpacklo_epi8(_mm_setzero_si128(), color));
2506                                 pix2 = _mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(Color_Diffusem, _mm_unpackhi_epi8(_mm_setzero_si128(), lightmap)), Color_Ambientm),
2507                                                                            _mm_unpackhi_epi8(_mm_setzero_si128(), color));
2508                                 _mm_storeu_si128((__m128i *)&pixel[x*4], _mm_packus_epi16(pix, pix2));
2509                                 x += 3;
2510                                 continue;
2511                         }
2512                         if(!pixelmask[x])
2513                                 continue;
2514                         color = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&buffer_texture_colorbgra8[x*4]));
2515                         lightmap = _mm_unpacklo_epi8(_mm_setzero_si128(), _mm_cvtsi32_si128(*(const int *)&buffer_texture_lightmapbgra8[x*4]));
2516                         pix = _mm_mulhi_epu16(_mm_add_epi16(_mm_mulhi_epu16(lightmap, Color_Diffusem), Color_Ambientm), color);
2517                         *(int *)&pixel[x*4] = _mm_cvtsi128_si32(_mm_packus_epi16(pix, pix));
2518                 }
2519         }
2520         if(pixel == buffer_FragColorbgra8)
2521                 DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2522 #endif
2523 }
2524
2525
2526
2527 void DPSOFTRAST_VertexShader_FakeLight(void)
2528 {
2529         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2530 }
2531
2532 void DPSOFTRAST_PixelShader_FakeLight(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2533 {
2534         // TODO: IMPLEMENT
2535         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2536         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2537         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2538         memset(buffer_FragColorbgra8, 0, span->length*4);
2539         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2540 }
2541
2542
2543
2544 void DPSOFTRAST_VertexShader_LightDirectionMap_ModelSpace(void)
2545 {
2546         DPSOFTRAST_VertexShader_Lightmap();
2547 }
2548
2549 void DPSOFTRAST_PixelShader_LightDirectionMap_ModelSpace(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2550 {
2551         DPSOFTRAST_PixelShader_Lightmap(span);
2552         // TODO: IMPLEMENT
2553 }
2554
2555
2556
2557 void DPSOFTRAST_VertexShader_LightDirectionMap_TangentSpace(void)
2558 {
2559         DPSOFTRAST_VertexShader_Lightmap();
2560 }
2561
2562 void DPSOFTRAST_PixelShader_LightDirectionMap_TangentSpace(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2563 {
2564         DPSOFTRAST_PixelShader_Lightmap(span);
2565         // TODO: IMPLEMENT
2566 }
2567
2568
2569
2570 void DPSOFTRAST_VertexShader_LightDirection(void)
2571 {
2572         int i;
2573         int numvertices = dpsoftrast.draw.numvertices;
2574         float LightDir[4];
2575         float LightVector[4];
2576         float EyePosition[4];
2577         float EyeVectorModelSpace[4];
2578         float EyeVector[4];
2579         float position[4];
2580         float svector[4];
2581         float tvector[4];
2582         float normal[4];
2583         LightDir[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightDir*4+0];
2584         LightDir[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightDir*4+1];
2585         LightDir[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightDir*4+2];
2586         LightDir[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightDir*4+3];
2587         EyePosition[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+0];
2588         EyePosition[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+1];
2589         EyePosition[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+2];
2590         EyePosition[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+3];
2591         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2592         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_TexMatrixM1);
2593         for (i = 0;i < numvertices;i++)
2594         {
2595                 position[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+0];
2596                 position[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+1];
2597                 position[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+2];
2598                 svector[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+0];
2599                 svector[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+1];
2600                 svector[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+2];
2601                 tvector[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+0];
2602                 tvector[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+1];
2603                 tvector[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+2];
2604                 normal[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+0];
2605                 normal[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+1];
2606                 normal[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+2];
2607                 LightVector[0] = svector[0] * LightDir[0] + svector[1] * LightDir[1] + svector[2] * LightDir[2];
2608                 LightVector[1] = tvector[0] * LightDir[0] + tvector[1] * LightDir[1] + tvector[2] * LightDir[2];
2609                 LightVector[2] = normal[0] * LightDir[0] + normal[1] * LightDir[1] + normal[2] * LightDir[2];
2610                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+0] = LightVector[0];
2611                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+1] = LightVector[1];
2612                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+2] = LightVector[2];
2613                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+3] = 0.0f;
2614                 EyeVectorModelSpace[0] = EyePosition[0] - position[0];
2615                 EyeVectorModelSpace[1] = EyePosition[1] - position[1];
2616                 EyeVectorModelSpace[2] = EyePosition[2] - position[2];
2617                 EyeVector[0] = svector[0] * EyeVectorModelSpace[0] + svector[1] * EyeVectorModelSpace[1] + svector[2] * EyeVectorModelSpace[2];
2618                 EyeVector[1] = tvector[0] * EyeVectorModelSpace[0] + tvector[1] * EyeVectorModelSpace[1] + tvector[2] * EyeVectorModelSpace[2];
2619                 EyeVector[2] = normal[0]  * EyeVectorModelSpace[0] + normal[1]  * EyeVectorModelSpace[1] + normal[2]  * EyeVectorModelSpace[2];
2620                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+0] = EyeVector[0];
2621                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+1] = EyeVector[1];
2622                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+2] = EyeVector[2];
2623                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+3] = 0.0f;
2624         }
2625 }
2626
2627 #define DPSOFTRAST_Min(a,b) ((a) < (b) ? (a) : (b))
2628 #define DPSOFTRAST_Max(a,b) ((a) > (b) ? (a) : (b))
2629 #define DPSOFTRAST_Vector3Dot(a,b) ((a)[0]*(b)[0]+(a)[1]*(b)[1]+(a)[2]*(b)[2])
2630 #define DPSOFTRAST_Vector3LengthSquared(v) (DPSOFTRAST_Vector3Dot((v),(v)))
2631 #define DPSOFTRAST_Vector3Length(v) (sqrt(DPSOFTRAST_Vector3LengthSquared(v)))
2632 #define DPSOFTRAST_Vector3Normalize(v)\
2633 do\
2634 {\
2635         float len = sqrt(DPSOFTRAST_Vector3Dot(v,v));\
2636         if (len)\
2637         {\
2638                 len = 1.0f / len;\
2639                 v[0] *= len;\
2640                 v[1] *= len;\
2641                 v[2] *= len;\
2642         }\
2643 }\
2644 while(0)
2645
2646 void DPSOFTRAST_PixelShader_LightDirection(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2647 {
2648         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2649         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2650         unsigned char buffer_texture_normalbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2651         unsigned char buffer_texture_glossbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2652         unsigned char buffer_texture_glowbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2653         unsigned char buffer_texture_pantsbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2654         unsigned char buffer_texture_shirtbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2655         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2656         int x, startx = span->startx, endx = span->endx;
2657         float Color_Ambient[4], Color_Diffuse[4], Color_Specular[4], Color_Glow[4], Color_Pants[4], Color_Shirt[4], LightColor[4];
2658         float LightVectordata[4];
2659         float LightVectorslope[4];
2660         float EyeVectordata[4];
2661         float EyeVectorslope[4];
2662         float z;
2663         float diffusetex[4];
2664         float glosstex[4];
2665         float surfacenormal[4];
2666         float lightnormal[4];
2667         float eyenormal[4];
2668         float specularnormal[4];
2669         float diffuse;
2670         float specular;
2671         float SpecularPower;
2672         int d[4];
2673         Color_Glow[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+0];
2674         Color_Glow[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+1];
2675         Color_Glow[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+2];
2676         Color_Glow[3] = 0.0f;
2677         Color_Ambient[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+0];
2678         Color_Ambient[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+1];
2679         Color_Ambient[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+2];
2680         Color_Ambient[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Alpha*4+0];
2681         Color_Pants[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+0];
2682         Color_Pants[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+1];
2683         Color_Pants[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+2];
2684         Color_Pants[3] = 0.0f;
2685         Color_Shirt[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+0];
2686         Color_Shirt[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+1];
2687         Color_Shirt[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+2];
2688         Color_Shirt[3] = 0.0f;
2689         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
2690         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_COLOR, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2691         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
2692         {
2693                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_pantsbgra8, GL20TU_PANTS, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2694                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_shirtbgra8, GL20TU_SHIRT, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2695         }
2696         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GLOW)
2697         {
2698                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_glowbgra8, GL20TU_GLOW, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2699         }
2700         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SPECULAR)
2701         {
2702                 Color_Diffuse[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+0];
2703                 Color_Diffuse[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+1];
2704                 Color_Diffuse[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+2];
2705                 Color_Diffuse[3] = 0.0f;
2706                 LightColor[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+0];
2707                 LightColor[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+1];
2708                 LightColor[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+2];
2709                 LightColor[3] = 0.0f;
2710                 LightVectordata[0]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][0];
2711                 LightVectordata[1]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][1];
2712                 LightVectordata[2]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][2];
2713                 LightVectordata[3]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][3];
2714                 LightVectorslope[0] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][0];
2715                 LightVectorslope[1] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][1];
2716                 LightVectorslope[2] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][2];
2717                 LightVectorslope[3] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][3];
2718                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_normalbgra8, GL20TU_NORMAL, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2719                 Color_Specular[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+0];
2720                 Color_Specular[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+1];
2721                 Color_Specular[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+2];
2722                 Color_Specular[3] = 0.0f;
2723                 SpecularPower = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_SpecularPower*4+0] * (1.0f / 255.0f);
2724                 EyeVectordata[0]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][0];
2725                 EyeVectordata[1]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][1];
2726                 EyeVectordata[2]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][2];
2727                 EyeVectordata[3]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][3];
2728                 EyeVectorslope[0]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][0];
2729                 EyeVectorslope[1]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][1];
2730                 EyeVectorslope[2]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][2];
2731                 EyeVectorslope[3]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][3];
2732                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_glossbgra8, GL20TU_GLOSS, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2733                 for (x = startx;x < endx;x++)
2734                 {
2735                         z = buffer_z[x];
2736                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
2737                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
2738                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
2739                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
2740                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
2741                         {
2742                                 diffusetex[0] += buffer_texture_pantsbgra8[x*4+0] * Color_Pants[0] + buffer_texture_shirtbgra8[x*4+0] * Color_Shirt[0];
2743                                 diffusetex[1] += buffer_texture_pantsbgra8[x*4+1] * Color_Pants[1] + buffer_texture_shirtbgra8[x*4+1] * Color_Shirt[1];
2744                                 diffusetex[2] += buffer_texture_pantsbgra8[x*4+2] * Color_Pants[2] + buffer_texture_shirtbgra8[x*4+2] * Color_Shirt[2];
2745                                 diffusetex[3] += buffer_texture_pantsbgra8[x*4+3] * Color_Pants[3] + buffer_texture_shirtbgra8[x*4+3] * Color_Shirt[3];
2746                         }
2747                         glosstex[0] = buffer_texture_glossbgra8[x*4+0];
2748                         glosstex[1] = buffer_texture_glossbgra8[x*4+1];
2749                         glosstex[2] = buffer_texture_glossbgra8[x*4+2];
2750                         glosstex[3] = buffer_texture_glossbgra8[x*4+3];
2751                         surfacenormal[0] = buffer_texture_normalbgra8[x*4+2] * (1.0f / 128.0f) - 1.0f;
2752                         surfacenormal[1] = buffer_texture_normalbgra8[x*4+1] * (1.0f / 128.0f) - 1.0f;
2753                         surfacenormal[2] = buffer_texture_normalbgra8[x*4+0] * (1.0f / 128.0f) - 1.0f;
2754                         DPSOFTRAST_Vector3Normalize(surfacenormal);
2755
2756                         lightnormal[0] = (LightVectordata[0] + LightVectorslope[0]*x) * z;
2757                         lightnormal[1] = (LightVectordata[1] + LightVectorslope[1]*x) * z;
2758                         lightnormal[2] = (LightVectordata[2] + LightVectorslope[2]*x) * z;
2759                         DPSOFTRAST_Vector3Normalize(lightnormal);
2760
2761                         eyenormal[0] = (EyeVectordata[0] + EyeVectorslope[0]*x) * z;
2762                         eyenormal[1] = (EyeVectordata[1] + EyeVectorslope[1]*x) * z;
2763                         eyenormal[2] = (EyeVectordata[2] + EyeVectorslope[2]*x) * z;
2764                         DPSOFTRAST_Vector3Normalize(eyenormal);
2765
2766                         specularnormal[0] = lightnormal[0] + eyenormal[0];
2767                         specularnormal[1] = lightnormal[1] + eyenormal[1];
2768                         specularnormal[2] = lightnormal[2] + eyenormal[2];
2769                         DPSOFTRAST_Vector3Normalize(specularnormal);
2770
2771                         diffuse = DPSOFTRAST_Vector3Dot(surfacenormal, lightnormal);if (diffuse < 0.0f) diffuse = 0.0f;
2772                         specular = DPSOFTRAST_Vector3Dot(surfacenormal, specularnormal);if (specular < 0.0f) specular = 0.0f;
2773                         specular = pow(specular, SpecularPower * glosstex[3]);
2774                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GLOW)
2775                         {
2776                                 d[0] = (int)(buffer_texture_glowbgra8[x*4+0] * Color_Glow[0] + diffusetex[0] * Color_Ambient[0] + (diffusetex[0] * Color_Diffuse[0] * diffuse + glosstex[0] * Color_Specular[0] * specular) * LightColor[0]);if (d[0] > 255) d[0] = 255;
2777                                 d[1] = (int)(buffer_texture_glowbgra8[x*4+1] * Color_Glow[1] + diffusetex[1] * Color_Ambient[1] + (diffusetex[1] * Color_Diffuse[1] * diffuse + glosstex[1] * Color_Specular[1] * specular) * LightColor[1]);if (d[1] > 255) d[1] = 255;
2778                                 d[2] = (int)(buffer_texture_glowbgra8[x*4+2] * Color_Glow[2] + diffusetex[2] * Color_Ambient[2] + (diffusetex[2] * Color_Diffuse[2] * diffuse + glosstex[2] * Color_Specular[2] * specular) * LightColor[2]);if (d[2] > 255) d[2] = 255;
2779                                 d[3] = (int)(                                                  diffusetex[3] * Color_Ambient[3]);if (d[3] > 255) d[3] = 255;
2780                         }
2781                         else
2782                         {
2783                                 d[0] = (int)(                                                  diffusetex[0] * Color_Ambient[0] + (diffusetex[0] * Color_Diffuse[0] * diffuse + glosstex[0] * Color_Specular[0] * specular) * LightColor[0]);if (d[0] > 255) d[0] = 255;
2784                                 d[1] = (int)(                                                  diffusetex[1] * Color_Ambient[1] + (diffusetex[1] * Color_Diffuse[1] * diffuse + glosstex[1] * Color_Specular[1] * specular) * LightColor[1]);if (d[1] > 255) d[1] = 255;
2785                                 d[2] = (int)(                                                  diffusetex[2] * Color_Ambient[2] + (diffusetex[2] * Color_Diffuse[2] * diffuse + glosstex[2] * Color_Specular[2] * specular) * LightColor[2]);if (d[2] > 255) d[2] = 255;
2786                                 d[3] = (int)(                                                  diffusetex[3] * Color_Ambient[3]);if (d[3] > 255) d[3] = 255;
2787                         }
2788                         buffer_FragColorbgra8[x*4+0] = d[0];
2789                         buffer_FragColorbgra8[x*4+1] = d[1];
2790                         buffer_FragColorbgra8[x*4+2] = d[2];
2791                         buffer_FragColorbgra8[x*4+3] = d[3];
2792                 }
2793         }
2794         else if (dpsoftrast.shader_permutation & SHADERPERMUTATION_DIFFUSE)
2795         {
2796                 Color_Diffuse[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+0];
2797                 Color_Diffuse[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+1];
2798                 Color_Diffuse[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+2];
2799                 Color_Diffuse[3] = 0.0f;
2800                 LightColor[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+0];
2801                 LightColor[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+1];
2802                 LightColor[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+2];
2803                 LightColor[3] = 0.0f;
2804                 LightVectordata[0]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][0];
2805                 LightVectordata[1]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][1];
2806                 LightVectordata[2]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][2];
2807                 LightVectordata[3]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][3];
2808                 LightVectorslope[0] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][0];
2809                 LightVectorslope[1] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][1];
2810                 LightVectorslope[2] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][2];
2811                 LightVectorslope[3] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][3];
2812                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_normalbgra8, GL20TU_NORMAL, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
2813                 for (x = startx;x < endx;x++)
2814                 {
2815                         z = buffer_z[x];
2816                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
2817                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
2818                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
2819                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
2820                         surfacenormal[0] = buffer_texture_normalbgra8[x*4+2] * (1.0f / 128.0f) - 1.0f;
2821                         surfacenormal[1] = buffer_texture_normalbgra8[x*4+1] * (1.0f / 128.0f) - 1.0f;
2822                         surfacenormal[2] = buffer_texture_normalbgra8[x*4+0] * (1.0f / 128.0f) - 1.0f;
2823                         DPSOFTRAST_Vector3Normalize(surfacenormal);
2824
2825                         lightnormal[0] = (LightVectordata[0] + LightVectorslope[0]*x) * z;
2826                         lightnormal[1] = (LightVectordata[1] + LightVectorslope[1]*x) * z;
2827                         lightnormal[2] = (LightVectordata[2] + LightVectorslope[2]*x) * z;
2828                         DPSOFTRAST_Vector3Normalize(lightnormal);
2829
2830                         diffuse = DPSOFTRAST_Vector3Dot(surfacenormal, lightnormal);if (diffuse < 0.0f) diffuse = 0.0f;
2831                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GLOW)
2832                         {
2833                                 d[0] = (int)(buffer_texture_glowbgra8[x*4+0] * Color_Glow[0] + diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse * LightColor[0]));if (d[0] > 255) d[0] = 255;
2834                                 d[1] = (int)(buffer_texture_glowbgra8[x*4+1] * Color_Glow[1] + diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse * LightColor[1]));if (d[1] > 255) d[1] = 255;
2835                                 d[2] = (int)(buffer_texture_glowbgra8[x*4+2] * Color_Glow[2] + diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse * LightColor[2]));if (d[2] > 255) d[2] = 255;
2836                                 d[3] = (int)(                                                  diffusetex[3] * (Color_Ambient[3]                                             ));if (d[3] > 255) d[3] = 255;
2837                         }
2838                         else
2839                         {
2840                                 d[0] = (int)(                                                + diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse * LightColor[0]));if (d[0] > 255) d[0] = 255;
2841                                 d[1] = (int)(                                                + diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse * LightColor[1]));if (d[1] > 255) d[1] = 255;
2842                                 d[2] = (int)(                                                + diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse * LightColor[2]));if (d[2] > 255) d[2] = 255;
2843                                 d[3] = (int)(                                                  diffusetex[3] * (Color_Ambient[3]                                             ));if (d[3] > 255) d[3] = 255;
2844                         }
2845                         buffer_FragColorbgra8[x*4+0] = d[0];
2846                         buffer_FragColorbgra8[x*4+1] = d[1];
2847                         buffer_FragColorbgra8[x*4+2] = d[2];
2848                         buffer_FragColorbgra8[x*4+3] = d[3];
2849                 }
2850         }
2851         else
2852         {
2853                 for (x = startx;x < endx;x++)
2854                 {
2855                         z = buffer_z[x];
2856                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
2857                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
2858                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
2859                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
2860
2861                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_GLOW)
2862                         {
2863                                 d[0] = (int)(buffer_texture_glowbgra8[x*4+0] * Color_Glow[0] + diffusetex[0] * Color_Ambient[0]);if (d[0] > 255) d[0] = 255;
2864                                 d[1] = (int)(buffer_texture_glowbgra8[x*4+1] * Color_Glow[1] + diffusetex[1] * Color_Ambient[1]);if (d[1] > 255) d[1] = 255;
2865                                 d[2] = (int)(buffer_texture_glowbgra8[x*4+2] * Color_Glow[2] + diffusetex[2] * Color_Ambient[2]);if (d[2] > 255) d[2] = 255;
2866                                 d[3] = (int)(                                                  diffusetex[3] * Color_Ambient[3]);if (d[3] > 255) d[3] = 255;
2867                         }
2868                         else
2869                         {
2870                                 d[0] = (int)(                                                  diffusetex[0] * Color_Ambient[0]);if (d[0] > 255) d[0] = 255;
2871                                 d[1] = (int)(                                                  diffusetex[1] * Color_Ambient[1]);if (d[1] > 255) d[1] = 255;
2872                                 d[2] = (int)(                                                  diffusetex[2] * Color_Ambient[2]);if (d[2] > 255) d[2] = 255;
2873                                 d[3] = (int)(                                                  diffusetex[3] * Color_Ambient[3]);if (d[3] > 255) d[3] = 255;
2874                         }
2875                         buffer_FragColorbgra8[x*4+0] = d[0];
2876                         buffer_FragColorbgra8[x*4+1] = d[1];
2877                         buffer_FragColorbgra8[x*4+2] = d[2];
2878                         buffer_FragColorbgra8[x*4+3] = d[3];
2879                 }
2880         }
2881         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
2882 }
2883
2884
2885
2886 void DPSOFTRAST_VertexShader_LightSource(void)
2887 {
2888         int i;
2889         int numvertices = dpsoftrast.draw.numvertices;
2890         float LightPosition[4];
2891         float LightVector[4];
2892         float LightVectorModelSpace[4];
2893         float EyePosition[4];
2894         float EyeVectorModelSpace[4];
2895         float EyeVector[4];
2896         float position[4];
2897         float svector[4];
2898         float tvector[4];
2899         float normal[4];
2900         LightPosition[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightPosition*4+0];
2901         LightPosition[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightPosition*4+1];
2902         LightPosition[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightPosition*4+2];
2903         LightPosition[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightPosition*4+3];
2904         EyePosition[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+0];
2905         EyePosition[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+1];
2906         EyePosition[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+2];
2907         EyePosition[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_EyePosition*4+3];
2908         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
2909         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD0], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_TexMatrixM1);
2910         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD3], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelToLightM1);
2911         DPSOFTRAST_Array_Copy(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD4], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD4], dpsoftrast.draw.numvertices);
2912         for (i = 0;i < numvertices;i++)
2913         {
2914                 position[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+0];
2915                 position[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+1];
2916                 position[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+2];
2917                 svector[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+0];
2918                 svector[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+1];
2919                 svector[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+2];
2920                 tvector[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+0];
2921                 tvector[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+1];
2922                 tvector[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+2];
2923                 normal[0] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+0];
2924                 normal[1] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+1];
2925                 normal[2] = dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_TEXCOORD3][i*4+2];
2926                 LightVectorModelSpace[0] = LightPosition[0] - position[0];
2927                 LightVectorModelSpace[1] = LightPosition[1] - position[1];
2928                 LightVectorModelSpace[2] = LightPosition[2] - position[2];
2929                 LightVector[0] = svector[0] * LightVectorModelSpace[0] + svector[1] * LightVectorModelSpace[1] + svector[2] * LightVectorModelSpace[2];
2930                 LightVector[1] = tvector[0] * LightVectorModelSpace[0] + tvector[1] * LightVectorModelSpace[1] + tvector[2] * LightVectorModelSpace[2];
2931                 LightVector[2] = normal[0]  * LightVectorModelSpace[0] + normal[1]  * LightVectorModelSpace[1] + normal[2]  * LightVectorModelSpace[2];
2932                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+0] = LightVector[0];
2933                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+1] = LightVector[1];
2934                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+2] = LightVector[2];
2935                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD1][i*4+3] = 0.0f;
2936                 EyeVectorModelSpace[0] = EyePosition[0] - position[0];
2937                 EyeVectorModelSpace[1] = EyePosition[1] - position[1];
2938                 EyeVectorModelSpace[2] = EyePosition[2] - position[2];
2939                 EyeVector[0] = svector[0] * EyeVectorModelSpace[0] + svector[1] * EyeVectorModelSpace[1] + svector[2] * EyeVectorModelSpace[2];
2940                 EyeVector[1] = tvector[0] * EyeVectorModelSpace[0] + tvector[1] * EyeVectorModelSpace[1] + tvector[2] * EyeVectorModelSpace[2];
2941                 EyeVector[2] = normal[0]  * EyeVectorModelSpace[0] + normal[1]  * EyeVectorModelSpace[1] + normal[2]  * EyeVectorModelSpace[2];
2942                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+0] = EyeVector[0];
2943                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+1] = EyeVector[1];
2944                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+2] = EyeVector[2];
2945                 dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_TEXCOORD2][i*4+3] = 0.0f;
2946         }
2947 }
2948
2949 void DPSOFTRAST_PixelShader_LightSource(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
2950 {
2951         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
2952         unsigned char buffer_texture_colorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2953         unsigned char buffer_texture_normalbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2954         unsigned char buffer_texture_glossbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2955         unsigned char buffer_texture_cubebgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2956         unsigned char buffer_texture_pantsbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2957         unsigned char buffer_texture_shirtbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2958         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
2959         int x, startx = span->startx, endx = span->endx;
2960         float Color_Ambient[4], Color_Diffuse[4], Color_Specular[4], Color_Glow[4], Color_Pants[4], Color_Shirt[4], LightColor[4];
2961         float CubeVectordata[4];
2962         float CubeVectorslope[4];
2963         float LightVectordata[4];
2964         float LightVectorslope[4];
2965         float EyeVectordata[4];
2966         float EyeVectorslope[4];
2967         float z;
2968         float diffusetex[4];
2969         float glosstex[4];
2970         float surfacenormal[4];
2971         float lightnormal[4];
2972         float eyenormal[4];
2973         float specularnormal[4];
2974         float diffuse;
2975         float specular;
2976         float SpecularPower;
2977         float CubeVector[4];
2978         float attenuation;
2979         int d[4];
2980         Color_Glow[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+0];
2981         Color_Glow[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+1];
2982         Color_Glow[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Glow*4+2];
2983         Color_Glow[3] = 0.0f;
2984         Color_Ambient[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+0];
2985         Color_Ambient[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+1];
2986         Color_Ambient[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Ambient*4+2];
2987         Color_Ambient[3] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Alpha*4+0];
2988         Color_Diffuse[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+0];
2989         Color_Diffuse[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+1];
2990         Color_Diffuse[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Diffuse*4+2];
2991         Color_Diffuse[3] = 0.0f;
2992         Color_Specular[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+0];
2993         Color_Specular[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+1];
2994         Color_Specular[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Specular*4+2];
2995         Color_Specular[3] = 0.0f;
2996         Color_Pants[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+0];
2997         Color_Pants[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+1];
2998         Color_Pants[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Pants*4+2];
2999         Color_Pants[3] = 0.0f;
3000         Color_Shirt[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+0];
3001         Color_Shirt[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+1];
3002         Color_Shirt[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_Color_Shirt*4+2];
3003         Color_Shirt[3] = 0.0f;
3004         LightColor[2] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+0];
3005         LightColor[1] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+1];
3006         LightColor[0] = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_LightColor*4+2];
3007         LightColor[3] = 0.0f;
3008         SpecularPower = dpsoftrast.uniform4f[DPSOFTRAST_UNIFORM_SpecularPower*4+0] * (1.0f / 255.0f);
3009         EyeVectordata[0]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][0];
3010         EyeVectordata[1]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][1];
3011         EyeVectordata[2]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][2];
3012         EyeVectordata[3]    = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD2][3];
3013         EyeVectorslope[0]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][0];
3014         EyeVectorslope[1]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][1];
3015         EyeVectorslope[2]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][2];
3016         EyeVectorslope[3]   = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD2][3];
3017         LightVectordata[0]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][0];
3018         LightVectordata[1]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][1];
3019         LightVectordata[2]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][2];
3020         LightVectordata[3]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD1][3];
3021         LightVectorslope[0] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][0];
3022         LightVectorslope[1] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][1];
3023         LightVectorslope[2] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][2];
3024         LightVectorslope[3] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD1][3];
3025         CubeVectordata[0]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD3][0];
3026         CubeVectordata[1]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD3][1];
3027         CubeVectordata[2]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD3][2];
3028         CubeVectordata[3]  = span->data[0][DPSOFTRAST_ARRAY_TEXCOORD3][3];
3029         CubeVectorslope[0] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD3][0];
3030         CubeVectorslope[1] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD3][1];
3031         CubeVectorslope[2] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD3][2];
3032         CubeVectorslope[3] = span->data[1][DPSOFTRAST_ARRAY_TEXCOORD3][3];
3033         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3034         memset(buffer_FragColorbgra8 + startx*4, 0, (endx-startx)*4); // clear first, because we skip writing black pixels, and there are a LOT of them...
3035         DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_colorbgra8, GL20TU_COLOR, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3036         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
3037         {
3038                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_pantsbgra8, GL20TU_PANTS, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3039                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_shirtbgra8, GL20TU_SHIRT, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3040         }
3041         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_CUBEFILTER)
3042                 DPSOFTRAST_Draw_Span_TextureCubeVaryingBGRA8(span, buffer_texture_cubebgra8, GL20TU_CUBE, DPSOFTRAST_ARRAY_TEXCOORD3, buffer_z);
3043         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SPECULAR)
3044         {
3045                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_normalbgra8, GL20TU_NORMAL, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3046                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_glossbgra8, GL20TU_GLOSS, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3047                 for (x = startx;x < endx;x++)
3048                 {
3049                         z = buffer_z[x];
3050                         CubeVector[0] = (CubeVectordata[0] + CubeVectorslope[0]*x) * z;
3051                         CubeVector[1] = (CubeVectordata[1] + CubeVectorslope[1]*x) * z;
3052                         CubeVector[2] = (CubeVectordata[2] + CubeVectorslope[2]*x) * z;
3053                         attenuation = 1.0f - DPSOFTRAST_Vector3LengthSquared(CubeVector);
3054                         if (attenuation < 0.01f)
3055                                 continue;
3056                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SHADOWMAP2D)
3057                         {
3058                                 attenuation *= DPSOFTRAST_SampleShadowmap(CubeVector);
3059                                 if (attenuation < 0.01f)
3060                                         continue;
3061                         }
3062
3063                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
3064                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
3065                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
3066                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
3067                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
3068                         {
3069                                 diffusetex[0] += buffer_texture_pantsbgra8[x*4+0] * Color_Pants[0] + buffer_texture_shirtbgra8[x*4+0] * Color_Shirt[0];
3070                                 diffusetex[1] += buffer_texture_pantsbgra8[x*4+1] * Color_Pants[1] + buffer_texture_shirtbgra8[x*4+1] * Color_Shirt[1];
3071                                 diffusetex[2] += buffer_texture_pantsbgra8[x*4+2] * Color_Pants[2] + buffer_texture_shirtbgra8[x*4+2] * Color_Shirt[2];
3072                                 diffusetex[3] += buffer_texture_pantsbgra8[x*4+3] * Color_Pants[3] + buffer_texture_shirtbgra8[x*4+3] * Color_Shirt[3];
3073                         }
3074                         glosstex[0] = buffer_texture_glossbgra8[x*4+0];
3075                         glosstex[1] = buffer_texture_glossbgra8[x*4+1];
3076                         glosstex[2] = buffer_texture_glossbgra8[x*4+2];
3077                         glosstex[3] = buffer_texture_glossbgra8[x*4+3];
3078                         surfacenormal[0] = buffer_texture_normalbgra8[x*4+2] * (1.0f / 128.0f) - 1.0f;
3079                         surfacenormal[1] = buffer_texture_normalbgra8[x*4+1] * (1.0f / 128.0f) - 1.0f;
3080                         surfacenormal[2] = buffer_texture_normalbgra8[x*4+0] * (1.0f / 128.0f) - 1.0f;
3081                         DPSOFTRAST_Vector3Normalize(surfacenormal);
3082
3083                         lightnormal[0] = (LightVectordata[0] + LightVectorslope[0]*x) * z;
3084                         lightnormal[1] = (LightVectordata[1] + LightVectorslope[1]*x) * z;
3085                         lightnormal[2] = (LightVectordata[2] + LightVectorslope[2]*x) * z;
3086                         DPSOFTRAST_Vector3Normalize(lightnormal);
3087
3088                         eyenormal[0] = (EyeVectordata[0] + EyeVectorslope[0]*x) * z;
3089                         eyenormal[1] = (EyeVectordata[1] + EyeVectorslope[1]*x) * z;
3090                         eyenormal[2] = (EyeVectordata[2] + EyeVectorslope[2]*x) * z;
3091                         DPSOFTRAST_Vector3Normalize(eyenormal);
3092
3093                         specularnormal[0] = lightnormal[0] + eyenormal[0];
3094                         specularnormal[1] = lightnormal[1] + eyenormal[1];
3095                         specularnormal[2] = lightnormal[2] + eyenormal[2];
3096                         DPSOFTRAST_Vector3Normalize(specularnormal);
3097
3098                         diffuse = DPSOFTRAST_Vector3Dot(surfacenormal, lightnormal);if (diffuse < 0.0f) diffuse = 0.0f;
3099                         specular = DPSOFTRAST_Vector3Dot(surfacenormal, specularnormal);if (specular < 0.0f) specular = 0.0f;
3100                         specular = pow(specular, SpecularPower * glosstex[3]);
3101                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_CUBEFILTER)
3102                         {
3103                                 // scale down the attenuation to account for the cubefilter multiplying everything by 255
3104                                 attenuation *= (1.0f / 255.0f);
3105                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse) + glosstex[0] * Color_Specular[0] * specular) * LightColor[0] * buffer_texture_cubebgra8[x*4+0] * attenuation);if (d[0] > 255) d[0] = 255;
3106                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse) + glosstex[1] * Color_Specular[1] * specular) * LightColor[1] * buffer_texture_cubebgra8[x*4+1] * attenuation);if (d[1] > 255) d[1] = 255;
3107                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse) + glosstex[2] * Color_Specular[2] * specular) * LightColor[2] * buffer_texture_cubebgra8[x*4+2] * attenuation);if (d[2] > 255) d[2] = 255;
3108                                 d[3] = (int)( diffusetex[3]                                                                                                                                                                );if (d[3] > 255) d[3] = 255;
3109                         }
3110                         else
3111                         {
3112                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse) + glosstex[0] * Color_Specular[0] * specular) * LightColor[0]                                   * attenuation);if (d[0] > 255) d[0] = 255;
3113                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse) + glosstex[1] * Color_Specular[1] * specular) * LightColor[1]                                   * attenuation);if (d[1] > 255) d[1] = 255;
3114                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse) + glosstex[2] * Color_Specular[2] * specular) * LightColor[2]                                   * attenuation);if (d[2] > 255) d[2] = 255;
3115                                 d[3] = (int)( diffusetex[3]                                                                                                                                                                );if (d[3] > 255) d[3] = 255;
3116                         }
3117                         buffer_FragColorbgra8[x*4+0] = d[0];
3118                         buffer_FragColorbgra8[x*4+1] = d[1];
3119                         buffer_FragColorbgra8[x*4+2] = d[2];
3120                         buffer_FragColorbgra8[x*4+3] = d[3];
3121                 }
3122         }
3123         else if (dpsoftrast.shader_permutation & SHADERPERMUTATION_DIFFUSE)
3124         {
3125                 DPSOFTRAST_Draw_Span_Texture2DVaryingBGRA8(span, buffer_texture_normalbgra8, GL20TU_NORMAL, DPSOFTRAST_ARRAY_TEXCOORD0, buffer_z);
3126                 for (x = startx;x < endx;x++)
3127                 {
3128                         z = buffer_z[x];
3129                         CubeVector[0] = (CubeVectordata[0] + CubeVectorslope[0]*x) * z;
3130                         CubeVector[1] = (CubeVectordata[1] + CubeVectorslope[1]*x) * z;
3131                         CubeVector[2] = (CubeVectordata[2] + CubeVectorslope[2]*x) * z;
3132                         attenuation = 1.0f - DPSOFTRAST_Vector3LengthSquared(CubeVector);
3133                         if (attenuation < 0.01f)
3134                                 continue;
3135                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SHADOWMAP2D)
3136                         {
3137                                 attenuation *= DPSOFTRAST_SampleShadowmap(CubeVector);
3138                                 if (attenuation < 0.01f)
3139                                         continue;
3140                         }
3141
3142                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
3143                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
3144                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
3145                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
3146                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
3147                         {
3148                                 diffusetex[0] += buffer_texture_pantsbgra8[x*4+0] * Color_Pants[0] + buffer_texture_shirtbgra8[x*4+0] * Color_Shirt[0];
3149                                 diffusetex[1] += buffer_texture_pantsbgra8[x*4+1] * Color_Pants[1] + buffer_texture_shirtbgra8[x*4+1] * Color_Shirt[1];
3150                                 diffusetex[2] += buffer_texture_pantsbgra8[x*4+2] * Color_Pants[2] + buffer_texture_shirtbgra8[x*4+2] * Color_Shirt[2];
3151                                 diffusetex[3] += buffer_texture_pantsbgra8[x*4+3] * Color_Pants[3] + buffer_texture_shirtbgra8[x*4+3] * Color_Shirt[3];
3152                         }
3153                         surfacenormal[0] = buffer_texture_normalbgra8[x*4+2] * (1.0f / 128.0f) - 1.0f;
3154                         surfacenormal[1] = buffer_texture_normalbgra8[x*4+1] * (1.0f / 128.0f) - 1.0f;
3155                         surfacenormal[2] = buffer_texture_normalbgra8[x*4+0] * (1.0f / 128.0f) - 1.0f;
3156                         DPSOFTRAST_Vector3Normalize(surfacenormal);
3157
3158                         lightnormal[0] = (LightVectordata[0] + LightVectorslope[0]*x) * z;
3159                         lightnormal[1] = (LightVectordata[1] + LightVectorslope[1]*x) * z;
3160                         lightnormal[2] = (LightVectordata[2] + LightVectorslope[2]*x) * z;
3161                         DPSOFTRAST_Vector3Normalize(lightnormal);
3162
3163                         diffuse = DPSOFTRAST_Vector3Dot(surfacenormal, lightnormal);if (diffuse < 0.0f) diffuse = 0.0f;
3164                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_CUBEFILTER)
3165                         {
3166                                 // scale down the attenuation to account for the cubefilter multiplying everything by 255
3167                                 attenuation *= (1.0f / 255.0f);
3168                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse)) * LightColor[0] * buffer_texture_cubebgra8[x*4+0] * attenuation);if (d[0] > 255) d[0] = 255;
3169                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse)) * LightColor[1] * buffer_texture_cubebgra8[x*4+1] * attenuation);if (d[1] > 255) d[1] = 255;
3170                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse)) * LightColor[2] * buffer_texture_cubebgra8[x*4+2] * attenuation);if (d[2] > 255) d[2] = 255;
3171                                 d[3] = (int)( diffusetex[3]                                                                                                                   );if (d[3] > 255) d[3] = 255;
3172                         }
3173                         else
3174                         {
3175                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0] + Color_Diffuse[0] * diffuse)) * LightColor[0]                                   * attenuation);if (d[0] > 255) d[0] = 255;
3176                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1] + Color_Diffuse[1] * diffuse)) * LightColor[1]                                   * attenuation);if (d[1] > 255) d[1] = 255;
3177                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2] + Color_Diffuse[2] * diffuse)) * LightColor[2]                                   * attenuation);if (d[2] > 255) d[2] = 255;
3178                                 d[3] = (int)( diffusetex[3]                                                                                                                                                                );if (d[3] > 255) d[3] = 255;
3179                         }
3180                         buffer_FragColorbgra8[x*4+0] = d[0];
3181                         buffer_FragColorbgra8[x*4+1] = d[1];
3182                         buffer_FragColorbgra8[x*4+2] = d[2];
3183                         buffer_FragColorbgra8[x*4+3] = d[3];
3184                 }
3185         }
3186         else
3187         {
3188                 for (x = startx;x < endx;x++)
3189                 {
3190                         z = buffer_z[x];
3191                         CubeVector[0] = (CubeVectordata[0] + CubeVectorslope[0]*x) * z;
3192                         CubeVector[1] = (CubeVectordata[1] + CubeVectorslope[1]*x) * z;
3193                         CubeVector[2] = (CubeVectordata[2] + CubeVectorslope[2]*x) * z;
3194                         attenuation = 1.0f - DPSOFTRAST_Vector3LengthSquared(CubeVector);
3195                         if (attenuation < 0.01f)
3196                                 continue;
3197                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_SHADOWMAP2D)
3198                         {
3199                                 attenuation *= DPSOFTRAST_SampleShadowmap(CubeVector);
3200                                 if (attenuation < 0.01f)
3201                                         continue;
3202                         }
3203
3204                         diffusetex[0] = buffer_texture_colorbgra8[x*4+0];
3205                         diffusetex[1] = buffer_texture_colorbgra8[x*4+1];
3206                         diffusetex[2] = buffer_texture_colorbgra8[x*4+2];
3207                         diffusetex[3] = buffer_texture_colorbgra8[x*4+3];
3208                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_COLORMAPPING)
3209                         {
3210                                 diffusetex[0] += buffer_texture_pantsbgra8[x*4+0] * Color_Pants[0] + buffer_texture_shirtbgra8[x*4+0] * Color_Shirt[0];
3211                                 diffusetex[1] += buffer_texture_pantsbgra8[x*4+1] * Color_Pants[1] + buffer_texture_shirtbgra8[x*4+1] * Color_Shirt[1];
3212                                 diffusetex[2] += buffer_texture_pantsbgra8[x*4+2] * Color_Pants[2] + buffer_texture_shirtbgra8[x*4+2] * Color_Shirt[2];
3213                                 diffusetex[3] += buffer_texture_pantsbgra8[x*4+3] * Color_Pants[3] + buffer_texture_shirtbgra8[x*4+3] * Color_Shirt[3];
3214                         }
3215                         if (dpsoftrast.shader_permutation & SHADERPERMUTATION_CUBEFILTER)
3216                         {
3217                                 // scale down the attenuation to account for the cubefilter multiplying everything by 255
3218                                 attenuation *= (1.0f / 255.0f);
3219                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0])) * LightColor[0] * buffer_texture_cubebgra8[x*4+0] * attenuation);if (d[0] > 255) d[0] = 255;
3220                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1])) * LightColor[1] * buffer_texture_cubebgra8[x*4+1] * attenuation);if (d[1] > 255) d[1] = 255;
3221                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2])) * LightColor[2] * buffer_texture_cubebgra8[x*4+2] * attenuation);if (d[2] > 255) d[2] = 255;
3222                                 d[3] = (int)( diffusetex[3]                                                                                      );if (d[3] > 255) d[3] = 255;
3223                         }
3224                         else
3225                         {
3226                                 d[0] = (int)((diffusetex[0] * (Color_Ambient[0])) * LightColor[0]                                   * attenuation);if (d[0] > 255) d[0] = 255;
3227                                 d[1] = (int)((diffusetex[1] * (Color_Ambient[1])) * LightColor[1]                                   * attenuation);if (d[1] > 255) d[1] = 255;
3228                                 d[2] = (int)((diffusetex[2] * (Color_Ambient[2])) * LightColor[2]                                   * attenuation);if (d[2] > 255) d[2] = 255;
3229                                 d[3] = (int)( diffusetex[3]                                                                                                                                                                );if (d[3] > 255) d[3] = 255;
3230                         }
3231                         buffer_FragColorbgra8[x*4+0] = d[0];
3232                         buffer_FragColorbgra8[x*4+1] = d[1];
3233                         buffer_FragColorbgra8[x*4+2] = d[2];
3234                         buffer_FragColorbgra8[x*4+3] = d[3];
3235                 }
3236         }
3237         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3238 }
3239
3240
3241
3242 void DPSOFTRAST_VertexShader_Refraction(void)
3243 {
3244         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
3245 }
3246
3247 void DPSOFTRAST_PixelShader_Refraction(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
3248 {
3249         // TODO: IMPLEMENT
3250         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3251         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
3252         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3253         memset(buffer_FragColorbgra8, 0, span->length*4);
3254         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3255 }
3256
3257
3258
3259 void DPSOFTRAST_VertexShader_Water(void)
3260 {
3261         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
3262 }
3263
3264
3265 void DPSOFTRAST_PixelShader_Water(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
3266 {
3267         // TODO: IMPLEMENT
3268         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3269         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
3270         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3271         memset(buffer_FragColorbgra8, 0, span->length*4);
3272         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3273 }
3274
3275
3276
3277 void DPSOFTRAST_VertexShader_ShowDepth(void)
3278 {
3279         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
3280 }
3281
3282 void DPSOFTRAST_PixelShader_ShowDepth(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
3283 {
3284         // TODO: IMPLEMENT
3285         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3286         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
3287         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3288         memset(buffer_FragColorbgra8, 0, span->length*4);
3289         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3290 }
3291
3292
3293
3294 void DPSOFTRAST_VertexShader_DeferredGeometry(void)
3295 {
3296         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
3297 }
3298
3299 void DPSOFTRAST_PixelShader_DeferredGeometry(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
3300 {
3301         // TODO: IMPLEMENT
3302         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3303         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
3304         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3305         memset(buffer_FragColorbgra8, 0, span->length*4);
3306         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3307 }
3308
3309
3310
3311 void DPSOFTRAST_VertexShader_DeferredLightSource(void)
3312 {
3313         DPSOFTRAST_Array_Transform(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.in_array4f[DPSOFTRAST_ARRAY_POSITION], dpsoftrast.draw.numvertices, dpsoftrast.uniform4f + 4*DPSOFTRAST_UNIFORM_ModelViewProjectionMatrixM1);
3314 }
3315
3316 void DPSOFTRAST_PixelShader_DeferredLightSource(const DPSOFTRAST_State_Draw_Span * RESTRICT span)
3317 {
3318         // TODO: IMPLEMENT
3319         float buffer_z[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3320         unsigned char buffer_FragColorbgra8[DPSOFTRAST_DRAW_MAXSPANLENGTH*4];
3321         DPSOFTRAST_Draw_Span_Begin(span, buffer_z);
3322         memset(buffer_FragColorbgra8, 0, span->length*4);
3323         DPSOFTRAST_Draw_Span_FinishBGRA8(span, buffer_FragColorbgra8);
3324 }
3325
3326
3327
3328 typedef struct DPSOFTRAST_ShaderModeInfo_s
3329 {
3330         int lodarrayindex;
3331         void (*Vertex)(void);
3332         void (*Span)(const DPSOFTRAST_State_Draw_Span * RESTRICT span);
3333 }
3334 DPSOFTRAST_ShaderModeInfo;
3335
3336 DPSOFTRAST_ShaderModeInfo DPSOFTRAST_ShaderModeTable[SHADERMODE_COUNT] =
3337 {
3338         {2, DPSOFTRAST_VertexShader_Generic,                        DPSOFTRAST_PixelShader_Generic,                      },
3339         {2, DPSOFTRAST_VertexShader_PostProcess,                    DPSOFTRAST_PixelShader_PostProcess,                  },
3340         {2, DPSOFTRAST_VertexShader_Depth_Or_Shadow,                DPSOFTRAST_PixelShader_Depth_Or_Shadow,              },
3341         {2, DPSOFTRAST_VertexShader_FlatColor,                      DPSOFTRAST_PixelShader_FlatColor,                    },
3342         {2, DPSOFTRAST_VertexShader_VertexColor,                    DPSOFTRAST_PixelShader_VertexColor,                  },
3343         {2, DPSOFTRAST_VertexShader_Lightmap,                       DPSOFTRAST_PixelShader_Lightmap,                     },
3344         {2, DPSOFTRAST_VertexShader_FakeLight,                      DPSOFTRAST_PixelShader_FakeLight,                    },
3345         {2, DPSOFTRAST_VertexShader_LightDirectionMap_ModelSpace,   DPSOFTRAST_PixelShader_LightDirectionMap_ModelSpace, },
3346         {2, DPSOFTRAST_VertexShader_LightDirectionMap_TangentSpace, DPSOFTRAST_PixelShader_LightDirectionMap_TangentSpace},
3347         {2, DPSOFTRAST_VertexShader_LightDirection,                 DPSOFTRAST_PixelShader_LightDirection,               },
3348         {2, DPSOFTRAST_VertexShader_LightSource,                    DPSOFTRAST_PixelShader_LightSource,                  },
3349         {2, DPSOFTRAST_VertexShader_Refraction,                     DPSOFTRAST_PixelShader_Refraction,                   },
3350         {2, DPSOFTRAST_VertexShader_Water,                          DPSOFTRAST_PixelShader_Water,                        },
3351         {2, DPSOFTRAST_VertexShader_ShowDepth,                      DPSOFTRAST_PixelShader_ShowDepth,                    },
3352         {2, DPSOFTRAST_VertexShader_DeferredGeometry,               DPSOFTRAST_PixelShader_DeferredGeometry,             },
3353         {2, DPSOFTRAST_VertexShader_DeferredLightSource,            DPSOFTRAST_PixelShader_DeferredLightSource,          }
3354 };
3355
3356
3357
3358 void DPSOFTRAST_Draw_ProcessSpans(void)
3359 {
3360         int i;
3361         int x;
3362         int startx;
3363         int endx;
3364         int numspans = dpsoftrast.draw.numspans;
3365 //      unsigned int c;
3366 //      unsigned int *colorpixel;
3367         unsigned int *depthpixel;
3368         float w;
3369         float wslope;
3370         int depth;
3371         int depthslope;
3372         unsigned int d;
3373         DPSOFTRAST_State_Draw_Span *span = dpsoftrast.draw.spanqueue;
3374         unsigned char pixelmask[DPSOFTRAST_DRAW_MAXSPANLENGTH];
3375         for (i = 0;i < numspans;i++, span++)
3376         {
3377                 w = span->data[0][DPSOFTRAST_ARRAY_TOTAL][3];
3378                 wslope = span->data[1][DPSOFTRAST_ARRAY_TOTAL][3];
3379                 if (dpsoftrast.user.depthtest && dpsoftrast.fb_depthpixels)
3380                 {
3381                         depthslope = (int)(wslope*DPSOFTRAST_DEPTHSCALE);
3382                         depth = (int)(w*DPSOFTRAST_DEPTHSCALE - DPSOFTRAST_DEPTHOFFSET*(dpsoftrast.user.polygonoffset[1] + fabs(wslope)*dpsoftrast.user.polygonoffset[0]));
3383                         depthpixel = dpsoftrast.fb_depthpixels + span->start;
3384                         switch(dpsoftrast.fb_depthfunc)
3385                         {
3386                         default:
3387                         case GL_ALWAYS:  for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = true; break;
3388                         case GL_LESS:    for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = depthpixel[x] < d; break;
3389                         case GL_LEQUAL:  for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = depthpixel[x] <= d; break;
3390                         case GL_EQUAL:   for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = depthpixel[x] == d; break;
3391                         case GL_GEQUAL:  for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = depthpixel[x] >= d; break;
3392                         case GL_GREATER: for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = depthpixel[x] > d; break;
3393                         case GL_NEVER:   for (x = 0, d = depth;x < span->length;x++, d += depthslope) pixelmask[x] = false; break;
3394                         }
3395                         //colorpixel = dpsoftrast.fb_colorpixels[0] + span->start;
3396                         //for (x = 0;x < span->length;x++)
3397                         //      colorpixel[x] = (depthpixel[x] & 0xFF000000) ? (0x00FF0000) : (depthpixel[x] & 0x00FF0000);
3398                         // if there is no color buffer, skip pixel shader
3399                         startx = 0;
3400                         endx = span->length;
3401                         while (startx < endx && !pixelmask[startx])
3402                                 startx++;
3403                         while (endx > startx && !pixelmask[endx-1])
3404                                 endx--;
3405                         if (startx >= endx)
3406                                 continue; // no pixels to fill
3407                         span->pixelmask = pixelmask;
3408                         span->startx = startx;
3409                         span->endx = endx;
3410                         // run pixel shader if appropriate
3411                         // do this before running depthmask code, to allow the pixelshader
3412                         // to clear pixelmask values for alpha testing
3413                         if (dpsoftrast.fb_colorpixels[0] && dpsoftrast.fb_colormask)
3414                                 DPSOFTRAST_ShaderModeTable[dpsoftrast.shader_mode].Span(span);
3415                         if (dpsoftrast.user.depthmask)
3416                                 for (x = 0, d = depth;x < span->length;x++, d += depthslope)
3417                                         if (pixelmask[x])
3418                                                 depthpixel[x] = d;
3419                 }
3420                 else
3421                 {
3422                         // no depth testing means we're just dealing with color...
3423                         // if there is no color buffer, skip pixel shader
3424                         if (dpsoftrast.fb_colorpixels[0] && dpsoftrast.fb_colormask)
3425                         {
3426                                 memset(pixelmask, 1, span->length);
3427                                 span->pixelmask = pixelmask;
3428                                 span->startx = 0;
3429                                 span->endx = span->length;
3430                                 DPSOFTRAST_ShaderModeTable[dpsoftrast.shader_mode].Span(span);
3431                         }
3432                 }
3433         }
3434 }
3435
3436 void DPSOFTRAST_Draw_ProcessTriangles(int firstvertex, int numtriangles, const int *element3i, const unsigned short *element3s, unsigned char *arraymask)
3437 {
3438 #ifdef SSE2_PRESENT
3439         int cullface = dpsoftrast.user.cullface;
3440         int width = dpsoftrast.fb_width;
3441         int height = dpsoftrast.fb_height;
3442         __m128i fbmax = _mm_sub_epi16(_mm_setr_epi16(width, height, width, height, width, height, width, height), _mm_set1_epi16(1));
3443         int i;
3444         int j;
3445         int k;
3446         int y;
3447         int e[3];
3448         __m128i screeny;
3449         int starty, endy;
3450         int numpoints;
3451         int edge0p;
3452         int edge0n;
3453         int edge1p;
3454         int edge1n;
3455         int startx;
3456         int endx;
3457         unsigned char mip[DPSOFTRAST_MAXTEXTUREUNITS];
3458         __m128 mipedgescale;
3459         float clipdist[4];
3460         __m128 clipped[DPSOFTRAST_ARRAY_TOTAL][4];
3461         __m128 screen[4];
3462         __m128 proj[DPSOFTRAST_ARRAY_TOTAL][4];
3463         DPSOFTRAST_Texture *texture;
3464         DPSOFTRAST_State_Draw_Span *span;
3465         DPSOFTRAST_State_Draw_Span *oldspan;
3466         for (i = 0;i < numtriangles;i++)
3467         {
3468                 // generate the 3 edges of this triangle
3469                 // generate spans for the triangle - switch based on left split or right split classification of triangle
3470                 if (element3i)
3471                 {
3472                         e[0] = element3i[i*3+0] - firstvertex;
3473                         e[1] = element3i[i*3+1] - firstvertex;
3474                         e[2] = element3i[i*3+2] - firstvertex;
3475                 }
3476                 else if (element3s)
3477                 {
3478                         e[0] = element3s[i*3+0] - firstvertex;
3479                         e[1] = element3s[i*3+1] - firstvertex;
3480                         e[2] = element3s[i*3+2] - firstvertex;
3481                 }
3482                 else
3483                 {
3484                         e[0] = i*3+0;
3485                         e[1] = i*3+1;
3486                         e[2] = i*3+2;
3487                 }
3488
3489 #define SKIPBACKFACE { \
3490                         __m128 triangleedge[2], triangleorigin, trianglenormal; \
3491                         triangleorigin = _mm_load_ps(&dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[1]*4]); \
3492                         triangleedge[0] = _mm_sub_ps(_mm_load_ps(&dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[0]*4]), triangleorigin); \
3493                         triangleedge[1] = _mm_sub_ps(_mm_load_ps(&dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[2]*4]), triangleorigin); \
3494                         /* store normal in 2, 0, 1 order instead of 0, 1, 2 as it requires fewer shuffles and leaves z component accessible as scalar */ \
3495                         trianglenormal = _mm_sub_ss(_mm_mul_ss(triangleedge[0], _mm_shuffle_ps(triangleedge[1], triangleedge[1], _MM_SHUFFLE(3, 0, 2, 1))), \
3496                                                                                 _mm_mul_ss(_mm_shuffle_ps(triangleedge[0], triangleedge[0], _MM_SHUFFLE(3, 0, 2, 1)), triangleedge[1])); \
3497                         /* apply current cullface mode (this culls many triangles) */ \
3498                         switch(cullface) \
3499                         { \
3500                         case GL_BACK: \
3501                                 if (_mm_ucomilt_ss(trianglenormal, _mm_setzero_ps())) \
3502                                         continue; \
3503                                 break; \
3504                         case GL_FRONT: \
3505                                 if (_mm_ucomigt_ss(trianglenormal, _mm_setzero_ps())) \
3506                                         continue; \
3507                                 break; \
3508                         } \
3509                 }
3510                         //trianglenormal = _mm_sub_ps(_mm_mul_ps(triangleedge[0], _mm_shuffle_ps(triangleedge[1], triangleedge[1], _MM_SHUFFLE(3, 0, 2, 1))),
3511                         //                                                _mm_mul_ps(_mm_shuffle_ps(triangleedge[0], triangleedge[0], _MM_SHUFFLE(3, 0, 2, 1)), triangleedge[1]));
3512                         //trianglenormal[2] = triangleedge[0][0] * triangleedge[1][1] - triangleedge[0][1] * triangleedge[1][0];
3513                         //trianglenormal[0] = triangleedge[0][1] * triangleedge[1][2] - triangleedge[0][2] * triangleedge[1][1];
3514                         //trianglenormal[1] = triangleedge[0][2] * triangleedge[1][0] - triangleedge[0][0] * triangleedge[1][2];
3515
3516                         // macros for clipping vertices
3517 #define CLIPPEDVERTEXLERP(k,p1,p2) { \
3518                         __m128 frac = _mm_set1_ps(clipdist[p1] / (clipdist[p1] - clipdist[p2]));\
3519                         for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL;j++)\
3520                         {\
3521                                 /*if (arraymask[j])*/\
3522                                 {\
3523                                         __m128 v1 = _mm_load_ps(&dpsoftrast.draw.post_array4f[j][e[p1]*4]), v2 = _mm_load_ps(&dpsoftrast.draw.post_array4f[j][e[p2]*4]); \
3524                                         clipped[j][k] = _mm_add_ps(v1, _mm_mul_ps(_mm_sub_ps(v2, v1), frac)); \
3525                                 }\
3526                         }\
3527                         screen[k] = DPSOFTRAST_Draw_ProjectVertex(clipped[DPSOFTRAST_ARRAY_POSITION][k]); \
3528                 }
3529 #define CLIPPEDVERTEXCOPY(k,p1) \
3530                         for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL;j++)\
3531                         {\
3532                                 /*if (arraymask[j])*/\
3533                                 {\
3534                                         clipped[j][k] = _mm_load_ps(&dpsoftrast.draw.post_array4f[j][e[p1]*4]); \
3535                                 }\
3536                         }\
3537                         screen[k] = _mm_load_ps(&dpsoftrast.draw.screencoord4f[e[p1]*4]);
3538
3539                 // calculate distance from nearplane
3540                 clipdist[0] = dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[0]*4+2] + 1.0f;
3541                 clipdist[1] = dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[1]*4+2] + 1.0f;
3542                 clipdist[2] = dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][e[2]*4+2] + 1.0f;
3543                 if (clipdist[0] >= 0.0f)
3544                 {
3545                         SKIPBACKFACE;
3546                         if (clipdist[1] >= 0.0f)
3547                         {
3548                                 CLIPPEDVERTEXCOPY(0,0);
3549                                 CLIPPEDVERTEXCOPY(1,1);
3550                                 if (clipdist[2] >= 0.0f)
3551                                 {
3552                                         // triangle is entirely in front of nearplane
3553                                         CLIPPEDVERTEXCOPY(2,2);
3554                                         numpoints = 3;
3555                                 }
3556                                 else
3557                                 {
3558                                         CLIPPEDVERTEXLERP(2,1,2);
3559                                         CLIPPEDVERTEXLERP(3,2,0);
3560                                         numpoints = 4;
3561                                 }
3562                         }
3563                         else
3564                         {
3565                                 CLIPPEDVERTEXCOPY(0,0);
3566                                 CLIPPEDVERTEXLERP(1,0,1);
3567                                 if (clipdist[2] >= 0.0f)
3568                                 {
3569                                         CLIPPEDVERTEXLERP(2,1,2);
3570                                         CLIPPEDVERTEXCOPY(3,2);
3571                                         numpoints = 4;
3572                                 }
3573                                 else
3574                                 {
3575                                         CLIPPEDVERTEXLERP(2,2,0);
3576                                         numpoints = 3;
3577                                 }
3578                         }
3579                 }
3580                 else if (clipdist[1] >= 0.0f)
3581                 {
3582                         SKIPBACKFACE;
3583                         CLIPPEDVERTEXLERP(0,0,1);
3584                         CLIPPEDVERTEXCOPY(1,1);
3585                         if (clipdist[2] >= 0.0f)
3586                         {
3587                                 CLIPPEDVERTEXCOPY(2,2);
3588                                 CLIPPEDVERTEXLERP(3,2,0);
3589                                 numpoints = 4;
3590                         }
3591                         else
3592                         {
3593                                 CLIPPEDVERTEXLERP(2,1,2);
3594                                 numpoints = 3;
3595                         }
3596                 }
3597                 else if (clipdist[2] >= 0.0f)
3598                 {
3599                         SKIPBACKFACE;
3600                         CLIPPEDVERTEXLERP(0,1,2);
3601                         CLIPPEDVERTEXCOPY(1,2);
3602                         CLIPPEDVERTEXLERP(2,2,0);
3603                         numpoints = 3;
3604                 }
3605                 else continue; // triangle is entirely behind nearpla
3606
3607                 {
3608                         // calculate integer y coords for triangle points
3609                         __m128i screeni = _mm_packs_epi32(_mm_cvttps_epi32(_mm_shuffle_ps(screen[0], screen[1], _MM_SHUFFLE(1, 0, 1, 0))),
3610                                                                                   _mm_cvttps_epi32(_mm_shuffle_ps(screen[2], numpoints <= 3 ? screen[2] : screen[3], _MM_SHUFFLE(1, 0, 1, 0)))),
3611                                         screenir = _mm_shuffle_epi32(screeni, _MM_SHUFFLE(1, 0, 3, 2)),
3612                                         screenmin = _mm_min_epi16(screeni, screenir),
3613                                         screenmax = _mm_max_epi16(screeni, screenir);
3614                         screenmin = _mm_min_epi16(screenmin, _mm_shufflelo_epi16(screenmin, _MM_SHUFFLE(1, 0, 3, 2)));
3615                         screenmax = _mm_max_epi16(screenmax, _mm_shufflelo_epi16(screenmax, _MM_SHUFFLE(1, 0, 3, 2)));
3616                         screenmin = _mm_max_epi16(screenmin, _mm_setzero_si128());
3617                         screenmax = _mm_min_epi16(screenmax, fbmax);
3618                         // skip offscreen triangles
3619                         if (_mm_cvtsi128_si32(_mm_cmplt_epi16(screenmax, screenmin)) == -1)
3620                                 continue;
3621                         starty = _mm_extract_epi16(screenmin, 1);
3622                         endy = _mm_extract_epi16(screenmax, 1)+1;
3623                         screeny = _mm_srai_epi32(screeni, 16);
3624                 }
3625
3626                 // okay, this triangle is going to produce spans, we'd better project
3627                 // the interpolants now (this is what gives perspective texturing),
3628                 // this consists of simply multiplying all arrays by the W coord
3629                 // (which is basically 1/Z), which will be undone per-pixel
3630                 // (multiplying by Z again) to get the perspective-correct array
3631                 // values
3632                 for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL;j++)
3633                 {
3634                         //if (arraymask[j])
3635                         {
3636                                 for (k = 0;k < numpoints;k++)
3637                                 {
3638                                         proj[j][k] = _mm_mul_ps(clipped[j][k], _mm_shuffle_ps(screen[k], screen[k], _MM_SHUFFLE(3, 3, 3, 3)));
3639                                 }
3640                         }
3641                 }
3642                 // adjust texture LOD by texture density, in the simplest way possible...
3643                 mipedgescale = _mm_sub_ps(_mm_shuffle_ps(screen[0], screen[2], _MM_SHUFFLE(1, 0, 1, 0)), _mm_shuffle_ps(screen[1], screen[1], _MM_SHUFFLE(1, 0, 1, 0)));
3644                 mipedgescale = _mm_mul_ps(mipedgescale, mipedgescale);
3645                 mipedgescale = _mm_rcp_ps(_mm_add_ps(mipedgescale, _mm_shuffle_ps(mipedgescale, mipedgescale, _MM_SHUFFLE(2, 3, 0, 1))));
3646                 for (j = 0;j < DPSOFTRAST_MAXTEXTUREUNITS;j++)
3647                 {
3648                         texture = dpsoftrast.texbound[j];
3649                         if (texture)
3650                         {
3651                                 __m128 mipedgetc;
3652                                 if (texture->filter <= DPSOFTRAST_TEXTURE_FILTER_LINEAR)
3653                                 {
3654                                         mip[j] = 0;
3655                                         continue;
3656                                 }
3657                                 k = DPSOFTRAST_ShaderModeTable[dpsoftrast.shader_mode].lodarrayindex;
3658                                 mipedgetc = _mm_sub_ps(_mm_shuffle_ps(clipped[k][0], clipped[k][2], _MM_SHUFFLE(1, 0, 1, 0)),
3659                                                                                 _mm_shuffle_ps(clipped[k][1], clipped[k][1], _MM_SHUFFLE(1, 0, 1, 0)));
3660                                 mipedgetc = _mm_mul_ps(mipedgetc, _mm_cvtepi32_ps(_mm_shuffle_epi32(_mm_loadl_epi64((const __m128i *)&texture->mipmap[0][2]), _MM_SHUFFLE(1, 0, 1, 0))));
3661                                 mipedgetc = _mm_mul_ps(mipedgetc, mipedgetc);
3662                                 mipedgetc = _mm_add_ps(mipedgetc, _mm_shuffle_ps(mipedgetc, mipedgetc, _MM_SHUFFLE(2, 3, 0, 1)));
3663                                 mipedgetc = _mm_mul_ps(mipedgetc, mipedgescale);
3664                                 mipedgetc = _mm_min_ss(mipedgetc, _mm_shuffle_ps(mipedgetc, mipedgetc, _MM_SHUFFLE(2, 2, 2, 2)));
3665                                 // this will be multiplied in the texturing routine by the texture resolution
3666                                 y = _mm_cvttss_si32(mipedgetc);
3667                                 if (y > 0)
3668                                 {
3669                                         y = (int)(log(y)/M_LN2);
3670                                         if (y > texture->mipmaps - 1)
3671                                                 y = texture->mipmaps - 1;
3672                                 }
3673                                 else y = 0;
3674                                 mip[j] = y;
3675                         }
3676                 }
3677                 // iterate potential spans
3678                 // TODO: optimize?  if we figured out the edge order beforehand, this
3679                 //         could do loops over the edges in the proper order rather than
3680                 //         selecting them for each span
3681                 // TODO: optimize?  the edges could have data slopes calculated
3682                 // TODO: optimize?  the data slopes could be calculated as a plane
3683                 //         (2D slopes) to avoid any interpolation along edges at all
3684                 for (y = starty+1;y < endy;)
3685                 {
3686                         int nexty = -1;
3687                         __m128 edge0lerp, edge1lerp, edge0scale, edge1scale;
3688                         __m128i screenycc = _mm_cmpgt_epi32(_mm_set1_epi32(y), screeny);
3689                         int screenymask = _mm_movemask_epi8(screenycc);
3690                         if (numpoints == 4)
3691                         {
3692                                 switch(screenymask)
3693                                 {
3694                                 default:
3695                                 case 0xFFFF: /*0000*/ y++; continue;
3696                                 case 0xFFF0: /*1000*/ edge0p = 3;edge0n = 0;edge1p = 1;edge1n = 0;break;
3697                                 case 0xFF0F: /*0100*/ edge0p = 0;edge0n = 1;edge1p = 2;edge1n = 1;break;
3698                                 case 0xFF00: /*1100*/ edge0p = 3;edge0n = 0;edge1p = 2;edge1n = 1;break;
3699                                 case 0xF0FF: /*0010*/ edge0p = 1;edge0n = 2;edge1p = 3;edge1n = 2;break;
3700                                 case 0xF0F0: /*1010*/ edge0p = 1;edge0n = 2;edge1p = 3;edge1n = 2;break; // concave - nonsense
3701                                 case 0xF00F: /*0110*/ edge0p = 0;edge0n = 1;edge1p = 3;edge1n = 2;break;
3702                                 case 0xF000: /*1110*/ edge0p = 3;edge0n = 0;edge1p = 3;edge1n = 2;break;
3703                                 case 0x0FFF: /*0001*/ edge0p = 2;edge0n = 3;edge1p = 0;edge1n = 3;break;
3704                                 case 0x0FF0: /*1001*/ edge0p = 2;edge0n = 3;edge1p = 1;edge1n = 0;break;
3705                                 case 0x0F0F: /*0101*/ edge0p = 2;edge0n = 3;edge1p = 1;edge1n = 2;break; // concave - nonsense
3706                                 case 0x0F00: /*1101*/ edge0p = 2;edge0n = 3;edge1p = 2;edge1n = 1;break;
3707                                 case 0x00FF: /*0011*/ edge0p = 1;edge0n = 2;edge1p = 0;edge1n = 3;break;
3708                                 case 0x00F0: /*1011*/ edge0p = 1;edge0n = 2;edge1p = 1;edge1n = 0;break;
3709                                 case 0x000F: /*0111*/ edge0p = 0;edge0n = 1;edge1p = 0;edge1n = 3;break;
3710                                 case 0x0000: /*1111*/ y++; continue;
3711                                 }
3712                         }
3713                         else
3714                         {
3715                                 switch(screenymask)
3716                                 {
3717                                 default:
3718                                 case 0xFFFF: /*000*/ y++; continue;
3719                                 case 0xFFF0: /*100*/ edge0p = 2;edge0n = 0;edge1p = 1;edge1n = 0;break;
3720                                 case 0xFF0F: /*010*/ edge0p = 0;edge0n = 1;edge1p = 2;edge1n = 1;break;
3721                                 case 0xFF00: /*110*/ edge0p = 2;edge0n = 0;edge1p = 2;edge1n = 1;break;
3722                                 case 0x00FF: /*001*/ edge0p = 1;edge0n = 2;edge1p = 0;edge1n = 2;break;
3723                                 case 0x00F0: /*101*/ edge0p = 1;edge0n = 2;edge1p = 1;edge1n = 0;break;
3724                                 case 0x000F: /*011*/ edge0p = 0;edge0n = 1;edge1p = 0;edge1n = 2;break;
3725                                 case 0x0000: /*111*/ y++; continue;
3726                                 }
3727                         }
3728                         screenycc = _mm_max_epi16(_mm_srli_epi16(screenycc, 1), screeny);
3729                         screenycc = _mm_min_epi16(screenycc, _mm_shuffle_epi32(screenycc, _MM_SHUFFLE(1, 0, 3, 2)));
3730                         screenycc = _mm_min_epi16(screenycc, _mm_shuffle_epi32(screenycc, _MM_SHUFFLE(2, 3, 0, 1)));
3731                         nexty = _mm_extract_epi16(screenycc, 0);
3732                         if(nexty >= endy) nexty = endy-1;
3733                         if (_mm_ucomigt_ss(_mm_max_ps(screen[edge0n], screen[edge0p]), _mm_min_ps(screen[edge1n], screen[edge1p])))
3734                         {
3735                                 int tmp = edge0n;
3736                                 edge0n = edge1n;
3737                                 edge1n = tmp;
3738                                 tmp = edge0p;
3739                                 edge0p = edge1p;
3740                                 edge1p = tmp;
3741                         }
3742                         edge0lerp = _mm_shuffle_ps(screen[edge0p], screen[edge0p], _MM_SHUFFLE(1, 1, 1, 1));
3743                         edge0scale = _mm_rcp_ss(_mm_sub_ss(_mm_shuffle_ps(screen[edge0n], screen[edge0n], _MM_SHUFFLE(1, 1, 1, 1)), edge0lerp));
3744                         edge0scale = _mm_shuffle_ps(edge0scale, edge0scale, _MM_SHUFFLE(0, 0, 0, 0));
3745                         edge0lerp = _mm_mul_ps(_mm_sub_ps(_mm_set1_ps(y), edge0lerp), edge0scale);
3746                         edge1lerp = _mm_shuffle_ps(screen[edge1p], screen[edge1p], _MM_SHUFFLE(1, 1, 1, 1));
3747                         edge1scale = _mm_rcp_ss(_mm_sub_ss(_mm_shuffle_ps(screen[edge1n], screen[edge1n], _MM_SHUFFLE(1, 1, 1, 1)), edge1lerp));
3748                         edge1scale = _mm_shuffle_ps(edge1scale, edge1scale, _MM_SHUFFLE(0, 0, 0, 0));
3749                         edge1lerp = _mm_mul_ps(_mm_sub_ps(_mm_set1_ps(y), edge1lerp), edge1scale);
3750                         for(; y <= nexty; y++, edge0lerp = _mm_add_ps(edge0lerp, edge0scale), edge1lerp = _mm_add_ps(edge1lerp, edge1scale))
3751                         {
3752                                 __m128 data0, data1, spanilength, startxlerp;
3753 #if 0
3754                                 if (_mm_ucomilt_ss(edge0lerp, _mm_setzero_ps()) || _mm_ucomilt_ss(edge1lerp, _mm_setzero_ps()) ||
3755                                         _mm_ucomigt_ss(edge0lerp, _mm_set1_ps(1)) || _mm_ucomigt_ss(edge1lerp, _mm_set1_ps(1)))
3756                                         continue;
3757 #endif
3758                                 data0 = _mm_add_ps(_mm_mul_ps(_mm_sub_ps(screen[edge0n], screen[edge0p]), edge0lerp), screen[edge0p]);
3759                                 data1 = _mm_add_ps(_mm_mul_ps(_mm_sub_ps(screen[edge1n], screen[edge1p]), edge1lerp), screen[edge1p]);
3760                                 startx = _mm_cvtss_si32(_mm_add_ss(data0, _mm_set1_ps(0.5f)));
3761                                 endx = _mm_cvtss_si32(_mm_add_ss(data1, _mm_set1_ps(0.5f)));
3762                                 if (startx < 0) startx = 0;
3763                                 if (endx > width) endx = width;
3764                                 if (startx >= endx) continue;
3765 #if 0
3766                                 _mm_store_ss(&startxf, data0);
3767                                 _mm_store_ss(&endxf, data1);
3768                                 if (startxf > startx || endxf < endx-1) { printf("%s:%i X wrong (%i to %i is outside %f to %f)\n", __FILE__, __LINE__, startx, endx, startxf, endxf); }
3769 #endif
3770                                 spanilength = _mm_rcp_ss(_mm_sub_ss(data1, data0));
3771                                 spanilength = _mm_shuffle_ps(spanilength, spanilength, _MM_SHUFFLE(0, 0, 0, 0));
3772                                 startxlerp = _mm_sub_ps(_mm_set1_ps(startx), _mm_shuffle_ps(data0, data0, _MM_SHUFFLE(0, 0, 0, 0)));
3773                                 span = &dpsoftrast.draw.spanqueue[dpsoftrast.draw.numspans++];
3774                                 memcpy(span->mip, mip, sizeof(span->mip));
3775                                 span->start = y * width + startx;
3776                                 span->length = endx - startx;
3777                                 j = DPSOFTRAST_ARRAY_TOTAL;
3778                                 data1 = _mm_mul_ps(_mm_sub_ps(data1, data0), spanilength);
3779                                 data0 = _mm_add_ps(data0, _mm_mul_ps(data1, startxlerp));
3780                                 _mm_store_ps(span->data[0][j], data0);
3781                                 _mm_store_ps(span->data[1][j], data1);
3782                                 for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL;j++)
3783                                 {
3784                                         //if (arraymask[j])
3785                                         {
3786                                                 data0 = _mm_add_ps(_mm_mul_ps(_mm_sub_ps(proj[j][edge0n], proj[j][edge0p]), edge0lerp), proj[j][edge0p]);
3787                                                 data1 = _mm_add_ps(_mm_mul_ps(_mm_sub_ps(proj[j][edge1n], proj[j][edge1p]), edge1lerp), proj[j][edge1p]);
3788                                                 data1 = _mm_mul_ps(_mm_sub_ps(data1, data0), spanilength);
3789                                                 data0 = _mm_add_ps(data0, _mm_mul_ps(data1, startxlerp));
3790                                                 _mm_store_ps(span->data[0][j], data0);
3791                                                 _mm_store_ps(span->data[1][j], data1);
3792                                         }
3793                                 }
3794                                 // to keep the shader routines from needing more than a small
3795                                 // buffer for pixel intermediate data, we split long spans...
3796                                 while (span->length > DPSOFTRAST_DRAW_MAXSPANLENGTH)
3797                                 {
3798                                         span->length = DPSOFTRAST_DRAW_MAXSPANLENGTH;
3799                                         if (dpsoftrast.draw.numspans >= DPSOFTRAST_DRAW_MAXSPANQUEUE)
3800                                         {
3801                                                 DPSOFTRAST_Draw_ProcessSpans();
3802                                                 dpsoftrast.draw.numspans = 0;
3803                                         }
3804                                         oldspan = span;
3805                                         span = &dpsoftrast.draw.spanqueue[dpsoftrast.draw.numspans++];
3806                                         *span = *oldspan;
3807                                         startx += DPSOFTRAST_DRAW_MAXSPANLENGTH;
3808                                         span->start = y * width + startx;
3809                                         span->length = endx - startx;
3810                                         j = DPSOFTRAST_ARRAY_TOTAL;
3811                                         _mm_store_ps(span->data[0][j], _mm_add_ps(_mm_load_ps(span->data[0][j]), _mm_mul_ps(_mm_load_ps(span->data[1][j]), _mm_set1_ps(DPSOFTRAST_DRAW_MAXSPANLENGTH))));
3812                                         for (j = 0;j < DPSOFTRAST_ARRAY_TOTAL;j++)
3813                                         {
3814                                                 //if (arraymask[j])
3815                                                 {
3816                                                          _mm_store_ps(span->data[0][j], _mm_add_ps(_mm_load_ps(span->data[0][j]), _mm_mul_ps(_mm_load_ps(span->data[1][j]), _mm_set1_ps(DPSOFTRAST_DRAW_MAXSPANLENGTH))));
3817                                                 }
3818                                         }
3819                                 }
3820                                 // after all that, we have a span suitable for the pixel shader...
3821                                 if (dpsoftrast.draw.numspans >= DPSOFTRAST_DRAW_MAXSPANQUEUE)
3822                                 {
3823                                         DPSOFTRAST_Draw_ProcessSpans();
3824                                         dpsoftrast.draw.numspans = 0;
3825                                 }
3826                         }
3827                 }
3828                 // draw outlines over triangle for debugging
3829         //      for (j = 0, k = numpoints-1;j < numpoints;k = j, j++)
3830         //              DPSOFTRAST_Draw_DebugEdgePoints(screen[k], screen[j]);
3831         }
3832         if (dpsoftrast.draw.numspans)
3833         {
3834                 DPSOFTRAST_Draw_ProcessSpans();
3835                 dpsoftrast.draw.numspans = 0;
3836         }
3837 #endif
3838 }
3839
3840 void DPSOFTRAST_Draw_DebugPoints(void)
3841 {
3842         int i;
3843         int x;
3844         int y;
3845         int numvertices = dpsoftrast.draw.numvertices;
3846         int w = dpsoftrast.fb_width;
3847         int bounds[4];
3848         unsigned int *pixels = dpsoftrast.fb_colorpixels[0];
3849         const float *c4f;
3850         bounds[0] = dpsoftrast.fb_viewportscissor[0];
3851         bounds[1] = dpsoftrast.fb_viewportscissor[1];
3852         bounds[2] = dpsoftrast.fb_viewportscissor[0] + dpsoftrast.fb_viewportscissor[2];
3853         bounds[3] = dpsoftrast.fb_viewportscissor[1] + dpsoftrast.fb_viewportscissor[3];
3854         for (i = 0;i < numvertices;i++)
3855         {
3856                 // check nearclip
3857                 //if (dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+3] != 1.0f)
3858                 //      continue;
3859                 x = (int)(dpsoftrast.draw.screencoord4f[i*4+0]);
3860                 y = (int)(dpsoftrast.draw.screencoord4f[i*4+1]);
3861                 //x = (int)(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+0] + 0.5f);
3862                 //y = (int)(dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+1] + 0.5f);
3863                 //x = (int)((dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+0] + 1.0f) * dpsoftrast.fb_width * 0.5f + 0.5f);
3864                 //y = (int)((dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION][i*4+1] + 1.0f) * dpsoftrast.fb_height * 0.5f + 0.5f);
3865                 if (x < bounds[0] || y < bounds[1] || x >= bounds[2] || y >= bounds[3])
3866                         continue;
3867                 c4f = dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_COLOR] + i*4;
3868                 pixels[y*w+x] = DPSOFTRAST_BGRA8_FROM_RGBA32F(c4f[0], c4f[1], c4f[2], c4f[3]);
3869         }
3870 }
3871
3872 void DPSOFTRAST_DrawTriangles(int firstvertex, int numvertices, int numtriangles, const int *element3i, const unsigned short *element3s)
3873 {
3874         unsigned char arraymask[DPSOFTRAST_ARRAY_TOTAL];
3875         arraymask[0] = true;
3876         arraymask[1] = dpsoftrast.fb_colorpixels[0] != NULL; // TODO: optimize (decide based on shadermode)
3877         arraymask[2] = dpsoftrast.pointer_texcoordf[0] != NULL;
3878         arraymask[3] = dpsoftrast.pointer_texcoordf[1] != NULL;
3879         arraymask[4] = dpsoftrast.pointer_texcoordf[2] != NULL;
3880         arraymask[5] = dpsoftrast.pointer_texcoordf[3] != NULL;
3881         arraymask[6] = dpsoftrast.pointer_texcoordf[4] != NULL;
3882         arraymask[7] = dpsoftrast.pointer_texcoordf[5] != NULL;
3883         arraymask[8] = dpsoftrast.pointer_texcoordf[6] != NULL;
3884         arraymask[9] = dpsoftrast.pointer_texcoordf[7] != NULL;
3885         DPSOFTRAST_Validate(DPSOFTRAST_VALIDATE_DRAW);
3886         DPSOFTRAST_Draw_LoadVertices(firstvertex, numvertices, true);
3887         DPSOFTRAST_ShaderModeTable[dpsoftrast.shader_mode].Vertex();
3888         DPSOFTRAST_Draw_ProjectVertices(dpsoftrast.draw.screencoord4f, dpsoftrast.draw.post_array4f[DPSOFTRAST_ARRAY_POSITION], numvertices);
3889         DPSOFTRAST_Draw_ProcessTriangles(firstvertex, numtriangles, element3i, element3s, arraymask);
3890 }
3891
3892 void DPSOFTRAST_Init(int width, int height, unsigned int *colorpixels, unsigned int *depthpixels)
3893 {
3894         union
3895         {
3896                 int i;
3897                 unsigned char b[4];
3898         }
3899         u;
3900         u.i = 1;
3901         memset(&dpsoftrast, 0, sizeof(dpsoftrast));
3902         dpsoftrast.bigendian = u.b[3];
3903         dpsoftrast.fb_width = width;
3904         dpsoftrast.fb_height = height;
3905         dpsoftrast.fb_depthpixels = depthpixels;
3906         dpsoftrast.fb_colorpixels[0] = colorpixels;
3907         dpsoftrast.fb_colorpixels[1] = NULL;
3908         dpsoftrast.fb_colorpixels[1] = NULL;
3909         dpsoftrast.fb_colorpixels[1] = NULL;
3910         dpsoftrast.texture_firstfree = 1;
3911         dpsoftrast.texture_end = 1;
3912         dpsoftrast.texture_max = 0;
3913         dpsoftrast.user.colormask[0] = 1;
3914         dpsoftrast.user.colormask[1] = 1;
3915         dpsoftrast.user.colormask[2] = 1;
3916         dpsoftrast.user.colormask[3] = 1;
3917         dpsoftrast.user.blendfunc[0] = GL_ONE;
3918         dpsoftrast.user.blendfunc[1] = GL_ZERO;
3919         dpsoftrast.user.depthmask = true;
3920         dpsoftrast.user.depthtest = true;
3921         dpsoftrast.user.depthfunc = GL_LEQUAL;
3922         dpsoftrast.user.scissortest = false;
3923         dpsoftrast.user.cullface = GL_BACK;
3924         dpsoftrast.user.alphatest = false;
3925         dpsoftrast.user.alphafunc = GL_GREATER;
3926         dpsoftrast.user.alphavalue = 0.5f;
3927         dpsoftrast.user.scissor[0] = 0;
3928         dpsoftrast.user.scissor[1] = 0;
3929         dpsoftrast.user.scissor[2] = dpsoftrast.fb_width;
3930         dpsoftrast.user.scissor[3] = dpsoftrast.fb_height;
3931         dpsoftrast.user.viewport[0] = 0;
3932         dpsoftrast.user.viewport[1] = 0;
3933         dpsoftrast.user.viewport[2] = dpsoftrast.fb_width;
3934         dpsoftrast.user.viewport[3] = dpsoftrast.fb_height;
3935         dpsoftrast.user.depthrange[0] = 0;
3936         dpsoftrast.user.depthrange[1] = 1;
3937         dpsoftrast.user.polygonoffset[0] = 0;
3938         dpsoftrast.user.polygonoffset[1] = 0;
3939         dpsoftrast.user.color[0] = 1;
3940         dpsoftrast.user.color[1] = 1;
3941         dpsoftrast.user.color[2] = 1;
3942         dpsoftrast.user.color[3] = 1;
3943         dpsoftrast.validate = -1;
3944         DPSOFTRAST_Validate(-1);
3945         dpsoftrast.validate = 0;
3946 }
3947
3948 void DPSOFTRAST_Shutdown(void)
3949 {
3950         int i;
3951         for (i = 0;i < dpsoftrast.texture_end;i++)
3952                 if (dpsoftrast.texture[i].bytes)
3953                         MM_FREE(dpsoftrast.texture[i].bytes);
3954         if (dpsoftrast.texture)
3955                 free(dpsoftrast.texture);
3956         memset(&dpsoftrast, 0, sizeof(dpsoftrast));
3957 }
3958