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

   1 program testdlasq4
   2 implicit none
   3
   4 ! Dlasq4 tests
   5 integer :: i0, n0, pp, ttype, iter, nDiv, n, n0in
   6 real(kind=8), dimension(401) :: z
   7 real(kind=8) :: dmin, desig, qmax, dmin1, dmin2, dn, dn1, dn2, g, tau, sigma
   8 logical :: ieee
   9
  10 z(1:401) = (/2.9266499599402414D+02, &
  11 2.9276699640004830D+02, &
  12 2.9210200041216615D+02, &
  13 2.9171973344922628D+02, &
  14 2.9254120660184350D+02, &
  15 2.9246923723528141D+02, &
  16 2.8739845025761718D+02, &
  17 2.8713781898215319D+02, &
  18 2.9216230331620721D+02, &
  19 2.9216609319427272D+02, &
  20 2.9037912623973870D+02, &
  21 2.9019341222762625D+02, &
  22 2.9184732870824223D+02, &
  23 2.9211597271231784D+02, &
  24 2.9133173969911888D+02, &
  25 2.9116703928303986D+02, &
  26 2.9158395177813247D+02, &
  27 2.9142053084921997D+02, &
  28 2.8936184155474774D+02, &
  29 2.9014388604362898D+02, &
  30 2.9216722380092682D+02, &
  31 2.9216917065163557D+02, &
  32 2.9033857172870131D+02, &
  33 2.9017428768626024D+02, &
  34 2.9187084534700494D+02, &
  35 2.9185392666684271D+02, &
  36 2.8851014988393416D+02, &
  37 2.8849290783118727D+02, &
  38 2.9182506564169046D+02, &
  39 2.9184679532940385D+02, &
  40 2.9022222656964715D+02, &
  41 2.9051820223150668D+02, &
  42 2.9219746119018231D+02, &
  43 2.9219997046052623D+02, &
  44 2.9076912988009974D+02, &
  45 2.9110108536012547D+02, &
  46 2.9226281758828725D+02, &
  47 2.9246883105892533D+02, &
  48 2.9221612201279271D+02, &
  49 2.9059157967819368D+02, &
  50 2.9112833071623868D+02, &
  51 2.9210680427275844D+02, &
  52 2.9199887004058536D+02, &
  53 2.9151030826460186D+02, &
  54 2.9154564759046622D+02, &
  55 2.9212462005167004D+02, &
  56 2.9212108612522559D+02, &
  57 2.9225720737270348D+02, &
  58 2.9226471402909704D+02, &
  59 2.9186949422820902D+02, &
  60 2.9179442772569342D+02, &
  61 2.9025839700680774D+02, &
  62 2.9028281327393216D+02, &
  63 2.9168280284451140D+02, &
  64 2.9168036127921908D+02, &
  65 2.9115346030764749D+02, &
  66 2.9115401102000044D+02, &
  67 2.9212282665919003D+02, &
  68 2.9212277159409672D+02, &
  69 2.9212554020600300D+02, &
  70 2.9212559651286460D+02, &
  71 2.9050907549403615D+02, &
  72 2.9050344542208029D+02, &
  73 2.9211004769628386D+02, &
  74 2.9211127836880860D+02, &
  75 2.9159257985266015D+02, &
  76 2.9158027318883302D+02, &
  77 2.9042716308833894D+02, &
  78 2.9043622253245718D+02, &
  79 2.9134221002346220D+02, &
  80 2.9134130414047053D+02, &
  81 2.9193179671375327D+02, &
  82 2.9193426987096007D+02, &
  83 2.9145366956406895D+02, &
  84 2.9145119646828226D+02, &
  85 2.9144876630819033D+02, &
  86 2.9145122608379586D+02, &
  87 2.9058395601500280D+02, &
  88 2.9055935887314877D+02, &
  89 2.9123982459393329D+02, &
  90 2.9125037059550806D+02, &
  91 2.9111624624209884D+02, &
  92 2.9110570030194424D+02, &
  93 2.9250004994137356D+02, &
  94 2.9254994096967755D+02, &
  95 2.9166071282413282D+02, &
  96 2.9116180260251298D+02, &
  97 2.8862845048087325D+02, &
  98 2.8925662460652461D+02, &
  99 2.9012609983518610D+02, &
 100 2.9010672249354360D+02, &
 101 2.9142999275190090D+02, &
 102 2.9150806548221971D+02, &
 103 2.9078477842534176D+02, &
 104 2.8840517363552317D+02, &
 105 2.8946267131307206D+02, &
 106 2.9189614533792309D+02, &
 107 2.9232617409967031D+02, &
 108 2.9223702224333306D+02, &
 109 2.9149978725880743D+02, &
 110 2.9168777367421643D+02, &
 111 2.9213296087125275D+02, &
 112 2.9211416223585388D+02, &
 113 2.9151431726536367D+02, &
 114 2.9159900781717892D+02, &
 115 2.9215362494632291D+02, &
 116 2.9214515589728342D+02, &
 117 2.9229900293632022D+02, &
 118 2.9231644811476662D+02, &
 119 2.9119079594655176D+02, &
 120 2.9016344223508082D+02, &
 121 2.9040052826233938D+02, &
 122 2.9146756071890462D+02, &
 123 2.9216592215429023D+02, &
 124 2.9212317137116514D+02, &
 125 2.9220217331561315D+02, &
 126 2.9227234439098311D+02, &
 127 2.9210931773273273D+02, &
 128 2.9139146663504778D+02, &
 129 2.9047426145779627D+02, &
 130 2.9113243832972501D+02, &
 131 2.9163901100086372D+02, &
 132 2.9155399887833846D+02, &
 133 2.9145039523229758D+02, &
 134 2.9151950918932886D+02, &
 135 2.9116352441156198D+02, &
 136 2.9094410515950784D+02, &
 137 2.9056083843332891D+02, &
 138 2.9097140423254604D+02, &
 139 2.9185385310744653D+02, &
 140 2.9181279658894493D+02, &
 141 2.9031697123096876D+02, &
 142 2.9033298229128292D+02, &
 143 2.9212507489852788D+02, &
 144 2.9212491479406674D+02, &
 145 2.9123248207180859D+02, &
 146 2.9123278004625678D+02, &
 147 2.9210890357464513D+02, &
 148 2.9210887378334229D+02, &
 149 2.9240010366769212D+02, &
 150 2.9240021314868147D+02, &
 151 2.9212046105624000D+02, &
 152 2.9212035158139264D+02, &
 153 2.9215242607939217D+02, &
 154 2.9215256473001551D+02, &
 155 2.9086373488934720D+02, &
 156 2.9084987044121846D+02, &
 157 2.9162085042968522D+02, &
 158 2.9163097873649599D+02, &
 159 2.9113228736441005D+02, &
 160 2.9112215911901944D+02, &
 161 2.9210393504686402D+02, &
 162 2.9211255232953414D+02, &
 163 2.9191445519642213D+02, &
 164 2.9182828243114074D+02, &
 165 2.9127013482394887D+02, &
 166 2.9129823908392501D+02, &
 167 2.9220880435635905D+02, &
 168 2.9220599393650343D+02, &
 169 2.9144676491963560D+02, &
 170 2.9145286023012017D+02, &
 171 2.9148058174349006D+02, &
 172 2.9147448649442555D+02, &
 173 2.9186619188590822D+02, &
 174 2.9187731897876392D+02, &
 175 2.9151148084374620D+02, &
 176 2.9150035381231066D+02, &
 177 2.9224292392669071D+02, &
 178 2.9224832614828313D+02, &
 179 2.9140236136264468D+02, &
 180 2.9134833920814032D+02, &
 181 2.9241531971951883D+02, &
 182 2.9247972954055604D+02, &
 183 2.9164613863843044D+02, &
 184 2.8920404894780614D+02, &
 185 2.8913485112658549D+02, &
 186 2.9142701775363003D+02, &
 187 2.9213449277295280D+02, &
 188 2.9191925854858397D+02, &
 189 2.9217029610406047D+02, &
 190 2.9224915292105180D+02, &
 191 2.9220325626287905D+02, &
 192 2.9212439945202976D+02, &
 193 2.9057154227612733D+02, &
 194 2.9093384291672965D+02, &
 195 2.9183595231931184D+02, &
 196 2.9179972231667176D+02, &
 197 2.9133916987034183D+02, &
 198 2.9135453536038978D+02, &
 199 2.9047085843212886D+02, &
 200 2.9031720414585078D+02, &
 201 2.9113106705127024D+02, &
 202 2.9119455617800747D+02, &
 203 2.9230127009862616D+02, &
 204 2.9229492119209442D+02, &
 205 2.9098170102485551D+02, &
 206 2.9110028345622749D+02, &
 207 2.9215236246918653D+02, &
 208 2.9215215113995436D+02, &
 209 2.9124650182241129D+02, &
 210 2.9124684419934567D+02, &
 211 2.9216244260160767D+02, &
 212 2.9216240837005626D+02, &
 213 2.9314769094184447D+02, &
 214 2.9314772207133655D+02, &
 215 2.9212215646682506D+02, &
 216 2.9212184517804621D+02, &
 217 2.9224739389310008D+02, &
 218 2.9224802593241526D+02, &
 219 2.8978611958992855D+02, &
 220 2.8915408641678846D+02, &
 221 2.9118533307571403D+02, &
 222 2.9194553410039072D+02, &
 223 2.9244688594134544D+02, &
 224 2.9237086584501975D+02, &
 225 2.9189894464487475D+02, &
 226 2.9210832103555418D+02, &
 227 2.9227685116407457D+02, &
 228 2.9225842459914986D+02, &
 229 2.9224821267514392D+02, &
 230 2.9226591109467802D+02, &
 231 2.9158625025645608D+02, &
 232 2.9140926612253548D+02, &
 233 2.9125298964526195D+02, &
 234 2.9136239316241648D+02, &
 235 2.9182040805585694D+02, &
 236 2.9171100460012258D+02, &
 237 2.9126693393315793D+02, &
 238 2.9130800749967943D+02, &
 239 2.9084722625682167D+02, &
 240 2.9043649120580767D+02, &
 241 2.9153791610078673D+02, &
 242 2.9210440912405335D+02, &
 243 2.9188686486416555D+02, &
 244 2.9138068978583891D+02, &
 245 2.9158124399039616D+02, &
 246 2.9213540811751852D+02, &
 247 2.9211624446790427D+02, &
 248 2.9138960755567365D+02, &
 249 2.9060415959100914D+02, &
 250 2.9118025885063923D+02, &
 251 2.9132902926505932D+02, &
 252 2.9120918643494110D+02, &
 253 2.9222955158764825D+02, &
 254 2.9236106160611405D+02, &
 255 2.9219941546106151D+02, &
 256 2.9167905448737707D+02, &
 257 2.8956733560159137D+02, &
 258 2.9016660738228080D+02, &
 259 2.9069372888915495D+02, &
 260 2.9058385568123452D+02, &
 261 2.9219861470697771D+02, &
 262 2.9223599112669444D+02, &
 263 2.9150234673408067D+02, &
 264 2.9112858259833354D+02, &
 265 2.9011809831967531D+02, &
 266 2.9046138673473865D+02, &
 267 2.9176526297220903D+02, &
 268 2.9173093419212285D+02, &
 269 2.9150067445343569D+02, &
 270 2.9152418894132262D+02, &
 271 2.9197021328944544D+02, &
 272 2.9194669886297868D+02, &
 273 2.9168012546172463D+02, &
 274 2.9169701864791432D+02, &
 275 2.9061916960844030D+02, &
 276 2.9045023836074438D+02, &
 277 2.9212529304266053D+02, &
 278 2.9217230349727680D+02, &
 279 2.9156357220873417D+02, &
 280 2.9093467723991887D+02, &
 281 2.9015918793416472D+02, &
 282 2.9095983823965622D+02, &
 283 2.9124303667219704D+02, &
 284 2.9116297170306802D+02, &
 285 2.9162557939634507D+02, &
 286 2.9193291492830070D+02, &
 287 2.9138731414547760D+02, &
 288 2.9079978674942083D+02, &
 289 2.9159254111122891D+02, &
 290 2.9228695093326087D+02, &
 291 2.9225423965335352D+02, &
 292 2.9126542837357596D+02, &
 293 2.9173807466948557D+02, &
 294 2.9270696227984581D+02, &
 295 2.9274075480094893D+02, &
 296 2.9210759999419366D+02, &
 297 2.9091701924155882D+02, &
 298 2.9163130710260498D+02, &
 299 2.9157158980248664D+02, &
 300 2.9031984685457667D+02, &
 301 2.8939239820067894D+02, &
 302 2.9070124438246870D+02, &
 303 2.9196350995833507D+02, &
 304 2.9189730956351508D+02, &
 305 2.9047434516138509D+02, &
 306 2.9050934070611282D+02, &
 307 2.9179387236503283D+02, &
 308 2.9179037287198020D+02, &
 309 2.9197909575321057D+02, &
 310 2.9198343084479546D+02, &
 311 2.9146450345905629D+02, &
 312 2.9146016842889151D+02, &
 313 2.9233889755268928D+02, &
 314 2.9234209014701219D+02, &
 315 2.9032302536354865D+02, &
 316 2.8837665454569191D+02, &
 317 2.8923259893670223D+02, &
 318 2.9119948081590582D+02, &
 319 2.9158422265927129D+02, &
 320 2.9138706789415261D+02, &
 321 2.9131767302238063D+02, &
 322 2.9147948119882238D+02, &
 323 2.9186937598758146D+02, &
 324 2.9170756787255982D+02, &
 325 2.9181460198203825D+02, &
 326 2.9210008868832023D+02, &
 327 2.9132088639794199D+02, &
 328 2.9113460155819826D+02, &
 329 2.9159448593680958D+02, &
 330 2.9214172380575934D+02, &
 331 2.9182660907964311D+02, &
 332 2.8938570202792107D+02, &
 333 2.8829979954353007D+02, &
 334 2.9064250900128923D+02, &
 335 2.9110347860364448D+02, &
 336 2.9039527564479215D+02, &
 337 2.9152783923893156D+02, &
 338 2.9213818222324863D+02, &
 339 2.9211881276563861D+02, &
 340 2.9133414472425187D+02, &
 341 2.9125375097531855D+02, &
 342 2.9190226937530139D+02, &
 343 2.9164903795105636D+02, &
 344 2.8851961249369930D+02, &
 345 2.8835729910731652D+02, &
 346 2.9144629542849515D+02, &
 347 2.9171485796658368D+02, &
 348 2.9126891989861593D+02, &
 349 2.9051052150188838D+02, &
 350 2.9113570969073521D+02, &
 351 2.9141652614094619D+02, &
 352 2.9133186866181995D+02, &
 353 2.8912484642922516D+02, &
 354 2.8915669391948876D+02, &
 355 2.9117142770169869D+02, &
 356 2.9117110928821620D+02, &
 357 2.9183490795481288D+02, &
 358 2.9183646161652371D+02, &
 359 2.9254379194824975D+02, &
 360 2.9254363658822064D+02, &
 361 2.9192127768678108D+02, &
 362 2.9192154096878102D+02, &
 363 2.9151073810022524D+02, &
 364 2.9151047487964541D+02, &
 365 2.9013993357169915D+02, &
 366 2.9014000572685774D+02, &
 367 2.9145714571416084D+02, &
 368 2.9145713856006512D+02, &
 369 2.9196091751473017D+02, &
 370 2.9196093255283131D+02, &
 371 2.9195731369267685D+02, &
 372 2.9195729871599588D+02, &
 373 2.9149460806187619D+02, &
 374 2.9149461579988633D+02, &
 375 2.9034687539457741D+02, &
 376 2.9034679862867739D+02, &
 377 2.9010765498373922D+02, &
 378 2.9010767881379803D+02, &
 379 2.8967492890608241D+02, &
 380 2.8967469674750663D+02, &
 381 2.9215189060671338D+02, &
 382 2.9215194287124120D+02, &
 383 2.8956132995364590D+02, &
 384 2.8950907156783751D+02, &
 385 2.9210783379750950D+02, &
 386 2.9211890340058585D+02, &
 387 2.9139896954309916D+02, &
 388 2.9128827357375536D+02, &
 389 2.9190012628763964D+02, &
 390 2.9212457714005234D+02, &
 391 2.9221336158434212D+02, &
 392 2.9217879707919576D+02, &
 393 2.9127472708852724D+02, &
 394 2.9132783649002283D+02, &
 395 2.9227509889951233D+02, &
 396 2.9226978796550475D+02, &
 397 2.9142078378447616D+02, &
 398 2.9142906715955786D+02, &
 399 2.9017517054475991D+02, &
 400 2.8992337408144391D+02, &
 401 2.8477362931531121D+02, &
 402 2.8514676291150971D+02, &
 403 2.8393423954726200D+02, &
 404 2.8313625395453965D+02, &
 405 2.8449999402409759D+02, &
 406 2.8534282615303522D+02, &
 407 2.8523140662657408D+02, &
 408 2.8739845025761718D+02, &
 409 2.8713781898215319D+02, &
 410 0.0000000000000000D+00/)
 411 i0 = 1
 412 n0 = 100
 413 pp = 1
 414 n0in = 100
 415 dmin = 7.9798559272231942D+282
 416 dmin1 = 7.9798559272231942D+282
 417 dmin2 = 2.3215857579138407D+285
 418 dn = 2.3140662657409930D+284
 419 dn1 = 7.9798559272231942D+282
 420 dn2 = 8.2833136615500372D+288
 421 g = 2.5000000000000000D-01
 422 tau = 6.1420124051338430D+283
 423 ttype = -4
 424
 425 ieee = .True.
 426
 427 print *, "Starting Dlasq4"
 428
 429 call Dlasq4(i0, n0, z, pp, n0in, dmin, dmin1, dmin2, dn, dn1, dn2, tau, ttype, g)
 430
 431 print *, "Done calling dlasq4"
 432 print *, "tau = ", tau
 433 print *, tau, ttype, g
 434
 435 end