LAPACK test segfault on zen/zen2/zen3 at bli_sgemmsup_rd_haswell_asm_1x16n about blis HOT 16 CLOSED

$ cat sgemmtest.f90
program sgemmtest
   IMPLICIT NONE

   REAL, ALLOCATABLE ::  Q(:, :), A(:, :), R(:, :)

   REAL ONE, ZERO
   PARAMETER(ONE=1.0, ZERO=0.0)

   INTEGER L, M, N

   INTRINSIC MAX, MIN

   M = 50
   N = 10

   L = MAX(M, N, 1)

   ALLOCATE (Q(L, L), A(M, N), R(M, L))

   CALL SLASET('A', M, N, ONE, ONE, A, M)
   CALL SLASET('A', L, L, ONE, ONE, Q, L)

   print *,'sgemmtest:'
   print *,'   M = ',M,' N = ',N,' L = ',L
   print *,'   A(1,1) is ',A(1,1),' Q(1,1) is ',Q(1,1)
   print *,' '
   print *,' R = Q**T * A, except Q is square but we use MxN of it'
   print *,' assert sum(A)==M*N*ONE is ', M*N*ONE == SUM(A)
   print *,' assert sum(Q)==L*L*ONE is ', L*L*ONE == SUM(Q)
   CALL SGEMM('T', 'N', M, N, M, ONE, Q, M, A, M, ZERO, R, M)
   print *,' r11 ',r(1,1), ' r211', r(2,1), ' rml ',r(M,L)
   print *,' matrix of MxN filled with M:'
   print *,' assert sum(R)==M*N*M is ', SUM(R) - M*N*M*ONE == ZERO

   print *,' done'
end

$ export GFORTRAN_UNBUFFERED_ALL=1
$ export MALLOC_OPTIONS=CFG
$ export BLIS_ARCH_DEBUG=1
$ ./tblis.x
 sgemmtest:
    M =           50  N =           10  L =           50
    A(1,1) is    1.00000000      Q(1,1) is    1.00000000

  R = Q**T * A, except Q is square but we use MxN of it
  assert sum(A)==M*N*ONE is  T
  assert sum(Q)==L*L*ONE is  T
libblis: selecting sub-configuration 'zen3'.
  r11    50.0000000      r211   50.0000000      rml    0.00000000
  matrix of MxN filled with M:
  assert sum(R)==M*N*M is  T
  done

$ egdb -q tblis.x
Reading symbols from tblis.x...
(gdb) run
Starting program: /home/jal/checkblis/tblis.x
 sgemmtest:
    M =           50  N =           10  L =           50
    A(1,1) is    1.00000000      Q(1,1) is    1.00000000

  R = Q**T * A, except Q is square but we use MxN of it
  assert sum(A)==M*N*ONE is  T
  assert sum(Q)==L*L*ONE is  T
libblis: selecting sub-configuration 'zen3'.

Program received signal SIGSEGV, Segmentation fault.
0x00000aad13098881 in bli_sgemmsup_rd_haswell_asm_1x16n ()
   from /usr/local/lib/libblis.so.0.0

(gdb) bt
#0  0x00000aad13098881 in bli_sgemmsup_rd_haswell_asm_1x16n ()
   from /usr/local/lib/libblis.so.0.0
#1  0x00000aad1309832a in bli_sgemmsup_rd_haswell_asm_6x16n ()
   from /usr/local/lib/libblis.so.0.0
#2  0x00000aad136fac87 in bli_gemmsup_ref_var1n ()
   from /usr/local/lib/libblis.so.0.0
#3  0x00000aad136f8a55 in bli_gemmsup_int ()
   from /usr/local/lib/libblis.so.0.0
#4  0x00000aad136f86aa in bli_l3_sup_thread_decorator_entry ()
   from /usr/local/lib/libblis.so.0.0
#5  0x00000aad136f85d7 in bli_l3_sup_thread_decorator ()
   from /usr/local/lib/libblis.so.0.0
