Numericalpatch

adding patch for numerical stability a sqrt calculation of tau 4-vector for derivatives.