void SimulationIsingModel::run(int nstep) { for (int i=0; i<nstep; i++) { for (int j=0; j<_nmovespercycle; j++) { int p = _rng.randInt()%_n; attemptFlip(p); } } }
void H2PolygonTriangulater::triangulate() { createSteinerPoints(); std::vector<uint> indices; uint i, N = fullPolygon.nbVertices(); indices.resize(N); for (i=0; i<N; i++) { indices[i] = i; } triangulateSubpolygon(indices); sortTriangles(); completeCutsAndSides(); uint counter = 0; while (attemptFlip() && counter!=maxNbFlips) {++counter;} sortTriangles(); completeCutsAndSides(); }