Lapack: [الإصدار 3.9.0] الحلقة اللانهائية في slasq3

تم إنشاؤها على ٢١ يناير ٢٠٢٠ · 7تعليقات · مصدر: Reference-LAPACK/lapack

عندما تتبع معلمات الإدخال إلى "slasq3" النمط التالي:
i0 = 3، n0 = 16، pp = 1، desig = 0.2f، qmax = 0.2f، nfail = 4، iter = 1، ndiv = 2، ieee = 0، ttype = -1
تقع هذه الوظيفة في حلقة لا نهائية بين سطر الكود 332 ~ 362.
100٪ تتكاثر.

مصدر

Peter9606

التعليق الأكثر فائدة

هل تصادف أن يكون لديك مُنسخ كامل لهذا ، خاصة مع القيم الموجودة في المصفوفة z ؟ لا يمكنني إعادة إنتاج المشكلة بالبيانات التي قدمتها و z العشوائية ، ~~وأتخيل حدوث أشياء سيئة إن~~ لم يكن حتى عندما يكون أي من عناصر المصفوفة هو NaN.

martin-frbg في ١ فبراير ٢٠٢٠

👍2

ال 7 كومينتر

هل تعتقد أن هذا كان يعمل في الإصدارات السابقة؟ AFAIK لم يتم تعديل أي من وظائف slasq منذ 3.7.0 على الأقل. من وصفك أنه يتكرر حيث أن TAU لا تصبح صغيرة بما يكفي؟

martin-frbg في ٢٦ يناير ٢٠٢٠

هل تعتقد أن هذا كان يعمل في الإصدارات السابقة؟ AFAIK لم يتم تعديل أي من وظائف slasq منذ 3.7.0 على الأقل. من وصفك أنه يتكرر حيث أن TAU لا تصبح صغيرة بما يكفي؟

عذرًا ، لقد قمت بالتحقق من الإصدار 3.9.0 فقط حتى الآن لأن هذا هو الإصدار الوحيد المتاح.

Peter9606 في ٣١ يناير ٢٠٢٠

martin-frbg في ١ فبراير ٢٠٢٠

👍2

إذا كانت dmin هي nan ، فيجب ضبطها على قيمة خاطئة شرعية ؛ سيمنع اللانهاية.

moe123 في ٢ فبراير ٢٠٢٠

فكرت في هذا أيضًا ، ولكن إذا كنت محقًا بشأن قواعد المقارنة مع NaN ، فربما لا تكون هذه هي المشكلة.

martin-frbg في ٢ فبراير ٢٠٢٠

لا توجد حلقة لانهائية باستخدام 4b3c7c249f6a2cfed76603d9987cf5164fa6f9e1 لبناء الكود:

      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 في ٢٢ أبريل ٢٠٢١

👍1

بناءً على مشاركة weslleyspereira @ ، دعنا نغلق هذه المشكلة بعد ذلك.
@ Peter9606 شكرًا على الإبلاغ عن المشكلة ، إذا عادت المشكلة ، يرجى الرجوع إلينا.
جوليان.