**********************************************************************
     *                                                                    *
     *                       check copy constructor                       *
     *                                                                    *
     **********************************************************************
Copy constructor works


     **********************************************************************
     *                                                                    *
     *                           does [] work?                            *
     *                                                                    *
     *                                                                    *
     **********************************************************************
It does work


     **********************************************************************
     *                                                                    *
     *        a followed by T(a) then by realmat ww = T(a) + T(a)         *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


     **********************************************************************
     *                                                                    *
     *                         a has elements i+j                         *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                       Check self assignment                        *
     *                           aa=a,a=a,a==aa                           *
     *                             aa+=aa+=aa                             *
     *                              bb=aa=a                               *
     *                                                                    *
     **********************************************************************

     (a == aa) = 1


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     8.00000    12.00000    16.00000    20.00000    24.00000
       Row  2    12.00000    16.00000    20.00000    24.00000    28.00000
       Row  3    16.00000    20.00000    24.00000    28.00000    32.00000
       Row  4    20.00000    24.00000    28.00000    32.00000    36.00000
       Row  5    24.00000    28.00000    32.00000    36.00000    40.00000
       Row  6    28.00000    32.00000    36.00000    40.00000    44.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                           Check division                           *
     *                             a/5.0==a/5                             *
     *                            (1.0/5.0)*a                             *
     *                               a/5.0                                *
     *                                                                    *
     **********************************************************************

     (a/5.0 == a/5) = 1


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     0.40000     0.60000     0.80000     1.00000     1.20000
       Row  2     0.60000     0.80000     1.00000     1.20000     1.40000
       Row  3     0.80000     1.00000     1.20000     1.40000     1.60000
       Row  4     1.00000     1.20000     1.40000     1.60000     1.80000
       Row  5     1.20000     1.40000     1.60000     1.80000     2.00000
       Row  6     1.40000     1.60000     1.80000     2.00000     2.20000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     0.40000     0.60000     0.80000     1.00000     1.20000
       Row  2     0.60000     0.80000     1.00000     1.20000     1.40000
       Row  3     0.80000     1.00000     1.20000     1.40000     1.60000
       Row  4     1.00000     1.20000     1.40000     1.60000     1.80000
       Row  5     1.20000     1.40000     1.60000     1.80000     2.00000
       Row  6     1.40000     1.60000     1.80000     2.00000     2.20000


     **********************************************************************
     *                                                                    *
     *                           Check division                           *
     *                          T(a)/5.0==T(a)/5                          *
     *                           (1.0/5.0)*T(a)                           *
     *                              T(a)/5.0                              *
     *                                                                    *
     **********************************************************************

     (T(a)/5.0==T(a)/5) = 1


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     0.40000     0.60000     0.80000     1.00000     1.20000     1.40000
 Row  2     0.60000     0.80000     1.00000     1.20000     1.40000     1.60000
 Row  3     0.80000     1.00000     1.20000     1.40000     1.60000     1.80000
 Row  4     1.00000     1.20000     1.40000     1.60000     1.80000     2.00000
 Row  5     1.20000     1.40000     1.60000     1.80000     2.00000     2.20000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     0.40000     0.60000     0.80000     1.00000     1.20000     1.40000
 Row  2     0.60000     0.80000     1.00000     1.20000     1.40000     1.60000
 Row  3     0.80000     1.00000     1.20000     1.40000     1.60000     1.80000
 Row  4     1.00000     1.20000     1.40000     1.60000     1.80000     2.00000
 Row  5     1.20000     1.40000     1.60000     1.80000     2.00000     2.20000


     **********************************************************************
     *                                                                    *
     *                               T(a)*a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


     **********************************************************************
     *                                                                    *
     *                           invpsd(T(a)*a)                           *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.89524    -1.58095         0.0         0.0         0.0
       Row  2    -1.58095     1.32381         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                            inv(T(a)*a)                             *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     0.70952         0.0         0.0         0.0    -0.39524
       Row  2    -0.58810         0.0         0.0         0.0     0.33095
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                      a*invpsd(T(a)*a)*T(a)*a                       *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                        a*inv(T(a)*a)*T(a)*a                        *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                       a*ginv(T(a)*a)*T(a)*a                        *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                               T(a)*a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


     **********************************************************************
     *                                                                    *
     *                          inv(T(a)*a,rank)                          *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     0.70952         0.0         0.0         0.0    -0.39524
       Row  2    -0.58810         0.0         0.0         0.0     0.33095
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0

     rank = 2


     **********************************************************************
     *                                                                    *
     *                         ginv(T(a)*a,rank)                          *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     0.13524    0.076190    0.017143   -0.041905    -0.10095
       Row  2    0.076190    0.042952  0.00971429   -0.023524   -0.056762
       Row  3    0.017143  0.00971429  0.00228571 -0.00514286   -0.012571
       Row  4   -0.041905   -0.023524 -0.00514286    0.013238    0.031619
       Row  5    -0.10095   -0.056762   -0.012571    0.031619    0.075810

     rank = 2


     **********************************************************************
     *                                                                    *
     *                        invpsd(T(a)*a,rank)                         *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.89524    -1.58095         0.0         0.0         0.0
       Row  2    -1.58095     1.32381         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0

     rank = 2


     **********************************************************************
     *                                                                    *
     *                            five*T(a)*a                             *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     695.000     830.000     965.000     1100.00     1235.00
       Row  2     830.000     995.000     1160.00     1325.00     1490.00
       Row  3     965.000     1160.00     1355.00     1550.00     1745.00
       Row  4     1100.00     1325.00     1550.00     1775.00     2000.00
       Row  5     1235.00     1490.00     1745.00     2000.00     2255.00


     **********************************************************************
     *                                                                    *
     *                            T(a)*a*five                             *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     695.000     830.000     965.000     1100.00     1235.00
       Row  2     830.000     995.000     1160.00     1325.00     1490.00
       Row  3     965.000     1160.00     1355.00     1550.00     1745.00
       Row  4     1100.00     1325.00     1550.00     1775.00     2000.00
       Row  5     1235.00     1490.00     1745.00     2000.00     2255.00


     **********************************************************************
     *                                                                    *
     *                 a*invpsd(five*T(a)*a)*T(a)*a*five                  *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                   a*inv(T(a)*a*five)*five*T(a)*a                   *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                           Inner Product                            *
     *                                                                    *
     **********************************************************************

     left 

                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000

     rite 

                                           Col  1

                               Row  1     2.00000
                               Row  2     3.00000
                               Row  3     4.00000
                               Row  4     5.00000
                               Row  5     6.00000

     left*rite

                                           Col  1

                               Row  1    90.00000

     left*T(left)

                                           Col  1

                               Row  1    90.00000

     T(rite)*rite

                                           Col  1

                               Row  1    90.00000


     **********************************************************************
     *                                                                    *
     *                      a transpose, three ways:                      *
     *                             cout<<T(a)                             *
     *                             T(a)(i,j)                              *
     *                              T(a)[i]                               *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     2.00000     3.00000     4.00000     5.00000     6.00000     7.00000
 Row  2     3.00000     4.00000     5.00000     6.00000     7.00000     8.00000
 Row  3     4.00000     5.00000     6.00000     7.00000     8.00000     9.00000
 Row  4     5.00000     6.00000     7.00000     8.00000     9.00000    10.00000
 Row  5     6.00000     7.00000     8.00000     9.00000    10.00000    11.00000

 2 3 4 5 6 7
 3 4 5 6 7 8
 4 5 6 7 8 9
 5 6 7 8 9 10
 6 7 8 9 10 11

 2
 3
 4
 5
 6

 3
 4
 5
 6
 7

 4
 5
 6
 7
 8

 5
 6
 7
 8
 9

 6
 7
 8
 9
 10

 7
 8
 9
 10
 11



     **********************************************************************
     *                                                                    *
     *                              Ta=T(a)                               *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     2.00000     3.00000     4.00000     5.00000     6.00000     7.00000
 Row  2     3.00000     4.00000     5.00000     6.00000     7.00000     8.00000
 Row  3     4.00000     5.00000     6.00000     7.00000     8.00000     9.00000
 Row  4     5.00000     6.00000     7.00000     8.00000     9.00000    10.00000
 Row  5     6.00000     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                                T(a)                                *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     2.00000     3.00000     4.00000     5.00000     6.00000     7.00000
 Row  2     3.00000     4.00000     5.00000     6.00000     7.00000     8.00000
 Row  3     4.00000     5.00000     6.00000     7.00000     8.00000     9.00000
 Row  4     5.00000     6.00000     7.00000     8.00000     9.00000    10.00000
 Row  5     6.00000     7.00000     8.00000     9.00000    10.00000    11.00000


     **********************************************************************
     *                                                                    *
     *                                Ta*a                                *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


     **********************************************************************
     *                                                                    *
     *                               T(a)*a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


     **********************************************************************
     *                                                                    *
     *                                a*Ta                                *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    90.00000     110.000     130.000     150.000     170.000     190.000
 Row  2     110.000     135.000     160.000     185.000     210.000     235.000
 Row  3     130.000     160.000     190.000     220.000     250.000     280.000
 Row  4     150.000     185.000     220.000     255.000     290.000     325.000
 Row  5     170.000     210.000     250.000     290.000     330.000     370.000
 Row  6     190.000     235.000     280.000     325.000     370.000     415.000


     **********************************************************************
     *                                                                    *
     *                               a*T(a)                               *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    90.00000     110.000     130.000     150.000     170.000     190.000
 Row  2     110.000     135.000     160.000     185.000     210.000     235.000
 Row  3     130.000     160.000     190.000     220.000     250.000     280.000
 Row  4     150.000     185.000     220.000     255.000     290.000     325.000
 Row  5     170.000     210.000     250.000     290.000     330.000     370.000
 Row  6     190.000     235.000     280.000     325.000     370.000     415.000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    90.00000     110.000     130.000     150.000     170.000     190.000
 Row  2     110.000     135.000     160.000     185.000     210.000     235.000
 Row  3     130.000     160.000     190.000     220.000     250.000     280.000
 Row  4     150.000     185.000     220.000     255.000     290.000     325.000
 Row  5     170.000     210.000     250.000     290.000     330.000     370.000
 Row  6     190.000     235.000     280.000     325.000     370.000     415.000


     **********************************************************************
     *                                                                    *
     *                     a = upper triangle of ones                     *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2         0.0     1.00000     1.00000     1.00000     1.00000
       Row  3         0.0         0.0     1.00000     1.00000     1.00000
       Row  4         0.0         0.0         0.0     1.00000     1.00000
       Row  5         0.0         0.0         0.0         0.0     1.00000
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                               T(a)*a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2     1.00000     2.00000     2.00000     2.00000     2.00000
       Row  3     1.00000     2.00000     3.00000     3.00000     3.00000
       Row  4     1.00000     2.00000     3.00000     4.00000     4.00000
       Row  5     1.00000     2.00000     3.00000     4.00000     5.00000


     **********************************************************************
     *                                                                    *
     *                           invpsd(T(a)*a)                           *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000    -1.00000         0.0         0.0         0.0
       Row  2    -1.00000     2.00000    -1.00000  1.4316e-15 -8.0358e-16
       Row  3         0.0    -1.00000     2.00000    -1.00000  1.8633e-15
       Row  4         0.0  1.4316e-15    -1.00000     2.00000    -1.00000
       Row  5         0.0 -8.0358e-16  1.8633e-15    -1.00000     1.00000


     **********************************************************************
     *                                                                    *
     *                       T(a)*a*invpsd(T(a)*a)                        *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000  4.0599e-16 -3.5717e-16  4.4409e-16  2.2204e-16
       Row  2         0.0     1.00000 -7.1433e-16  8.8818e-16  4.4409e-16
       Row  3         0.0  9.9593e-16     1.00000         0.0  1.3323e-15
       Row  4         0.0  1.6240e-15 -2.3168e-15     1.00000         0.0
       Row  5         0.0  8.2039e-16 -4.5357e-16  1.7764e-15     1.00000


     **********************************************************************
     *                                                                    *
     *                              b=2.0*a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     2.00000     2.00000     2.00000     2.00000
       Row  2         0.0     2.00000     2.00000     2.00000     2.00000
       Row  3         0.0         0.0     2.00000     2.00000     2.00000
       Row  4         0.0         0.0         0.0     2.00000     2.00000
       Row  5         0.0         0.0         0.0         0.0     2.00000
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                              c = a+b                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     3.00000     3.00000     3.00000     3.00000     3.00000
       Row  2         0.0     3.00000     3.00000     3.00000     3.00000
       Row  3         0.0         0.0     3.00000     3.00000     3.00000
       Row  4         0.0         0.0         0.0     3.00000     3.00000
       Row  5         0.0         0.0         0.0         0.0     3.00000
       Row  6         0.0         0.0         0.0         0.0         0.0

	c has 6 rows and 5 columns.


     **********************************************************************
     *                                                                    *
     *                                T(c)                                *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     3.00000         0.0         0.0         0.0         0.0         0.0
 Row  2     3.00000     3.00000         0.0         0.0         0.0         0.0
 Row  3     3.00000     3.00000     3.00000         0.0         0.0         0.0
 Row  4     3.00000     3.00000     3.00000     3.00000         0.0         0.0
 Row  5     3.00000     3.00000     3.00000     3.00000     3.00000         0.0

	T(c) has 5 rows and 6 columns.

	T(c) has size 30