#6  0x00000aad136f9f7b in bli_gemmsup_ref ()
   from /usr/local/lib/libblis.so.0.0
#7  0x00000aad136f82db in bli_gemmsup () from /usr/local/lib/libblis.so.0.0
#8  0x00000aad136f6a1b in bli_gemm_ex () from /usr/local/lib/libblis.so.0.0
#9  0x00000aad137464e3 in sgemm_ () from /usr/local/lib/libblis.so.0.0
#10 0x00000aab12ff8ca5 in sgemmtest () at sgemmtest.f90:30
#11 0x00000aab12ff9059 in main (argc=1, argv=0x7571737f3dd0)
    at sgemmtest.f90:36
#12 0x00000aab12ff7e7b in _start ()

o75snap$ egfortran --version
GNU Fortran (GCC) 8.4.0
Copyright (C) 2018 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

o75snap$ cc --version
OpenBSD clang version 16.0.6
Target: amd64-unknown-openbsd7.5
Thread model: posix
InstalledDir: /usr/bin

$ gmake showconfig
configuration family:       x86_64
sub-configurations:         skx knl haswell sandybridge penryn zen3 zen2 zen excavator steamroller piledriver bulldozer generic
requisite kernels sets:     skx knl sandybridge penryn zen3 zen2 haswell zen piledriver bulldozer generic
kernel-to-config map:       bulldozer:bulldozer generic:generic haswell:haswell knl:knl penryn:penryn piledriver:piledriver sandybridge:sandybridge skx:skx zen:zen zen2:zen2 zen3:zen3
-------------------------
BLIS version string:        1.0
.so major version:          0
.so minor.build vers:       0
install libdir:             /usr/local/lib
install includedir:         /usr/local/include
install sharedir:           /usr/local/share
debugging status:           off
enable AddressSanitizer?    no
enabled threading model(s): single
enable BLAS API?            yes
enable CBLAS API?           yes
build static library?       yes
build shared library?       yes
ARG_MAX hack enabled?       no

24         print *,'malloctest:'
25         write (*,'(A,Z16)') '  location A is ', loc(A)
26         write (*,'(A,Z16)') '  location Q is ', loc(Q)
27         write (*,'(A,Z16)') '  location R is ', loc(R)
28         CALL SGEMM('T', 'N', M, N, M, ONE, Q, M, A, M, ZERO, R, M)

(gdb) break *bli_sgemmsup_rd_haswell_asm_1x16n+721
Breakpoint 3 at 0xd826ecae9a1
(gdb) c
Continuing.
 malloctest:
  location A is      D82E68897D0
  location Q is      D8214FBA000
  location R is      D822BE43000
libblis: selecting sub-configuration 'zen3'.

Breakpoint 3, 0x00000d826ecae9a1 in bli_sgemmsup_rd_haswell_asm_1x16n ()
   from /usr/local/lib/libblis.so.0.0
(gdb) x/i $pc
=> 0xd826ecae9a1 <bli_sgemmsup_rd_haswell_asm_1x16n+721>:
    vmovss (%rax,%r8,1),%xmm1
(gdb) info reg r8 rax rsp
r8             0xc8                200
rax            0xd82e6889f98       14855864623000
rsp            0x734c53fa1830      0x734c53fa1830
(gdb) p 0xd82e6889f98 - 0xD82E68897D0
$1 = 1992

(gdb) c
Continuing.

Program received signal SIGSEGV, Segmentation fault.
0x00000d826ecae9a1 in bli_sgemmsup_rd_haswell_asm_1x16n ()
   from /usr/local/lib/libblis.so.0.0
(gdb) p/x *(0xD82E68897D0)
$2 = 0x3f800000
(gdb) p/x *($rax)
$3 = 0x3f800000

(gdb) p/x *($rax+$r8)
Cannot access memory at address 0xd82e688a060