github.com/gonum/lapack@v0.0.0-20181123203213-e4cdc5a0bff9/internal/testdata/dlasqtest/testdlasq4.f90 (about) 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