public void MakeConstrains() { int u; int v; constrains = new Constrained[0]; // Structural for (v = 0; v < NrV; v++) { for (u = 0; u < NrU; u++) { constrains.Push(new Constrained(particles[index(u, v)], particles[index(u, v + 1)])); constrains.Push(new Constrained(particles[index(u, v)], particles[index(u + 1, v)])); } } for (v = 0; v < NrV; v++) { u = NrU; constrains.Push(new Constrained(particles[index(u, v)], particles[index(u, v + 1)])); } for (u = 0; u < NrU; u++) { v = NrV; constrains.Push(new Constrained(particles[index(u, v)], particles[index(u + 1, v)])); } }