Функция, которая считает момент совпадения двух точек при равноускоренном движении

double coordinate(const double *x0,const double *a,const double *v,const double *t){
    return (*a)*pow(*t,2)/2+(*v)*(*t)+*x0;
}

double overlap(const double *x0_1,const double *a1,const double *v1,const double *x0_2,const double *a2,const double *v2){
    double valueT=0;
    double k=0; //шаг
    double value1,value2; 
    double D=pow((*v1-*v2),2)-2*(*a1-*a2)*(*x0_1-*x0_2); //для проверки условия совпадения двух парабол
    do{
        if (D<0){ //если такой дискриминант < 0, то графики, а соответственно и точки, не пересекутся
            return nan("");
        }
        if(x0_1==x0_2){ 
            return 0;
        }
        value1 = coordinate(x0_1, a1, v1, &k); // рассчет координаты 1-й точки в момент времени k
        value2 = coordinate(x0_2, a2, v2, &k); // рассчет координаты 2-й точки в момент времени k
        if(fabs(value1-value2)<=DBL_MIN){ //если разница бесконечно мала, то это и есть момент пересечения
            valueT=k; 
            break;
        }
/* И вот здесь кроется проблема, т.к. с таким малым шагом цикл крутится бесконечно.
Пробовал шагать по 0.1(откровенный костыль) - для маленьких значений x0,v,a все в порядке, для больших снова бесконечный цикл, т.к. точки если и встречаются, то где-то в бесконечности.
Как вариант, сравнить производные для двух функций, но как?*/
        k+=DBL_MIN; 
    }while(fabs(value1-value2)>=DBL_MIN); /*Та же проблема*/
    valueT=k;
    return valueT;
}