Algoritmo DDA para líneas

DDA.cpp

void DDA (GLint x1,GLint y1,GLint x2,GLint y2)
{
   float ax,ay,x,y,res;
   GLint i;
   if(abs(x2-x1)>=abs(y2-y1))
    //si la variacion en x es mayor o igual que la variacion en y
    res=abs(x2-x1);
        else
    //si la variacion en y es mayor que la variacion en x
    res=abs(y2-y1);

   ax=(x2-x1)/res;//el valor a aumentar en x
   ay=(y2-y1)/res;//el valor a aumentar en y

   //se realiza casteo a float porque los resultados de la división es un real

   x=(float)x1;
   y=(float)y1;

   i=0;
   while(i<=res)
   {
        setPixel(roundf(x),roundf(y));
        x=x+ax;
        y=y+ay;
        i++;
   }
}