410     f32 mag;  in C_QUATNormalize()  local
415     mag = (src->x * src->x) + (src->y * src->y) + (src->z * src->z) + (src->w * src->w);  in C_QUATNormalize()
417     if ( mag >= QUAT_EPSILON )  in C_QUATNormalize()
419         mag = 1.0F / sqrtf(mag);  in C_QUATNormalize()
421         unit->x = src->x * mag;  in C_QUATNormalize()
422         unit->y = src->y * mag;  in C_QUATNormalize()
423         unit->z = src->z * mag;  in C_QUATNormalize()
424         unit->w = src->w * mag;  in C_QUATNormalize()
441     register f32    mag, rsqmag, diff, c_zero;  in PSQUATNormalize()  local
452         ps_mul      mag, sxy, sxy  in PSQUATNormalize()
459         ps_madd     mag, szw, szw, mag  in PSQUATNormalize()
461         ps_sum0     mag, mag, mag, mag  in PSQUATNormalize()
464         frsqrte     rsqmag, mag  in PSQUATNormalize()
466         ps_sub      diff, mag, epsilon  in PSQUATNormalize()
470         fnmsub      nwork0, nwork0, mag, c_three  in PSQUATNormalize()
501     f32 mag, norminv;  in C_QUATInverse()  local
506     mag = ( src->x*src->x + src->y*src->y + src->z*src->z + src->w*src->w );  in C_QUATInverse()
508     if ( mag == 0.0f )  in C_QUATInverse()
510         mag = 1.0f;  in C_QUATInverse()
514     norminv = 1.0f / mag;  in C_QUATInverse()
530     register f32    mag, nmag, norminv, nninv, nwork0, c_two, c_zero;  in PSQUATInverse()  local
539         ps_mul      mag, sxy, sxy  in PSQUATInverse()
547         ps_madd     mag, szw, szw, mag  in PSQUATInverse()
551         ps_sum0     mag, mag, mag, mag  in PSQUATInverse()
554         fcmpu       cr0, mag, c_zero  in PSQUATInverse()
558         fres        norminv, mag  in PSQUATInverse()
560         ps_neg      nmag, mag  in PSQUATInverse()
562         ps_nmsub    nwork0, mag, norminv, c_two  in PSQUATInverse()
689         f32 mag;  in C_QUATLogN()  local
690         mag = scale + q->z*q->z;  in C_QUATLogN()
691         if ( mag < 1.0F - QUAT_EPSILON || mag > 1.0F + QUAT_EPSILON )  in C_QUATLogN()