/*
Copyright (C) 1996-1997 Id Software, Inc.

This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License
as published by the Free Software Foundation; either version 2
of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  

See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.

*/
// d_part.c: software driver module for drawing particles

#include "quakedef.h"
#include "d_local.h"
#include "r_local.h"

vec3_t			r_pright, r_pup, r_ppn;

//Spike: Particles are depth sorted. So why depth write? They are the last to be drawn anyway.
#define PARTICLEFACTOR 0x8000 // Change DP_Partfac in ASM to match this

/*
==============
D_EndParticles
==============
*/
void D_EndParticles (void)
{
// not used by software driver
}


/*
==============
D_StartParticles
==============
*/
void D_StartParticles (void)
{
// not used by software driver
}

/*
==============
D_DrawParticle
==============
*/
#if	!id386

void D_DrawParticle8S (vec3_t porg, float palpha, float pscale, unsigned int pcolour)
{
	vec3_t	local, transformed;
	float	zi;
	qbyte	*pdest;
	short	*pz;
	int		i, izi, pix, count, u, v;

// transform point
	VectorSubtract (porg, r_origin, local);

	transformed[0] = DotProduct(local, r_pright);
	transformed[1] = DotProduct(local, r_pup);
	transformed[2] = DotProduct(local, r_ppn);		

	if (transformed[2] < PARTICLE_Z_CLIP)
		return;

// project the point
// FIXME: preadjust xcenter and ycenter
	zi = 1.0 / transformed[2];
	u = (int)(xcenter + zi * transformed[0] + 0.5);
	v = (int)(ycenter - zi * transformed[1] + 0.5);

	if ((v > d_vrectbottom_particle) || 
		(u > d_vrectright_particle) ||
		(v < d_vrecty) ||
		(u < d_vrectx))
	{
		return;
	}

	pz = d_pzbuffer + (d_zwidth * v) + u;
	pdest = d_viewbuffer + d_scantable[v] + u;
	izi = (int)(zi * PARTICLEFACTOR);

	pix = izi >> d_pix_shift;
	pix *= pscale;

	if (pix < d_pix_min)
		pix = d_pix_min;
	else if (pix > d_pix_max)
		pix = d_pix_max;

	switch (pix)
	{
	case 1:
		count = 1 << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			if (pz[0] <= izi)
			{
//				pz[0] = izi;
				pdest[0] = pcolour;
			}
		}
		break;

	case 2:
		count = 2 << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			if (pz[0] <= izi)
			{
//				pz[0] = izi;
				pdest[0] = pcolour;
			}

			if (pz[1] <= izi)
			{
//				pz[1] = izi;
				pdest[1] = pcolour;
			}
		}
		break;

	case 3:
		count = 3 << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			if (pz[0] <= izi)
			{
//				pz[0] = izi;
				pdest[0] = pcolour;
			}

			if (pz[1] <= izi)
			{
//				pz[1] = izi;
				pdest[1] = pcolour;
			}

			if (pz[2] <= izi)
			{
//				pz[2] = izi;
				pdest[2] = pcolour;
			}
		}
		break;

	case 4:
		count = 4 << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			if (pz[0] <= izi)
			{
//				pz[0] = izi;
				pdest[0] = pcolour;
			}

			if (pz[1] <= izi)
			{
//				pz[1] = izi;
				pdest[1] = pcolour;
			}

			if (pz[2] <= izi)
			{
//				pz[2] = izi;
				pdest[2] = pcolour;
			}

			if (pz[3] <= izi)
			{
//				pz[3] = izi;
				pdest[3] = pcolour;
			}
		}
		break;

	default:
		count = pix << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			for (i=0 ; i<pix ; i++)
			{
				if (pz[i] <= izi)
				{
//					pz[i] = izi;
					pdest[i] = pcolour;
				}
			}
		}
		break;
	}
}

#endif	// !id386

