Adaptive Step Size

double NumericalIntegrator::getDTau(fourVec Accl, double min, double max)
{
    double dtauAdpt=min;
    fourVec dTau(0,0,0,0);
    for(int alpha=0;alpha<4;alpha++){
        if(Accl[alpha]!=0)
            dTau[alpha] = 1.0/Accl[alpha];
        else
            dTau[alpha] = max;
    }
    for(int alpha=0;alpha<4;alpha++){
        dtauAdpt = qMax(dTau[alpha],dtauAdpt);
    }
    if(max!=0)
        return qMin(dtauAdpt,max);
    return dtauAdpt;
}