Lapack: DSYEVR은 비직교 벡터를 반환합니다.

Miro님, BUG 126을 GitHub 문제로 추가해주셔서 감사합니다. 이것은 훌륭합니다. Osni( @oamarques )는 #100 작업을 하고 있습니다. 테스트 코드를 공유해 주셔서 감사합니다. 네, LAPACK을 위한 새로운 테스트 스위트를 작업할 계획이 있으며 어떤 불쾌한 매트릭스도 환영합니다. Osni에는 이러한 매트릭스가 많이 있으며 목표는 테스트를 위한 인프라와 함께 이러한 매트릭스 컬렉션을 성장시키는 것입니다. 어쨌든 BUG 126을 GitHub 문제로 추가해주셔서 감사합니다. 이것은 훌륭합니다. 건배, 줄리앙.

이 버그는 해당 소스 파일의 606행에서 시작하는 DSTEMR에서 사용되는 분할 기준과 관련이 있습니다(LAPACK 3.7.0에서와 같이). 입력 변수 TRYRAC를 기반으로 서브루틴 DLARRR은 INFO=1을 반환합니다(이는 DLARRR의 예상 동작과 모순되는 것으로 보입니다). 이는 결과적으로 입력으로 전달된 (특정) 삼대각 행렬의 유해한 분할을 트리거합니다. _DLARRR이 하는 일을 무시하고 "분할 기준 설정" 직전에 IINFO=0을 설정하면 반환된 벡터가 직교합니다. 새 버전의 DLARRF를 사용할 때 문제 100을 참조하십시오. 그러나 이것은 이 버그에 대한 적절한 수정이 아닙니다. DLARRR의 올바른 동작이 무엇인지 파악해야 합니다._)

안녕하세요, 이 버그 해결에 관한 소식이 있습니까?

@oamarques 기본 DSTEMR 대신 DSTEGR을 사용하면 이 문제가 완화됩니까?

@uihsnv : 좋은 질문입니다. STEGR은 더 이상 유지되지 않습니다. 문제가 있습니다. 우리는 STEMR에 대해서만 작업하고 유지하고 있습니다. 그러나 예, STEGR이 이 문제에 대해 무엇을 하는지 보는 것이 좋습니다.

@langou 오! 나는 아무 생각이 없었다!! 내 코드에서 몇 가지 STEGR 인스턴스를 사용하고 있습니다. STEMR을 사용하도록 전환해야 합니까?

또한 @miroi 가 공유한 테스트에 일부 수정을 가했는데, 이를 수락하면 차이점을 볼 수 있습니다.

내 코드에서 몇 가지 STEGR 인스턴스를 사용하고 있습니다. STEMR을 사용하도록 전환해야 합니까?

그렇게 생각 해요. @oamarques 가 곧 이에 대한 후속 조치를 취하기를 바랍니다.

또한 @miroi 가 공유한 테스트에 일부 수정을 가했는데, 이를 수락하면 차이점을 볼 수 있습니다.

수정된 테스트를 @oamarques 와 공유하는 것이 좋습니다.

안녕하세요, 이 버그에 대한 소식이 있습니까? 그리고 왜 "Documentation"이라고 이름을 붙였습니까?

그래도 신선한 lapack은 이 오류를 제공합니다(gfortran 4.9.2로 테스트됨, https://github.com/miroi/lapack-dsyevr-test ).

~~~
DIRAC의 대각화 루틴 테스트 프로그램에 오신 것을 환영합니다!
읽기 N= 9 NZ= 1
대각화 행렬을 읽었습니다.
RS 서브루틴에 대한 배열이 할당됨
Eispack의 RS 고유값:
1 -1.5000000000000062
2 -1.5000000000000002
3 -1.4999999999999982
4 0.49999999999999772
5 0.49999999999999856
6 0.49999999999999994
7 0.50000000000000033
8 0.50000000000000122
9 2.4999999999999969
* 에스팩 RS *
U^{+} A U - eps ?= 0> norm/diag:0.1912D-15 norm/offdiag:0.1008D-15
U^{+} U - I ?= 0> norm/diag:0.2344D-15 norm/offdiag:0.7239D-16U U^{+} - I ?= 0> norm/diag:0.9869D-16 norm/offdiag:0.9320D-16
LAPACK DSYEVR 고유값:
1 -1.5000000000000056
2 -1.5000000000000009
3 -1.4999999999999973
4 0.49999999999999523
5 0.49999999999999867
6 0.49999999999999939
7 0.49999999999999950
8 0.50000000000000100
9 2.4999999999999862
* LAPACK DSYEVR *
U^{+} A U - eps ?= 0> norm/diag:0.1807D-14 norm/offdiag:0.9012D-07
U^{+} U - I ?= 0> norm/diag:0.2591D-15 norm/offdiag:0.1802D-06U U^{+} - I ?= 0> norm/diag:0.7543D-06 norm/offdiag:0.4951D-06
LAPACK DSTEGR 고유값:
1 -1.5000000000000056
2 -1.5000000000000009
3 -1.4999999999999973
4 0.49999999999999523
5 0.49999999999999867
6 0.49999999999999939
7 0.49999999999999950
8 0.50000000000000100
9 2.4999999999999862
* LAPACK DSYTRD+DSTEGR *
U^{+} A U - eps ?= 0> norm/diag:0.1617D+01 norm/offdiag:0.1644D+00
U^{+} U - I ?= 0> norm/diag:0.1234D-15 norm/offdiag:0.1802D-06U U^{+} - I ?= 0> normal/ diag:0.5411D- 06norm/ offdiag:0.2301D-06

~~~

이것을 다시 읽을 때 무한 루프를 수정하려는 시도로 인해 발생한 #100 수정에 달려 있는 것으로 보이며(최근에 내가 저사양 솔루션을 전문적으로 제안하고 있다고 생각함) 모든 것이 무한 루프에 빠졌습니다. 이상적인 솔루션이 무엇인지에 대한 의구심. 3.5.0으로 돌아가서 #100에서 언급한 플래그 변수를 변경하기 전에 이 경우가 어떻게 작동했는지 확인하는 사람이 있습니까?