Padé replacement function

gistfile1.c

// - - - - pade BEGIN - - - -
float coefficient(int iDimension, int iCoefficient)
{
   int iRow, iColumn;
   
   iRow = iDimension*2 + (iCoefficient/4);
   iColumn = iCoefficient % 4;
   
   return tV[iRow][iColumn];

}

float4 PadePosP(float4 PosW)
{

	//KCNote: This bit's changed
	//so that it reads back
	//perspective shift as distance
	//from camera.
	
	
    //perform Pade
    float4 PosVP = float4(0,0,0,1);
    float Pos2D[2];
    float w[2];
    for (int i=0; i<2; i++)
	{
		w[i] = (PosW.x*coefficient(i,4) +
                    PosW.y*coefficient(i,5) +
                    PosW.z*coefficient(i,6) +
                    coefficient(i,7));

        Pos2D[i] = (PosW.x*coefficient(i,0) +
                    PosW.y*coefficient(i,1) +
                    PosW.z*coefficient(i,2) +
                    coefficient(i,3))
                   / w[i];
	}


    return float4(Pos2D[0],Pos2D[1],(1-1/(w[0]+w[1]))/1000.0,1);

}
// - - - - pade END - - - -