Lapack: [Versi 3.9.0] Loop tak terbatas di slasq3

Dibuat pada 21 Jan 2020 · 7Komentar · Sumber: Reference-LAPACK/lapack

Saat memasukkan parameter ke "slasq3" mengikuti pola di bawah ini:
i0 =3, n0=16, pp=1,desain=0.2f, qmax=0.2f, nfail=4,iter=1, ndiv = 2, ieee = 0, ttype = -1
fungsi ini jatuh ke loop tak terbatas di antara baris kode 332 ~ 362.
100% reproduksi.

Sumber

Peter9606

Komentar yang paling membantu

Apakah Anda kebetulan memiliki reproduksi lengkap untuk ini, khususnya dengan nilai-nilai yang terkandung dalam array z ? Saya tidak dapat mereproduksi masalah dengan data yang Anda berikan dan z acak, ~~dan saya membayangkan hal-hal buruk akan terjadi jika~~ tidak bahkan ketika salah satu elemen array adalah NaN.

martin-frbg pada 1 Feb 2020

👍2

Semua 7 komentar

Apakah menurut Anda ini berfungsi di versi sebelumnya? AFAIK tidak ada fungsi slasq yang dimodifikasi setidaknya sejak 3.7.0. Dari deskripsi Anda, itu berulang karena TAU tidak pernah cukup kecil?

martin-frbg pada 26 Jan 2020

Apakah menurut Anda ini berfungsi di versi sebelumnya? AFAIK tidak ada fungsi slasq yang dimodifikasi setidaknya sejak 3.7.0. Dari deskripsi Anda, itu berulang karena TAU tidak pernah cukup kecil?

Maaf, saya hanya memverifikasi pada 3.9.0 sejauh ini karena ini adalah satu-satunya versi yang ada.

Peter9606 pada 31 Jan 2020

martin-frbg pada 1 Feb 2020

👍2

jika dmin adalah nan, itu harus disetel ke nilai salah yang sah; akan mencegah tak terhingga.

moe123 pada 2 Feb 2020

Saya juga memikirkan hal ini, tetapi jika saya benar tentang aturan untuk membandingkan dengan NaN, ini mungkin bukan masalahnya.

martin-frbg pada 2 Feb 2020

Tidak ada loop tak terbatas menggunakan 4b3c7c249f6a2cfed76603d9987cf5164fa6f9e1 untuk membuat kode:

      PROGRAM testSLASQ3
      IMPLICIT NONE

*     .. Scalar Arguments ..
      LOGICAL            IEEE
      INTEGER            I0, ITER, N0, NDIV, NFAIL, PP, TTYPE
      REAL               DESIG, DMIN, DMIN1, DMIN2, DN, DN1, DN2, G,
     &                   QMAX, SIGMA, TAU
*     .. Array Arguments ..
      REAL               Z( 64 )
*     .. Subroutines ..
      EXTERNAL           SLASQ3

      PARAMETER ( IEEE=.FALSE. )
      PARAMETER ( I0 = 3 )
      PARAMETER ( QMAX = 1.0 )

      ITER=1
      N0=16
      NDIV=2
      NFAIL=4
      PP=1
      TTYPE=-1
      DESIG=0.2
      DMIN1=0.1
      DMIN2=0.1
      DN=1.0
      DN1=1.0
      DN2=1.0
      G=1.0
      TAU=1.0

      Z = (/0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,
     &      0.1,0.2,0.3,0.4,0.5,0.6,
     &      0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,
     &      0.1,0.2,0.3,0.4,0.5,0.6,
     &      0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,
     &      0.1,0.2,0.3,0.4,0.5,0.6,
     &      0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,
     &      0.1,0.2,0.3,0.4,0.5,0.6/)

      CALL SLASQ3( I0, N0, Z, PP, DMIN, SIGMA, DESIG, QMAX, NFAIL,
     &                   ITER, NDIV, IEEE, TTYPE, DMIN1, DMIN2, DN, DN1,
     &                   DN2, G, TAU )

      END PROGRAM testSLASQ3

weslleyspereira pada 22 Apr 2021

👍1

Berdasarkan postingan @weslleyspereira , mari kita tutup masalah ini.
@ Peter9606 terima kasih telah melaporkan masalah ini, jika masalah muncul kembali, silakan kembali kepada kami.
Julien.