PrincetonUniversity/athena

Array indexing errors in llf_mhd_rel_no_transform.cpp

Closed this issue · 2 comments

The prim_l and prim_r arrays in hydro/rsolvers/mhd/llf_mhd_rel_no_transform.cpp are indexed incorrectly.

Instead of prim_l(IDN,i) they are indexed prim_l(IDN,k,j,i). This causes undefined behavior for the primitive values and in my case nans. The other llf files seem to be correctly indexed.

@smressle can you test my fix in a9cc89e ?