Got AssertionError exception when running solver:

[Yuxiang Lin]

Hi,

I am running an unmodified demo case of dam break.inp but met the assertion error:

=== EXCEPTION ===== EXCEPTION ===== EXCEPTION ===== EXCEPTION ===== EXCEPTION == Traceback:

File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/run.py", line 92, in run_simulation simulation.solver.run() File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solvers/ipcs_algebraic.py", line 634, in run self.slope_limit_velocities() File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/utils/timer.py", line 27, in wrapper ret = f(args, *kwds) File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solvers/ipcs_algebraic.py", line 542, in slope_limit_velocities self.slope_limiter.run() File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solver_parts/slope_limiter_velocity/componentwise_limiter.py", line 97, in run lim.run() File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solver_parts/slope_limiter/hierarchical_taylor.py", line 128, in run boundary_dof_value, File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solver_parts/slope_limiter/hierarchical_taylor.py", line 210, in _run_cpp boundary_dof_type, boundary_dof_value, self.enforce_boundary_conditions File "/workspace/Yuxiang/src/miniconda2/envs/fenicsproject/lib/python3.6/site-packages/ocellaris/solver_parts/slope_limiter/limiter_cpp_utils.py", line 75, in set_boundary_values assert boundary_dof_type.min() >= 0 and boundary_dof_type.max() <= 3

Got AssertionError exception when running solver:

Any idea on how to fix it?

[Tormod Landet]

I somehow did not see this bug report before now. Sorry about that!

Please let me know if you are still looking at running Ocellaris. I have made some changes to this part of the code lately, so it might have been fixed, but the code is still in a bit of a flux since I am working towards finalizing a release (and a Ph.D ...). The demos should all work, but since running them takes some time, only a few are actually tested with every change that is made to the code, the rest are only tested when I happen to want to test something, so things unfortunately break sometimes :-(

[Tormod Landet]

Fixed in master (commit 1a58c1c)