<div dir="ltr"><div>This bug was correctly retitled by the maintainer as related to the AVX-512 kernel.</div><div>However, it is filed against the libopenblas0-pthread debian package while I have found another case that gives wrong results even using one thread.</div><div>The script solves a linear elasticity problem using the finite element method in octave, the files are shared here: <a href="https://drive.google.com/file/d/1xkw-1YfX3W-mfqD92I1fNYnQM78yA9Sy/view?usp=sharing">https://drive.google.com/file/d/1xkw-1YfX3W-mfqD92I1fNYnQM78yA9Sy/view?usp=sharing</a></div><div>After building the system of linear equations, it is solved using the octave operator "\", and by "lu" factorization, the results are checked by computing the residual vectors (should have a norm around the truncation error).</div><div><br></div><div>These are the results of the tests:</div><div>-Using defaults:<br></div><div>$ octave fem_lame2d.m <br>Residualbackslash = 17.947<br>ResidualLUPQbackslash = 7.2444<br></div><div><br></div><div>-Forcing sequential mode (1 thread):</div><div>$ OMP_NUM_THREADS=1 octave fem_lame2d.m <br>Residualbackslash = 17.947<br>ResidualLUPQbackslash = 7.2444</div><div><br></div><div>-Avoiding the use of the openblas AVX-512 kernel (working workaround !!):</div><div>$ OPENBLAS_CORETYPE=Haswell octave fem_lame2d.m <br>Residualbackslash = 1.4547e-13<br>ResidualLUPQbackslash = 1.4879e-13</div><div><br></div><div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><span>Enzo A. Dari<br></span></div><div>Profesor Titular</div><div><a href="http://www.ib.edu.ar/" target="_blank">Instituto Balseiro</a><br></div></div></div></div></div>