Lapack: [рд╕рдВрд╕реНрдХрд░рдг 3.9.0] slasq3 рдореЗрдВ рдЕрдирдВрдд рд▓реВрдк

рдХреЛ рдирд┐рд░реНрдорд┐рдд 21 рдЬрдире░ 2020  ┬╖  7рдЯрд┐рдкреНрдкрдгрд┐рдпрд╛рдБ  ┬╖  рд╕реНрд░реЛрдд: Reference-LAPACK/lapack

рдЬрдм "slasq3" рдХреЗ рд▓рд┐рдП рдЗрдирдкреБрдЯ рдкреИрд░рд╛рдореАрдЯрд░ рдиреАрдЪреЗ рджрд┐рдП рдЧрдП рдкреИрдЯрд░реНрди рдХрд╛ рдЕрдиреБрд╕рд░рдг рдХрд░рддреЗ рд╣реИрдВ:
i0 =3, n0=16, pp=1,desig=0.2f, qmax=0.2f, nfail=4,iter=1, ndiv = 2, рдпрд╛рдиреА = 0, ttype = -1
рдпрд╣ рдлрд╝рдВрдХреНрд╢рди рдХреЛрдб рд▓рд╛рдЗрди 332 ~ 362 рдХреЗ рдмреАрдЪ рдЕрдирдВрдд рд▓реВрдк рдореЗрдВ рдЖрддрд╛ рд╣реИред
100% рдкреБрдирд░реБрддреНрдкрд╛рджрдиред

рд╕рдмрд╕реЗ рдЙрдкрдпреЛрдЧреА рдЯрд┐рдкреНрдкрдгреА

рдХреНрдпрд╛ рдЖрдкрдХреЗ рдкрд╛рд╕ рдЗрд╕рдХреЗ рд▓рд┐рдП рдПрдХ рдкреВрд░реНрдг рдкреБрдирд░реБрддреНрдкрд╛рджрдХ рд╣реЛрдЧрд╛, рд╡рд┐рд╢реЗрд╖ рд░реВрдк рд╕реЗ рд╕рд░рдгреА z рдореЗрдВ рдирд┐рд╣рд┐рдд рдорд╛рдиреЛрдВ рдХреЗ рд╕рд╛рде? рдореИрдВ рдЖрдкрдХреЗ рджреНрд╡рд╛рд░рд╛ рдкреНрд░рджрд╛рди рдХрд┐рдП рдЧрдП рдбреЗрдЯрд╛ рдФрд░ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ z рдХреЗ рд╕рд╛рде рд╕рдорд╕реНрдпрд╛ рдХреЛ рдкреБрди: рдЙрддреНрдкрдиреНрди рдирд╣реАрдВ рдХрд░ рд╕рдХрддрд╛, рдФрд░ рдореБрдЭреЗ рд▓рдЧрддрд╛ рд╣реИ рдХрд┐ рдмреБрд░реА рдЪреАрдЬреЗрдВ рддрдм рднреА

рд╕рднреА 7 рдЯрд┐рдкреНрдкрдгрд┐рдпрд╛рдБ

рдХреНрдпрд╛ рдЖрдкрдХреЛ рд▓рдЧрддрд╛ рд╣реИ рдХрд┐ рдпрд╣ рдкрд╣рд▓реЗ рдХреЗ рд╕рдВрд╕реНрдХрд░рдгреЛрдВ рдореЗрдВ рдХрд╛рдо рдХрд░рддрд╛ рдерд╛? AFAIK рдХрдо рд╕реЗ рдХрдо 3.7.0 рдХреЗ рдмрд╛рдж рд╕реЗ рдХреЛрдИ рднреА рд╕реНрд▓реИрд╕ рдлрд╝рдВрдХреНрд╢рди рд╕рдВрд╢реЛрдзрд┐рдд рдирд╣реАрдВ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдерд╛ред рдЖрдкрдХреЗ рд╡рд┐рд╡рд░рдг рд╕реЗ рдпрд╣ рд▓реВрдкрд┐рдВрдЧ рд╣реИ рдХреНрдпреЛрдВрдХрд┐ рдЯреАрдПрдпреВ рдХрднреА рднреА рдЫреЛрдЯрд╛ рдирд╣реАрдВ рд╣реЛрддрд╛ рд╣реИ?

