LAPACK  3.6.1 LAPACK: Linear Algebra PACKage
 subroutine check2 ( double precision SFAC )

Definition at line 337 of file dblat1.f.

337 * .. Parameters ..
338  INTEGER nout
339  parameter (nout=6)
340 * .. Scalar Arguments ..
341  DOUBLE PRECISION sfac
342 * .. Scalars in Common ..
343  INTEGER icase, incx, incy, n
344  LOGICAL pass
345 * .. Local Scalars ..
346  DOUBLE PRECISION sa
347  INTEGER i, j, ki, kn, kni, kpar, ksize, lenx, leny,
348  \$ mx, my
349 * .. Local Arrays ..
350  DOUBLE PRECISION dt10x(7,4,4), dt10y(7,4,4), dt7(4,4),
351  \$ dt8(7,4,4), dx1(7),
352  \$ dy1(7), ssize1(4), ssize2(14,2), ssize(7),
353  \$ stx(7), sty(7), sx(7), sy(7),
354  \$ dpar(5,4), dt19x(7,4,16),dt19xa(7,4,4),
355  \$ dt19xb(7,4,4), dt19xc(7,4,4),dt19xd(7,4,4),
356  \$ dt19y(7,4,16), dt19ya(7,4,4),dt19yb(7,4,4),
357  \$ dt19yc(7,4,4), dt19yd(7,4,4), dtemp(5)
358  INTEGER incxs(4), incys(4), lens(4,2), ns(4)
359 * .. External Functions ..
360  DOUBLE PRECISION ddot, dsdot
361  EXTERNAL ddot, dsdot
362 * .. External Subroutines ..
363  EXTERNAL daxpy, dcopy, drotm, dswap, stest, stest1
364 * .. Intrinsic Functions ..
365  INTRINSIC abs, min
366 * .. Common blocks ..
367  COMMON /combla/icase, n, incx, incy, pass
368 * .. Data statements ..
369  equivalence(dt19x(1,1,1),dt19xa(1,1,1)),(dt19x(1,1,5),
370  a dt19xb(1,1,1)),(dt19x(1,1,9),dt19xc(1,1,1)),
371  b (dt19x(1,1,13),dt19xd(1,1,1))
372  equivalence(dt19y(1,1,1),dt19ya(1,1,1)),(dt19y(1,1,5),
373  a dt19yb(1,1,1)),(dt19y(1,1,9),dt19yc(1,1,1)),
374  b (dt19y(1,1,13),dt19yd(1,1,1))
375
376  DATA sa/0.3d0/
377  DATA incxs/1, 2, -2, -1/
378  DATA incys/1, -2, 1, -2/
379  DATA lens/1, 1, 2, 4, 1, 1, 3, 7/
380  DATA ns/0, 1, 2, 4/
381  DATA dx1/0.6d0, 0.1d0, -0.5d0, 0.8d0, 0.9d0, -0.3d0,
382  + -0.4d0/
383  DATA dy1/0.5d0, -0.9d0, 0.3d0, 0.7d0, -0.6d0, 0.2d0,
384  + 0.8d0/
385  DATA dt7/0.0d0, 0.30d0, 0.21d0, 0.62d0, 0.0d0,
386  + 0.30d0, -0.07d0, 0.85d0, 0.0d0, 0.30d0, -0.79d0,
387  + -0.74d0, 0.0d0, 0.30d0, 0.33d0, 1.27d0/
388  DATA dt8/0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
389  + 0.0d0, 0.68d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
390  + 0.0d0, 0.0d0, 0.68d0, -0.87d0, 0.0d0, 0.0d0,
391  + 0.0d0, 0.0d0, 0.0d0, 0.68d0, -0.87d0, 0.15d0,
392  + 0.94d0, 0.0d0, 0.0d0, 0.0d0, 0.5d0, 0.0d0,
393  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.68d0,
394  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
395  + 0.35d0, -0.9d0, 0.48d0, 0.0d0, 0.0d0, 0.0d0,
396  + 0.0d0, 0.38d0, -0.9d0, 0.57d0, 0.7d0, -0.75d0,
397  + 0.2d0, 0.98d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0,
398  + 0.0d0, 0.0d0, 0.0d0, 0.68d0, 0.0d0, 0.0d0,
399  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.35d0, -0.72d0,
400  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.38d0,
401  + -0.63d0, 0.15d0, 0.88d0, 0.0d0, 0.0d0, 0.0d0,
402  + 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
403  + 0.68d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
404  + 0.0d0, 0.68d0, -0.9d0, 0.33d0, 0.0d0, 0.0d0,
405  + 0.0d0, 0.0d0, 0.68d0, -0.9d0, 0.33d0, 0.7d0,
406  + -0.75d0, 0.2d0, 1.04d0/
407  DATA dt10x/0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
408  + 0.0d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
409  + 0.0d0, 0.5d0, -0.9d0, 0.0d0, 0.0d0, 0.0d0,
410  + 0.0d0, 0.0d0, 0.5d0, -0.9d0, 0.3d0, 0.7d0,
411  + 0.0d0, 0.0d0, 0.0d0, 0.6d0, 0.0d0, 0.0d0, 0.0d0,
412  + 0.0d0, 0.0d0, 0.0d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0,
413  + 0.0d0, 0.0d0, 0.0d0, 0.3d0, 0.1d0, 0.5d0, 0.0d0,
414  + 0.0d0, 0.0d0, 0.0d0, 0.8d0, 0.1d0, -0.6d0,
415  + 0.8d0, 0.3d0, -0.3d0, 0.5d0, 0.6d0, 0.0d0,
416  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.5d0, 0.0d0,
417  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, -0.9d0,
418  + 0.1d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.7d0,
419  + 0.1d0, 0.3d0, 0.8d0, -0.9d0, -0.3d0, 0.5d0,
420  + 0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
421  + 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
422  + 0.5d0, 0.3d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
423  + 0.5d0, 0.3d0, -0.6d0, 0.8d0, 0.0d0, 0.0d0,
424  + 0.0d0/
425  DATA dt10y/0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
426  + 0.0d0, 0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
427  + 0.0d0, 0.6d0, 0.1d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
428  + 0.0d0, 0.6d0, 0.1d0, -0.5d0, 0.8d0, 0.0d0,
429  + 0.0d0, 0.0d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
430  + 0.0d0, 0.0d0, 0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
431  + 0.0d0, 0.0d0, -0.5d0, -0.9d0, 0.6d0, 0.0d0,
432  + 0.0d0, 0.0d0, 0.0d0, -0.4d0, -0.9d0, 0.9d0,
433  + 0.7d0, -0.5d0, 0.2d0, 0.6d0, 0.5d0, 0.0d0,
434  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.6d0, 0.0d0,
435  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, -0.5d0,
436  + 0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
437  + -0.4d0, 0.9d0, -0.5d0, 0.6d0, 0.0d0, 0.0d0,
438  + 0.0d0, 0.5d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
439  + 0.0d0, 0.6d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
440  + 0.0d0, 0.6d0, -0.9d0, 0.1d0, 0.0d0, 0.0d0,
441  + 0.0d0, 0.0d0, 0.6d0, -0.9d0, 0.1d0, 0.7d0,
442  + -0.5d0, 0.2d0, 0.8d0/
443  DATA ssize1/0.0d0, 0.3d0, 1.6d0, 3.2d0/
444  DATA ssize2/0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
445  + 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0, 0.0d0,
446  + 0.0d0, 1.17d0, 1.17d0, 1.17d0, 1.17d0, 1.17d0,
447  + 1.17d0, 1.17d0, 1.17d0, 1.17d0, 1.17d0, 1.17d0,
448  + 1.17d0, 1.17d0, 1.17d0/
449 *
450 * FOR DROTM
451 *
452  DATA dpar/-2.d0, 0.d0,0.d0,0.d0,0.d0,
453  a -1.d0, 2.d0, -3.d0, -4.d0, 5.d0,
454  b 0.d0, 0.d0, 2.d0, -3.d0, 0.d0,
455  c 1.d0, 5.d0, 2.d0, 0.d0, -4.d0/
456 * TRUE X RESULTS F0R ROTATIONS DROTM
457  DATA dt19xa/.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
458  a .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
459  b .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
460  c .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
461  d .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
462  e -.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
463  f -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
464  g 3.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
465  h .6d0, .1d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
466  i -.8d0, 3.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
467  j -.9d0, 2.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
468  k 3.5d0, -.4d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
469  l .6d0, .1d0, -.5d0, .8d0, 0.d0,0.d0,0.d0,
470  m -.8d0, 3.8d0, -2.2d0, -1.2d0, 0.d0,0.d0,0.d0,
471  n -.9d0, 2.8d0, -1.4d0, -1.3d0, 0.d0,0.d0,0.d0,
472  o 3.5d0, -.4d0, -2.2d0, 4.7d0, 0.d0,0.d0,0.d0/
473 *
474  DATA dt19xb/.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
475  a .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
476  b .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
477  c .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
478  d .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
479  e -.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
480  f -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
481  g 3.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
482  h .6d0, .1d0, -.5d0, 0.d0,0.d0,0.d0,0.d0,
483  i 0.d0, .1d0, -3.0d0, 0.d0,0.d0,0.d0,0.d0,
484  j -.3d0, .1d0, -2.0d0, 0.d0,0.d0,0.d0,0.d0,
485  k 3.3d0, .1d0, -2.0d0, 0.d0,0.d0,0.d0,0.d0,
486  l .6d0, .1d0, -.5d0, .8d0, .9d0, -.3d0, -.4d0,
487  m -2.0d0, .1d0, 1.4d0, .8d0, .6d0, -.3d0, -2.8d0,
488  n -1.8d0, .1d0, 1.3d0, .8d0, 0.d0, -.3d0, -1.9d0,
489  o 3.8d0, .1d0, -3.1d0, .8d0, 4.8d0, -.3d0, -1.5d0 /
490 *
491  DATA dt19xc/.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
492  a .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
493  b .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
494  c .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
495  d .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
496  e -.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
497  f -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
498  g 3.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
499  h .6d0, .1d0, -.5d0, 0.d0,0.d0,0.d0,0.d0,
500  i 4.8d0, .1d0, -3.0d0, 0.d0,0.d0,0.d0,0.d0,
501  j 3.3d0, .1d0, -2.0d0, 0.d0,0.d0,0.d0,0.d0,
502  k 2.1d0, .1d0, -2.0d0, 0.d0,0.d0,0.d0,0.d0,
503  l .6d0, .1d0, -.5d0, .8d0, .9d0, -.3d0, -.4d0,
504  m -1.6d0, .1d0, -2.2d0, .8d0, 5.4d0, -.3d0, -2.8d0,
505  n -1.5d0, .1d0, -1.4d0, .8d0, 3.6d0, -.3d0, -1.9d0,
506  o 3.7d0, .1d0, -2.2d0, .8d0, 3.6d0, -.3d0, -1.5d0 /
507 *
508  DATA dt19xd/.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
509  a .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
510  b .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
511  c .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
512  d .6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
513  e -.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
514  f -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
515  g 3.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
516  h .6d0, .1d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
517  i -.8d0, -1.0d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
518  j -.9d0, -.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
519  k 3.5d0, .8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
520  l .6d0, .1d0, -.5d0, .8d0, 0.d0,0.d0,0.d0,
521  m -.8d0, -1.0d0, 1.4d0, -1.6d0, 0.d0,0.d0,0.d0,
522  n -.9d0, -.8d0, 1.3d0, -1.6d0, 0.d0,0.d0,0.d0,
523  o 3.5d0, .8d0, -3.1d0, 4.8d0, 0.d0,0.d0,0.d0/
524 * TRUE Y RESULTS FOR ROTATIONS DROTM
525  DATA dt19ya/.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
526  a .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
527  b .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
528  c .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
529  d .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
530  e .7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
531  f 1.7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
532  g -2.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
533  h .5d0, -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
534  i .7d0, -4.8d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
535  j 1.7d0, -.7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
536  k -2.6d0, 3.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
537  l .5d0, -.9d0, .3d0, .7d0, 0.d0,0.d0,0.d0,
538  m .7d0, -4.8d0, 3.0d0, 1.1d0, 0.d0,0.d0,0.d0,
539  n 1.7d0, -.7d0, -.7d0, 2.3d0, 0.d0,0.d0,0.d0,
540  o -2.6d0, 3.5d0, -.7d0, -3.6d0, 0.d0,0.d0,0.d0/
541 *
542  DATA dt19yb/.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
543  a .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
544  b .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
545  c .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
546  d .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
547  e .7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
548  f 1.7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
549  g -2.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
550  h .5d0, -.9d0, .3d0, 0.d0,0.d0,0.d0,0.d0,
551  i 4.0d0, -.9d0, -.3d0, 0.d0,0.d0,0.d0,0.d0,
552  j -.5d0, -.9d0, 1.5d0, 0.d0,0.d0,0.d0,0.d0,
553  k -1.5d0, -.9d0, -1.8d0, 0.d0,0.d0,0.d0,0.d0,
554  l .5d0, -.9d0, .3d0, .7d0, -.6d0, .2d0, .8d0,
555  m 3.7d0, -.9d0, -1.2d0, .7d0, -1.5d0, .2d0, 2.2d0,
556  n -.3d0, -.9d0, 2.1d0, .7d0, -1.6d0, .2d0, 2.0d0,
557  o -1.6d0, -.9d0, -2.1d0, .7d0, 2.9d0, .2d0, -3.8d0 /
558 *
559  DATA dt19yc/.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
560  a .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
561  b .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
562  c .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
563  d .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
564  e .7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
565  f 1.7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
566  g -2.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
567  h .5d0, -.9d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
568  i 4.0d0, -6.3d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
569  j -.5d0, .3d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
570  k -1.5d0, 3.0d0, 0.d0,0.d0,0.d0,0.d0,0.d0,
571  l .5d0, -.9d0, .3d0, .7d0, 0.d0,0.d0,0.d0,
572  m 3.7d0, -7.2d0, 3.0d0, 1.7d0, 0.d0,0.d0,0.d0,
573  n -.3d0, .9d0, -.7d0, 1.9d0, 0.d0,0.d0,0.d0,
574  o -1.6d0, 2.7d0, -.7d0, -3.4d0, 0.d0,0.d0,0.d0/
575 *
576  DATA dt19yd/.5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
577  a .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
578  b .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
579  c .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
580  d .5d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
581  e .7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
582  f 1.7d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
583  g -2.6d0, 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0,
584  h .5d0, -.9d0, .3d0, 0.d0,0.d0,0.d0,0.d0,
585  i .7d0, -.9d0, 1.2d0, 0.d0,0.d0,0.d0,0.d0,
586  j 1.7d0, -.9d0, .5d0, 0.d0,0.d0,0.d0,0.d0,
587  k -2.6d0, -.9d0, -1.3d0, 0.d0,0.d0,0.d0,0.d0,
588  l .5d0, -.9d0, .3d0, .7d0, -.6d0, .2d0, .8d0,
589  m .7d0, -.9d0, 1.2d0, .7d0, -1.5d0, .2d0, 1.6d0,
590  n 1.7d0, -.9d0, .5d0, .7d0, -1.6d0, .2d0, 2.4d0,
591  o -2.6d0, -.9d0, -1.3d0, .7d0, 2.9d0, .2d0, -4.0d0 /
592 *
593 * .. Executable Statements ..
594 *
595  DO 120 ki = 1, 4
596  incx = incxs(ki)
597  incy = incys(ki)
598  mx = abs(incx)
599  my = abs(incy)
600 *
601  DO 100 kn = 1, 4
602  n = ns(kn)
603  ksize = min(2,kn)
604  lenx = lens(kn,mx)
605  leny = lens(kn,my)
606 * .. Initialize all argument arrays ..
607  DO 20 i = 1, 7
608  sx(i) = dx1(i)
609  sy(i) = dy1(i)
610  20 CONTINUE
611 *
612  IF (icase.EQ.1) THEN
613 * .. DDOT ..
614  CALL stest1(ddot(n,sx,incx,sy,incy),dt7(kn,ki),ssize1(kn)
615  + ,sfac)
616  ELSE IF (icase.EQ.2) THEN
617 * .. DAXPY ..
618  CALL daxpy(n,sa,sx,incx,sy,incy)
619  DO 40 j = 1, leny
620  sty(j) = dt8(j,kn,ki)
621  40 CONTINUE
622  CALL stest(leny,sy,sty,ssize2(1,ksize),sfac)
623  ELSE IF (icase.EQ.5) THEN
624 * .. DCOPY ..
625  DO 60 i = 1, 7
626  sty(i) = dt10y(i,kn,ki)
627  60 CONTINUE
628  CALL dcopy(n,sx,incx,sy,incy)
629  CALL stest(leny,sy,sty,ssize2(1,1),1.0d0)
630  ELSE IF (icase.EQ.6) THEN
631 * .. DSWAP ..
632  CALL dswap(n,sx,incx,sy,incy)
633  DO 80 i = 1, 7
634  stx(i) = dt10x(i,kn,ki)
635  sty(i) = dt10y(i,kn,ki)
636  80 CONTINUE
637  CALL stest(lenx,sx,stx,ssize2(1,1),1.0d0)
638  CALL stest(leny,sy,sty,ssize2(1,1),1.0d0)
639  ELSE IF (icase.EQ.12) THEN
640 * .. DROTM ..
641  kni=kn+4*(ki-1)
642  DO kpar=1,4
643  DO i=1,7
644  sx(i) = dx1(i)
645  sy(i) = dy1(i)
646  stx(i)= dt19x(i,kpar,kni)
647  sty(i)= dt19y(i,kpar,kni)
648  END DO
649 *
650  DO i=1,5
651  dtemp(i) = dpar(i,kpar)
652  END DO
653 *
654  DO i=1,lenx
655  ssize(i)=stx(i)
656  END DO
657 * SEE REMARK ABOVE ABOUT DT11X(1,2,7)
658 * AND DT11X(5,3,8).
659  IF ((kpar .EQ. 2) .AND. (kni .EQ. 7))
660  \$ ssize(1) = 2.4d0
661  IF ((kpar .EQ. 3) .AND. (kni .EQ. 8))
662  \$ ssize(5) = 1.8d0
663 *
664  CALL drotm(n,sx,incx,sy,incy,dtemp)
665  CALL stest(lenx,sx,stx,ssize,sfac)
666  CALL stest(leny,sy,sty,sty,sfac)
667  END DO
668  ELSE IF (icase.EQ.13) THEN
669 * .. DSDOT ..
670  CALL testdsdot(REAL(DSDOT(N,REAL(SX),INCX,REAL(SY),INCY)),
671  \$ REAL(DT7(KN,KI)),REAL(SSIZE1(KN)), .3125e-1)
672  ELSE
673  WRITE (nout,*) ' Shouldn''t be here in CHECK2'
674  stop
675  END IF
676  100 CONTINUE
677  120 CONTINUE
678  RETURN
subroutine dcopy(N, DX, INCX, DY, INCY)
DCOPY
Definition: dcopy.f:53
double precision function ddot(N, DX, INCX, DY, INCY)
DDOT
Definition: ddot.f:53
subroutine testdsdot(SCOMP, STRUE, SSIZE, SFAC)
Definition: dblat1.f:943
subroutine daxpy(N, DA, DX, INCX, DY, INCY)
DAXPY
Definition: daxpy.f:54
subroutine dswap(N, DX, INCX, DY, INCY)
DSWAP
Definition: dswap.f:53
subroutine drotm(N, DX, INCX, DY, INCY, DPARAM)
DROTM
Definition: drotm.f:100
subroutine stest1(SCOMP1, STRUE1, SSIZE, SFAC)
Definition: cblat1.f:620
subroutine stest(LEN, SCOMP, STRUE, SSIZE, SFAC)
Definition: cblat1.f:564
double precision function dsdot(N, SX, INCX, SY, INCY)
DSDOT
Definition: dsdot.f:121

Here is the call graph for this function: