411     f32 mag;  in C_QUATNormalize()  local
416     mag = (src->x * src->x) + (src->y * src->y) + (src->z * src->z) + (src->w * src->w);  in C_QUATNormalize()
418     if ( mag >= QUAT_EPSILON )  in C_QUATNormalize()
420         mag = 1.0F / sqrtf(mag);  in C_QUATNormalize()
422         unit->x = src->x * mag;  in C_QUATNormalize()
423         unit->y = src->y * mag;  in C_QUATNormalize()
424         unit->z = src->z * mag;  in C_QUATNormalize()
425         unit->w = src->w * mag;  in C_QUATNormalize()
442     register f32    mag, rsqmag, diff, c_zero;  in PSQUATNormalize()  local
453         ps_mul      mag, sxy, sxy  in PSQUATNormalize()
460         ps_madd     mag, szw, szw, mag  in PSQUATNormalize()
462         ps_sum0     mag, mag, mag, mag  in PSQUATNormalize()
465         frsqrte     rsqmag, mag  in PSQUATNormalize()
467         ps_sub      diff, mag, epsilon  in PSQUATNormalize()
471         fnmsub      nwork0, nwork0, mag, c_three  in PSQUATNormalize()
502     f32 mag, norminv;  in C_QUATInverse()  local
507     mag = ( src->x*src->x + src->y*src->y + src->z*src->z + src->w*src->w );  in C_QUATInverse()
509     if ( mag == 0.0f )  in C_QUATInverse()
511         mag = 1.0f;  in C_QUATInverse()
515     norminv = 1.0f / mag;  in C_QUATInverse()
531     register f32    mag, nmag, norminv, nninv, nwork0, c_two, c_zero;  in PSQUATInverse()  local
540         ps_mul      mag, sxy, sxy  in PSQUATInverse()
548         ps_madd     mag, szw, szw, mag  in PSQUATInverse()
552         ps_sum0     mag, mag, mag, mag  in PSQUATInverse()
555         fcmpu       cr0, mag, c_zero  in PSQUATInverse()
559         fres        norminv, mag  in PSQUATInverse()
561         ps_neg      nmag, mag  in PSQUATInverse()
563         ps_nmsub    nwork0, mag, norminv, c_two  in PSQUATInverse()
690         f32 mag;  in C_QUATLogN()  local
691         mag = scale + q->z*q->z;  in C_QUATLogN()
692         if ( mag < 1.0F - QUAT_EPSILON || mag > 1.0F + QUAT_EPSILON )  in C_QUATLogN()