.
	T(c) has 5 rows and 6 columns.

	T(c) has size 30
.

     **********************************************************************
     *                                                                    *
     *                               c=c-c                                *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                              d(6,6,2)                              *
     *                                d*a                                 *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     4.00000     6.00000     8.00000    10.00000
       Row  2     2.00000     4.00000     6.00000     8.00000    10.00000
       Row  3     2.00000     4.00000     6.00000     8.00000    10.00000
       Row  4     2.00000     4.00000     6.00000     8.00000    10.00000
       Row  5     2.00000     4.00000     6.00000     8.00000    10.00000
       Row  6     2.00000     4.00000     6.00000     8.00000    10.00000


     **********************************************************************
     *                                                                    *
     *                             e(6,6,-1)                              *
     *                                e*a                                 *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000
       Row  2    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000
       Row  3    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000
       Row  4    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000
       Row  5    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000
       Row  6    -1.00000    -2.00000    -3.00000    -4.00000    -5.00000


     **********************************************************************
     *                                                                    *
     *                             -e and +e                              *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000
 Row  2     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000
 Row  3     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000
 Row  4     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000
 Row  5     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000
 Row  6     1.00000     1.00000     1.00000     1.00000     1.00000     1.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  2    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  3    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  4    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  5    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  6    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000


     **********************************************************************
     *                                                                    *
     *                               a += a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     2.00000     2.00000     2.00000     2.00000
       Row  2         0.0     2.00000     2.00000     2.00000     2.00000
       Row  3         0.0         0.0     2.00000     2.00000     2.00000
       Row  4         0.0         0.0         0.0     2.00000     2.00000
       Row  5         0.0         0.0         0.0         0.0     2.00000
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                               a -= a                               *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                           ++a  then --a                            *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  3     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  4     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  5     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  6     1.00000     1.00000     1.00000     1.00000     1.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                                                                    *
     *                               line 1                               *
     *                               line 2                               *
     *                               line 3                               *
     *                                                                    *
     *                                                                    *
     **********************************************************************


     **********************************************************************
     *                                                                    *
     *                                                                    *
     *                               line 1                               *
     *                               line 2                               *
     *                               line 3                               *
     *                                                                    *
     *                                                                    *
     **********************************************************************


     **********************************************************************
     *                                                                    *
     *                        zz=z where z is null                        *
     *                         before and after =                         *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0

                                  Null matrix


     **********************************************************************
     *                                                                    *
     *                    Try to invert a null matrix                     *
     *                                                                    *
     **********************************************************************
