vendor/gonum.org/v1/gonum/lapack/internal/testdata/dlasqtest/testdlasq3.f90

   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