isnt this equal to:



_acos = -c*log(dotprods[j]+c*sqrt(one-(dotprods[j]*dotprods[j])));