Error, realmat, invpsd, null matrix

     **********************************************************************
     *                                                                    *
     *                     Try to multiply two null's                     *
     *                                                                    *
     **********************************************************************
Error, realmat, operator *, null matrix

     **********************************************************************
     *                                                                    *
     *              Try to multiply nonconformable matrices               *
     *                                                                    *
     **********************************************************************
Error, realmat, operator *, matrices not conformable

     **********************************************************************
     *                                                                    *
     *                        test unary operators                        *
     *             e = -e; e = (-2)*e; e = (0.5)*e; e = +e;               *
     *                         e before and after                         *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  2    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  3    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  4    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  5    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  6    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  2    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  3    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  4    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  5    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  6    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000


     **********************************************************************
     *                                                                    *
     *                              e and a                               *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  2    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  3    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  4    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  5    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000
 Row  6    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000    -1.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *         e == e, e != e, a == e, a != e, nullmat == nullmat         *
     *                                                                    *
     **********************************************************************
 1 0 0 1 1

     **********************************************************************
     *                                                                    *
     *                   e <= e, e < e, a <= e, a < e,                    *
     *               nullmat <= nullmat, nullmat < nullmat                *
     *                                                                    *
     **********************************************************************
 1 0 0 0 1 0

     **********************************************************************
     *                                                                    *
     *                   e >= e, e > e, a >= e, a > e,                    *
     *               nullmat >= nullmat, nullmat > nullmat                *
     *                                                                    *
     **********************************************************************
 1 0 0 0 1 0

     **********************************************************************
     *                                                                    *
     *                  e <= -e, e < -e, e >= -e, e > -e                  *
     *                                                                    *
     **********************************************************************
 1 1 0 0

     **********************************************************************
     *                                                                    *
     *                   Try to increment a null matrix                   *
     *                                                                    *
     **********************************************************************
