SimplePETSc

This benchmark is available as a Jupyter notebook: Notebooks/SimplePETSc.ipynb.

The following shows running a simple steady-state diffusion benchmark running on 2 cores.

import os

prj_name = "square_1e1_neumann"
data_dir = os.environ.get('OGS_DATA_DIR', '../../../Data')
prj_file = f"{data_dir}/EllipticPETSc/{prj_name}.prj"
out_dir = os.environ.get('OGS_TESTRUNNER_OUT_DIR', '_out')

if not os.path.exists(out_dir):
    os.makedirs(out_dir)
os.chdir(out_dir)

print(f"mpirun -np 2 ogs {prj_file} > out.txt")
! mpirun -np 2 ogs {prj_file} > out.txt

from datetime import datetime
print(datetime.now())
mpirun -np 2 ogs /home/jovyan/work/Tests/Data/EllipticPETSc/square_1e1_neumann.prj > out.txt
2022-05-19 16:54:01.262938
import vtuIO

pvdfile = vtuIO.PVDIO(f"{prj_name}.pvd", dim=2)
time = pvdfile.timesteps
points={'pt0': (0.3,0.5,0.0), 'pt1': (0.24,0.21,0.0)}
pressure_linear = pvdfile.read_time_series("pressure", points)

import matplotlib.pyplot as plt
plt.plot(time, pressure_linear["pt0"], "b-", label="pt0 linear interpolated")
plt.plot(time, pressure_linear["pt1"], "r-", label="pt1 linear interpolated")
plt.legend()
plt.xlabel("t")
plt.ylabel("p")
Text(0, 0.5, 'p')

png


This article was written by Lars Bilke. If you are missing something or you find an error please let us know. Generated with Hugo 0.96.0. Last revision: November 9, 2021