github.com/gonum/lapack@v0.0.0-20181123203213-e4cdc5a0bff9/internal/testdata/dlasqtest/testdlasq3.f90 (about)

     1  program testdlasq3
     2  implicit none
     3  
     4  ! Dlasq3 tests
     5  integer :: i0, n0, pp,ttype, nFail,iter,nDiv,n
     6  real(kind=8), dimension(84) :: z
     7  real(kind=8) :: dmin,desig,qmax,dmin1,dmin2,dn,dn1,dn2,g,tau,sigma
     8  logical :: ieee
     9  
    10  n = 84
    11  
    12  ieee = .True.
    13  
    14  z(1:84) = (/1.9143652461808209D+00, &
    15  1.5643842977038900D+00, &
    16  2.4933891621438989D-01, &
    17  3.4998094847693051D-01, &
    18  1.3159965131315450D+00, &
    19  1.3638621124906269D+00, &
    20  9.8984666119707590D-02, &
    21  2.0147331685530781D-01, &
    22  6.0239739795872871D-01, &
    23  6.4655447927417942D-01, &
    24  2.2100334106387809D-03, &
    25  5.4827584804256833D-02, &
    26  9.8618572336789667D-02, &
    27  2.4281908107454921D-02, &
    28  4.7563214844548191D-01, &
    29  7.6546697639973529D-02, &
    30  2.5887481436771148D-01, &
    31  6.1277840695087704D-01, &
    32  1.0786113766900041D-01, &
    33  1.2172855586231640D-01, &
    34  6.4428964922552456D-01, &
    35  2.2938358048981550D-01, &
    36  6.2032304866397048D-01, &
    37  5.2276720640470942D-01, &
    38  3.6956606786075852D-01, &
    39  7.6452331847458654D-01, &
    40  5.3788380542522651D-01, &
    41  2.2536579805014259D-01, &
    42  3.5625331812646233D-01, &
    43  8.8204867223354833D-01, &
    44  2.2221324964361449D-02, &
    45  1.2088451318140349D-02, &
    46  1.2750943030216850D+00, &
    47  6.5487468521633574D-01, &
    48  1.6473243548212180D-01, &
    49  6.4244094276971109D-01, &
    50  1.0075305765438660D+00, &
    51  3.2695517365467008D-01, &
    52  3.4538816017831181D-01, &
    53  8.4530783837131718D-01, &
    54  2.6793917191534039D-01, &
    55  4.1167148387782809D-01, &
    56  7.3286777366837230D-01, &
    57  2.0165584821582411D-01, &
    58  8.3608281383074101D-01, &
    59  9.7375794521953263D-01, &
    60  4.8136607095928219D-01, &
    61  5.9519264227958080D-01, &
    62  6.4953705136764595D-01, &
    63  6.7618762481481709D-01, &
    64  2.3254758802226480D-01, &
    65  4.5471549751211121D-01, &
    66  1.9936248028938069D-01, &
    67  3.3218193673422552D-01, &
    68  3.7823189169112570D-01, &
    69  9.9728131577419957D-02, &
    70  9.8304494035037460D-01, &
    71  7.5610809968448423D-01, &
    72  4.4297338640403672D-01, &
    73  6.0516873235701607D-01, &
    74  1.1732795506024030D+00, &
    75  7.1957244803166864D-01, &
    76  5.0355240691445868D-01, &
    77  8.9668048897477137D-01, &
    78  3.0589803950585209D-01, &
    79  6.5888323539286620D-01, &
    80  3.0146344334154529D-01, &
    81  1.5056721102744461D-01, &
    82  1.2894222375675779D+00, &
    83  6.1246453109936005D-01, &
    84  7.5833643057994404D-01, &
    85  9.7842114980976291D-01, &
    86  4.9778147794615712D-01, &
    87  9.9938135774918690D-01, &
    88  2.8414688478625982D-01, &
    89  2.5673655077691432D-01, &
    90  9.2575397942057647D-01, &
    91  5.5092683856146663D-01, &
    92  5.2313556054509897D-05, &
    93  6.5897402564536967D-01, &
    94  2.1178692213810329D-05, &
    95  7.3492248268320243D-05, &
    96  0.0000000000000000D+00, &
    97  0.0000000000000000D+00/)
    98  i0 = 1
    99  n0 = 21
   100  pp = 0
   101  dmin = -0.0000000000000000D+00
   102  sigma = 0.0000000000000000D+00
   103  desig = 0.0000000000000000D+00
   104  qmax = 2.1637041623952107D+00
   105  nFail = 0
   106  iter = 2
   107  nDiv = 40
   108  ttype = 0
   109  dmin1 = 0.0000000000000000D+00
   110  dmin2 = 0.0000000000000000D+00
   111  dn = 0.0000000000000000D+00
   112  dn1 = 0.0000000000000000D+00
   113  dn2 = 0.0000000000000000D+00
   114  g = 0.0000000000000000D+00
   115  tau = 0.0000000000000000D+00
   116  
   117  print *, "Starting Dlasq3"
   118  open(unit = 4, file = "gen4tests.txt")
   119  open(unit = 5, file = "gen5tests.txt")
   120  call Dlasq3(i0, n0, z, pp, dmin, sigma, desig, qmax, nFail, iter, nDiv, ieee, ttype, dmin1, dmin2, dn, dn1, dn2, g, tau)
   121  close(4)
   122  close(5)
   123  
   124  print *, "Done calling dlasq3"
   125  print *, Z(1:n)
   126  
   127  print *, I0, N0, PP
   128  print *, DMIN, DESIG, SIGMA, QMAX
   129  print *, NFAIL, ITER, NDIV, TTYPE
   130  print *, DMIN1, DMIN2, DN, DN1, DN2, G, TAU
   131  
   132  end