Error, realmat, operator ++, can't increment a null matrix

     **********************************************************************
     *                                                                    *
     *                   Try to decrement a null matrix                   *
     *                                                                    *
     **********************************************************************
Error, realmat, operator --, can't decrement a null matrix

     **********************************************************************
     *                                                                    *
     *               Check precedence ++a - a then a++ - a                *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  3     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  4     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  5     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  6     1.00000     1.00000     1.00000     1.00000     1.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  3     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  4     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  5     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  6     1.00000     1.00000     1.00000     1.00000     1.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  2     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  3     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  4     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  5     1.00000     1.00000     1.00000     1.00000     1.00000
       Row  6     1.00000     1.00000     1.00000     1.00000     1.00000



     **********************************************************************
     *                                                                    *
     *                        Check fill and get_x                        *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  2     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  3     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  4     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  5     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  6     5.00000     5.00000     5.00000     5.00000     5.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  2     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  3     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  4     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  5     5.00000     5.00000     5.00000     5.00000     5.00000
       Row  6     5.00000     5.00000     5.00000     5.00000     5.00000


     **********************************************************************
     *                                                                    *
     *                The matrix d and submatrix s(sr,sc)                 *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000

     sr = 1,2,4:6  sc = 1:3,5:6,1,1,1

             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    25.00000    31.00000     1.00000
 Row  2     2.00000     8.00000    14.00000    26.00000    32.00000     2.00000
 Row  3     4.00000    10.00000    16.00000    28.00000    34.00000     4.00000
 Row  4     5.00000    11.00000    17.00000    29.00000    35.00000     5.00000
 Row  5     6.00000    12.00000    18.00000    30.00000    36.00000     6.00000


             Col  7      Col  8

 Row  1     1.00000     1.00000
 Row  2     2.00000     2.00000
 Row  3     4.00000     4.00000
 Row  4     5.00000     5.00000
 Row  5     6.00000     6.00000



     **********************************************************************
     *                                                                    *
     *                          Try a bad string                          *
     *                The matrix d and submatrix s(sr,sc)                 *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0     5.00000    10.00000    15.00000    25.00000    30.00000
 Row  2         0.0     5.00000    10.00000    20.00000    25.00000    30.00000
 Row  3         0.0     5.00000    15.00000    20.00000    25.00000    30.00000
 Row  4         0.0    10.00000    15.00000    20.00000    25.00000    30.00000
 Row  5     5.00000    10.00000    15.00000    20.00000    25.00000    35.00000
 Row  6     5.00000    10.00000    15.00000    20.00000    30.00000    35.00000

     sr0 = 1,2,4:6  sc0 = 1:3,5:6-,8
