/**
  * @return true if all cells in m have the same value as the corresponding cell in bs. The method
  *     assumes that m.dimensions==b.dimensions.
  */
 private static final boolean equivalent(BooleanMatrix m, BooleanMatrix b) {
   if (!(equivalent(m.dimensions(), b.dimensions()) && m.density() == b.density())) return false;
   int max = m.dimensions().capacity();
   for (int i = 0; i < max; i++) {
     if (m.get(i) != b.get(i)) return false;
   }
   return true;
 }
  @Test
  public final void testCrossProduct() {
    final BooleanMatrix mT43 = mF43.not();

    fill(mT43, range(0, dim43.capacity() - 1));
    fill(mF4, range(0, dim4.capacity() - 1));

    assertCrossProductFalse(mF324, mT43);
    assertCrossProductFalse(mF324, mF4);

    fill(mT324, range(0, dim324.capacity() - 1));

    BooleanValue[] result = new BooleanValue[dim324.cross(dim43).capacity()];
    Arrays.fill(result, TRUE);

    final int c324 = dim324.capacity(), c43 = dim43.capacity();
    final int c32443 = c324 * c43;

    for (int i = 0; i < c32443; i++) {
      result[i] = f.and(mT324.get(i / c43), mT43.get(i % c43));
    }
    assertTrue(equivalent(mT324.cross(mT43), result));

    mT324.set(1, TRUE);
    for (int i = c43; i < c43 * 2; i++) {
      result[i] = mT43.get(i % c43);
    }

    assertTrue(equivalent(mT324.cross(mT43), result));

    mT43.set(5, TRUE);
    for (int i = 0; i < c324; i++) {
      result[i * c43 + 5] = mT324.get(i);
    }

    assertTrue(equivalent(mT324.cross(mT43), result));

    mT324.set(10, FALSE);
    for (int i = c43 * 10; i < c43 * 11; i++) {
      result[i] = FALSE;
    }

    assertTrue(equivalent(mT324.cross(mT43), result));
    // System.out.println(Arrays.asList(result));
    // System.out.println(mT324.crossProduct(mT43));

  }
 /**
  * @return true if all cells in m have the same value as the corresponding cell in the given array
  *     of formulas. The method assumes that m.dimensions.capacity==cells.length.
  */
 private static final boolean equivalent(BooleanMatrix m, BooleanValue[] cells) {
   for (int i = 0; i < cells.length; i++) {
     if (m.get(i) != cells[i]) {
       // System.out.println(i + ": " + m.get(i) + " vs. " + cells[i]);
       return false;
     }
   }
   return true;
 }
  @Test
  public final void testDotProduct() {
    fill(mF43, range(0, dim43.capacity() - 1));
    fill(mF4, range(0, dim4.capacity() - 1));

    assertDotProductFalse(mF324, mF43);
    assertDotProductFalse(mF324, mF4);

    fill(mF324, range(0, dim324.capacity() - 1));

    BooleanValue[] result = new BooleanValue[dim324.dot(dim43).capacity()];
    for (int i = 0; i < result.length; i++) {
      result[i] = BooleanAccumulator.treeGate(Operator.Nary.OR);
    }
    int rows43 = dim324.dimension(dim324.numDimensions() - 1);
    int rows324 = dim324.capacity() / rows43;
    int cols43 = dim43.capacity() / rows43;
    for (int i = 0; i < rows324; i++) {
      for (int j = 0; j < cols43; j++) {
        for (int k = 0; k < rows43; k++) {
          int index324 = i * rows43 + k;
          int index43 = j + k * cols43;
          int indexR = cols43 * i + j;
          ((BooleanAccumulator) result[indexR]).add(f.and(mF324.get(index324), mF43.get(index43)));
        }
      }
    }
    for (int i = 0; i < result.length; i++) {
      result[i] = f.accumulate((BooleanAccumulator) result[i]);
    }

    assertTrue(equivalent(mF324.dot(mF43), result));

    for (int i = 0; i < dim324.capacity(); i += 2) {
      mF324.set(i, FALSE);
    }
    for (int i = 1; i < dim4.capacity(); i += 2) {
      mF4.set(i, FALSE);
    }

    assertDotProductFalse(mF324, mF4);
  }
  @Test
  public final void testOverride() {
    assertTrue(equivalent(mT324.override(mT324), mT324));
    assertTrue(equivalent(mT324.override(mF324), mT324));
    assertTrue(equivalent(mF324.override(mT324), mT324));
    assertTrue(equivalent(mF324.override(mF324), mF324));

    final BooleanMatrix mF324c = mF324.clone(), mT324c = mT324.clone();

    mF324.set(3, vars[3]);
    mF324.set(17, vars[17]);
    mF324.set(22, vars[22]);
    assertTrue(equivalent(mF324.override(mF324c), mF324));
    assertTrue(equivalent(mF324.override(mT324), mT324));

    mF324c.set(9, vars[9]);
    assertTrue(equivalent(mF324.override(mF324c), mF324.or(mF324c)));

    mT324.set(0, BooleanConstant.FALSE);
    assertTrue(equivalent(mF324.override(mT324), mT324));
    assertTrue(equivalent(mT324.override(mT324c), mT324c));
    assertTrue(equivalent(mT324c.override(mT324), mT324));

    final BooleanMatrix mFoF = f.matrix(dim324);
    mF324.set(10, vars[10]);
    mF324c.set(3, vars[4]);
    mF324c.set(20, vars[20]);
    mF324c.set(19, vars[19]);

    mFoF.set(3, f.or(vars[4], f.and(vars[3], f.not(vars[4]))));
    mFoF.set(9, vars[9]);
    mFoF.set(10, f.and(vars[10], f.not(vars[9])));
    mFoF.set(17, f.and(vars[17], f.and(f.not(vars[19]), f.not(vars[20]))));
    mFoF.set(22, f.and(vars[22], f.and(f.not(vars[19]), f.not(vars[20]))));
    mFoF.set(20, vars[20]);
    mFoF.set(19, vars[19]);

    assertTrue(equivalent(mF324.override(mF324c), mFoF));

    mT324.set(3, vars[4]);
    mT324.set(11, vars[11]);
    for (int i = 16; i < 24; i++) mT324.set(i, vars[i - 16]);

    final BooleanMatrix mFoT = f.matrix(dim324);
    for (int i = 0; i < 16; i++) mFoT.set(i, mT324.get(i));

    final BooleanAccumulator g = BooleanAccumulator.treeGate(Operator.AND);
    for (int i = 0; i < 8; i++) g.add(f.not(vars[i]));
    final BooleanValue v3 = f.accumulate(g);

    for (int i = 16; i < 24; i++) mFoT.set(i, f.or(f.and(v3, mF324.get(i)), mT324.get(i)));

    assertTrue(equivalent(mF324.override(mT324), mFoT));

    final BooleanMatrix mToF = f.matrix(dim324);
    for (int i = 0; i < 8; i++)
      mToF.set(i, f.or(mF324.get(i), f.and(mT324.get(i), f.not(vars[3]))));
    for (int i = 8; i < 16; i++)
      mToF.set(i, f.or(mF324.get(i), f.and(mT324.get(i), f.not(vars[10]))));
    for (int i = 16; i < 24; i++)
      mToF.set(i, f.or(mF324.get(i), f.and(mT324.get(i), f.and(f.not(vars[17]), f.not(vars[22])))));

    assertTrue(equivalent(mT324.override(mF324), mToF));

    final BooleanMatrix mToT = f.matrix(dim324).not();
    mT324c.set(11, vars[12]);
    mT324c.set(12, vars[13]);
    mT324c.set(18, vars[18]);
    for (int i = 0; i < 16; i++) mToT.set(i, mT324.get(i));
    for (int i = 16; i < 24; i++) mToT.set(i, f.or(mT324.get(i), f.and(mT324c.get(i), v3)));

    assertTrue(equivalent(mT324c.override(mT324), mToT));
  }