рдХреНрдпрд╛ рдЖрдкрдХреЛ рд▓рдЧрддрд╛ рд╣реИ рдХрд┐ рдпрд╣ рдкрд╣рд▓реЗ рдХреЗ рд╕рдВрд╕реНрдХрд░рдгреЛрдВ рдореЗрдВ рдХрд╛рдо рдХрд░рддрд╛ рдерд╛? AFAIK рдХрдо рд╕реЗ рдХрдо 3.7.0 рдХреЗ рдмрд╛рдж рд╕реЗ рдХреЛрдИ рднреА рд╕реНрд▓реИрд╕ рдлрд╝рдВрдХреНрд╢рди рд╕рдВрд╢реЛрдзрд┐рдд рдирд╣реАрдВ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдерд╛ред рдЖрдкрдХреЗ рд╡рд┐рд╡рд░рдг рд╕реЗ рдпрд╣ рд▓реВрдкрд┐рдВрдЧ рд╣реИ рдХреНрдпреЛрдВрдХрд┐ рдЯреАрдПрдпреВ рдХрднреА рднреА рдЫреЛрдЯрд╛ рдирд╣реАрдВ рд╣реЛрддрд╛ рд╣реИ?

рдХреНрд╖рдорд╛ рдХрд░реЗрдВ, рдореИрдВрдиреЗ рдЕрднреА рддрдХ рдХреЗрд╡рд▓ 3.9.0 рдкрд░ рд╕рддреНрдпрд╛рдкрд┐рдд рдХрд┐рдпрд╛ рд╣реИ рдХреНрдпреЛрдВрдХрд┐ рдпрд╣ рдПрдХрдорд╛рддреНрд░ рд╕рдВрд╕реНрдХрд░рдг рд╣реИред

рдХреНрдпрд╛ рдЖрдкрдХреЗ рдкрд╛рд╕ рдЗрд╕рдХреЗ рд▓рд┐рдП рдПрдХ рдкреВрд░реНрдг рдкреБрдирд░реБрддреНрдкрд╛рджрдХ рд╣реЛрдЧрд╛, рд╡рд┐рд╢реЗрд╖ рд░реВрдк рд╕реЗ рд╕рд░рдгреА z рдореЗрдВ рдирд┐рд╣рд┐рдд рдорд╛рдиреЛрдВ рдХреЗ рд╕рд╛рде? рдореИрдВ рдЖрдкрдХреЗ рджреНрд╡рд╛рд░рд╛ рдкреНрд░рджрд╛рди рдХрд┐рдП рдЧрдП рдбреЗрдЯрд╛ рдФрд░ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ z рдХреЗ рд╕рд╛рде рд╕рдорд╕реНрдпрд╛ рдХреЛ рдкреБрди: рдЙрддреНрдкрдиреНрди рдирд╣реАрдВ рдХрд░ рд╕рдХрддрд╛, рдФрд░ рдореБрдЭреЗ рд▓рдЧрддрд╛ рд╣реИ рдХрд┐ рдмреБрд░реА рдЪреАрдЬреЗрдВ рддрдм рднреА

рдпрджрд┐ 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 рдкреЛрд╕реНрдЯ рдХреЗ рдЖрдзрд╛рд░ рдкрд░, рдЖрдЗрдП рд╣рдо рдЗрд╕ рдореБрджреНрджреЗ рдХреЛ рдмрдВрдж рдХрд░рддреЗ рд╣реИрдВред
@ рдкреАрдЯрд░9606 рд╕рдорд╕реНрдпрд╛ рдХреА рд░рд┐рдкреЛрд░реНрдЯ рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП рдзрдиреНрдпрд╡рд╛рдж, рдпрджрд┐ рд╕рдорд╕реНрдпрд╛ рд╡рд╛рдкрд╕ рдЖрддреА рд╣реИ, рддреЛ рдХреГрдкрдпрд╛ рд╣рдорд╛рд░реЗ рдкрд╛рд╕ рд╡рд╛рдкрд╕ рдЖрдПрдВред
рдЬреВрд▓рд┐рдпрдиред

рдХреНрдпрд╛ рдпрд╣ рдкреГрд╖реНрда рдЙрдкрдпреЛрдЧреА рдерд╛?
0 / 5 - 0 рд░реЗрдЯрд┐рдВрдЧреНрд╕