Différences
Ci-dessous, les différences entre deux révisions de la page.
Les deux révisions précédentes Révision précédente | |||
teaching:progappchim:matplotlib_gallery:potentiel_energy_surface [2016/12/14 11:57] – [Références] villersd | teaching:progappchim:matplotlib_gallery:potentiel_energy_surface [2020/12/07 17:22] (Version actuelle) – [Programme] villersd | ||
---|---|---|---|
Ligne 32: | Ligne 32: | ||
==== Programme ==== | ==== Programme ==== | ||
- | <sxh python; title : PES-contour-01.py> | + | <code python PES-contour-01.py> |
# | # | ||
# -*- coding: utf-8 -*- | # -*- coding: utf-8 -*- | ||
Ligne 69: | Ligne 69: | ||
X, Y = np.meshgrid(x_1d, | X, Y = np.meshgrid(x_1d, | ||
print X.shape, X.dtype, X.ndim, Y.shape, Y.dtype, Y.ndim | print X.shape, X.dtype, X.ndim, Y.shape, Y.dtype, Y.ndim | ||
- | E=(Q(X)+Q(Y)+Q(X+Y)-np.sqrt(2.*((a(X)-a(Y))**2.+(a(Y)-a(X+Y))**2.+(a(X+Y)-a(X))**2.)))/ | + | E=(Q(X)+Q(Y)+Q(X+Y)-np.sqrt(2.*( (a(X)-a(Y) )**2.+(a(Y)-a(X+Y) )**2.+(a(X+Y)-a(X) )**2.) ))/(1.+k) |
print np.min(E) | print np.min(E) | ||
Ligne 86: | Ligne 86: | ||
plt.ylabel(' | plt.ylabel(' | ||
- | fig = plt.figure(2, | + | fig = plt.figure(2, |
ax = Axes3D(fig) | ax = Axes3D(fig) | ||
ax.plot_surface(X, | ax.plot_surface(X, | ||
Ligne 93: | Ligne 93: | ||
ax.set_zlabel(' | ax.set_zlabel(' | ||
plt.show() | plt.show() | ||
- | </sxh> | + | </code> |
Avec les paramètres essayés, la valeur minimale de E est environ -1.603 | Avec les paramètres essayés, la valeur minimale de E est environ -1.603 |