Lapack: [버전 3.9.0] slasq3의 무한 루프

이것이 이전 버전에서 작동했다고 생각하십니까? AFAIK 적어도 3.7.0 이후로 slasq 함수는 수정되지 않았습니다. 귀하의 설명에서 TAU가 충분히 작아지지 않기 때문에 반복되고 있습니까?

이것이 이전 버전에서 작동했다고 생각하십니까? AFAIK 적어도 3.7.0 이후로 slasq 함수는 수정되지 않았습니다. 귀하의 설명에서 TAU가 충분히 작아지지 않기 때문에 반복되고 있습니까?

죄송합니다. 지금까지 3.9.0에서만 확인했습니다. 이 버전이 현재 유일한 버전이기 때문입니다.

특히 배열 z 에 포함된 값을 사용하여 이에 대한 완전한 재현자가 ​​있습니까? 나는 당신이 제공한 데이터와 임의의 z로 문제를 재현할 수 없으며 배열 요소 중 하나라도 NaN이 아니더라도 나쁜 일이 발생할 것이라고 상상 합니다.

dmin이 nan이면 올바른 잘못된 값으로 설정해야 합니다. 무한을 막을 것입니다.

나도 이것에 대해 생각했지만, NaN과 비교하는 규칙이 맞다면 이것은 아마도 문제가 아닐 것입니다.

코드를 빌드하기 위해 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 게시물을 기반으로 이 문제를 종료하겠습니다.
@Peter9606 문제를 신고해 주셔서 감사합니다. 문제가 다시 발생하면 다시 문의해 주세요.
줄리앙.