In [1]:

    

import matplotlib.pyplot as plt
import numpy as np


    

In [2]:

    

from floodfill.vreco import Vreco

v = Vreco('data/colvar_mtd',
          'data/parvar_mtd')
pot, grid_vecs = v.pot(4000)


    

In [20]:

    

extents = [(gv[0], gv[-1]) for gv in grid_vecs]


    

In [45]:

    

fig = plt.figure(figsize=(8, 6), dpi=80)
ax = fig.add_axes([0., 0., 1, 1])
ax.imshow(-pot.T, origin='lower', aspect=0.75,
          interpolation='bilinear')
ax.grid(True)
plt.show()