added some more tests. Something goes wrong in the updating of the neighbors, that is, some reactions comes out of the initializations active yet illegal. Was stuck for a while due to mess with matrix indices.
merged some work.