128 LOGICAL LSAME, LSAMEN
143 IF( .NOT.( sord .OR. corz ) )
146 IF(
lsamen( 2, p2,
'GE' ) )
THEN
150 WRITE( iounit, fmt = 9999 )path
151 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
152 WRITE( iounit, fmt = 9979 )
153 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
154 WRITE( iounit, fmt = 9962 )1
155 WRITE( iounit, fmt = 9961 )2
156 WRITE( iounit, fmt = 9960 )3
157 WRITE( iounit, fmt = 9959 )4
158 WRITE( iounit, fmt = 9958 )5
159 WRITE( iounit, fmt = 9957 )6
160 WRITE( iounit, fmt = 9956 )7
161 WRITE( iounit, fmt = 9955 )8
162 WRITE( iounit, fmt =
'( '' Messages:'' )' )
164 ELSE IF(
lsamen( 2, p2,
'GB' ) )
THEN
168 WRITE( iounit, fmt = 9998 )path
169 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
170 WRITE( iounit, fmt = 9978 )
171 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
172 WRITE( iounit, fmt = 9962 )1
173 WRITE( iounit, fmt = 9960 )2
174 WRITE( iounit, fmt = 9959 )3
175 WRITE( iounit, fmt = 9958 )4
176 WRITE( iounit, fmt = 9957 )5
177 WRITE( iounit, fmt = 9956 )6
178 WRITE( iounit, fmt = 9955 )7
179 WRITE( iounit, fmt =
'( '' Messages:'' )' )
181 ELSE IF(
lsamen( 2, p2,
'GT' ) )
THEN
185 WRITE( iounit, fmt = 9997 )path
186 WRITE( iounit, fmt = 9977 )
187 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
188 WRITE( iounit, fmt = 9962 )1
189 WRITE( iounit, fmt = 9960 )2
190 WRITE( iounit, fmt = 9959 )3
191 WRITE( iounit, fmt = 9958 )4
192 WRITE( iounit, fmt = 9957 )5
193 WRITE( iounit, fmt = 9956 )6
194 WRITE( iounit, fmt = 9955 )7
195 WRITE( iounit, fmt =
'( '' Messages:'' )' )
197 ELSE IF(
lsamen( 2, p2,
'PO' ) .OR.
lsamen( 2, p2,
'PP' ) )
THEN
207 IF(
lsame( c3,
'O' ) )
THEN
208 WRITE( iounit, fmt = 9996 )path, sym
210 WRITE( iounit, fmt = 9995 )path, sym
212 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
213 WRITE( iounit, fmt = 9975 )path
214 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
215 WRITE( iounit, fmt = 9954 )1
216 WRITE( iounit, fmt = 9961 )2
217 WRITE( iounit, fmt = 9960 )3
218 WRITE( iounit, fmt = 9959 )4
219 WRITE( iounit, fmt = 9958 )5
220 WRITE( iounit, fmt = 9957 )6
221 WRITE( iounit, fmt = 9956 )7
222 WRITE( iounit, fmt = 9955 )8
223 WRITE( iounit, fmt =
'( '' Messages:'' )' )
225 ELSE IF(
lsamen( 2, p2,
'PS' ) )
THEN
234 IF(
lsame( c1,
'S' ) .OR.
lsame( c1,
'C' ) )
THEN
239 WRITE( iounit, fmt = 9995 )path, sym
240 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
241 WRITE( iounit, fmt = 8973 )eigcnm, eigcnm, eigcnm
242 WRITE( iounit, fmt =
'( '' Difference:'' )' )
243 WRITE( iounit, fmt = 8972 )c1
244 WRITE( iounit, fmt =
'( '' Test ratio:'' )' )
245 WRITE( iounit, fmt = 8950 )
246 WRITE( iounit, fmt =
'( '' Messages:'' )' )
247 ELSE IF(
lsamen( 2, p2,
'PB' ) )
THEN
252 WRITE( iounit, fmt = 9994 )path,
'Symmetric'
254 WRITE( iounit, fmt = 9994 )path,
'Hermitian'
256 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
257 WRITE( iounit, fmt = 9973 )path
258 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
259 WRITE( iounit, fmt = 9954 )1
260 WRITE( iounit, fmt = 9960 )2
261 WRITE( iounit, fmt = 9959 )3
262 WRITE( iounit, fmt = 9958 )4
263 WRITE( iounit, fmt = 9957 )5
264 WRITE( iounit, fmt = 9956 )6
265 WRITE( iounit, fmt = 9955 )7
266 WRITE( iounit, fmt =
'( '' Messages:'' )' )
268 ELSE IF(
lsamen( 2, p2,
'PT' ) )
THEN
273 WRITE( iounit, fmt = 9993 )path,
'Symmetric'
275 WRITE( iounit, fmt = 9993 )path,
'Hermitian'
277 WRITE( iounit, fmt = 9976 )
278 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
279 WRITE( iounit, fmt = 9952 )1
280 WRITE( iounit, fmt = 9960 )2
281 WRITE( iounit, fmt = 9959 )3
282 WRITE( iounit, fmt = 9958 )4
283 WRITE( iounit, fmt = 9957 )5
284 WRITE( iounit, fmt = 9956 )6
285 WRITE( iounit, fmt = 9955 )7
286 WRITE( iounit, fmt =
'( '' Messages:'' )' )
288 ELSE IF(
lsamen( 2, p2,
'SY' ) )
THEN
293 IF(
lsame( c3,
'Y' ) )
THEN
294 WRITE( iounit, fmt = 9992 )path,
'Symmetric'
296 WRITE( iounit, fmt = 9991 )path,
'Symmetric'
298 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
300 WRITE( iounit, fmt = 9972 )
302 WRITE( iounit, fmt = 9971 )
304 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
305 WRITE( iounit, fmt = 9953 )1
306 WRITE( iounit, fmt = 9961 )2
307 WRITE( iounit, fmt = 9960 )3
308 WRITE( iounit, fmt = 9960 )4
309 WRITE( iounit, fmt = 9959 )5
310 WRITE( iounit, fmt = 9958 )6
311 WRITE( iounit, fmt = 9956 )7
312 WRITE( iounit, fmt = 9957 )8
313 WRITE( iounit, fmt = 9955 )9
314 WRITE( iounit, fmt =
'( '' Messages:'' )' )
316 ELSE IF(
lsamen( 2, p2,
'SR' ) .OR.
lsamen( 2, p2,
'SK') )
THEN
327 WRITE( iounit, fmt = 9892 )path,
'Symmetric'
329 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
331 WRITE( iounit, fmt = 9972 )
333 WRITE( iounit, fmt = 9971 )
336 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
337 WRITE( iounit, fmt = 9953 )1
338 WRITE( iounit, fmt = 9961 )2
339 WRITE( iounit, fmt = 9927 )3
340 WRITE( iounit, fmt = 9928 )
341 WRITE( iounit, fmt = 9926 )4
342 WRITE( iounit, fmt = 9928 )
343 WRITE( iounit, fmt = 9960 )5
344 WRITE( iounit, fmt = 9959 )6
345 WRITE( iounit, fmt = 9955 )7
346 WRITE( iounit, fmt =
'( '' Messages:'' )' )
348 ELSE IF(
lsamen( 2, p2,
'SP' ) )
THEN
353 IF(
lsame( c3,
'Y' ) )
THEN
354 WRITE( iounit, fmt = 9992 )path,
'Symmetric'
356 WRITE( iounit, fmt = 9991 )path,
'Symmetric'
358 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
360 WRITE( iounit, fmt = 9972 )
362 WRITE( iounit, fmt = 9971 )
364 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
365 WRITE( iounit, fmt = 9953 )1
366 WRITE( iounit, fmt = 9961 )2
367 WRITE( iounit, fmt = 9960 )3
368 WRITE( iounit, fmt = 9959 )4
369 WRITE( iounit, fmt = 9958 )5
370 WRITE( iounit, fmt = 9956 )6
371 WRITE( iounit, fmt = 9957 )7
372 WRITE( iounit, fmt = 9955 )8
373 WRITE( iounit, fmt =
'( '' Messages:'' )' )
375 ELSE IF(
lsamen( 2, p2,
'HA' ) )
THEN
380 WRITE( iounit, fmt = 9992 )path,
'Hermitian'
382 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
383 WRITE( iounit, fmt = 9972 )
385 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
386 WRITE( iounit, fmt = 9953 )1
387 WRITE( iounit, fmt = 9961 )2
388 WRITE( iounit, fmt = 9960 )3
389 WRITE( iounit, fmt = 9960 )4
390 WRITE( iounit, fmt = 9959 )5
391 WRITE( iounit, fmt = 9958 )6
392 WRITE( iounit, fmt = 9956 )7
393 WRITE( iounit, fmt = 9957 )8
394 WRITE( iounit, fmt = 9955 )9
395 WRITE( iounit, fmt =
'( '' Messages:'' )' )
397 ELSE IF(
lsamen( 2, p2,
'HE' ) )
THEN
402 WRITE( iounit, fmt = 9992 )path,
'Hermitian'
404 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
405 WRITE( iounit, fmt = 9972 )
407 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
408 WRITE( iounit, fmt = 9953 )1
409 WRITE( iounit, fmt = 9961 )2
410 WRITE( iounit, fmt = 9960 )3
411 WRITE( iounit, fmt = 9960 )4
412 WRITE( iounit, fmt = 9959 )5
413 WRITE( iounit, fmt = 9958 )6
414 WRITE( iounit, fmt = 9956 )7
415 WRITE( iounit, fmt = 9957 )8
416 WRITE( iounit, fmt = 9955 )9
417 WRITE( iounit, fmt =
'( '' Messages:'' )' )
419 ELSE IF(
lsamen( 2, p2,
'HR' ) .OR.
lsamen( 2, p2,
'HR' ) )
THEN
430 WRITE( iounit, fmt = 9892 )path,
'Hermitian'
432 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
433 WRITE( iounit, fmt = 9972 )
435 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
436 WRITE( iounit, fmt = 9953 )1
437 WRITE( iounit, fmt = 9961 )2
438 WRITE( iounit, fmt = 9927 )3
439 WRITE( iounit, fmt = 9928 )
440 WRITE( iounit, fmt = 9926 )4
441 WRITE( iounit, fmt = 9928 )
442 WRITE( iounit, fmt = 9960 )5
443 WRITE( iounit, fmt = 9959 )6
444 WRITE( iounit, fmt = 9955 )7
445 WRITE( iounit, fmt =
'( '' Messages:'' )' )
447 ELSE IF(
lsamen( 2, p2,
'HP' ) )
THEN
452 IF(
lsame( c3,
'E' ) )
THEN
453 WRITE( iounit, fmt = 9992 )path,
'Hermitian'
455 WRITE( iounit, fmt = 9991 )path,
'Hermitian'
457 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
458 WRITE( iounit, fmt = 9972 )
459 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
460 WRITE( iounit, fmt = 9953 )1
461 WRITE( iounit, fmt = 9961 )2
462 WRITE( iounit, fmt = 9960 )3
463 WRITE( iounit, fmt = 9959 )4
464 WRITE( iounit, fmt = 9958 )5
465 WRITE( iounit, fmt = 9956 )6
466 WRITE( iounit, fmt = 9957 )7
467 WRITE( iounit, fmt = 9955 )8
468 WRITE( iounit, fmt =
'( '' Messages:'' )' )
470 ELSE IF(
lsamen( 2, p2,
'TR' ) .OR.
lsamen( 2, p2,
'TP' ) )
THEN
475 IF(
lsame( c3,
'R' ) )
THEN
476 WRITE( iounit, fmt = 9990 )path
477 subnam = path( 1: 1 ) //
'LATRS'
479 WRITE( iounit, fmt = 9989 )path
480 subnam = path( 1: 1 ) //
'LATPS'
482 WRITE( iounit, fmt = 9966 )path
483 WRITE( iounit, fmt = 9965 )subnam(1:len_trim( subnam ))
484 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
485 WRITE( iounit, fmt = 9961 )1
486 WRITE( iounit, fmt = 9960 )2
487 WRITE( iounit, fmt = 9959 )3
488 WRITE( iounit, fmt = 9958 )4
489 WRITE( iounit, fmt = 9957 )5
490 WRITE( iounit, fmt = 9956 )6
491 WRITE( iounit, fmt = 9955 )7
492 WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 8
493 WRITE( iounit, fmt =
'( '' Messages:'' )' )
495 ELSE IF(
lsamen( 2, p2,
'TB' ) )
THEN
499 WRITE( iounit, fmt = 9988 )path
500 subnam = path( 1: 1 ) //
'LATBS'
501 WRITE( iounit, fmt = 9964 )path
502 WRITE( iounit, fmt = 9963 )subnam(1:len_trim( subnam ))
503 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
504 WRITE( iounit, fmt = 9960 )1
505 WRITE( iounit, fmt = 9959 )2
506 WRITE( iounit, fmt = 9958 )3
507 WRITE( iounit, fmt = 9957 )4
508 WRITE( iounit, fmt = 9956 )5
509 WRITE( iounit, fmt = 9955 )6
510 WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 7
511 WRITE( iounit, fmt =
'( '' Messages:'' )' )
513 ELSE IF(
lsamen( 2, p2,
'QR' ) )
THEN
517 WRITE( iounit, fmt = 9987 )path,
'QR'
518 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
519 WRITE( iounit, fmt = 9970 )
520 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
521 WRITE( iounit, fmt = 9950 )1
522 WRITE( iounit, fmt = 6950 )8
523 WRITE( iounit, fmt = 9946 )2
524 WRITE( iounit, fmt = 9944 )3,
'M'
525 WRITE( iounit, fmt = 9943 )4,
'M'
526 WRITE( iounit, fmt = 9942 )5,
'M'
527 WRITE( iounit, fmt = 9941 )6,
'M'
528 WRITE( iounit, fmt = 9960 )7
529 WRITE( iounit, fmt = 6660 )9
530 WRITE( iounit, fmt =
'( '' Messages:'' )' )
532 ELSE IF(
lsamen( 2, p2,
'LQ' ) )
THEN
536 WRITE( iounit, fmt = 9987 )path,
'LQ'
537 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
538 WRITE( iounit, fmt = 9970 )
539 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
540 WRITE( iounit, fmt = 9949 )1
541 WRITE( iounit, fmt = 9945 )2
542 WRITE( iounit, fmt = 9944 )3,
'N'
543 WRITE( iounit, fmt = 9943 )4,
'N'
544 WRITE( iounit, fmt = 9942 )5,
'N'
545 WRITE( iounit, fmt = 9941 )6,
'N'
546 WRITE( iounit, fmt = 9960 )7
547 WRITE( iounit, fmt =
'( '' Messages:'' )' )
549 ELSE IF(
lsamen( 2, p2,
'QL' ) )
THEN
553 WRITE( iounit, fmt = 9987 )path,
'QL'
554 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
555 WRITE( iounit, fmt = 9970 )
556 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
557 WRITE( iounit, fmt = 9948 )1
558 WRITE( iounit, fmt = 9946 )2
559 WRITE( iounit, fmt = 9944 )3,
'M'
560 WRITE( iounit, fmt = 9943 )4,
'M'
561 WRITE( iounit, fmt = 9942 )5,
'M'
562 WRITE( iounit, fmt = 9941 )6,
'M'
563 WRITE( iounit, fmt = 9960 )7
564 WRITE( iounit, fmt =
'( '' Messages:'' )' )
566 ELSE IF(
lsamen( 2, p2,
'RQ' ) )
THEN
570 WRITE( iounit, fmt = 9987 )path,
'RQ'
571 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
572 WRITE( iounit, fmt = 9970 )
573 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
574 WRITE( iounit, fmt = 9947 )1
575 WRITE( iounit, fmt = 9945 )2
576 WRITE( iounit, fmt = 9944 )3,
'N'
577 WRITE( iounit, fmt = 9943 )4,
'N'
578 WRITE( iounit, fmt = 9942 )5,
'N'
579 WRITE( iounit, fmt = 9941 )6,
'N'
580 WRITE( iounit, fmt = 9960 )7
581 WRITE( iounit, fmt =
'( '' Messages:'' )' )
583 ELSE IF(
lsamen( 2, p2,
'QP' ) )
THEN
587 WRITE( iounit, fmt = 9986 )path
588 WRITE( iounit, fmt = 9969 )
589 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
590 WRITE( iounit, fmt = 9940 )1
591 WRITE( iounit, fmt = 9939 )2
592 WRITE( iounit, fmt = 9938 )3
593 WRITE( iounit, fmt =
'( '' Messages:'' )' )
595 ELSE IF(
lsamen( 2, p2,
'TZ' ) )
THEN
599 WRITE( iounit, fmt = 9985 )path
600 WRITE( iounit, fmt = 9968 )
601 WRITE( iounit, fmt = 9929 )c1
602 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
603 WRITE( iounit, fmt = 9940 )1
604 WRITE( iounit, fmt = 9937 )2
605 WRITE( iounit, fmt = 9938 )3
606 WRITE( iounit, fmt =
'( '' Messages:'' )' )
608 ELSE IF(
lsamen( 2, p2,
'LS' ) )
THEN
613 WRITE( iounit, fmt = 9984 )path
614 WRITE( iounit, fmt = 9967 )
615 WRITE( iounit, fmt = 9921 )c1, c1, c1, c1
616 WRITE( iounit, fmt = 9935 )1
617 WRITE( iounit, fmt = 9931 )2
618 WRITE( iounit, fmt = 9933 )3
619 WRITE( iounit, fmt = 9935 )4
620 WRITE( iounit, fmt = 9934 )5
621 WRITE( iounit, fmt = 9932 )6
622 WRITE( iounit, fmt = 9920 )
623 WRITE( iounit, fmt =
'( '' Messages:'' )' )
625 ELSE IF(
lsamen( 2, p2,
'LU' ) )
THEN
629 WRITE( iounit, fmt = 9983 )path
630 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
631 WRITE( iounit, fmt = 9979 )
632 WRITE( iounit, fmt =
'( '' Test ratio:'' )' )
633 WRITE( iounit, fmt = 9962 )1
634 WRITE( iounit, fmt =
'( '' Messages:'' )' )
636 ELSE IF(
lsamen( 2, p2,
'CH' ) )
THEN
640 WRITE( iounit, fmt = 9982 )path
641 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
642 WRITE( iounit, fmt = 9974 )
643 WRITE( iounit, fmt =
'( '' Test ratio:'' )' )
644 WRITE( iounit, fmt = 9954 )1
645 WRITE( iounit, fmt =
'( '' Messages:'' )' )
647 ELSE IF(
lsamen( 2, p2,
'QS' ) )
THEN
651 WRITE( iounit, fmt = 9981 )path
652 WRITE( iounit, fmt =
'( '' Matrix types:'' )' )
653 WRITE( iounit, fmt = 9970 )
654 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
656 ELSE IF(
lsamen( 2, p2,
'QT' ) )
THEN
660 WRITE( iounit, fmt = 8000 ) path
661 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
662 WRITE( iounit, fmt = 8011 ) 1
663 WRITE( iounit, fmt = 8012 ) 2
664 WRITE( iounit, fmt = 8013 ) 3
665 WRITE( iounit, fmt = 8014 ) 4
666 WRITE( iounit, fmt = 8015 ) 5
667 WRITE( iounit, fmt = 8016 ) 6
669 ELSE IF(
lsamen( 2, p2,
'QX' ) )
THEN
673 WRITE( iounit, fmt = 8001 ) path
674 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
675 WRITE( iounit, fmt = 8017 ) 1
676 WRITE( iounit, fmt = 8018 ) 2
677 WRITE( iounit, fmt = 8019 ) 3
678 WRITE( iounit, fmt = 8020 ) 4
679 WRITE( iounit, fmt = 8021 ) 5
680 WRITE( iounit, fmt = 8022 ) 6
682 ELSE IF(
lsamen( 2, p2,
'TQ' ) )
THEN
686 WRITE( iounit, fmt = 8002 ) path
687 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
688 WRITE( iounit, fmt = 8023 ) 1
689 WRITE( iounit, fmt = 8024 ) 2
690 WRITE( iounit, fmt = 8025 ) 3
691 WRITE( iounit, fmt = 8026 ) 4
692 WRITE( iounit, fmt = 8027 ) 5
693 WRITE( iounit, fmt = 8028 ) 6
695 ELSE IF(
lsamen( 2, p2,
'XQ' ) )
THEN
699 WRITE( iounit, fmt = 8003 ) path
700 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
701 WRITE( iounit, fmt = 8029 ) 1
702 WRITE( iounit, fmt = 8030 ) 2
703 WRITE( iounit, fmt = 8031 ) 3
704 WRITE( iounit, fmt = 8032 ) 4
705 WRITE( iounit, fmt = 8033 ) 5
706 WRITE( iounit, fmt = 8034 ) 6
708 ELSE IF(
lsamen( 2, p2,
'TS' ) )
THEN
712 WRITE( iounit, fmt = 8004 ) path
713 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
714 WRITE( iounit, fmt = 8035 ) 1
715 WRITE( iounit, fmt = 8036 ) 2
716 WRITE( iounit, fmt = 8037 ) 3
717 WRITE( iounit, fmt = 8038 ) 4
718 WRITE( iounit, fmt = 8039 ) 5
719 WRITE( iounit, fmt = 8040 ) 6
721 ELSE IF(
lsamen( 2, p2,
'HH' ) )
THEN
725 WRITE( iounit, fmt = 8005 ) path
726 WRITE( iounit, fmt =
'( '' Test ratios:'' )' )
727 WRITE( iounit, fmt = 8050 ) 1
728 WRITE( iounit, fmt = 8051 ) 2
729 WRITE( iounit, fmt = 8052 ) 3
730 WRITE( iounit, fmt = 8053 ) 4
731 WRITE( iounit, fmt = 8054 ) 5
732 WRITE( iounit, fmt = 8055 ) 6
738 WRITE( iounit, fmt = 9980 )path
743 9999
FORMAT( / 1x, a3,
': General dense matrices' )
744 9998
FORMAT( / 1x, a3,
': General band matrices' )
745 9997
FORMAT( / 1x, a3,
': General tridiagonal' )
746 9996
FORMAT( / 1x, a3,
': ', a9,
' positive definite matrices' )
747 9995
FORMAT( / 1x, a3,
': ', a9,
' positive definite packed matrices'
749 9994
FORMAT( / 1x, a3,
': ', a9,
' positive definite band matrices' )
750 9993
FORMAT( / 1x, a3,
': ', a9,
' positive definite tridiagonal' )
751 9992
FORMAT( / 1x, a3,
': ', a9,
' indefinite matrices',
752 $
', partial (Bunch-Kaufman) pivoting' )
753 9991
FORMAT( / 1x, a3,
': ', a9,
' indefinite packed matrices',
754 $
', partial (Bunch-Kaufman) pivoting' )
755 9892
FORMAT( / 1x, a3,
': ', a9,
' indefinite matrices',
756 $
', "rook" (bounded Bunch-Kaufman) pivoting' )
757 9891
FORMAT( / 1x, a3,
': ', a9,
' indefinite packed matrices',
758 $
', "rook" (bounded Bunch-Kaufman) pivoting' )
759 9990
FORMAT( / 1x, a3,
': Triangular matrices' )
760 9989
FORMAT( / 1x, a3,
': Triangular packed matrices' )
761 9988
FORMAT( / 1x, a3,
': Triangular band matrices' )
762 9987
FORMAT( / 1x, a3,
': ', a2,
' factorization of general matrices'
764 9986
FORMAT( / 1x, a3,
': QR factorization with column pivoting' )
765 9985
FORMAT( / 1x, a3,
': RQ factorization of trapezoidal matrix' )
766 9984
FORMAT( / 1x, a3,
': Least squares driver routines' )
767 9983
FORMAT( / 1x, a3,
': LU factorization variants' )
768 9982
FORMAT( / 1x, a3,
': Cholesky factorization variants' )
769 9981
FORMAT( / 1x, a3,
': QR factorization variants' )
770 9980
FORMAT( / 1x, a3,
': No header available' )
771 8000
FORMAT( / 1x, a3,
': QRT factorization for general matrices' )
772 8001
FORMAT( / 1x, a3,
': QRT factorization for ',
773 $
'triangular-pentagonal matrices' )
774 8002
FORMAT( / 1x, a3,
': LQT factorization for general matrices' )
775 8003
FORMAT( / 1x, a3,
': LQT factorization for ',
776 $
'triangular-pentagonal matrices' )
777 8004
FORMAT( / 1x, a3,
': TS factorization for ',
778 $
'tall-skinny or short-wide matrices' )
779 8005
FORMAT( / 1x, a3,
': Householder recostruction from TSQR',
780 $
' factorization output ', /,
' for tall-skinny matrices.' )
784 9979
FORMAT( 4x,
'1. Diagonal', 24x,
'7. Last n/2 columns zero', / 4x,
785 $
'2. Upper triangular', 16x,
786 $
'8. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
787 $
'3. Lower triangular', 16x,
'9. Random, CNDNUM = 0.1/EPS',
788 $ / 4x,
'4. Random, CNDNUM = 2', 13x,
789 $
'10. Scaled near underflow', / 4x,
'5. First column zero',
790 $ 14x,
'11. Scaled near overflow', / 4x,
791 $
'6. Last column zero' )
795 9978
FORMAT( 4x,
'1. Random, CNDNUM = 2', 14x,
796 $
'5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
797 $
'2. First column zero', 15x,
'6. Random, CNDNUM = .01/EPS',
798 $ / 4x,
'3. Last column zero', 16x,
799 $
'7. Scaled near underflow', / 4x,
800 $
'4. Last n/2 columns zero', 11x,
'8. Scaled near overflow' )
804 9977
FORMAT(
' Matrix types (1-6 have specified condition numbers):',
805 $ / 4x,
'1. Diagonal', 24x,
'7. Random, unspecified CNDNUM',
806 $ / 4x,
'2. Random, CNDNUM = 2', 14x,
'8. First column zero',
807 $ / 4x,
'3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
808 $
'9. Last column zero', / 4x,
'4. Random, CNDNUM = 0.1/EPS',
809 $ 7x,
'10. Last n/2 columns zero', / 4x,
810 $
'5. Scaled near underflow', 10x,
811 $
'11. Scaled near underflow', / 4x,
812 $
'6. Scaled near overflow', 11x,
'12. Scaled near overflow' )
816 9976
FORMAT(
' Matrix types (1-6 have specified condition numbers):',
817 $ / 4x,
'1. Diagonal', 24x,
'7. Random, unspecified CNDNUM',
818 $ / 4x,
'2. Random, CNDNUM = 2', 14x,
819 $
'8. First row and column zero', / 4x,
820 $
'3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
821 $
'9. Last row and column zero', / 4x,
822 $
'4. Random, CNDNUM = 0.1/EPS', 7x,
823 $
'10. Middle row and column zero', / 4x,
824 $
'5. Scaled near underflow', 10x,
825 $
'11. Scaled near underflow', / 4x,
826 $
'6. Scaled near overflow', 11x,
'12. Scaled near overflow' )
830 9975
FORMAT( 4x,
'1. Diagonal', 24x,
831 $
'6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
832 $
'2. Random, CNDNUM = 2', 14x,
'7. Random, CNDNUM = 0.1/EPS',
833 $ / 3x,
'*3. First row and column zero', 7x,
834 $
'8. Scaled near underflow', / 3x,
835 $
'*4. Last row and column zero', 8x,
836 $
'9. Scaled near overflow', / 3x,
837 $
'*5. Middle row and column zero', / 3x,
838 $
'(* - tests error exits from ', a3,
839 $
'TRF, no test ratios are computed)' )
843 9974
FORMAT( 4x,
'1. Diagonal', 24x,
844 $
'6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
845 $
'2. Random, CNDNUM = 2', 14x,
'7. Random, CNDNUM = 0.1/EPS',
846 $ / 3x,
'*3. First row and column zero', 7x,
847 $
'8. Scaled near underflow', / 3x,
848 $
'*4. Last row and column zero', 8x,
849 $
'9. Scaled near overflow', / 3x,
850 $
'*5. Middle row and column zero', / 3x,
851 $
'(* - tests error exits, no test ratios are computed)' )
855 8973
FORMAT( 4x,
'1. Diagonal', / 4x,
'2. Random, CNDNUM = 2', 14x,
856 $ / 3x,
'*3. Nonzero eigenvalues of: D(1:RANK-1)=1 and ',
857 $
'D(RANK) = 1.0/', a4, / 3x,
858 $
'*4. Nonzero eigenvalues of: D(1)=1 and ',
859 $
' D(2:RANK) = 1.0/', a4, / 3x,
860 $
'*5. Nonzero eigenvalues of: D(I) = ', a4,
861 $
'**(-(I-1)/(RANK-1)) ',
' I=1:RANK', / 4x,
862 $
'6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
863 $
'7. Random, CNDNUM = 0.1/EPS', / 4x,
864 $
'8. Scaled near underflow', / 4x,
'9. Scaled near overflow',
865 $ / 3x,
'(* - Semi-definite tests )' )
866 8972
FORMAT( 3x,
'RANK minus computed rank, returned by ', a,
'PSTRF' )
870 9973
FORMAT( 4x,
'1. Random, CNDNUM = 2', 14x,
871 $
'5. Random, CNDNUM = sqrt(0.1/EPS)', / 3x,
872 $
'*2. First row and column zero', 7x,
873 $
'6. Random, CNDNUM = 0.1/EPS', / 3x,
874 $
'*3. Last row and column zero', 8x,
875 $
'7. Scaled near underflow', / 3x,
876 $
'*4. Middle row and column zero', 6x,
877 $
'8. Scaled near overflow', / 3x,
878 $
'(* - tests error exits from ', a3,
879 $
'TRF, no test ratios are computed)' )
883 9972
FORMAT( 4x,
'1. Diagonal', 24x,
884 $
'6. Last n/2 rows and columns zero', / 4x,
885 $
'2. Random, CNDNUM = 2', 14x,
886 $
'7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
887 $
'3. First row and column zero', 7x,
888 $
'8. Random, CNDNUM = 0.1/EPS', / 4x,
889 $
'4. Last row and column zero', 8x,
890 $
'9. Scaled near underflow', / 4x,
891 $
'5. Middle row and column zero', 5x,
892 $
'10. Scaled near overflow' )
896 9971
FORMAT( 4x,
'1. Diagonal', 24x,
897 $
'7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
898 $
'2. Random, CNDNUM = 2', 14x,
'8. Random, CNDNUM = 0.1/EPS',
899 $ / 4x,
'3. First row and column zero', 7x,
900 $
'9. Scaled near underflow', / 4x,
901 $
'4. Last row and column zero', 7x,
902 $
'10. Scaled near overflow', / 4x,
903 $
'5. Middle row and column zero', 5x,
904 $
'11. Block diagonal matrix', / 4x,
905 $
'6. Last n/2 rows and columns zero' )
909 9970
FORMAT( 4x,
'1. Diagonal', 24x,
910 $
'5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
911 $
'2. Upper triangular', 16x,
'6. Random, CNDNUM = 0.1/EPS',
912 $ / 4x,
'3. Lower triangular', 16x,
913 $
'7. Scaled near underflow', / 4x,
'4. Random, CNDNUM = 2',
914 $ 14x,
'8. Scaled near overflow' )
918 9969
FORMAT(
' Matrix types (2-6 have condition 1/EPS):', / 4x,
919 $
'1. Zero matrix', 21x,
'4. First n/2 columns fixed', / 4x,
920 $
'2. One small eigenvalue', 12x,
'5. Last n/2 columns fixed',
921 $ / 4x,
'3. Geometric distribution', 10x,
922 $
'6. Every second column fixed' )
926 9968
FORMAT(
' Matrix types (2-3 have condition 1/EPS):', / 4x,
927 $
'1. Zero matrix', / 4x,
'2. One small eigenvalue', / 4x,
928 $
'3. Geometric distribution' )
932 9967
FORMAT(
' Matrix types (1-3: full rank, 4-6: rank deficient):',
933 $ / 4x,
'1 and 4. Normal scaling', / 4x,
934 $
'2 and 5. Scaled near overflow', / 4x,
935 $
'3 and 6. Scaled near underflow' )
939 9966
FORMAT(
' Matrix types for ', a3,
' routines:', / 4x,
940 $
'1. Diagonal', 24x,
'6. Scaled near overflow', / 4x,
941 $
'2. Random, CNDNUM = 2', 14x,
'7. Identity', / 4x,
942 $
'3. Random, CNDNUM = sqrt(0.1/EPS) ',
943 $
'8. Unit triangular, CNDNUM = 2', / 4x,
944 $
'4. Random, CNDNUM = 0.1/EPS', 8x,
945 $
'9. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
946 $
'5. Scaled near underflow', 10x,
947 $
'10. Unit, CNDNUM = 0.1/EPS' )
948 9965
FORMAT(
' Special types for testing ', a,
':', / 3x,
949 $
'11. Matrix elements are O(1), large right hand side', / 3x,
950 $
'12. First diagonal causes overflow,',
951 $
' offdiagonal column norms < 1', / 3x,
952 $
'13. First diagonal causes overflow,',
953 $
' offdiagonal column norms > 1', / 3x,
954 $
'14. Growth factor underflows, solution does not overflow',
955 $ / 3x,
'15. Small diagonal causes gradual overflow', / 3x,
956 $
'16. One zero diagonal element', / 3x,
957 $
'17. Large offdiagonals cause overflow when adding a column'
958 $ , / 3x,
'18. Unit triangular with large right hand side' )
962 9964
FORMAT(
' Matrix types for ', a3,
' routines:', / 4x,
963 $
'1. Random, CNDNUM = 2', 14x,
'6. Identity', / 4x,
964 $
'2. Random, CNDNUM = sqrt(0.1/EPS) ',
965 $
'7. Unit triangular, CNDNUM = 2', / 4x,
966 $
'3. Random, CNDNUM = 0.1/EPS', 8x,
967 $
'8. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
968 $
'4. Scaled near underflow', 11x,
969 $
'9. Unit, CNDNUM = 0.1/EPS', / 4x,
970 $
'5. Scaled near overflow' )
971 9963
FORMAT(
' Special types for testing ', a,
':', / 3x,
972 $
'10. Matrix elements are O(1), large right hand side', / 3x,
973 $
'11. First diagonal causes overflow,',
974 $
' offdiagonal column norms < 1', / 3x,
975 $
'12. First diagonal causes overflow,',
976 $
' offdiagonal column norms > 1', / 3x,
977 $
'13. Growth factor underflows, solution does not overflow',
978 $ / 3x,
'14. Small diagonal causes gradual overflow', / 3x,
979 $
'15. One zero diagonal element', / 3x,
980 $
'16. Large offdiagonals cause overflow when adding a column'
981 $ , / 3x,
'17. Unit triangular with large right hand side' )
985 9962
FORMAT( 3x, i2,
': norm( L * U - A ) / ( N * norm(A) * EPS )' )
986 9961
FORMAT( 3x, i2,
': norm( I - A*AINV ) / ',
987 $
'( N * norm(A) * norm(AINV) * EPS )' )
988 9960
FORMAT( 3x, i2,
': norm( B - A * X ) / ',
989 $
'( norm(A) * norm(X) * EPS )' )
990 6660
FORMAT( 3x, i2,
': diagonal is not non-negative')
991 9959
FORMAT( 3x, i2,
': norm( X - XACT ) / ',
992 $
'( norm(XACT) * CNDNUM * EPS )' )
993 9958
FORMAT( 3x, i2,
': norm( X - XACT ) / ',
994 $
'( norm(XACT) * CNDNUM * EPS ), refined' )
995 9957
FORMAT( 3x, i2,
': norm( X - XACT ) / ',
996 $
'( norm(XACT) * (error bound) )' )
997 9956
FORMAT( 3x, i2,
': (backward error) / EPS' )
998 9955
FORMAT( 3x, i2,
': RCOND * CNDNUM - 1.0' )
999 9954
FORMAT( 3x, i2,
': norm( U'' * U - A ) / ( N * norm(A) * EPS )',
1000 $
', or', / 7x,
'norm( L * L'' - A ) / ( N * norm(A) * EPS )'
1003 $
'norm( P * U'' * U * P'' - A ) / ( N * norm(A) * EPS )',
1005 $
'norm( P * L * L'' * P'' - A ) / ( N * norm(A) * EPS )' )
1006 9953
FORMAT( 3x, i2,
': norm( U*D*U'' - A ) / ( N * norm(A) * EPS )',
1007 $
', or', / 7x,
'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
1009 9952
FORMAT( 3x, i2,
': norm( U''*D*U - A ) / ( N * norm(A) * EPS )',
1010 $
', or', / 7x,
'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
1012 9951
FORMAT(
' Test ratio for ', a,
':', / 3x, i2,
1013 $
': norm( s*b - A*x ) / ( norm(A) * norm(x) * EPS )' )
1014 9950
FORMAT( 3x, i2,
': norm( R - Q'' * A ) / ( M * norm(A) * EPS )' )
1015 6950
FORMAT( 3x, i2,
': norm( R - Q'' * A ) / ( M * norm(A) * EPS )
1017 9949
FORMAT( 3x, i2,
': norm( L - A * Q'' ) / ( N * norm(A) * EPS )' )
1018 9948
FORMAT( 3x, i2,
': norm( L - Q'' * A ) / ( M * norm(A) * EPS )' )
1019 9947
FORMAT( 3x, i2,
': norm( R - A * Q'' ) / ( N * norm(A) * EPS )' )
1020 9946
FORMAT( 3x, i2,
': norm( I - Q''*Q ) / ( M * EPS )' )
1021 9945
FORMAT( 3x, i2,
': norm( I - Q*Q'' ) / ( N * EPS )' )
1022 9944
FORMAT( 3x, i2,
': norm( Q*C - Q*C ) / ',
'( ', a1,
1023 $
' * norm(C) * EPS )' )
1024 9943
FORMAT( 3x, i2,
': norm( C*Q - C*Q ) / ',
'( ', a1,
1025 $
' * norm(C) * EPS )' )
1026 9942
FORMAT( 3x, i2,
': norm( Q''*C - Q''*C )/ ',
'( ', a1,
1027 $
' * norm(C) * EPS )' )
1028 9941
FORMAT( 3x, i2,
': norm( C*Q'' - C*Q'' )/ ',
'( ', a1,
1029 $
' * norm(C) * EPS )' )
1030 9940
FORMAT( 3x, i2,
': norm(svd(A) - svd(R)) / ',
1031 $
'( M * norm(svd(R)) * EPS )' )
1032 9939
FORMAT( 3x, i2,
': norm( A*P - Q*R ) / ( M * norm(A) * EPS )'
1034 9938
FORMAT( 3x, i2,
': norm( I - Q''*Q ) / ( M * EPS )' )
1035 9937
FORMAT( 3x, i2,
': norm( A - R*Q ) / ( M * norm(A) * EPS )'
1037 9935
FORMAT( 3x, i2,
': norm( B - A * X ) / ',
1038 $
'( max(M,N) * norm(A) * norm(X) * EPS )' )
1039 9934
FORMAT( 3x, i2,
': norm( (A*X-B)'' *A ) / ',
1040 $
'( max(M,N,NRHS) * norm(A) * norm(B) * EPS )' )
1041 9933
FORMAT( 3x, i2,
': norm(svd(A)-svd(R)) / ',
1042 $
'( min(M,N) * norm(svd(R)) * EPS )' )
1043 9932
FORMAT( 3x, i2,
': Check if X is in the row space of A or A''' )
1044 9931
FORMAT( 3x, i2,
': norm( (A*X-B)'' *A ) / ',
1045 $
'( max(M,N,NRHS) * norm(A) * norm(B) * EPS )', / 7x,
1046 $
'if TRANS=''N'.GE.
' and MN or TRANS=''T'.LT.
' and MN, ',
1047 $
'otherwise', / 7x,
1048 $
'check if X is in the row space of A or A'' ',
1049 $
'(overdetermined case)' )
1050 9929
FORMAT(
' Test ratios (1-3: ', a1,
'TZRZF):' )
1051 9920
FORMAT( 3x,
' 7-10: same as 3-6', 3x,
' 11-14: same as 3-6' )
1052 9921
FORMAT(
' Test ratios:', /
' (1-2: ', a1,
'GELS, 3-6: ', a1,
1053 $
'GELSY, 7-10: ', a1,
'GELSS, 11-14: ', a1,
'GELSD, 15-16: ',
1055 9928
FORMAT( 7x,
'where ALPHA = ( 1 + SQRT( 17 ) ) / 8' )
1056 9927
FORMAT( 3x, i2,
': ABS( Largest element in L )', / 12x,
1057 $
' - ( 1 / ( 1 - ALPHA ) ) + THRESH' )
1058 9926
FORMAT( 3x, i2,
': Largest 2-Norm of 2-by-2 pivots', / 12x,
1059 $
' - ( ( 1 + ALPHA ) / ( 1 - ALPHA ) ) + THRESH' )
1060 8011
FORMAT(3x,i2,
': norm( R - Q''*A ) / ( M * norm(A) * EPS )' )
1061 8012
FORMAT(3x,i2,
': norm( I - Q''*Q ) / ( M * EPS )' )
1062 8013
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( M * norm(C) * EPS )' )
1063 8014
FORMAT(3x,i2,
': norm( Q''*C - Q''*C ) / ( M * norm(C) * EPS )')
1064 8015
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( M * norm(C) * EPS )' )
1065 8016
FORMAT(3x,i2,
': norm( C*Q'' - C*Q'' ) / ( M * norm(C) * EPS )')
1066 8017
FORMAT(3x,i2,
': norm( R - Q''*A ) / ( (M+N) * norm(A) * EPS )' )
1067 8018
FORMAT(3x,i2,
': norm( I - Q''*Q ) / ( (M+N) * EPS )' )
1068 8019
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1070 $
': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1071 8021
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1073 $
': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1074 8023
FORMAT(3x,i2,
': norm( L - A*Q'' ) / ( (M+N) * norm(A) * EPS )' )
1075 8024
FORMAT(3x,i2,
': norm( I - Q*Q'' ) / ( (M+N) * EPS )' )
1076 8025
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1078 $
': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1079 8027
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1081 $
': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1082 8029
FORMAT(3x,i2,
': norm( L - A*Q'' ) / ( (M+N) * norm(A) * EPS )' )
1083 8030
FORMAT(3x,i2,
': norm( I - Q*Q'' ) / ( (M+N) * EPS )' )
1084 8031
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1086 $
': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1087 8033
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1089 $
': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1090 8035
FORMAT(3x,i2,
': norm( R - Q''*A ) / ( (M+N) * norm(A) * EPS )' )
1091 8036
FORMAT(3x,i2,
': norm( I - Q''*Q ) / ( (M+N) * EPS )' )
1092 8037
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1094 $
': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1095 8039
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1097 $
': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1099 8050
FORMAT(3x,i2,
': norm( R - Q''*A ) / ( M * norm(A) * EPS )' )
1100 8051
FORMAT(3x,i2,
': norm( I - Q''*Q ) / ( M * EPS )' )
1101 8052
FORMAT(3x,i2,
': norm( Q*C - Q*C ) / ( M * norm(C) * EPS )' )
1102 8053
FORMAT(3x,i2,
': norm( Q''*C - Q''*C ) / ( M * norm(C) * EPS )')
1103 8054
FORMAT(3x,i2,
': norm( C*Q - C*Q ) / ( M * norm(C) * EPS )' )
1104 8055
FORMAT(3x,i2,
': norm( C*Q'' - C*Q'' ) / ( M * norm(C) * EPS )')
logical function lsamen(N, CA, CB)
LSAMEN
logical function lsame(CA, CB)
LSAME