Usinf ert from mater commit

Author: daavid00

docs/_sources/installation.rst.txt

@@ -7,6 +7,14 @@ While using package managers such as Anaconda, Miniforge, or Mamba might work, t
 The supported Python versions are 3.11 to 3.13. We will update the documentation when Python3.14 is supported
 (e.g., the ert Python package is not yet available via pip install in Python 3.14).
 
+.. note::
+
+    In Ubuntu, one also needs to install freeglut3-dev:
+
+    .. code-block:: bash
+        
+        sudo apt-get install freeglut3-dev
+
 .. _vpycopm:
 
 Python package

docs/installation.html

@@ -95,6 +95,13 @@ <h1>Installation<a class="headerlink" href="#installation" title="Link to this h
 While using package managers such as Anaconda, Miniforge, or Mamba might work, these are not tested.
 The supported Python versions are 3.11 to 3.13. We will update the documentation when Python3.14 is supported
 (e.g., the ert Python package is not yet available via pip install in Python 3.14).</p>
+<div class="admonition note">
+<p class="admonition-title">Note</p>
+<p>In Ubuntu, one also needs to install freeglut3-dev:</p>
+<div class="highlight-bash notranslate"><div class="highlight"><pre><span></span>sudo<span class="w"> </span>apt-get<span class="w"> </span>install<span class="w"> </span>freeglut3-dev
+</pre></div>
+</div>
+</div>
 <section id="python-package">
 <span id="vpycopm"></span><h2>Python package<a class="headerlink" href="#python-package" title="Link to this heading"></a></h2>
 <p>To install the <strong>pycopm</strong> executable from the development version:</p>

docs/searchindex.js

@@ -7,6 +7,14 @@ While using package managers such as Anaconda, Miniforge, or Mamba might work, t
 The supported Python versions are 3.11 to 3.13. We will update the documentation when Python3.14 is supported
 (e.g., the ert Python package is not yet available via pip install in Python 3.14).
 
+.. note::
+
+    In Ubuntu, one also needs to install freeglut3-dev:
+
+    .. code-block:: bash
+        
+        sudo apt-get install freeglut3-dev
+
 .. _vpycopm:
 
 Python package

docs/text/installation.rst

@@ -24,7 +24,7 @@ classifiers = [
 ]
 dependencies = [
     "alive_progress<=3.3.0",
-    "ert==16.0.8",
+    "ert@git+https://github.com/equinor/ert#egg=d79d7c9", # Using master since the latest available version (16.0.8) has a few issues 
     "h5py<=3.14.0",
     "mako<=1.3.10",
     "opm<=2025.10"

pyproject.toml

@@ -80,7 +80,7 @@ def pycopm():
             )
     if (
         not isinstance(dic["rcijk"], str)
-        or sum(len(dic[f"{i}ref"]) for i in ["x", "y", "z"]) > 0
+        or np.sum(len(dic[f"{i}ref"]) for i in ["x", "y", "z"]) > 0
     ):
         dic["refinement"] = True
     else:

src/pycopm/core/pycopm.py