void D_DrawParticle16S (vec3_t porg, float palpha, float pscale, unsigned int pcolour)
{
	vec3_t	local, transformed;
	float	zi;
	unsigned short	*pdest;
	int a;
	short	*pz;
	int		i, izi, pix, count, u, v;

	if (palpha <= 0.2)
		return;

// transform point
	VectorSubtract (porg, r_origin, local);

	transformed[0] = DotProduct(local, r_pright);
	transformed[1] = DotProduct(local, r_pup);
	transformed[2] = DotProduct(local, r_ppn);		

	if (transformed[2] < PARTICLE_Z_CLIP)
		return;

// project the point
// FIXME: preadjust xcenter and ycenter
	zi = 1.0 / transformed[2];
	u = (int)(xcenter + zi * transformed[0] + 0.5);
	v = (int)(ycenter - zi * transformed[1] + 0.5);

	if ((v > d_vrectbottom_particle) || 
		(u > d_vrectright_particle) ||
		(v < d_vrecty) ||
		(u < d_vrectx))
	{
		return;
	}

	pz = d_pzbuffer + (d_zwidth * v) + u;	
	izi = (int)(zi * PARTICLEFACTOR);

	pix = ((int)(izi*pscale));	

	if (pix < d_pix_min)
		pix = d_pix_min;
	else if (pix > d_pix_max)
		pix = d_pix_max;

	u -= pix/2;
	v -= pix/2;
	if (u < 0) u = 0;
	if (v < 0) v = 0;
	pdest = (unsigned short	*)d_viewbuffer + ((d_scantable[v] + u));

	a = 255*palpha;

	switch (pix)
	{
	default:
		count = pix << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
		{
			for (i=0 ; i<pix ; i++)
			{
				if (pz[i] <= izi)
				{
//					pz[i] = izi;
					pdest[i] = d_8to16table[pcolour];
				}
			}
		}
		break;
	}
}

#pragma message("fixme: D_DrawParticle16T is not implemented")
#define D_DrawParticle16T D_DrawParticle16S


void D_DrawParticle32T (vec3_t porg, float palpha, float pscale, unsigned int pcolour)
{
	vec3_t	local, transformed;
	float	zi;
	qbyte	*pdest;
	qbyte	*pal;
	int a;
	short	*pz;
	int		i, izi, pix, count, u, v;

	if (palpha <= 0.0)
		return;

// transform point
	VectorSubtract (porg, r_origin, local);

	transformed[0] = DotProduct(local, r_pright);
	transformed[1] = DotProduct(local, r_pup);
	transformed[2] = DotProduct(local, r_ppn);		

	if (transformed[2] < PARTICLE_Z_CLIP)
		return;

// project the point
// FIXME: preadjust xcenter and ycenter
	zi = 1.0 / transformed[2];
	u = (int)(xcenter + zi * transformed[0] + 0.5);
	v = (int)(ycenter - zi * transformed[1] + 0.5);

	if ((v > d_vrectbottom_particle) || 
		(u > d_vrectright_particle) ||
		(v < d_vrecty) ||
		(u < d_vrectx))
	{
		return;
	}

	pz = d_pzbuffer + (d_zwidth * v) + u;	
	izi = (int)(zi * PARTICLEFACTOR);

	pix = ((int)(izi*pscale)) >> d_pix_shift;	

	if (pix < d_pix_min)
		pix = d_pix_min;
	else if (pix > d_pix_max)
		pix = d_pix_max;

	u -= pix/2;
	v -= pix/2;
	if (u < 0) u = 0;
	if (v < 0) v = 0;
	pdest = d_viewbuffer + ((d_scantable[v] + u)<<2);
	pal = (qbyte *)&d_8to32table[pcolour];

	a = 255*palpha;

	switch (pix)
	{
	default:
		count = pix << d_y_aspect_shift;

		for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth<<2)
		{
			for (i=0 ; i<pix ; i++)
			{
				if (pz[i] <= izi)
				{
//					pz[i] = izi;
					pdest[(i<<2)+0] = (pdest[(i<<2)+0]*(255-a) + pal[0]*a) / 255;
					pdest[(i<<2)+1] = (pdest[(i<<2)+1]*(255-a) + pal[1]*a) / 255;
					pdest[(i<<2)+2] = (pdest[(i<<2)+2]*(255-a) + pal[2]*a) / 255;
				}
			}
		}
		break;
	}
}

