visualize-heat-capacity-generational

#!/usr/bin/env python

import numpy as np
import matplotlib.pyplot as plt
import matplotlib
from os import path, makedirs

# TODO: make these scripts take these as params
loadfile = 'sim-20191114-000009_server'
folder = 'save/' + loadfile
# iter_gen = np.arange(0, 2000, 250)
# iter_gen = np.append(iter_gen, 1999)
# iter_gen = [0, 252, 504, 756, 1002, 1254, 1500, 1752, 1998]
# iter_gen = [0, 1, 2, 5, 10, 20, 50, 100, 250, 1000, 1999,
#             2250, 2500, 2750, 3000, 3250, 3500, 3750, 3999]
# iter_gen = [0, 250, 500, 750, 1000, 1250, 1500, 1750, 1999,
#            2250, 2500, 2750, 3000, 3250, 3500, 3750, 3999]
iter_gen = [500, 1999]

R = 10
Nbetas = 101
betas = 10 ** np.linspace(-1, 1, Nbetas)
numAgents = 150
size = 10

C = np.zeros((R, numAgents, Nbetas, len(iter_gen)))

print('Loading data...')
for ii, iter in enumerate(iter_gen):
    #for bind in np.arange(0, 100):
    for bind in np.arange(1, 100):
        filename = folder + '/C/C_' + str(iter) + '/C-size_' + str(size) + '-Nbetas_' + \
                   str(Nbetas) + '-bind_' + str(bind) + '.npy'
        C[:, :, bind, ii] = np.load(filename)
print('Done.')

plt.rc('text', usetex=True)
font = {'family': 'serif', 'size': 28, 'serif': ['computer modern roman']}
plt.rc('font', **font)
plt.rc('legend', **{'fontsize': 20})

b = 0.8
alpha = 0.3

print('Generating figures...')
for ii, iter in enumerate(iter_gen):

    fig, ax = plt.subplots(1, 1, figsize=(11, 10), sharex=True)
    fig.text(0.51, 0.035, r'$\beta$', ha='center', fontsize=28)
    fig.text(0.005, 0.5, r'$C/N$', va='center', rotation='vertical', fontsize=28)
    title = 'Specific Heat of Foraging Community\n Generation: ' + str(iter)
    fig.suptitle(title)

    # CHANGE THIS TO CUSTOMIZE HEIGHT OF PLOT
    upperbound = 1.5 * np.max(np.mean(np.mean(C[:, :, :-40, :], axis=0), axis=0))
    # upperbound = np.max(np.mean(np.mean(C, axis=0)), axis=0)
    upperbound = 0.4

    label = iter

    for numOrg in range(numAgents):
        c = np.dot(np.random.random(), [1, 1, 1])
        ax.scatter(betas, np.mean(C[:, numOrg, :, ii], axis=0),
                   color=[0, 0, 0], s=30, alpha=alpha, marker='x', label=label)

    xticks = [0.1, 0.5, 1, 2, 4, 10]
    ax.set_xscale("log", nonposx='clip')
    ax.set_xticks(xticks)
    ax.get_xaxis().set_major_formatter(matplotlib.ticker.ScalarFormatter())

    plt.axis([0.1, 10, 0, upperbound])

    # leg = plt.legend(loc=2, title='Generation')
    #
    # for lh in leg.legendHandles:
    #     lh.set_alpha(1)
    #     lh.set_sizes(30)

    savefolder = folder + '/figs/C/'
    savefilename = savefolder + 'C-size_' + str(size) + '-Nbetas_' + \
                   str(Nbetas) + '-gen_' + str(iter) + '.png'
    if not path.exists(savefolder):
        makedirs(savefolder)

    plt.savefig(savefilename, bbox_inches='tight')
    plt.close()
    # plt.clf()
    savemsg = 'Saving ' + savefilename
    print(savemsg)
    # plt.show()
    # plt.pause(0.1)