Yury Zabegaev

According to the line profiler, 90% of time spent in `Operator.evaluate` (value without Jacobian) is in `MergedOperator.parse`, see [profile.md](https://github.com/pmgbergen/porepy/files/13512635/value.md.txt). Inside, the majority of time is spent in `scipy.sparse.block_diag`. A more...

Renaming of `bc_type_` methods and raising errors on old names

1

Recently, methods for setting a type of boundary condition have been renamed, e.g. `bc_type_mobrho` and `bc_type_darcy` have evolved to `bc_type_fluid_flux` and `bc_type_darcy_flux` respectively. Currently, these methods are named after the...

In the `TimeManager`, there are some public methods and variables that are in risk of being used incorrectly due to their names. Specifically: * `dt_history` and `time_history` - they are...

Make `DenseArray` and `SparseArray` immutable

5

During the discussion with @keileg, the idea was considered to make the array and matrix AD operators immutable. Effectively, this can be achieved by setting the certain flag for the...

Example config

6

As we discussed with @IvarStefansson, it would be nice to see the up-to-date example of the `params` which are used for porepy models, with docstrings explaining the keys and possible...

`PlaneFracture.is_convex()` is slow

3

This performance bottleneck occurs when importing large fracture networks from csv with the `network_3d_from_csv` function. For instance, the case 4 of the 3d benchmark. I can share the code to...