The R Project SVN R

Rev

Rev 4117 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4117 Rev 7185
Line 579... Line 579...
579 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 579 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
580 
C
 580 
C
581 
C     THIS VERSION DATED AUGUST 1983.
 581 
C     THIS VERSION DATED AUGUST 1983.
582 
C
 582 
C
583 
C     ------------------------------------------------------------------
 583 
C     ------------------------------------------------------------------
- 
 
 584 
c
- 
 
 585 
c     unnecessary initialization of L to keep g77 -Wall happy
- 
 
 586 
c
- 
 
 587 
      L = 0
584 
C
 588 
C
585 
      IERR = 0
 589 
      IERR = 0
586 
      IF (LOW .EQ. IGH) GO TO 180
 590 
      IF (LOW .EQ. IGH) GO TO 180
587 
C     .......... CREATE REAL SUBDIAGONAL ELEMENTS ..........
 591 
C     .......... CREATE REAL SUBDIAGONAL ELEMENTS ..........
588 
      L = LOW + 1
 592 
      L = LOW + 1
Line 817... Line 821...
817 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 821 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
818 
C
 822 
C
819 
C     THIS VERSION DATED AUGUST 1983.
 823 
C     THIS VERSION DATED AUGUST 1983.
820 
C
 824 
C
821 
C     ------------------------------------------------------------------
 825 
C     ------------------------------------------------------------------
- 
 
 826 
c
- 
 
 827 
c     unnecessary initialization of L to keep g77 -Wall happy
- 
 
 828 
c
- 
 
 829 
      L = 0
822 
C
 830 
C
823 
      IERR = 0
 831 
      IERR = 0
824 
C     .......... INITIALIZE EIGENVECTOR MATRIX ..........
 832 
C     .......... INITIALIZE EIGENVECTOR MATRIX ..........
825 
      DO 101 J = 1, N
 833 
      DO 101 J = 1, N
826 
C
 834 
C
Line 1567... Line 1575...
1567 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 1575 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
1568 
C
 1576 
C
1569 
C     THIS VERSION DATED AUGUST 1983.
 1577 
C     THIS VERSION DATED AUGUST 1983.
1570 
C
 1578 
C
1571 
C     ------------------------------------------------------------------
 1579 
C     ------------------------------------------------------------------
- 
 
 1580 
c
- 
 
 1581 
c     unnecessary initialization of L M P Q R to keep g77 -Wall happy
- 
 
 1582 
c
- 
 
 1583 
      L = 0
- 
 
 1584 
      M = 0
- 
 
 1585 
      P = 0.0D0
- 
 
 1586 
      Q = 0.0D0
- 
 
 1587 
      R = 0.0D0
1572 
C
 1588 
C
1573 
      IERR = 0
 1589 
      IERR = 0
1574 
      NORM = 0.0D0
 1590 
      NORM = 0.0D0
1575 
      K = 1
 1591 
      K = 1
1576 
C     .......... STORE ROOTS ISOLATED BY BALANC
 1592 
C     .......... STORE ROOTS ISOLATED BY BALANC
Line 1817... Line 1833...
1817 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 1833 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
1818 
C
 1834 
C
1819 
C     THIS VERSION DATED AUGUST 1983.
 1835 
C     THIS VERSION DATED AUGUST 1983.
1820 
C
 1836 
C
1821 
C     ------------------------------------------------------------------
 1837 
C     ------------------------------------------------------------------
- 
 
 1838 
c
- 
 
 1839 
c     unnecessary initialization of L M P R S to keep g77 -Wall happy
- 
 
 1840 
c
- 
 
 1841 
      L = 0
- 
 
 1842 
      M = 0
- 
 
 1843 
      P = 0.0D0
- 
 
 1844 
      R = 0.0D0
- 
 
 1845 
      S = 0.0D0
1822 
C
 1846 
C
1823 
      IERR = 0
 1847 
      IERR = 0
1824 
      NORM = 0.0D0
 1848 
      NORM = 0.0D0
1825 
      K = 1
 1849 
      K = 1
1826 
C     .......... STORE ROOTS ISOLATED BY BALANC
 1850 
C     .......... STORE ROOTS ISOLATED BY BALANC
Line 2504... Line 2528...
2504 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 2528 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
2505 
C
 2529 
C
2506 
C     THIS VERSION DATED AUGUST 1983.
 2530 
C     THIS VERSION DATED AUGUST 1983.
2507 
C
 2531 
C
2508 
C     ------------------------------------------------------------------
 2532 
C     ------------------------------------------------------------------
- 
 
 2533 
c
- 
 
 2534 
c     unnecessary initialization of C3 and S2 to keep g77 -Wall happy
- 
 
 2535 
c
- 
 
 2536 
      C3 = 0.0D0
- 
 
 2537 
      S2 = 0.0D0
2509 
C
 2538 
C
2510 
      IERR = 0
 2539 
      IERR = 0
2511 
      IF (N .EQ. 1) GO TO 1001
 2540 
      IF (N .EQ. 1) GO TO 1001
2512 
C
 2541 
C
2513 
      DO 100 I = 2, N
 2542 
      DO 100 I = 2, N
Line 2655... Line 2684...
2655 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 2684 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
2656 
C
 2685 
C
2657 
C     THIS VERSION DATED AUGUST 1983.
 2686 
C     THIS VERSION DATED AUGUST 1983.
2658 
C
 2687 
C
2659 
C     ------------------------------------------------------------------
 2688 
C     ------------------------------------------------------------------
- 
 
 2689 
c
- 
 
 2690 
c     unnecessary initialization of C3 and S2 to keep g77 -Wall happy
- 
 
 2691 
c
- 
 
 2692 
      C3 = 0.0D0
- 
 
 2693 
      S2 = 0.0D0
2660 
C
 2694 
C
2661 
      IERR = 0
 2695 
      IERR = 0
2662 
      IF (N .EQ. 1) GO TO 1001
 2696 
      IF (N .EQ. 1) GO TO 1001
2663 
C
 2697 
C
2664 
      DO 100 I = 2, N
 2698 
      DO 100 I = 2, N
Line 2807... Line 2841...
2807 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 2841 
C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
2808 
C
 2842 
C
2809 
C     THIS VERSION DATED AUGUST 1983.
 2843 
C     THIS VERSION DATED AUGUST 1983.
2810 
C
 2844 
C
2811 
C     ------------------------------------------------------------------
 2845 
C     ------------------------------------------------------------------
- 
 
 2846 
c
- 
 
 2847 
c     unnecessary initialization of B and C to keep g77 -Wall happy
- 
 
 2848 
c
- 
 
 2849 
      B = 0.0D0
- 
 
 2850 
      C = 0.0D0
2812 
C
 2851 
C
2813 
      IERR = 0
 2852 
      IERR = 0
2814 
      IF (N .EQ. 1) GO TO 1001
 2853 
      IF (N .EQ. 1) GO TO 1001
2815 
C
 2854 
C
2816 
      DO 100 I = 2, N
 2855 
      DO 100 I = 2, N