fixed slow implementation by vectorizing

orchidas · Oct 23, 2021 · 226cd3b · 226cd3b
1 parent 05e9eea
commit 226cd3b
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/Source/CoupledMixingMatrix.cpp b/Source/CoupledMixingMatrix.cpp
@@ -64,15 +64,17 @@ void CoupledMixingMatrix::updateCouplingFilters(){
 Eigen::VectorXf CoupledMixingMatrix::process(Eigen::VectorXf delayLineOutput){
 
     Eigen::VectorXf delayLineInput(delayLineOutput.size());
-    delayLineInput.setZero();
+    //delayLineInput.setZero();
+    delayLineInput = (M_block.cwiseProduct(couplingScalars)) * delayLineOutput;
 
-    for(int i = 0; i < nDelayLines; i++){
+    //doing this in a loop is too slow and breaks things
+    /*for(int i = 0; i < nDelayLines; i++){
         for(int j = 0; j < nDelayLines; j++){
             if (!isFilter)
-                delayLineInput(i) += M_block(i,j)*couplingScalars(i,j)*delayLineOutput(i);
+                delayLineInput(i) += M_block(i,j)*couplingScalars(i,j)*delayLineOutput(j);
             
         }
-    }
+    }*/
 
     return delayLineInput;
 }