487     f32 mag;  in C_QUATNormalize()  local
492     mag = (src->x * src->x) + (src->y * src->y) + (src->z * src->z) + (src->w * src->w);  in C_QUATNormalize()
494     if ( mag >= QUAT_EPSILON )  in C_QUATNormalize()
496         mag = 1.0F / sqrtf(mag);  in C_QUATNormalize()
498         unit->x = src->x * mag;  in C_QUATNormalize()
499         unit->y = src->y * mag;  in C_QUATNormalize()
500         unit->z = src->z * mag;  in C_QUATNormalize()
501         unit->w = src->w * mag;  in C_QUATNormalize()
519     f32x2    mag, rsqmag, diff;  in PSQUATNormalize()  local
529     mag = __PS_MUL(sxy, sxy);  in PSQUATNormalize()
540     mag = __PS_MADD(szw, szw, mag);  in PSQUATNormalize()
544     mag = __PS_SUM0(mag, mag, mag);  in PSQUATNormalize()
548     rsqmag = __PS_RSQRTE(mag);  in PSQUATNormalize()
552     diff = __PS_SUB(mag, epsilon);  in PSQUATNormalize()
562     nwork0 = __PS_NMSUB(nwork0, mag, c33);  in PSQUATNormalize()
606     f32 mag, norminv;  in C_QUATInverse()  local
611     mag = ( src->x*src->x + src->y*src->y + src->z*src->z + src->w*src->w );  in C_QUATInverse()
613     if ( mag == 0.0f )  in C_QUATInverse()
615         mag = 1.0f;  in C_QUATInverse()
619     norminv = 1.0f / mag;  in C_QUATInverse()
636     f32x2    mag, norminv, nninv, nwork0;  in PSQUATInverse()  local
646     mag = __PS_MUL(sxy, sxy);  in PSQUATInverse()
659     mag = __PS_MADD(szw, szw, mag);  in PSQUATInverse()
666     mag = __PS_SUM0(mag, mag, mag);  in PSQUATInverse()
669     if ( mag[0] == 0.0f )  in PSQUATInverse()
671         mag[0] = mag[1] = 1.0f;  in PSQUATInverse()
676     norminv = __PS_RES(mag);  in PSQUATInverse()
680     nwork0 = __PS_NMSUB(mag, norminv, c22);  in PSQUATInverse()
817         f32 mag;  in C_QUATLogN()  local
818         mag = scale + q->z*q->z;  in C_QUATLogN()
819         if ( mag < 1.0F - QUAT_EPSILON || mag > 1.0F + QUAT_EPSILON )  in C_QUATLogN()