Our implementation of the BANMF algorithm cycles between updating the NMF factors \(W\) , \(H\) , and the corresponding auxiliary matrix \(Y\) . We use the multiplicative update method introduced by Lee and Seung [1]} for the \(W\)  and \(H\)  updates, but any other alternating NMF algorithms are also suitable. For the update of the auxiliary variable \(Y\) , we project \(Y\)  onto the feasible set defined by the constraints. alg:NMC is an outline of a procedure to solve the BANMF problem (REF ). To find the appropriate threshold parameters \(\delta _W\)  and \(\delta _H\) , and the corresponding Boolean matrices, \(\hat{W}\)  and \(\hat{H}\) , we employ alg:booleanization.