@@ -1,10 +1,13 @@
+<%!
+import numpy as np
+%>\
 {
 %for j in range(len(dic["actnum_c"])):
 %if dic["actnum_c"][j] == 1:
 %if j == last:
-"${dic["rock"][i][0]}${sum(dic["actnum_c"][0:j])}": <${dic["rock"][i][0]}${sum(dic["actnum_c"][0:j])}>
+"${dic["rock"][i][0]}${np.sum(dic["actnum_c"][0:j])}": <${dic["rock"][i][0]}${np.sum(dic["actnum_c"][0:j])}>
 %else:
-"${dic["rock"][i][0]}${sum(dic["actnum_c"][0:j])}": <${dic["rock"][i][0]}${sum(dic["actnum_c"][0:j])}>,
+"${dic["rock"][i][0]}${np.sum(dic["actnum_c"][0:j])}": <${dic["rock"][i][0]}${np.sum(dic["actnum_c"][0:j])}>,
 %endif
 %endif
 %endfor

src/pycopm/template_scripts/common/perm.mako

@@ -1,3 +1,6 @@
+<%!
+import numpy as np
+%>\
 #!/usr/bin/env python3
 # SPDX-FileCopyrightText: 2024-2025 NORCE Research AS
 # SPDX-License-Identifier: GPL-3.0
@@ -23,7 +26,7 @@ def perm_evaluation():
     % for k in range(len(dic['actnum_c'])):
     % if dic['rock'][i][1] > 0 and dic['mode'] in ["files","ert"]:
     % if dic['actnum_c'][k] == 1:
-    P[${k}] = max(1e1,min(C["${dic['rock'][i][0]}${sum(dic['actnum_c'][0:k])}"],1e9))
+    P[${k}] = max(1e1,min(C["${dic['rock'][i][0]}${np.sum(dic['actnum_c'][0:k])}"],1e9))
     % endif
     % else:
     P[${k}] = ${k_c[k]}

src/pycopm/template_scripts/common/perm_eval.mako

@@ -1,5 +1,8 @@
+<%!
+import numpy as np
+%>\
 %for j in range(len(dic['actnum_c'])):
 %if dic['actnum_c'][j] == 1:
-${dic["rock"][i][0]}${sum(dic['actnum_c'][0:j])} UNIFORM ${min_max[j][0]} ${min_max[j][1]}
+${dic["rock"][i][0]}${np.sum(dic['actnum_c'][0:j])} UNIFORM ${min_max[j][0]} ${min_max[j][1]}
 %endif
 %endfor

src/pycopm/template_scripts/common/perm_priors.mako

@@ -188,7 +188,7 @@ def visualizeData():
     for type, ftype, dens in zip(["WWPR","WOPR","WGPR"],["FWPT","FOPT","FGPT"],[999.04100, 852.95669, 0.90358]):
         fmass[ftype] = dens*smspec[ftype+"H"][-1]
         for well in nwells:
-            if sum(smspec[type+"H:"+well])>0:
+            if np.sum(smspec[type+"H:"+well])>0:
                 wells[type].append(type+":"+well)
                 j += 1
     linei = [[ ] for _ in range(j)]
@@ -417,14 +417,14 @@ def visualizeData():
             data = smspec[wells[type][i]]
             datah = smspec[wells[type][i][:4]+"H"+wells[type][i][4:]]
             axs[j].plot(smsp_dates, data, color='m', label = 'opm-tests')
-            if sum(datah>0):
+            if np.sum(datah>0):
                 axs[j].errorbar(smsp_dates, datah, yerr= [max(minerr[type],var[type]*d_2) for d_2 in datah], color=[128 / 255.0, 128 / 255.0, 128 / 255.0], markersize='.5', elinewidth=.5, fmt='o', linestyle='', label = 'Data')
             axs[j].set_ylabel(f'{wells[type][i]} [SM3/day]', fontsize=12)
             axs[j].set_xlabel('Time [years]', fontsize=12)
             axs[j].xaxis.set_tick_params(size=6, rotation=45)
             axs[j].legend()
             axs[j].set_ylim(bottom=0)
-            if sum(datah>0):
+            if np.sum(datah>0):
                 figs[j].savefig(f"{output_folder}/postprocessing/wells/HISTO_DATA_{wells[type][i][:4]}_{wells[type][i][5:]}.png", bbox_inches="tight")
             else:
                 figs[j].savefig(f"{output_folder}/postprocessing/wells/{wells[type][i][:4]}_{wells[type][i][5:]}.png", bbox_inches="tight")
@@ -455,7 +455,7 @@ def visualizeData():
     figs, ax = plt.subplots()
     tab20s = matplotlib.colormaps['tab20']
     for i in range(I):
-        ax.plot(i,sum(error_hist[i])/(len(error_hist[i])), markersize='10', marker='o', label=r"$N_{ens}=$"+f"{n_e[i]}", c=tab20s.colors[i])
+        ax.plot(i,np.sum(error_hist[i])/(len(error_hist[i])), markersize='10', marker='o', label=r"$N_{ens}=$"+f"{n_e[i]}", c=tab20s.colors[i])
     ax.axhline(y=error_standard, color="black", ls="--", lw=1, label=f'opm-tests (#Active cells:{numCellsDefault})')
     ax.set_title(r"$O_i=\frac{1}{N_{ens}}\sum_{j}^{N_{e}}O_{i,j}$, "+f"#HM parameters: {len(csvData)}, #Active cells: {numCellsHm}")
     ax.legend()
@@ -578,7 +578,7 @@ def visualizeData():
             f.write(f'Missmatch (standard simulation from opm-test deck): {error_standard : .4e}\n')
             for i in range(I):
                 f.write(f'Iteration {i}; Number of ensembles {int(n_e[i])}\n')
-                f.write(f'Missmatch (mean): {sum(error_hist[i])/(len(error_hist[i])) : .4e} \n')
+                f.write(f'Missmatch (mean): {np.sum(error_hist[i])/(len(error_hist[i])) : .4e} \n')
                 f.write(f'Missmatch (closest realization to all obs): {error_ens[i][eobs[i][0][0]]: .4e}\n')
         print(f'\nClosest final realization to all obs:{eobs[-1][0][0]}')
         print(f'Number of parameters to HM: {len(csvData)}')
@@ -587,7 +587,7 @@ def visualizeData():
         print(f'Missmatch (standard simulation from opm-test deck): {error_standard : .4e}')
         for i in range(I):
             print(f'Iteration {i}; Number of ensembles {int(n_e[i])}')
-            print(f'Missmatch (mean): {sum(error_hist[i])/(len(error_hist[i])) : .4e}')
+            print(f'Missmatch (mean): {np.sum(error_hist[i])/(len(error_hist[i])) : .4e}')
             print(f'Missmatch (closest realization to all obs): {error_ens[i][eobs[i][0][0]] : .4e}')
     else:
         with open(f"{output_folder}/postprocessing/errors.txt", 'w') as f: