rounding through intrinsics

Author: -

src/cubeb_mixer.cpp

@@ -377,22 +377,24 @@ MixerContext::init()
 
       for (uint32_t j = 0; j < _in_ch_count; j++) {
         double target = _matrix[i][j] * 32768 + rem;
-        int value = lrintf(target);
+        int value = _mm_cvtsd_si32(_mm_set_sd(target)); // round double to int
         rem += target - value;
         sum += std::abs(value);
       }
       maxsum = std::max(maxsum, sum);
     }
-    if (maxsum > 32768) {
+    if (maxsum > 32768) { //> or >=?
       _clipping = true;
     }
   }
 
   // FIXME quantize for integers
+  // June 23, 2022: I made a change but have no idea what was wanted
   for (uint32_t i = 0; i < CHANNELS_MAX; i++) {
     int ch_in = 0;
     for (uint32_t j = 0; j < CHANNELS_MAX; j++) {
-      _matrix32[i][j] = lrintf(_matrix[i][j] * 32768);
+      // round double to int
+      _matrix32[i][j] = _mm_cvtsd_si32(_mm_set_sd(_matrix[i][j] * 32768));
       if (_matrix[i][j]) {
         _matrix_ch[i][++ch_in] = j;
       }

src/cubeb_sndio.c

@@ -106,7 +106,8 @@ float_to_s16(void * ptr, long nsamp, float volume)
   int s;
 
   while (nsamp-- > 0) {
-    s = lrintf(*(src++) * mult);
+    // round float to int
+    s = _mm_cvt_ss2si(_mm_set_ss(*(src++) * mult));
     if (s < -32768)
       s = -32768;
     else if (s > 32767)

src/cubeb_wasapi.cpp

@@ -938,7 +938,8 @@ refill(cubeb_stream * stm, void * input_buffer, long input_frames_count,
       case CUBEB_SAMPLE_S16NE: {
         short * buf = static_cast<short *>(dest);
         for (long i = 0; i < out_samples; ++i) {
-          buf[i] = static_cast<short>(static_cast<float>(buf[i]) * volume);
+          // round float to int
+          buf[i] = _mm_cvt_ss2si(_mm_set_ss(buf[i] * volume));
         }
         break;
       }