outfile << '{';

    bool outputComma = false;

        for (float tau : tauValues) {

            DegreeSequence ds(numVertices);

            // For a single n,tau take samples over several instances of

            // the degree distribution.

            // 500 samples seems to give reasonable results

                // Generate a graph

                // might require multiple tries

                for (int i = 1; ; ++i) {

                std::cout << "Running n = " << numVertices << ", tau = " << tau

                          << ". \t" << std::flush;

                constexpr int measurements = 50;

                constexpr int measureSkip =

                    200; // Take a sample every ... steps

graphdata_exponent_mixN.m

    output: {{n,tau},avgTriangles}

    n from 200 to 2000 with step 200

    initial ErdosGallai

    mixingTime = N * (32.0f - 26.0f*(tau - 2.0f)) * n

    measurements = 50

graphdata_initialtris2.m

    same as above but simulation was stopped during (n,tau)=(1400,2.6)

graphdata_gcm_partial.m

    ??

    output: {{n,tau},triangleseq,ds, {...??...} }