#define draw(x, y) x=Trans((qbyte)x,(qbyte)y)
#define addblend(x, y) x=AddBlend((qbyte)x,(qbyte)y)
void D_DrawParticleTrans (vec3_t porg, float palpha, float pscale, unsigned int pcolour, blendmode_t blendmode)
{
	vec3_t	local, transformed;
	float	zi;
	qbyte	*pdest;
	short	*pz;
	int		i, izi, pix, count, u, v;

	if (r_pixbytes == 4)
	{
		D_DrawParticle32T(porg, palpha, pscale, pcolour);
		return;
	}
	if (r_pixbytes == 2)
	{
		D_DrawParticle16T(porg, palpha, pscale, pcolour);
		return;
	}

	if (palpha < TRANS_LOWER_CAP)
		return;

	if (palpha > TRANS_UPPER_CAP && blendmode == BM_BLEND)
	{
		D_DrawParticle8S(porg, palpha, pscale, pcolour);
		return;
	}

	D_SetTransLevel(palpha, blendmode);

// transform point
	VectorSubtract (porg, r_origin, local);

	transformed[0] = DotProduct(local, r_pright);
	transformed[1] = DotProduct(local, r_pup);
	transformed[2] = DotProduct(local, r_ppn);		

	if (transformed[2] < PARTICLE_Z_CLIP)
		return;

// project the point
// FIXME: preadjust xcenter and ycenter
	zi = 1.0 / transformed[2];
	u = (int)(xcenter + zi * transformed[0] + 0.5);
	v = (int)(ycenter - zi * transformed[1] + 0.5);

	if ((v > d_vrectbottom_particle) || 
		(u > d_vrectright_particle) ||
		(v < d_vrecty) ||
		(u < d_vrectx))
	{
		return;
	}

	pz = d_pzbuffer + (d_zwidth * v) + u;	
	izi = (int)(zi * PARTICLEFACTOR);

	pix = ((int)(izi*pscale)) >> d_pix_shift;	

	if (pix < d_pix_min)
		pix = d_pix_min;
	else if (pix > d_pix_max)
		pix = d_pix_max;

	u -= pix/2;
	v -= pix/2;
	if (u < 0) u = 0;
	if (v < 0) v = 0;
	pdest = d_viewbuffer + d_scantable[v] + u;

	if (blendmode == BM_ADD) // additive drawing
	{
		switch (pix)
		{
		case 1:
			count = 1 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					addblend(pdest[0], pcolour);
				}
			}
			break;

		case 2:
			count = 2 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					addblend(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					addblend(pdest[1], pcolour);
				}
			}
			break;

		case 3:
			count = 3 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					addblend(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					addblend(pdest[1], pcolour);
				}

				if (pz[2] <= izi)
				{
	//				pz[2] = izi;
					addblend(pdest[2], pcolour);
				}
			}
			break;

		case 4:
			count = 4 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					addblend(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					addblend(pdest[1], pcolour);
				}

				if (pz[2] <= izi)
				{
	//				pz[2] = izi;
					addblend(pdest[2], pcolour);
				}

				if (pz[3] <= izi)
				{
	//				pz[3] = izi;
					addblend(pdest[3], pcolour);
				}
			}
			break;

		default:
			count = pix << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				for (i=0 ; i<pix ; i++)
				{
					if (pz[i] <= izi)
					{
	//					pz[i] = izi;
						addblend(pdest[i], pcolour);
					}
				}
			}
			break;
		}
	}
	else // merge drawing
	{
		switch (pix)
		{
		case 1:
			count = 1 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					draw(pdest[0], pcolour);
				}
			}
			break;

		case 2:
			count = 2 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					draw(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					draw(pdest[1], pcolour);
				}
			}
			break;

		case 3:
			count = 3 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					draw(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					draw(pdest[1], pcolour);
				}

				if (pz[2] <= izi)
				{
	//				pz[2] = izi;
					draw(pdest[2], pcolour);
				}
			}
			break;

		case 4:
			count = 4 << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				if (pz[0] <= izi)
				{
	//				pz[0] = izi;
					draw(pdest[0], pcolour);
				}

				if (pz[1] <= izi)
				{
	//				pz[1] = izi;
					draw(pdest[1], pcolour);
				}

				if (pz[2] <= izi)
				{
	//				pz[2] = izi;
					draw(pdest[2], pcolour);
				}

				if (pz[3] <= izi)
				{
	//				pz[3] = izi;
					draw(pdest[3], pcolour);
				}
			}
			break;

		default:
			count = pix << d_y_aspect_shift;

			for ( ; count ; count--, pz += d_zwidth, pdest += screenwidth)
			{
				for (i=0 ; i<pix ; i++)
				{
					if (pz[i] <= izi)
					{
	//					pz[i] = izi;
						draw(pdest[i], pcolour);
					}
				}
			}
			break;
		}
	}
}