Error, intvec, intvec, bad input string: 1:3,5:6-,8
                                  Null matrix



     **********************************************************************
     *                                                                    *
     *                                 d                                  *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000
s0, s1, s2, s4, s5 = 1,3,4 1:2   2:1   2:4   5     


     **********************************************************************
     *                                                                    *
     *                      cbind(d(s0,s1),d(s0,s2))                      *
     *                                                                    *
     **********************************************************************


                         Col  1      Col  2      Col  3      Col  4

             Row  1     1.00000     7.00000     7.00000     1.00000
             Row  2     3.00000     9.00000     9.00000     3.00000
             Row  3     4.00000    10.00000    10.00000     4.00000


     **********************************************************************
     *                                                                    *
     *                      rbind(d(s4,s0),d(s5,s0))                      *
     *                                                                    *
     **********************************************************************


                               Col  1      Col  2      Col  3

                   Row  1     2.00000    14.00000    20.00000
                   Row  2     3.00000    15.00000    21.00000
                   Row  3     4.00000    16.00000    22.00000
                   Row  4     5.00000    17.00000    23.00000


     **********************************************************************
     *                                                                    *
     *                            How's this?                             *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     1.00000     1.00000         0.0         0.0         0.0
 Row  2     1.00000     1.00000     1.00000         0.0         0.0         0.0
 Row  3     1.00000     1.00000     1.00000         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0     1.00000     1.00000     1.00000
 Row  5         0.0         0.0         0.0     1.00000     1.00000     1.00000
 Row  6         0.0         0.0         0.0     1.00000     1.00000     1.00000


                         Col  1      Col  2      Col  3      Col  4

             Row  1     1.00000     1.00000         0.0         0.0
             Row  2     1.00000     1.00000         0.0         0.0
             Row  3         0.0         0.0     1.00000     1.00000
             Row  4         0.0         0.0     1.00000     1.00000


     **********************************************************************
     *                                                                    *
     *                         Check subscripting                         *
     *                                                                    *
     **********************************************************************
