random_1_series_parallel.py

import sys
import os
import subprocess

DIR='random_series_parallel'

try:
  os.mkdir(DIR)
except:
  pass

CMR_BUILD = './build'
CMR_GENERATE_SERIES_PARALLEL = CMR_BUILD + '/cmr-generate-series-parallel'
CMR_PERTURB = CMR_BUILD + '/cmr-perturb-random'
CMR_SERIES_PARALLEL = CMR_BUILD + '/cmr-series-parallel'
CMR_GRAPHIC = CMR_BUILD + '/cmr-graphic'
CMR_NETWORK = CMR_BUILD + '/cmr-network'

N = 10000
#N = 2000 # for debugging
repetitions = 100
#repetitions = 10 # for debugging

for kary in ['binary', 'ternary']:
    for beta in [0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]:
        alpha = 0.0
        PROBABILITY = 1
        BASE = max(1, int(alpha*N))
        UNIT = min(N - BASE, int(beta*N))
        gamma = 1.0 - alpha- beta
        COPY = N - BASE - UNIT
        assert COPY >= 0
        delta = 0.0
        makeTernary = '-t' if kary == 'ternary' else ''
        testBinary = '-b' if kary == 'binary' else ''
        PREFIX=f'{DIR}/{kary}-{N}-{alpha:.2f}-{beta:.2f}-{gamma:.2f}-{delta:.2f}-{PROBABILITY}'
        open(f'{PREFIX}.sp', 'w').close()
        open(f'{PREFIX}.gra', 'w').close()
        for i in range(repetitions):
            subprocess.call(f'{CMR_GENERATE_SERIES_PARALLEL} {BASE} {BASE} -p {PROBABILITY} -u {UNIT} {UNIT} -c {COPY} {COPY} {makeTernary} > {PREFIX}-{i}.sparse', shell=True)
            subprocess.call(f'{CMR_SERIES_PARALLEL} -i sparse {PREFIX}-{i}.sparse {testBinary} -R - -N - -s >> {PREFIX}.sp', shell=True)
            if kary == 'binary':
                subprocess.call(f'{CMR_GRAPHIC} -i sparse {PREFIX}-{i}.sparse -s >> {PREFIX}.gra', shell=True)
            os.unlink(f'{PREFIX}-{i}.sparse')