void D_2dPos(vec3_t pos, int *u, int *v, int *z)
{
	float	zi;
	vec3_t	local, transformed;

	// transform point
	VectorSubtract (pos, r_origin, local);

	transformed[2] = DotProduct(local, r_ppn);		

	if (transformed[2] < PARTICLE_Z_CLIP)	//near clip
	{
		*u = -1;	//send it off the side intentionally.
		return;
	}

	transformed[0] = DotProduct(local, r_pright);
	transformed[1] = DotProduct(local, r_pup);

// project the point
	zi = 1.0 / transformed[2];
	*u = (int)(xcenter + zi * transformed[0] + 0.5);
	*v = (int)(ycenter - zi * transformed[1] + 0.5);

	*z = (int)(zi * 0x8000);
}
vec_t VI2Length(int x, int y)
{
	float	length;
	length = (float)x*x + (float)y*y;	
	length = sqrt (length);
	return length;
}
void D_DrawSpark32T (vec3_t src, vec3_t dest, float palpha, unsigned int pcolour)	//draw a line in 3d space
{
	/*
	Finds 2d coords for the points, then draws a line between them with an appropriate alpha
	*/
	unsigned char *pdest;
	unsigned char *pal;
	short	*pz;
	int		count, u1, v1, z1, a1, a, ia;
	int u2, v2, z2;

	int du, dv, dz, da;

	if (palpha <= 0.0)
		return;

	D_2dPos(src, &u1, &v1, &z1);
	D_2dPos(dest, &u2, &v2, &z2);

	if ((v1 > d_vrectbottom_particle) || 
		(u1 > d_vrectright_particle) ||
		(v1 < d_vrecty) ||
		(u1 < d_vrectx))
	{
		return;
	}

	if ((v2 > d_vrectbottom_particle) || 
		(u2 > d_vrectright_particle) ||
		(v2 < d_vrecty) ||
		(u2 < d_vrectx))
	{
		return;
	}	
	pal = (qbyte *)(d_8to32table + pcolour);
	a1 = 255 * palpha;

	du = u2 - u1;
	dv = v2 - v1;
	dz = z2 - z1;
	da = 0 - a1;

	if (!du && !dv)
		count = 1;
	else
	{
		count = VI2Length(du, dv);
		if (!count)
			count = 1;
	}

	du *= 256*256;
	dv *= 256*256;
	dz *= 256*256;
	da *= 256*256;
	u1 = u1<<16;
	v1 = v1<<16;
	z1 = z1<<16;
	a1 = a1<<16;
		{
			du /= count;
			dv /= count;
			dz /= count;
			da /= count;
		}
	do
	{		
		pz = d_pzbuffer + (d_zwidth * (v1>>16)) + (u1>>16);		

		if (*pz <= z1>>16)
		{
//			*pz = z1>>16;

			a = a1>>16;
			ia = 255-a;
			pdest = (qbyte *)((unsigned int *)d_viewbuffer + ((d_scantable[v1>>16] + (u1>>16))));
			pdest[0] = (pdest[0]*((ia)) + pal[0]*(a))/255;
			pdest[1] = (pdest[1]*((ia)) + pal[1]*(a))/255;
			pdest[2] = (pdest[2]*((ia)) + pal[2]*(a))/255;
		}

		u1 += du;
		v1 += dv;
		z1 += dz;
		a1 += da;
	} while (count--);
}