1 1 1
1 1
2 2 2
7 7
3 3 3
13 13
4 4 4
19 19
5 5 5
25 25
6 6 6
31 31
7 7 7
2 2
8 8 8
8 8
9 9 9
14 14
10 10 10
20 20
11 11 11
26 26
12 12 12
32 32
13 13 13
3 3
14 14 14
9 9
15 15 15
15 15
16 16 16
21 21
17 17 17
27 27
18 18 18
33 33
19 19 19
4 4
20 20 20
10 10
21 21 21
16 16
22 22 22
22 22
23 23 23
28 28
24 24 24
34 34
25 25 25
5 5
26 26 26
11 11
27 27 27
17 17
28 28 28
23 23
29 29 29
29 29
30 30 30
35 35
31 31 31
6 6
32 32 32
12 12
33 33 33
18 18
34 34 34
24 24
35 35 35
30 30
36 36 36
36 36
1 1 1
1 1
7 7 7
2 2
13 13 13
3 3
19 19 19
4 4
25 25 25
5 5
31 31 31
6 6
2 2 2
7 7
8 8 8
8 8
14 14 14
9 9
20 20 20
10 10
26 26 26
11 11
32 32 32
12 12
3 3 3
13 13
9 9 9
14 14
15 15 15
15 15
21 21 21
16 16
27 27 27
17 17
33 33 33
18 18
4 4 4
19 19
10 10 10
20 20
16 16 16
21 21
22 22 22
22 22
28 28 28
23 23
34 34 34
24 24
5 5 5
25 25
11 11 11
26 26
17 17 17
27 27
23 23 23
28 28
29 29 29
29 29
35 35 35
30 30
6 6 6
31 31
12 12 12
32 32
18 18 18
33 33
24 24 24
34 34
30 30 30
35 35
36 36 36
36 36


     **********************************************************************
     *                                                                    *
     *              Here is d once again, then jvec and ivec              *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000

  2  3  4
  -5  -4  -3  -2  -1  0  1  2  3  4  5  6  7


     **********************************************************************
     *                                                                    *
     *                       Here is d(ivec,jvec)                         *
     *                                                                    *
     **********************************************************************


                               Col  1      Col  2      Col  3

                   Row  1     7.00000    13.00000    19.00000
                   Row  2     8.00000    14.00000    20.00000
                   Row  3     9.00000    15.00000    21.00000
                   Row  4    10.00000    16.00000    22.00000
                   Row  5    11.00000    17.00000    23.00000
                   Row  6    12.00000    18.00000    24.00000


     **********************************************************************
     *                                                                    *
     *              Here is d once again, then jvec and ivec              *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000

  -2  -3  -4
  -5  -4  -3  -2  -1  0  1  2  3  4  5  6  7


     **********************************************************************
     *                                                                    *
     *                       Here is d(ivec,jvec)                         *
     *                                                                    *
     **********************************************************************

                                  Null matrix


     **********************************************************************
     *                                                                    *
     *              Here is d once again, then jvec and ivec              *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000

  2  2  1  1
  -5  -4  -3  -2  -1  0  1  2  3  4  5  6  7


     **********************************************************************
     *                                                                    *
     *                       Here is d(ivec,jvec)                         *
     *                                                                    *
     **********************************************************************


                         Col  1      Col  2      Col  3      Col  4

             Row  1     7.00000     7.00000     1.00000     1.00000
             Row  2     8.00000     8.00000     2.00000     2.00000
             Row  3     9.00000     9.00000     3.00000     3.00000
             Row  4    10.00000    10.00000     4.00000     4.00000
             Row  5    11.00000    11.00000     5.00000     5.00000
             Row  6    12.00000    12.00000     6.00000     6.00000


     **********************************************************************
     *                                                                    *
     *                 First d, then d's rows and columns                 *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000


                                           Col  1

                               Row  1     1.00000
                               Row  2     2.00000
                               Row  3     3.00000
                               Row  4     4.00000
                               Row  5     5.00000
                               Row  6     6.00000


                                           Col  1

                               Row  1     7.00000
                               Row  2     8.00000
                               Row  3     9.00000
                               Row  4    10.00000
                               Row  5    11.00000
                               Row  6    12.00000


                                           Col  1

                               Row  1    13.00000
                               Row  2    14.00000
                               Row  3    15.00000
                               Row  4    16.00000
                               Row  5    17.00000
                               Row  6    18.00000


                                           Col  1

                               Row  1    19.00000
                               Row  2    20.00000
                               Row  3    21.00000
                               Row  4    22.00000
                               Row  5    23.00000
                               Row  6    24.00000


                                           Col  1

                               Row  1    25.00000
                               Row  2    26.00000
                               Row  3    27.00000
                               Row  4    28.00000
                               Row  5    29.00000
                               Row  6    30.00000


                                           Col  1

                               Row  1    31.00000
                               Row  2    32.00000
                               Row  3    33.00000
                               Row  4    34.00000
                               Row  5    35.00000
                               Row  6    36.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     7.00000    13.00000    19.00000    25.00000    31.00000
 Row  2     2.00000     8.00000    14.00000    20.00000    26.00000    32.00000
 Row  3     3.00000     9.00000    15.00000    21.00000    27.00000    33.00000
 Row  4     4.00000    10.00000    16.00000    22.00000    28.00000    34.00000
 Row  5     5.00000    11.00000    17.00000    23.00000    29.00000    35.00000
 Row  6     6.00000    12.00000    18.00000    24.00000    30.00000    36.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     2.00000     3.00000     4.00000     5.00000     6.00000
 Row  2     7.00000     8.00000     9.00000    10.00000    11.00000    12.00000
 Row  3    13.00000    14.00000    15.00000    16.00000    17.00000    18.00000
 Row  4    19.00000    20.00000    21.00000    22.00000    23.00000    24.00000
 Row  5    25.00000    26.00000    27.00000    28.00000    29.00000    30.00000
 Row  6    31.00000    32.00000    33.00000    34.00000    35.00000    36.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0
 Row  6         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                       check a*realmat(T(b))                        *
     *                                                                    *
     **********************************************************************


                                     Col  1      Col  2

                         Row  1     2166.00     2262.00
                         Row  2     2262.00     2364.00
                         Row  3     2358.00     2466.00
                         Row  4     2454.00     2568.00
                         Row  5     2550.00     2670.00
                         Row  6     2646.00     2772.00


                                     Col  1      Col  2

                         Row  1     2166.00     2262.00
                         Row  2     2262.00     2364.00
                         Row  3     2358.00     2466.00
                         Row  4     2454.00     2568.00
                         Row  5     2550.00     2670.00
                         Row  6     2646.00     2772.00


     **********************************************************************
     *                                                                    *
     *                       check realmat(T(a))*b                        *
     *                                                                    *
     **********************************************************************


                                     Col  1      Col  2

                         Row  1    91.00000     217.000
                         Row  2     217.000     559.000
                         Row  3     343.000     901.000
                         Row  4     469.000     1243.00
                         Row  5     595.000     1585.00
                         Row  6     721.000     1927.00


                                     Col  1      Col  2

                         Row  1    91.00000     217.000
                         Row  2     217.000     559.000
                         Row  3     343.000     901.000
                         Row  4     469.000     1243.00
                         Row  5     595.000     1585.00
                         Row  6     721.000     1927.00


     **********************************************************************
     *                                                                    *
     *                         check realmat_cmp                          *
     *                                                                    *
     **********************************************************************
