bug fix in invertLsMtx2D()

leomccormack · Aug 20, 2021 · e99bf9a · e99bf9a
1 parent e147322
commit e99bf9a
Showing 1 changed file with 3 additions and 5 deletions.
diff --git a/framework/modules/saf_vbap/saf_vbap.c b/framework/modules/saf_vbap/saf_vbap.c
@@ -946,15 +946,13 @@ void invertLsMtx2D
         /* get the unit vectors for the current group */
         for(i=0; i<2; i++)
             for(j=0; j<2; j++)
-                tempGroup[j*2+i] = U_spkr[ls_pairs[n*2+i]*2 + j];
+                tempGroup[j*2+i] = U_spkr[ls_pairs[n*2+i]*2 + j]; /* ^T */
 
         /* get inverse of current group */
         utility_sinv(hSinv, tempGroup, tempInv, 2);
 
-        /* store the vectorized inverse as a row the output */
-        for(i=0; i<2; i++)
-            for(j=0; j<2; j++)
-                (*layoutInvMtx)[n*4+(i*2+j)] = tempInv[j*2+i];
+        /* store the vectorised inverse as a row in the output */
+        cblas_scopy(4, tempInv, 1, (*layoutInvMtx) + n*4, 1);
     }
 
     utility_sinv_destroy(&hSinv);