void D_DrawSpark16S (vec3_t src, vec3_t dest, float palpha, unsigned int pcolour)	//draw a line in 3d space, 8bpp
{
	unsigned short	*pdest;
	short	*pz;
	int		count, u1, v1, z1;
	int u2, v2, z2;

	int du, dv, dz;

	if (palpha <= 0.0)
		return;

	D_2dPos(src, &u1, &v1, &z1);
	D_2dPos(dest, &u2, &v2, &z2);

	if ((v1 > d_vrectbottom_particle) || 
		(u1 > d_vrectright_particle) ||
		(v1 < d_vrecty) ||
		(u1 < d_vrectx))
	{
		return;
	}

	if ((v2 > d_vrectbottom_particle) || 
		(u2 > d_vrectright_particle) ||
		(v2 < d_vrecty) ||
		(u2 < d_vrectx))
	{
		return;
	}	

	du = u2 - u1;
	dv = v2 - v1;
	dz = z2 - z1;

	if (!du && !dv)
		count = 1;
	else
	{
		count = VI2Length(du, dv);
		if (!count)
			count = 1;
	}

	du *= 256*256;
	dv *= 256*256;
	dz *= 256*256;
	u1 = u1<<16;
	v1 = v1<<16;
	z1 = z1<<16;
		{
			du /= count;
			dv /= count;
			dz /= count;
		}
	do
	{		
		pz = d_pzbuffer + (d_zwidth * (v1>>16)) + (u1>>16);

		if (*pz <= z1>>16)
		{
//			*pz = z1>>16;

			pdest = (unsigned short*)d_viewbuffer + d_scantable[v1>>16] + (u1>>16);
			*pdest = d_8to16table[pcolour];
		}

		u1 += du;
		v1 += dv;
		z1 += dz;
	} while (count--);
}
#pragma message("fixme: D_DrawSpark16T is not implemented")
#define D_DrawSpark16T D_DrawSpark16S

void D_DrawSparkTrans (vec3_t src, vec3_t dest, float palpha, unsigned int pcolour, blendmode_t blendmode)	//draw a line in 3d space, 8bpp
{
	qbyte	*pdest;
	short	*pz;
	int		count, u1, v1, z1;
	int u2, v2, z2;

	int du, dv, dz;

	if (r_pixbytes == 4)
	{
		D_DrawSpark32T(src, dest, palpha, pcolour);
		return;
	}
	if (r_pixbytes == 2)
	{
		D_DrawSpark16T(src, dest, palpha, pcolour);
		return;
	}

	D_SetTransLevel(palpha, blendmode);

	D_2dPos(src, &u1, &v1, &z1);
	D_2dPos(dest, &u2, &v2, &z2);

	if ((v1 > d_vrectbottom_particle) || 
		(u1 > d_vrectright_particle) ||
		(v1 < d_vrecty) ||
		(u1 < d_vrectx))
	{
		return;
	}

	if ((v2 > d_vrectbottom_particle) || 
		(u2 > d_vrectright_particle) ||
		(v2 < d_vrecty) ||
		(u2 < d_vrectx))
	{
		return;
	}	

	du = u2 - u1;
	dv = v2 - v1;
	dz = z2 - z1;

	if (!du && !dv)
		count = 1;
	else
	{
		count = VI2Length(du, dv);
		if (!count)
			count = 1;
	}

	du *= 256*256;
	dv *= 256*256;
	dz *= 256*256;
	u1 = u1<<16;
	v1 = v1<<16;
	z1 = z1<<16;
	du /= count;
	dv /= count;
	dz /= count;

	if (blendmode == BM_ADD) // additive
	{
		do
		{		
			pz = d_pzbuffer + (d_zwidth * (v1>>16)) + (u1>>16);

			if (*pz <= z1>>16)
			{
	//				*pz = z1>>16;
				pdest = d_viewbuffer + d_scantable[v1>>16] + (u1>>16);
				addblend(*pdest, (qbyte)pcolour);
			}

			u1 += du;
			v1 += dv;
			z1 += dz;
		} while (count--);
	}
	else // merge blend
	{
		do
		{		
			pz = d_pzbuffer + (d_zwidth * (v1>>16)) + (u1>>16);

			if (*pz <= z1>>16)
			{
	//				*pz = z1>>16;
				pdest = d_viewbuffer + d_scantable[v1>>16] + (u1>>16);
				draw(*pdest, (qbyte)pcolour);
			}

			u1 += du;
			v1 += dv;
			z1 += dz;
		} while (count--);
	}
}