计算几何卡精度

求某些函数时,例如acos(),需要对边界情况进行特判,若超过边界,函数会返回nan,从而导致计算错误

db ran(P a,P b){ 
	db pp=cran(a,b);
	if(sign(pp-1)==0) return 0;
	if(sign(pp+1)==0) return acos(-1);
	return acos(cran(a,b)); 
}