search

BjornNyberg / NetworkGT

The NetworkGT (Network Geometry and Topology) Toolbox is a set of tools designed for the geometrical and topological analysis of fracture networks.
GNU General Public License v3.0
30 stars 3 forks source link

1D fracture network tracer error #39

Open BjornNyberg opened 11 months ago

BjornNyberg commented 11 months ago

Running 1D Flow receives the following error

'integer' not defined from the setmembership script. Changing that line to 'if issubclass(mat.dtype.type, int) and tol < 0.5:' causes a new problem unfortunately,

ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (3,) + inhomogeneous part. from the hyperbolic interface laws.

New errors seem to be in part caused by updated numpy version. As a result, I have updated the porepy module to be compatiable with Numpy 2.0+ in the latest push release found here. This edits seems to make the 2D flow work. However, I have noticed that the tracer seems to go beyond the 0 to 1 value (sometimes 500+) which im not sure the cause of at this time.

alessiofumagalli commented 11 months ago

Dear @BjornNyberg, can I have a minimal example to be able to reproduce the error? Similar for the tracer transport. Thanks