false
true
false


     **********************************************************************
     *                                                                    *
     *                    check resize to null matrix                     *
     *                                                                    *
     **********************************************************************

                                  Null matrix


     **********************************************************************
     *                                                                    *
     *                         check push_back()                          *
     *                                                                    *
     **********************************************************************


                                           Col  1

                               Row  1     1.00000
                               Row  2     2.00000
                               Row  3     3.00000
                               Row  4     4.00000
                               Row  5     5.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     2.00000     3.00000     4.00000     5.00000     6.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     1.00000     2.00000     3.00000     4.00000     5.00000     6.00000


             Col  7      Col  8      Col  9      Col 10      Col 11      Col 12

 Row  1     1.00000     2.00000     3.00000     4.00000     5.00000     6.00000


     **********************************************************************
     *                                                                    *
     *                         check T(a) + T(a)                          *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                           check T(a) + b                           *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                           check a + T(b)                           *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     4.00000     6.00000     8.00000    10.00000    12.00000
       Row  2     6.00000     8.00000    10.00000    12.00000    14.00000
       Row  3     8.00000    10.00000    12.00000    14.00000    16.00000
       Row  4    10.00000    12.00000    14.00000    16.00000    18.00000
       Row  5    12.00000    14.00000    16.00000    18.00000    20.00000
       Row  6    14.00000    16.00000    18.00000    20.00000    22.00000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1     4.00000     6.00000     8.00000    10.00000    12.00000    14.00000
 Row  2     6.00000     8.00000    10.00000    12.00000    14.00000    16.00000
 Row  3     8.00000    10.00000    12.00000    14.00000    16.00000    18.00000
 Row  4    10.00000    12.00000    14.00000    16.00000    18.00000    20.00000
 Row  5    12.00000    14.00000    16.00000    18.00000    20.00000    22.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                            check  +T(b)                            *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     2.00000     3.00000     4.00000     5.00000     6.00000
       Row  2     3.00000     4.00000     5.00000     6.00000     7.00000
       Row  3     4.00000     5.00000     6.00000     7.00000     8.00000
       Row  4     5.00000     6.00000     7.00000     8.00000     9.00000
       Row  5     6.00000     7.00000     8.00000     9.00000    10.00000
       Row  6     7.00000     8.00000     9.00000    10.00000    11.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                         check T(a) - T(a)                          *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                           check T(a) - b                           *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                           check a - T(b)                           *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                            check  -T(b)                            *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1    -2.00000    -3.00000    -4.00000    -5.00000    -6.00000
       Row  2    -3.00000    -4.00000    -5.00000    -6.00000    -7.00000
       Row  3    -4.00000    -5.00000    -6.00000    -7.00000    -8.00000
       Row  4    -5.00000    -6.00000    -7.00000    -8.00000    -9.00000
       Row  5    -6.00000    -7.00000    -8.00000    -9.00000   -10.00000
       Row  6    -7.00000    -8.00000    -9.00000   -10.00000   -11.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1    -2.00000    -3.00000    -4.00000    -5.00000    -6.00000
       Row  2    -3.00000    -4.00000    -5.00000    -6.00000    -7.00000
       Row  3    -4.00000    -5.00000    -6.00000    -7.00000    -8.00000
       Row  4    -5.00000    -6.00000    -7.00000    -8.00000    -9.00000
       Row  5    -6.00000    -7.00000    -8.00000    -9.00000   -10.00000
       Row  6    -7.00000    -8.00000    -9.00000   -10.00000   -11.00000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0
       Row  6         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                          check T(a)*T(b)                           *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                            check T(a)*b                            *
     *                                                                    *
     **********************************************************************


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1     139.000     166.000     193.000     220.000     247.000
       Row  2     166.000     199.000     232.000     265.000     298.000
       Row  3     193.000     232.000     271.000     310.000     349.000
       Row  4     220.000     265.000     310.000     355.000     400.000
       Row  5     247.000     298.000     349.000     400.000     451.000


                   Col  1      Col  2      Col  3      Col  4      Col  5

       Row  1         0.0         0.0         0.0         0.0         0.0
       Row  2         0.0         0.0         0.0         0.0         0.0
       Row  3         0.0         0.0         0.0         0.0         0.0
       Row  4         0.0         0.0         0.0         0.0         0.0
       Row  5         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                            check a*T(b)                            *
     *                                                                    *
     **********************************************************************


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    90.00000     110.000     130.000     150.000     170.000     190.000
 Row  2     110.000     135.000     160.000     185.000     210.000     235.000
 Row  3     130.000     160.000     190.000     220.000     250.000     280.000
 Row  4     150.000     185.000     220.000     255.000     290.000     325.000
 Row  5     170.000     210.000     250.000     290.000     330.000     370.000
 Row  6     190.000     235.000     280.000     325.000     370.000     415.000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1    90.00000     110.000     130.000     150.000     170.000     190.000
 Row  2     110.000     135.000     160.000     185.000     210.000     235.000
 Row  3     130.000     160.000     190.000     220.000     250.000     280.000
 Row  4     150.000     185.000     220.000     255.000     290.000     325.000
 Row  5     170.000     210.000     250.000     290.000     330.000     370.000
 Row  6     190.000     235.000     280.000     325.000     370.000     415.000


             Col  1      Col  2      Col  3      Col  4      Col  5      Col  6

 Row  1         0.0         0.0         0.0         0.0         0.0         0.0
 Row  2         0.0         0.0         0.0         0.0         0.0         0.0
 Row  3         0.0         0.0         0.0         0.0         0.0         0.0
 Row  4         0.0         0.0         0.0         0.0         0.0         0.0
 Row  5         0.0         0.0         0.0         0.0         0.0         0.0
 Row  6         0.0         0.0         0.0         0.0         0.0         0.0


     **********************************************************************
     *                                                                    *
     *                          check iterators                           *
     *                                                                    *
     **********************************************************************

	 iterators work


     **********************************************************************
     *                                                                    *
     *           checks on matrices over and under cblas sizes            *
     *                                                                    *
     **********************************************************************
	1
	true
	true
	true
	2
	true
	true
	true
	3
	true
	true
	true
	4
	true
	true
	true
	5
	true
	true
	true
	6
	true
	true
	true
	7
	true
	true
	true
	true
	8
	true
	true
	true
	true