When simulating with negative w values, we would expect to get more and more straight chains the more we grow w to the negative direction. Intuitively, a negative w should push a chemical species to prefer to be neighbour of different species instead of other cells of the same one.
However with current code we don't get this.
@lhilbert suggested this could be caused by wrong initialization of deltaEmin = -10 * omega instead of deltaEmin = -10 * fabs(omega).
When simulating with negative w values, we would expect to get more and more straight chains the more we grow w to the negative direction. Intuitively, a negative w should push a chemical species to prefer to be neighbour of different species instead of other cells of the same one. However with current code we don't get this. @lhilbert suggested this could be caused by wrong initialization of
deltaEmin = -10 * omega
instead ofdeltaEmin = -10 * fabs(omega)
.