import numpy as np
import matplotlib.pyplot as plt
from time import time

import odeg as od

plt.rcParams['figure.dpi'] = 300
plt.rcParams['savefig.dpi'] = 300


if __name__ == '__main__':
    
    od.utils.blockevs.output_time = True
    
    cfg = {
        'r0': 1,
        'rs': 2,
        'p_max': 4,
        'pos_charge': 10,
        'theta': .5,
        'mu': -.5,
        'load_evs': False,
        'parallel': False,
        }

    print("\n\nconfiguration:\n", cfg, "\n\n")
    
    start = time()
    h_exp, n_exp, n2_exp = od.fctExp(od.energyFct, od.numberFct, od.number2Fct, **cfg)
    print("t[fctexp] =", time() - start)
    # free_energy = od.freeEnergy(**cfg)
    # entropy = od.entropy(**cfg)