       a1[i *  nr_rowsA + k] += (a2[i *  nr_rowsA + k] + a3[i *  nr_rowsA + k]); 
     }
   }
-  matrixMultSymm_cuBlas(handle, tmpC2, a1, d_aS, d_bS,  tmpTransferLongToGMP,
+  matrixMultSymm_cuBlas(handle, tmpC2, a1, d_aS,  tmpTransferLongToGMP,
 			tmp, sizeMatrix, realExpMatrix, signMatrix, d_prodRes, d_res,
 			nr_rowsA, nr_colsA, prec);
  

      a1[i *  nr_rowsA + k] -= (2 * a2[i *  nr_rowsA + k]); 
    }
  }
  matrixMultSymm_cuBlas(handle, tmpC3, a1, d_aS, d_bS,  tmpTransferLongToGMP,
			tmp, sizeMatrix, realExpMatrix, signMatrix, d_prodRes, d_res,
			nr_rowsA, nr_colsA, prec);
  

      
      gettimeofday(&t1, NULL);
      vecAdd__wSign<<<gridSize, blockSize>>>(d_prodRes, d_res,  nr_rowsA * nr_colsB);
      cudaThreadSynchronize();
      gettimeofday(&t2, NULL);
      etAdd += (((t2.tv_sec*uS_PER_SEC)+t2.tv_usec) - ((t1.tv_sec*uS_PER_SEC)+t1.tv_usec))/(float)uS_PER_mS;
      // This is safe with overflow until there are 1024 of our own limbs that need to be summed up
