Error in Zhegvd

Open discussion for MAGMA

Error in Zhegvd

Postby fdahm » Fri Sep 02, 2011 4:50 am

First of all, thanks again for all the work done in magma.
I try the brand new zhegvd in my code, and it works well and fast for almost all my case.
However i got incorrect results with one case of very small size (N=26).
I submit the matrix who produce this incorrect result if it can helps
Code: Select all
matrix a:
a( 1, 1)=  3.0055E-01 -3.0466E-17
a( 2, 1)= -9.7361E-17 -2.1250E-17
a( 3, 1)= -1.3195E-16  1.0126E-16
a( 4, 1)=  1.2967E-16  1.1314E-16
a( 5, 1)= -5.2042E-18 -6.1800E-17
a( 6, 1)= -4.4398E-17  3.3285E-17
a( 7, 1)= -1.7347E-16  9.0531E-18
a( 8, 1)= -1.5613E-16 -2.2985E-17
a( 9, 1)= -1.6892E-16  3.1252E-17
a(10, 1)= -2.3852E-18  4.2501E-17
a(11, 1)=  2.3809E-16  6.1366E-17
a(12, 1)=  1.3227E-17  4.9873E-17
a(13, 1)=  3.0358E-17 -3.4694E-18
a(14, 1)=  0.0000E+00  0.0000E+00
a(15, 1)=  0.0000E+00  0.0000E+00
a(16, 1)=  0.0000E+00  0.0000E+00
a(17, 1)=  0.0000E+00  0.0000E+00
a(18, 1)=  0.0000E+00  0.0000E+00
a(19, 1)=  0.0000E+00  0.0000E+00
a(20, 1)=  0.0000E+00  0.0000E+00
a(21, 1)=  0.0000E+00  0.0000E+00
a(22, 1)=  0.0000E+00  0.0000E+00
a(23, 1)=  0.0000E+00  0.0000E+00
a(24, 1)=  0.0000E+00  0.0000E+00
a(25, 1)=  0.0000E+00  0.0000E+00
a(26, 1)=  0.0000E+00  0.0000E+00
a( 1, 2)= -8.8905E-17  3.0791E-17
a( 2, 2)=  3.6382E-01 -1.5884E-17
a( 3, 2)= -1.3878E-16  2.2551E-17
a( 4, 2)= -4.9440E-17 -1.7130E-17
a( 5, 2)= -1.0040E-16 -7.8063E-17
a( 6, 2)=  1.6263E-16  6.4727E-17
a( 7, 2)=  9.9747E-17  8.1749E-17
a( 8, 2)= -9.5843E-17 -8.3700E-17
a( 9, 2)=  5.5511E-17 -9.9747E-17
a(10, 2)= -1.0582E-16  2.4286E-17
a(11, 2)= -4.6187E-17 -9.7145E-17
a(12, 2)=  1.4745E-17  6.4618E-17
a(13, 2)=  2.2118E-17  6.5486E-17
a(14, 2)=  0.0000E+00  0.0000E+00
a(15, 2)=  0.0000E+00  0.0000E+00
a(16, 2)=  0.0000E+00  0.0000E+00
a(17, 2)=  0.0000E+00  0.0000E+00
a(18, 2)=  0.0000E+00  0.0000E+00
a(19, 2)=  0.0000E+00  0.0000E+00
a(20, 2)=  0.0000E+00  0.0000E+00
a(21, 2)=  0.0000E+00  0.0000E+00
a(22, 2)=  0.0000E+00  0.0000E+00
a(23, 2)=  0.0000E+00  0.0000E+00
a(24, 2)=  0.0000E+00  0.0000E+00
a(25, 2)=  0.0000E+00  0.0000E+00
a(26, 2)=  0.0000E+00  0.0000E+00
a( 1, 3)= -1.5282E-16 -1.5179E-18
a( 2, 3)= -1.6480E-16 -3.4044E-17
a( 3, 3)=  3.8240E-01  1.7076E-17
a( 4, 3)= -1.9932E-16  3.7513E-17
a( 5, 3)=  7.0690E-17 -6.4835E-17
a( 6, 3)=  1.0582E-16 -4.5862E-17
a( 7, 3)= -1.4745E-16 -8.9338E-17
a( 8, 3)= -6.0715E-18  4.3802E-17
a( 9, 3)= -6.2884E-17  1.3065E-16
a(10, 3)=  7.9364E-17 -9.5410E-18
a(11, 3)= -8.4568E-17  2.0600E-17
a(12, 3)= -1.1536E-16 -1.7347E-17
a(13, 3)=  1.7260E-16 -8.3917E-17
a(14, 3)=  0.0000E+00  0.0000E+00
a(15, 3)=  0.0000E+00  0.0000E+00
a(16, 3)=  0.0000E+00  0.0000E+00
a(17, 3)=  0.0000E+00  0.0000E+00
a(18, 3)=  0.0000E+00  0.0000E+00
a(19, 3)=  0.0000E+00  0.0000E+00
a(20, 3)=  0.0000E+00  0.0000E+00
a(21, 3)=  0.0000E+00  0.0000E+00
a(22, 3)=  0.0000E+00  0.0000E+00
a(23, 3)=  0.0000E+00  0.0000E+00
a(24, 3)=  0.0000E+00  0.0000E+00
a(25, 3)=  0.0000E+00  0.0000E+00
a(26, 3)=  0.0000E+00  0.0000E+00
a( 1, 4)=  4.5970E-17 -7.1991E-17
a( 2, 4)= -8.3971E-17  4.9873E-17
a( 3, 4)= -2.3072E-16 -1.8648E-17
a( 4, 4)=  4.1597E-01 -1.1493E-17
a( 5, 4)= -1.2837E-16 -8.4568E-18
a( 6, 4)=  7.7629E-17  6.6570E-17
a( 7, 4)= -8.0556E-17  8.2833E-17
a( 8, 4)=  6.5066E-17  7.8063E-18
a( 9, 4)= -2.1684E-19 -1.1102E-16
a(10, 4)=  7.9147E-17  3.9899E-17
a(11, 4)= -3.3122E-17 -3.2092E-17
a(12, 4)= -4.8139E-17  1.0408E-17
a(13, 4)=  3.4694E-17 -1.2935E-16
a(14, 4)=  0.0000E+00  0.0000E+00
a(15, 4)=  0.0000E+00  0.0000E+00
a(16, 4)=  0.0000E+00  0.0000E+00
a(17, 4)=  0.0000E+00  0.0000E+00
a(18, 4)=  0.0000E+00  0.0000E+00
a(19, 4)=  0.0000E+00  0.0000E+00
a(20, 4)=  0.0000E+00  0.0000E+00
a(21, 4)=  0.0000E+00  0.0000E+00
a(22, 4)=  0.0000E+00  0.0000E+00
a(23, 4)=  0.0000E+00  0.0000E+00
a(24, 4)=  0.0000E+00  0.0000E+00
a(25, 4)=  0.0000E+00  0.0000E+00
a(26, 4)=  0.0000E+00  0.0000E+00
a( 1, 5)=  9.2808E-17  3.4044E-17
a( 2, 5)= -1.1970E-16  8.2183E-17
a( 3, 5)=  2.5587E-17  7.4376E-17
a( 4, 5)= -5.1608E-17 -2.1034E-17
a( 5, 5)=  4.4218E-01 -3.1876E-17
a( 6, 5)=  3.1117E-17 -1.1905E-16
a( 7, 5)= -8.1966E-17  1.0148E-16
a( 8, 5)= -5.6704E-17  1.7239E-17
a( 9, 5)= -6.1366E-17 -1.1514E-16
a(10, 5)=  2.4026E-16  5.2042E-17
a(11, 5)=  6.2233E-17  1.7347E-17
a(12, 5)= -8.4134E-17  1.9787E-16
a(13, 5)= -3.0054E-16 -3.5649E-16
a(14, 5)=  0.0000E+00  0.0000E+00
a(15, 5)=  0.0000E+00  0.0000E+00
a(16, 5)=  0.0000E+00  0.0000E+00
a(17, 5)=  0.0000E+00  0.0000E+00
a(18, 5)=  0.0000E+00  0.0000E+00
a(19, 5)=  0.0000E+00  0.0000E+00
a(20, 5)=  0.0000E+00  0.0000E+00
a(21, 5)=  0.0000E+00  0.0000E+00
a(22, 5)=  0.0000E+00  0.0000E+00
a(23, 5)=  0.0000E+00  0.0000E+00
a(24, 5)=  0.0000E+00  0.0000E+00
a(25, 5)=  0.0000E+00  0.0000E+00
a(26, 5)=  0.0000E+00  0.0000E+00
a( 1, 6)= -1.2685E-17 -5.9089E-17
a( 2, 6)=  1.4875E-16 -8.9609E-17
a( 3, 6)=  7.8063E-17  6.5919E-17
a( 4, 6)=  1.0365E-16 -1.9732E-17
a( 5, 6)=  1.4745E-17  7.1557E-18
a( 6, 6)=  4.7056E-01 -5.4210E-18
a( 7, 6)= -1.2186E-16  1.0387E-16
a( 8, 6)=  3.3827E-17  8.8905E-17
a( 9, 6)=  1.2186E-16 -9.1073E-17
a(10, 6)=  1.1970E-16  6.1583E-17
a(11, 6)=  2.0817E-16  9.9367E-17
a(12, 6)=  1.4051E-16 -2.2855E-16
a(13, 6)= -3.7297E-16  8.7170E-17
a(14, 6)=  0.0000E+00  0.0000E+00
a(15, 6)=  0.0000E+00  0.0000E+00
a(16, 6)=  0.0000E+00  0.0000E+00
a(17, 6)=  0.0000E+00  0.0000E+00
a(18, 6)=  0.0000E+00  0.0000E+00
a(19, 6)=  0.0000E+00  0.0000E+00
a(20, 6)=  0.0000E+00  0.0000E+00
a(21, 6)=  0.0000E+00  0.0000E+00
a(22, 6)=  0.0000E+00  0.0000E+00
a(23, 6)=  0.0000E+00  0.0000E+00
a(24, 6)=  0.0000E+00  0.0000E+00
a(25, 6)=  0.0000E+00  0.0000E+00
a(26, 6)=  0.0000E+00  0.0000E+00
a( 1, 7)= -8.3267E-17  3.4694E-18
a( 2, 7)=  1.1796E-16 -7.7412E-17
a( 3, 7)= -1.2967E-16  7.1991E-17
a( 4, 7)= -4.5726E-17 -3.4044E-17
a( 5, 7)= -3.6429E-17 -7.4051E-17
a( 6, 7)= -1.4995E-16 -1.1709E-16
a( 7, 7)=  5.1871E-01  7.5894E-19
a( 8, 7)=  3.1258E-16 -3.9899E-17
a( 9, 7)= -1.4843E-16 -3.6646E-17
a(10, 7)= -1.7391E-16 -9.3241E-17
a(11, 7)=  1.1146E-16  6.5919E-17
a(12, 7)= -1.7174E-16 -7.8063E-18
a(13, 7)=  1.5179E-16 -2.2378E-16
a(14, 7)=  0.0000E+00  0.0000E+00
a(15, 7)=  0.0000E+00  0.0000E+00
a(16, 7)=  0.0000E+00  0.0000E+00
a(17, 7)=  0.0000E+00  0.0000E+00
a(18, 7)=  0.0000E+00  0.0000E+00
a(19, 7)=  0.0000E+00  0.0000E+00
a(20, 7)=  0.0000E+00  0.0000E+00
a(21, 7)=  0.0000E+00  0.0000E+00
a(22, 7)=  0.0000E+00  0.0000E+00
a(23, 7)=  0.0000E+00  0.0000E+00
a(24, 7)=  0.0000E+00  0.0000E+00
a(25, 7)=  0.0000E+00  0.0000E+00
a(26, 7)=  0.0000E+00  0.0000E+00
a( 1, 8)= -8.1749E-17  8.2833E-17
a( 2, 8)= -1.3010E-16  1.0408E-16
a( 3, 8)=  2.1684E-18 -7.5894E-17
a( 4, 8)=  8.2616E-17  2.2985E-17
a( 5, 8)= -8.3700E-17  6.1800E-17
a( 6, 8)=  9.7578E-18 -1.1102E-16
a( 7, 8)=  3.3545E-16 -1.7347E-17
a( 8, 8)=  5.4370E-01  3.3258E-17
a( 9, 8)=  4.2696E-16 -1.3878E-17
a(10, 8)=  3.9899E-17  1.3563E-16
a(11, 8)= -3.1225E-17  3.6917E-17
a(12, 8)= -2.0318E-16  3.9465E-17
a(13, 8)=  1.1709E-17 -1.8757E-16
a(14, 8)=  0.0000E+00  0.0000E+00
a(15, 8)=  0.0000E+00  0.0000E+00
a(16, 8)=  0.0000E+00  0.0000E+00
a(17, 8)=  0.0000E+00  0.0000E+00
a(18, 8)=  0.0000E+00  0.0000E+00
a(19, 8)=  0.0000E+00  0.0000E+00
a(20, 8)=  0.0000E+00  0.0000E+00
a(21, 8)=  0.0000E+00  0.0000E+00
a(22, 8)=  0.0000E+00  0.0000E+00
a(23, 8)=  0.0000E+00  0.0000E+00
a(24, 8)=  0.0000E+00  0.0000E+00
a(25, 8)=  0.0000E+00  0.0000E+00
a(26, 8)=  0.0000E+00  0.0000E+00
a( 1, 9)= -1.7748E-16 -4.4019E-17
a( 2, 9)=  4.8139E-17  6.2667E-17
a( 3, 9)= -9.2808E-17 -1.5106E-16
a( 4, 9)= -7.9635E-17  7.9797E-17
a( 5, 9)= -1.0777E-16  4.6295E-17
a( 6, 9)=  5.6812E-17  4.0766E-17
a( 7, 9)= -1.1753E-16 -2.3310E-17
a( 8, 9)=  3.7080E-16 -6.1149E-17
a( 9, 9)=  5.8379E-01  4.0820E-17
a(10, 9)=  2.3614E-16  1.2197E-16
a(11, 9)= -1.5840E-16  7.8063E-18
a(12, 9)= -1.7521E-16  3.2092E-17
a(13, 9)=  5.9848E-17 -2.4796E-16
a(14, 9)=  0.0000E+00  0.0000E+00
a(15, 9)=  0.0000E+00  0.0000E+00
a(16, 9)=  0.0000E+00  0.0000E+00
a(17, 9)=  0.0000E+00  0.0000E+00
a(18, 9)=  0.0000E+00  0.0000E+00
a(19, 9)=  0.0000E+00  0.0000E+00
a(20, 9)=  0.0000E+00  0.0000E+00
a(21, 9)=  0.0000E+00  0.0000E+00
a(22, 9)=  0.0000E+00  0.0000E+00
a(23, 9)=  0.0000E+00  0.0000E+00
a(24, 9)=  0.0000E+00  0.0000E+00
a(25, 9)=  0.0000E+00  0.0000E+00
a(26, 9)=  0.0000E+00  0.0000E+00
a( 1,10)= -5.8222E-17 -6.0715E-17
a( 2,10)= -1.1579E-16 -2.3419E-17
a( 3,10)=  6.8088E-17 -1.1276E-17
a( 4,10)=  8.6953E-17 -5.3776E-17
a( 5,10)=  1.3227E-16  8.2399E-18
a( 6,10)=  1.4919E-16 -1.4138E-16
a( 7,10)= -1.9592E-16  5.5294E-17
a( 8,10)=  4.3368E-18 -1.4919E-16
a( 9,10)=  2.3332E-16 -5.6433E-17
a(10,10)=  6.6843E-01  4.9114E-17
a(11,10)= -1.9299E-16  2.9761E-17
a(12,10)= -1.9299E-16  5.4102E-17
a(13,10)=  1.1883E-16 -3.7513E-17
a(14,10)=  0.0000E+00  0.0000E+00
a(15,10)=  0.0000E+00  0.0000E+00
a(16,10)=  0.0000E+00  0.0000E+00
a(17,10)=  0.0000E+00  0.0000E+00
a(18,10)=  0.0000E+00  0.0000E+00
a(19,10)=  0.0000E+00  0.0000E+00
a(20,10)=  0.0000E+00  0.0000E+00
a(21,10)=  0.0000E+00  0.0000E+00
a(22,10)=  0.0000E+00  0.0000E+00
a(23,10)=  0.0000E+00  0.0000E+00
a(24,10)=  0.0000E+00  0.0000E+00
a(25,10)=  0.0000E+00  0.0000E+00
a(26,10)=  0.0000E+00  0.0000E+00
a( 1,11)=  1.4854E-16 -9.1832E-17
a( 2,11)= -2.7972E-17  7.7629E-17
a( 3,11)= -8.3267E-17 -6.4510E-17
a( 4,11)= -4.9656E-17  9.4976E-17
a( 5,11)=  5.9848E-17 -2.4286E-17
a( 6,11)=  1.7000E-16 -1.8459E-16
a( 7,11)=  2.6129E-17 -4.9440E-17
a( 8,11)= -1.8648E-17 -5.8113E-17
a( 9,11)= -1.4138E-16  1.0408E-17
a(10,11)= -1.6989E-16 -2.3852E-18
a(11,11)=  8.2960E-01  4.1986E-17
a(12,11)=  1.5179E-16  8.2399E-17
a(13,11)= -1.1243E-16  6.3751E-17
a(14,11)=  0.0000E+00  0.0000E+00
a(15,11)=  0.0000E+00  0.0000E+00
a(16,11)=  0.0000E+00  0.0000E+00
a(17,11)=  0.0000E+00  0.0000E+00
a(18,11)=  0.0000E+00  0.0000E+00
a(19,11)=  0.0000E+00  0.0000E+00
a(20,11)=  0.0000E+00  0.0000E+00
a(21,11)=  0.0000E+00  0.0000E+00
a(22,11)=  0.0000E+00  0.0000E+00
a(23,11)=  0.0000E+00  0.0000E+00
a(24,11)=  0.0000E+00  0.0000E+00
a(25,11)=  0.0000E+00  0.0000E+00
a(26,11)=  0.0000E+00  0.0000E+00
a( 1,12)=  9.5952E-17 -2.7322E-17
a( 2,12)=  9.5410E-18 -5.5077E-17
a( 3,12)= -8.5652E-17  2.9490E-17
a( 4,12)= -2.2985E-17  3.4694E-18
a( 5,12)= -1.5959E-16 -7.9797E-17
a( 6,12)=  3.2960E-17  1.4566E-16
a( 7,12)= -1.3921E-16 -1.1276E-17
a( 8,12)= -2.4048E-16  3.9031E-18
a( 9,12)= -1.3791E-16 -1.5613E-17
a(10,12)= -1.4171E-16 -3.6646E-17
a(11,12)=  1.8041E-16 -1.2620E-16
a(12,12)=  9.6320E-01  1.1265E-16
a(13,12)= -7.0256E-17 -1.6632E-16
a(14,12)=  0.0000E+00  0.0000E+00
a(15,12)=  0.0000E+00  0.0000E+00
a(16,12)=  0.0000E+00  0.0000E+00
a(17,12)=  0.0000E+00  0.0000E+00
a(18,12)=  0.0000E+00  0.0000E+00
a(19,12)=  0.0000E+00  0.0000E+00
a(20,12)=  0.0000E+00  0.0000E+00
a(21,12)=  0.0000E+00  0.0000E+00
a(22,12)=  0.0000E+00  0.0000E+00
a(23,12)=  0.0000E+00  0.0000E+00
a(24,12)=  0.0000E+00  0.0000E+00
a(25,12)=  0.0000E+00  0.0000E+00
a(26,12)=  0.0000E+00  0.0000E+00
a( 1,13)=  5.6433E-17 -2.4286E-17
a( 2,13)= -3.8164E-17 -1.4615E-16
a( 3,13)=  3.0358E-17 -1.0755E-16
a( 4,13)=  4.2067E-17 -5.1933E-17
a( 5,13)=  5.9631E-19 -9.0856E-17
a( 6,13)=  1.0669E-16 -5.2042E-17
a( 7,13)=  1.2837E-16 -1.2143E-17
a( 8,13)=  1.2143E-17 -2.4720E-17
a( 9,13)= -1.7044E-16  2.1174E-16
a(10,13)=  6.3101E-17  1.0148E-16
a(11,13)= -3.8858E-16 -9.9963E-17
a(12,13)= -8.9338E-17 -5.8547E-18
a(13,13)=  1.1121E+00  9.7958E-17
a(14,13)=  0.0000E+00  0.0000E+00
a(15,13)=  0.0000E+00  0.0000E+00
a(16,13)=  0.0000E+00  0.0000E+00
a(17,13)=  0.0000E+00  0.0000E+00
a(18,13)=  0.0000E+00  0.0000E+00
a(19,13)=  0.0000E+00  0.0000E+00
a(20,13)=  0.0000E+00  0.0000E+00
a(21,13)=  0.0000E+00  0.0000E+00
a(22,13)=  0.0000E+00  0.0000E+00
a(23,13)=  0.0000E+00  0.0000E+00
a(24,13)=  0.0000E+00  0.0000E+00
a(25,13)=  0.0000E+00  0.0000E+00
a(26,13)=  0.0000E+00  0.0000E+00
a( 1,14)=  4.5617E-01  2.7918E-18
a( 2,14)= -1.7008E-03  7.1380E-03
a( 3,14)= -6.4922E-03  1.5873E-02
a( 4,14)=  2.0612E-02  1.7606E-02
a( 5,14)= -1.6200E-02 -5.9392E-03
a( 6,14)= -3.7463E-03 -3.7213E-04
a( 7,14)=  4.6655E-03 -1.0786E-03
a( 8,14)=  1.5898E-02  1.0385E-02
a( 9,14)=  1.3939E-02 -1.5992E-02
a(10,14)= -1.3338E-02 -1.6825E-03
a(11,14)=  9.7969E-03 -8.5816E-04
a(12,14)= -5.9013E-03  3.1637E-03
a(13,14)= -1.6774E-02  7.5404E-02
a(14,14)=  1.2188E+00  6.1854E-17
a(15,14)=  1.9952E-02 -1.4445E-02
a(16,14)= -1.9068E-02 -5.0621E-03
a(17,14)=  2.3616E-02 -2.2899E-03
a(18,14)= -7.5032E-03 -2.8950E-02
a(19,14)= -5.1700E-02 -2.9616E-02
a(20,14)=  3.5164E-02 -3.9827E-02
a(21,14)= -4.5942E-02  5.8231E-02
a(22,14)= -1.3699E-02 -1.0679E-01
a(23,14)= -2.2654E-02 -3.3386E-02
a(24,14)=  3.6913E-03 -3.4083E-02
a(25,14)= -2.1019E-03 -2.9005E-02
a(26,14)=  6.1151E-03 -6.6322E-02
a( 1,15)= -5.4984E-04  8.3210E-03
a( 2,15)=  4.6203E-01  1.0327E-04
a( 3,15)=  9.0982E-03 -1.6336E-02
a( 4,15)=  1.8395E-02  9.9966E-03
a( 5,15)=  1.5555E-03 -3.0044E-02
a( 6,15)= -4.7462E-03 -7.9789E-03
a( 7,15)= -2.6738E-02  9.1075E-03
a( 8,15)= -1.9243E-02 -2.6127E-02
a( 9,15)=  9.7271E-03  7.1269E-03
a(10,15)=  1.7626E-02 -4.6594E-03
a(11,15)=  5.3992E-02 -4.8575E-03
a(12,15)= -1.1041E-02 -8.9671E-03
a(13,15)= -4.1520E-02 -9.1484E-03
a(14,15)=  1.9952E-02  1.4445E-02
a(15,15)=  9.0270E-01  3.0791E-17
a(16,15)=  4.7900E-03 -3.1708E-02
a(17,15)=  4.8457E-02  2.6428E-02
a(18,15)= -9.7077E-03 -4.0872E-02
a(19,15)= -1.3520E-02 -3.1460E-02
a(20,15)= -3.2383E-02 -2.7021E-02
a(21,15)= -1.5540E-02 -1.1614E-02
a(22,15)= -5.8758E-03  1.4040E-02
a(23,15)=  4.0432E-02 -4.1111E-02
a(24,15)=  3.5519E-02 -1.9683E-02
a(25,15)= -3.3795E-03  1.2219E-03
a(26,15)= -6.0544E-02 -6.7730E-03
a( 1,16)= -1.9201E-03  3.1147E-03
a( 2,16)=  7.5483E-04  9.4789E-03
a( 3,16)=  5.1282E-01 -5.0212E-05
a( 4,16)= -3.0382E-03 -1.5462E-02
a( 5,16)=  1.1870E-03  4.4019E-03
a( 6,16)=  1.3485E-03 -2.6720E-02
a( 7,16)=  1.3382E-02 -1.4379E-02
a( 8,16)= -9.8592E-04 -9.4790E-03
a( 9,16)= -2.1937E-02 -6.0945E-03
a(10,16)= -1.4658E-03  3.5846E-03
a(11,16)=  1.8008E-02  1.4990E-03
a(12,16)=  1.7253E-02  1.8179E-04
a(13,16)= -1.7847E-02  3.5067E-02
a(14,16)= -1.9068E-02  5.0621E-03
a(15,16)=  4.7900E-03  3.1708E-02
a(16,16)=  1.2015E+00 -1.3010E-17
a(17,16)=  1.7324E-02 -3.3080E-02
a(18,16)= -5.5480E-03  1.0559E-02
a(19,16)=  9.1335E-03 -4.7573E-02
a(20,16)=  2.5542E-03  4.7500E-04
a(21,16)=  1.1249E-02 -5.0906E-02
a(22,16)=  2.8881E-03  4.2306E-03
a(23,16)=  1.1931E-02 -7.3736E-03
a(24,16)=  1.3198E-02  2.8620E-03
a(25,16)=  1.3512E-02 -2.6299E-02
a(26,16)=  5.0408E-03  4.3310E-02
a( 1,17)= -4.0710E-04 -7.7040E-04
a( 2,17)=  7.6684E-03 -5.1995E-03
a( 3,17)=  4.6679E-03  1.0754E-03
a( 4,17)=  5.1899E-01 -3.0816E-05
a( 5,17)= -1.2630E-03  1.9415E-03
a( 6,17)= -1.4095E-03 -1.2932E-02
a( 7,17)=  2.7275E-03 -7.2216E-03
a( 8,17)=  1.3867E-02 -6.3510E-03
a( 9,17)= -1.9169E-02 -8.3415E-03
a(10,17)= -1.5984E-03  2.9873E-03
a(11,17)= -1.1375E-02  1.5405E-02
a(12,17)=  1.2604E-02  1.4038E-02
a(13,17)= -9.5623E-02 -3.4365E-02
a(14,17)=  2.3616E-02  2.2899E-03
a(15,17)=  4.8457E-02 -2.6428E-02
a(16,17)=  1.7324E-02  3.3080E-02
a(17,17)=  1.1860E+00 -1.8648E-17
a(18,17)= -1.9645E-02 -1.9618E-02
a(19,17)= -3.0540E-04 -3.5618E-02
a(20,17)=  1.3485E-02 -8.3844E-03
a(21,17)= -5.9828E-03  2.6750E-03
a(22,17)= -8.5402E-03  1.2914E-02
a(23,17)=  1.7710E-02 -2.0182E-02
a(24,17)= -3.0196E-03 -1.1503E-02
a(25,17)=  3.8478E-02  2.3090E-03
a(26,17)= -9.6095E-02 -1.9265E-02
a( 1,18)= -2.0064E-03  9.8480E-04
a( 2,18)= -7.7889E-04  5.2859E-04
a( 3,18)= -3.0272E-04  2.0560E-03
a( 4,18)= -5.8374E-03  1.3442E-03
a( 5,18)=  5.5018E-01 -6.5721E-05
a( 6,18)= -2.2987E-02  5.4516E-03
a( 7,18)= -7.1588E-03  1.4484E-02
a( 8,18)=  2.4296E-02  1.5248E-02
a( 9,18)= -3.8024E-02  2.8488E-02
a(10,18)= -8.0991E-03 -2.2841E-02
a(11,18)= -3.3315E-02 -1.8671E-02
a(12,18)=  3.8216E-02 -9.3058E-03
a(13,18)= -2.9590E-02  6.6754E-02
a(14,18)= -7.5032E-03  2.8950E-02
a(15,18)= -9.7077E-03  4.0872E-02
a(16,18)= -5.5480E-03 -1.0559E-02
a(17,18)= -1.9645E-02  1.9618E-02
a(18,18)=  1.2397E+00  3.6212E-17
a(19,18)= -3.0620E-02  1.7459E-02
a(20,18)=  1.7058E-03  1.6742E-03
a(21,18)=  3.2858E-03  1.2264E-02
a(22,18)= -7.6113E-02  9.1631E-03
a(23,18)=  2.8322E-03 -1.5462E-02
a(24,18)= -1.9898E-03 -1.6409E-02
a(25,18)=  4.5580E-02  4.5170E-03
a(26,18)= -2.2874E-02  6.0332E-02
a( 1,19)= -9.1928E-03  6.2711E-03
a( 2,19)= -6.7156E-03  7.5104E-03
a( 3,19)=  4.9334E-03  5.2277E-03
a( 4,19)= -3.2588E-03  1.1759E-02
a( 5,19)= -9.6643E-03  3.9155E-03
a( 6,19)=  5.6422E-01 -3.9358E-04
a( 7,19)= -1.4274E-02 -1.5717E-02
a( 8,19)=  2.7374E-02  5.7493E-03
a( 9,19)= -3.1008E-02  2.1516E-02
a(10,19)=  2.1692E-02 -8.5180E-03
a(11,19)=  4.5824E-03 -1.5251E-02
a(12,19)=  1.8155E-02 -2.7292E-02
a(13,19)=  1.5438E-02  1.7907E-02
a(14,19)= -5.1700E-02  2.9616E-02
a(15,19)= -1.3520E-02  3.1460E-02
a(16,19)=  9.1335E-03  4.7573E-02
a(17,19)= -3.0540E-04  3.5618E-02
a(18,19)= -3.0620E-02 -1.7459E-02
a(19,19)=  1.2121E+00 -1.0679E-17
a(20,19)= -4.3192E-03 -6.2658E-02
a(21,19)=  5.6437E-02  3.5819E-02
a(22,19)= -9.3580E-02  5.4438E-02
a(23,19)=  3.4426E-02 -1.8101E-02
a(24,19)=  3.4166E-02  9.2572E-03
a(25,19)= -1.7671E-02 -6.1010E-02
a(26,19)=  1.5552E-02  3.8042E-02
a( 1,20)=  3.4777E-03  7.3032E-03
a( 2,20)= -8.6693E-03  3.3129E-04
a( 3,20)=  7.4785E-04 -3.2239E-03
a( 4,20)=  3.8178E-03 -4.5019E-03
a( 5,20)= -2.3323E-03 -1.2696E-03
a( 6,20)=  1.9859E-03  1.3853E-02
a( 7,20)=  5.5545E-01 -7.1685E-04
a( 8,20)=  2.2899E-02  2.4200E-02
a( 9,20)= -2.7483E-02  5.5986E-02
a(10,20)= -3.4813E-02 -4.9999E-02
a(11,20)=  2.6442E-02 -1.6320E-02
a(12,20)=  4.1716E-02  1.7916E-03
a(13,20)=  1.4683E-03  4.2168E-02
a(14,20)=  3.5164E-02  3.9827E-02
a(15,20)= -3.2383E-02  2.7021E-02
a(16,20)=  2.5542E-03 -4.7500E-04
a(17,20)=  1.3485E-02  8.3844E-03
a(18,20)=  1.7058E-03 -1.6742E-03
a(19,20)= -4.3192E-03  6.2658E-02
a(20,20)=  1.0996E+00  9.6494E-18
a(21,20)=  4.7974E-02  2.0378E-02
a(22,20)= -3.9933E-02  8.0533E-02
a(23,20)= -2.7113E-02 -2.5206E-02
a(24,20)=  3.5971E-02 -1.6749E-02
a(25,20)=  1.9266E-02 -9.0752E-03
a(26,20)=  1.8178E-02  9.6824E-02
a( 1,21)= -6.2550E-03 -5.5394E-03
a( 2,21)= -1.2592E-03 -2.9399E-03
a( 3,21)=  3.7377E-03  9.6945E-03
a( 4,21)= -1.2108E-03  1.6273E-03
a( 5,21)= -8.7528E-04 -4.4239E-03
a( 6,21)=  1.1089E-02 -2.0136E-03
a( 7,21)=  1.2203E-02 -4.9330E-03
a( 8,21)=  5.7689E-01 -4.6329E-04
a( 9,21)=  4.3673E-02 -2.4188E-02
a(10,21)=  9.6791E-03  2.1015E-02
a(11,21)= -2.5556E-02  3.7100E-02
a(12,21)= -3.0207E-02  6.1290E-03
a(13,21)=  2.7292E-02  3.2171E-03
a(14,21)= -4.5942E-02 -5.8231E-02
a(15,21)= -1.5540E-02  1.1614E-02
a(16,21)=  1.1249E-02  5.0906E-02
a(17,21)= -5.9828E-03 -2.6750E-03
a(18,21)=  3.2858E-03 -1.2264E-02
a(19,21)=  5.6437E-02 -3.5819E-02
a(20,21)=  4.7974E-02 -2.0378E-02
a(21,21)=  1.1644E+00 -2.2958E-17
a(22,21)=  7.8565E-02 -9.3241E-03
a(23,21)= -8.0908E-03  3.9720E-02
a(24,21)= -1.7437E-02  4.6827E-02
a(25,21)=  1.6223E-02 -1.7327E-02
a(26,21)= -1.7430E-02 -4.3046E-03
a( 1,22)=  9.3013E-04  8.6659E-03
a( 2,22)= -1.2239E-03 -2.1618E-03
a( 3,22)= -6.7991E-04 -1.6131E-03
a( 4,22)=  2.2812E-04  1.8380E-03
a( 5,22)= -1.0391E-02  7.8498E-04
a( 6,22)= -1.3395E-02 -1.2231E-02
a( 7,22)= -9.1809E-03 -1.6046E-02
a( 8,22)=  7.5760E-03 -1.4479E-03
a( 9,22)=  5.6393E-01 -2.4030E-04
a(10,22)=  1.3328E-02 -3.9003E-03
a(11,22)=  8.9142E-02 -5.9044E-02
a(12,22)=  5.2077E-02  3.9445E-02
a(13,22)= -1.4021E-02 -8.6058E-02
a(14,22)= -1.3699E-02  1.0679E-01
a(15,22)= -5.8758E-03 -1.4040E-02
a(16,22)=  2.8881E-03 -4.2306E-03
a(17,22)= -8.5402E-03 -1.2914E-02
a(18,22)= -7.6113E-02 -9.1631E-03
a(19,22)= -9.3580E-02 -5.4438E-02
a(20,22)= -3.9933E-02 -8.0533E-02
a(21,22)=  7.8565E-02  9.3241E-03
a(22,22)=  1.0602E+00 -1.6588E-17
a(23,22)=  2.4115E-02  2.3906E-02
a(24,22)=  1.0594E-01 -7.4607E-02
a(25,22)= -1.1312E-02  3.9615E-02
a(26,22)=  1.7289E-02 -9.1186E-02
a( 1,23)= -4.6053E-03  5.6771E-03
a( 2,23)=  1.0524E-02  6.1212E-03
a( 3,23)=  1.7687E-03 -3.5173E-03
a( 4,23)=  4.1253E-03  1.3669E-03
a( 5,23)= -8.7845E-04 -5.3192E-03
a( 6,23)=  4.1660E-03  3.8328E-03
a( 7,23)= -3.0289E-03  3.1960E-03
a( 8,23)= -2.3322E-03 -1.3564E-02
a( 9,23)= -1.6397E-04 -9.6081E-03
a(10,23)=  5.1643E-01 -2.5787E-04
a(11,23)=  1.6704E-01  6.4501E-02
a(12,23)= -2.8497E-03 -3.7774E-02
a(13,23)= -6.9129E-02 -9.4374E-02
a(14,23)= -2.2654E-02  3.3386E-02
a(15,23)=  4.0432E-02  4.1111E-02
a(16,23)=  1.1931E-02  7.3736E-03
a(17,23)=  1.7710E-02  2.0182E-02
a(18,23)=  2.8322E-03  1.5462E-02
a(19,23)=  3.4426E-02  1.8101E-02
a(20,23)= -2.7113E-02  2.5206E-02
a(21,23)= -8.0908E-03 -3.9720E-02
a(22,23)=  2.4115E-02 -2.3906E-02
a(23,23)=  8.3232E-01 -1.7889E-18
a(24,23)=  1.5152E-01  7.5009E-02
a(25,23)=  1.2402E-02 -8.1728E-04
a(26,23)= -1.1331E-01 -1.3227E-01
a( 1,24)=  1.7391E-03  1.0230E-02
a( 2,24)=  6.3773E-03  2.1783E-03
a( 3,24)=  3.0560E-03 -3.2533E-03
a( 4,24)=  5.6551E-03  1.2798E-03
a( 5,24)= -2.5013E-03 -8.0859E-04
a( 6,24)=  9.1191E-03 -5.7405E-04
a( 7,24)=  1.0564E-02  6.6717E-03
a( 8,24)= -3.8295E-03 -1.0804E-02
a( 9,24)=  1.7039E-02  1.3536E-02
a(10,24)=  3.8544E-02 -1.8050E-02
a(11,24)=  7.2107E-01 -4.9290E-05
a(12,24)= -1.1055E-01 -8.2524E-02
a(13,24)= -7.3233E-02 -3.5741E-02
a(14,24)=  3.6913E-03  3.4083E-02
a(15,24)=  3.5519E-02  1.9683E-02
a(16,24)=  1.3198E-02 -2.8620E-03
a(17,24)= -3.0196E-03  1.1503E-02
a(18,24)= -1.9898E-03  1.6409E-02
a(19,24)=  3.4166E-02 -9.2572E-03
a(20,24)=  3.5971E-02  1.6749E-02
a(21,24)= -1.7437E-02 -4.6827E-02
a(22,24)=  1.0594E-01  7.4607E-02
a(23,24)=  1.5152E-01 -7.5009E-02
a(24,24)=  1.1239E+00 -9.9963E-17
a(25,24)= -1.3526E-01 -8.0904E-02
a(26,24)= -1.0544E-01 -5.0978E-02
a( 1,25)=  2.0218E-03  7.8751E-03
a( 2,25)= -6.7830E-03  5.5796E-03
a( 3,25)=  2.8604E-03  2.8982E-03
a( 4,25)= -5.3904E-04 -2.9263E-03
a( 5,25)=  1.0779E-02 -4.6812E-03
a( 6,25)= -4.5575E-03  1.1417E-02
a( 7,25)=  3.6607E-03 -6.1027E-03
a( 8,25)=  1.1735E-02  5.5774E-03
a( 9,25)= -3.8767E-03  5.0341E-03
a(10,25)= -1.2353E-03  2.1755E-03
a(11,25)= -3.9966E-02  2.3196E-02
a(12,25)=  6.9050E-01  5.9057E-04
a(13,25)= -2.5329E-01  1.5086E-03
a(14,25)= -2.1019E-03  2.9005E-02
a(15,25)= -3.3795E-03 -1.2219E-03
a(16,25)=  1.3512E-02  2.6299E-02
a(17,25)=  3.8478E-02 -2.3090E-03
a(18,25)=  4.5580E-02 -4.5170E-03
a(19,25)= -1.7671E-02  6.1010E-02
a(20,25)=  1.9266E-02  9.0752E-03
a(21,25)=  1.6223E-02  1.7327E-02
a(22,25)= -1.1312E-02 -3.9615E-02
a(23,25)=  1.2402E-02  8.1728E-04
a(24,25)= -1.3526E-01  8.0904E-02
a(25,25)=  9.8294E-01 -2.7972E-17
a(26,25)= -2.8731E-01 -1.0378E-03
a( 1,26)=  3.6838E-03  2.3940E-02
a( 2,26)= -7.4490E-03  1.8829E-03
a( 3,26)=  7.3235E-04 -1.6659E-03
a( 4,26)= -6.0464E-03 -1.4578E-03
a( 5,26)= -3.6460E-03  1.8876E-04
a( 6,26)= -1.9354E-03 -1.6676E-03
a( 7,26)=  1.9822E-03 -1.6124E-02
a( 8,26)= -4.8205E-03  1.6864E-04
a( 9,26)=  3.5479E-04  5.5712E-03
a(10,26)= -1.2980E-02  1.4776E-02
a(11,26)= -2.4663E-02  1.2851E-02
a(12,26)= -2.0165E-02  2.8520E-03
a(13,26)=  9.6147E-01  3.3491E-04
a(14,26)=  6.1151E-03  6.6322E-02
a(15,26)= -6.0544E-02  6.7730E-03
a(16,26)=  5.0408E-03 -4.3310E-02
a(17,26)= -9.6095E-02  1.9265E-02
a(18,26)= -2.2874E-02 -6.0332E-02
a(19,26)=  1.5552E-02 -3.8042E-02
a(20,26)=  1.8178E-02 -9.6824E-02
a(21,26)= -1.7430E-02  4.3046E-03
a(22,26)=  1.7289E-02  9.1186E-02
a(23,26)= -1.1331E-01  1.3227E-01
a(24,26)= -1.0544E-01  5.0978E-02
a(25,26)= -2.8731E-01  1.0378E-03
a(26,26)=  1.4175E+00  4.2501E-16

Code: Select all
matrix b:
( 1, 1)=  1.0000E+00  0.0000E+00
b( 2, 1)=  0.0000E+00  0.0000E+00
b( 3, 1)=  0.0000E+00  0.0000E+00
b( 4, 1)=  0.0000E+00  0.0000E+00
b( 5, 1)=  0.0000E+00  0.0000E+00
b( 6, 1)=  0.0000E+00  0.0000E+00
b( 7, 1)=  0.0000E+00  0.0000E+00
b( 8, 1)=  0.0000E+00  0.0000E+00
b( 9, 1)=  0.0000E+00  0.0000E+00
b(10, 1)=  0.0000E+00  0.0000E+00
b(11, 1)=  0.0000E+00  0.0000E+00
b(12, 1)=  0.0000E+00  0.0000E+00
b(13, 1)=  0.0000E+00  0.0000E+00
b(14, 1)=  0.0000E+00  0.0000E+00
b(15, 1)=  0.0000E+00  0.0000E+00
b(16, 1)=  0.0000E+00  0.0000E+00
b(17, 1)=  0.0000E+00  0.0000E+00
b(18, 1)=  0.0000E+00  0.0000E+00
b(19, 1)=  0.0000E+00  0.0000E+00
b(20, 1)=  0.0000E+00  0.0000E+00
b(21, 1)=  0.0000E+00  0.0000E+00
b(22, 1)=  0.0000E+00  0.0000E+00
b(23, 1)=  0.0000E+00  0.0000E+00
b(24, 1)=  0.0000E+00  0.0000E+00
b(25, 1)=  0.0000E+00  0.0000E+00
b(26, 1)=  0.0000E+00  0.0000E+00
b( 1, 2)=  0.0000E+00  0.0000E+00
b( 2, 2)=  1.0000E+00  0.0000E+00
b( 3, 2)=  0.0000E+00  0.0000E+00
b( 4, 2)=  0.0000E+00  0.0000E+00
b( 5, 2)=  0.0000E+00  0.0000E+00
b( 6, 2)=  0.0000E+00  0.0000E+00
b( 7, 2)=  0.0000E+00  0.0000E+00
b( 8, 2)=  0.0000E+00  0.0000E+00
b( 9, 2)=  0.0000E+00  0.0000E+00
b(10, 2)=  0.0000E+00  0.0000E+00
b(11, 2)=  0.0000E+00  0.0000E+00
b(12, 2)=  0.0000E+00  0.0000E+00
b(13, 2)=  0.0000E+00  0.0000E+00
b(14, 2)=  0.0000E+00  0.0000E+00
b(15, 2)=  0.0000E+00  0.0000E+00
b(16, 2)=  0.0000E+00  0.0000E+00
b(17, 2)=  0.0000E+00  0.0000E+00
b(18, 2)=  0.0000E+00  0.0000E+00
b(19, 2)=  0.0000E+00  0.0000E+00
b(20, 2)=  0.0000E+00  0.0000E+00
b(21, 2)=  0.0000E+00  0.0000E+00
b(22, 2)=  0.0000E+00  0.0000E+00
b(23, 2)=  0.0000E+00  0.0000E+00
b(24, 2)=  0.0000E+00  0.0000E+00
b(25, 2)=  0.0000E+00  0.0000E+00
b(26, 2)=  0.0000E+00  0.0000E+00
b( 1, 3)=  0.0000E+00  0.0000E+00
b( 2, 3)=  0.0000E+00  0.0000E+00
b( 3, 3)=  1.0000E+00  0.0000E+00
b( 4, 3)=  0.0000E+00  0.0000E+00
b( 5, 3)=  0.0000E+00  0.0000E+00
b( 6, 3)=  0.0000E+00  0.0000E+00
b( 7, 3)=  0.0000E+00  0.0000E+00
b( 8, 3)=  0.0000E+00  0.0000E+00
b( 9, 3)=  0.0000E+00  0.0000E+00
b(10, 3)=  0.0000E+00  0.0000E+00
b(11, 3)=  0.0000E+00  0.0000E+00
b(12, 3)=  0.0000E+00  0.0000E+00
b(13, 3)=  0.0000E+00  0.0000E+00
b(14, 3)=  0.0000E+00  0.0000E+00
b(15, 3)=  0.0000E+00  0.0000E+00
b(16, 3)=  0.0000E+00  0.0000E+00
b(17, 3)=  0.0000E+00  0.0000E+00
b(18, 3)=  0.0000E+00  0.0000E+00
b(19, 3)=  0.0000E+00  0.0000E+00
b(20, 3)=  0.0000E+00  0.0000E+00
b(21, 3)=  0.0000E+00  0.0000E+00
b(22, 3)=  0.0000E+00  0.0000E+00
b(23, 3)=  0.0000E+00  0.0000E+00
b(24, 3)=  0.0000E+00  0.0000E+00
b(25, 3)=  0.0000E+00  0.0000E+00
b(26, 3)=  0.0000E+00  0.0000E+00
b( 1, 4)=  0.0000E+00  0.0000E+00
b( 2, 4)=  0.0000E+00  0.0000E+00
b( 3, 4)=  0.0000E+00  0.0000E+00
b( 4, 4)=  1.0000E+00  0.0000E+00
b( 5, 4)=  0.0000E+00  0.0000E+00
b( 6, 4)=  0.0000E+00  0.0000E+00
b( 7, 4)=  0.0000E+00  0.0000E+00
b( 8, 4)=  0.0000E+00  0.0000E+00
b( 9, 4)=  0.0000E+00  0.0000E+00
b(10, 4)=  0.0000E+00  0.0000E+00
b(11, 4)=  0.0000E+00  0.0000E+00
b(12, 4)=  0.0000E+00  0.0000E+00
b(13, 4)=  0.0000E+00  0.0000E+00
b(14, 4)=  0.0000E+00  0.0000E+00
b(15, 4)=  0.0000E+00  0.0000E+00
b(16, 4)=  0.0000E+00  0.0000E+00
b(17, 4)=  0.0000E+00  0.0000E+00
b(18, 4)=  0.0000E+00  0.0000E+00
b(19, 4)=  0.0000E+00  0.0000E+00
b(20, 4)=  0.0000E+00  0.0000E+00
b(21, 4)=  0.0000E+00  0.0000E+00
b(22, 4)=  0.0000E+00  0.0000E+00
b(23, 4)=  0.0000E+00  0.0000E+00
b(24, 4)=  0.0000E+00  0.0000E+00
b(25, 4)=  0.0000E+00  0.0000E+00
b(26, 4)=  0.0000E+00  0.0000E+00
b( 1, 5)=  0.0000E+00  0.0000E+00
b( 2, 5)=  0.0000E+00  0.0000E+00
b( 3, 5)=  0.0000E+00  0.0000E+00
b( 4, 5)=  0.0000E+00  0.0000E+00
b( 5, 5)=  1.0000E+00  0.0000E+00
b( 6, 5)=  0.0000E+00  0.0000E+00
b( 7, 5)=  0.0000E+00  0.0000E+00
b( 8, 5)=  0.0000E+00  0.0000E+00
b( 9, 5)=  0.0000E+00  0.0000E+00
b(10, 5)=  0.0000E+00  0.0000E+00
b(11, 5)=  0.0000E+00  0.0000E+00
b(12, 5)=  0.0000E+00  0.0000E+00
b(13, 5)=  0.0000E+00  0.0000E+00
b(14, 5)=  0.0000E+00  0.0000E+00
b(15, 5)=  0.0000E+00  0.0000E+00
b(16, 5)=  0.0000E+00  0.0000E+00
b(17, 5)=  0.0000E+00  0.0000E+00
b(18, 5)=  0.0000E+00  0.0000E+00
b(19, 5)=  0.0000E+00  0.0000E+00
b(20, 5)=  0.0000E+00  0.0000E+00
b(21, 5)=  0.0000E+00  0.0000E+00
b(22, 5)=  0.0000E+00  0.0000E+00
b(23, 5)=  0.0000E+00  0.0000E+00
b(24, 5)=  0.0000E+00  0.0000E+00
b(25, 5)=  0.0000E+00  0.0000E+00
b(26, 5)=  0.0000E+00  0.0000E+00
b( 1, 6)=  0.0000E+00  0.0000E+00
b( 2, 6)=  0.0000E+00  0.0000E+00
b( 3, 6)=  0.0000E+00  0.0000E+00
b( 4, 6)=  0.0000E+00  0.0000E+00
b( 5, 6)=  0.0000E+00  0.0000E+00
b( 6, 6)=  1.0000E+00  0.0000E+00
b( 7, 6)=  0.0000E+00  0.0000E+00
b( 8, 6)=  0.0000E+00  0.0000E+00
b( 9, 6)=  0.0000E+00  0.0000E+00
b(10, 6)=  0.0000E+00  0.0000E+00
b(11, 6)=  0.0000E+00  0.0000E+00
b(12, 6)=  0.0000E+00  0.0000E+00
b(13, 6)=  0.0000E+00  0.0000E+00
b(14, 6)=  0.0000E+00  0.0000E+00
b(15, 6)=  0.0000E+00  0.0000E+00
b(16, 6)=  0.0000E+00  0.0000E+00
b(17, 6)=  0.0000E+00  0.0000E+00
b(18, 6)=  0.0000E+00  0.0000E+00
b(19, 6)=  0.0000E+00  0.0000E+00
b(20, 6)=  0.0000E+00  0.0000E+00
b(21, 6)=  0.0000E+00  0.0000E+00
b(22, 6)=  0.0000E+00  0.0000E+00
b(23, 6)=  0.0000E+00  0.0000E+00
b(24, 6)=  0.0000E+00  0.0000E+00
b(25, 6)=  0.0000E+00  0.0000E+00
b(26, 6)=  0.0000E+00  0.0000E+00
b( 1, 7)=  0.0000E+00  0.0000E+00
b( 2, 7)=  0.0000E+00  0.0000E+00
b( 3, 7)=  0.0000E+00  0.0000E+00
b( 4, 7)=  0.0000E+00  0.0000E+00
b( 5, 7)=  0.0000E+00  0.0000E+00
b( 6, 7)=  0.0000E+00  0.0000E+00
b( 7, 7)=  1.0000E+00  0.0000E+00
b( 8, 7)=  0.0000E+00  0.0000E+00
b( 9, 7)=  0.0000E+00  0.0000E+00
b(10, 7)=  0.0000E+00  0.0000E+00
b(11, 7)=  0.0000E+00  0.0000E+00
b(12, 7)=  0.0000E+00  0.0000E+00
b(13, 7)=  0.0000E+00  0.0000E+00
b(14, 7)=  0.0000E+00  0.0000E+00
b(15, 7)=  0.0000E+00  0.0000E+00
b(16, 7)=  0.0000E+00  0.0000E+00
b(17, 7)=  0.0000E+00  0.0000E+00
b(18, 7)=  0.0000E+00  0.0000E+00
b(19, 7)=  0.0000E+00  0.0000E+00
b(20, 7)=  0.0000E+00  0.0000E+00
b(21, 7)=  0.0000E+00  0.0000E+00
b(22, 7)=  0.0000E+00  0.0000E+00
b(23, 7)=  0.0000E+00  0.0000E+00
b(24, 7)=  0.0000E+00  0.0000E+00
b(25, 7)=  0.0000E+00  0.0000E+00
b(26, 7)=  0.0000E+00  0.0000E+00
b( 1, 8)=  0.0000E+00  0.0000E+00
b( 2, 8)=  0.0000E+00  0.0000E+00
b( 3, 8)=  0.0000E+00  0.0000E+00
b( 4, 8)=  0.0000E+00  0.0000E+00
b( 5, 8)=  0.0000E+00  0.0000E+00
b( 6, 8)=  0.0000E+00  0.0000E+00
b( 7, 8)=  0.0000E+00  0.0000E+00
b( 8, 8)=  1.0000E+00  0.0000E+00
b( 9, 8)=  0.0000E+00  0.0000E+00
b(10, 8)=  0.0000E+00  0.0000E+00
b(11, 8)=  0.0000E+00  0.0000E+00
b(12, 8)=  0.0000E+00  0.0000E+00
b(13, 8)=  0.0000E+00  0.0000E+00
b(14, 8)=  0.0000E+00  0.0000E+00
b(15, 8)=  0.0000E+00  0.0000E+00
b(16, 8)=  0.0000E+00  0.0000E+00
b(17, 8)=  0.0000E+00  0.0000E+00
b(18, 8)=  0.0000E+00  0.0000E+00
b(19, 8)=  0.0000E+00  0.0000E+00
b(20, 8)=  0.0000E+00  0.0000E+00
b(21, 8)=  0.0000E+00  0.0000E+00
b(22, 8)=  0.0000E+00  0.0000E+00
b(23, 8)=  0.0000E+00  0.0000E+00
b(24, 8)=  0.0000E+00  0.0000E+00
b(25, 8)=  0.0000E+00  0.0000E+00
b(26, 8)=  0.0000E+00  0.0000E+00
b( 1, 9)=  0.0000E+00  0.0000E+00
b( 2, 9)=  0.0000E+00  0.0000E+00
b( 3, 9)=  0.0000E+00  0.0000E+00
b( 4, 9)=  0.0000E+00  0.0000E+00
b( 5, 9)=  0.0000E+00  0.0000E+00
b( 6, 9)=  0.0000E+00  0.0000E+00
b( 7, 9)=  0.0000E+00  0.0000E+00
b( 8, 9)=  0.0000E+00  0.0000E+00
b( 9, 9)=  1.0000E+00  0.0000E+00
b(10, 9)=  0.0000E+00  0.0000E+00
b(11, 9)=  0.0000E+00  0.0000E+00
b(12, 9)=  0.0000E+00  0.0000E+00
b(13, 9)=  0.0000E+00  0.0000E+00
b(14, 9)=  0.0000E+00  0.0000E+00
b(15, 9)=  0.0000E+00  0.0000E+00
b(16, 9)=  0.0000E+00  0.0000E+00
b(17, 9)=  0.0000E+00  0.0000E+00
b(18, 9)=  0.0000E+00  0.0000E+00
b(19, 9)=  0.0000E+00  0.0000E+00
b(20, 9)=  0.0000E+00  0.0000E+00
b(21, 9)=  0.0000E+00  0.0000E+00
b(22, 9)=  0.0000E+00  0.0000E+00
b(23, 9)=  0.0000E+00  0.0000E+00
b(24, 9)=  0.0000E+00  0.0000E+00
b(25, 9)=  0.0000E+00  0.0000E+00
b(26, 9)=  0.0000E+00  0.0000E+00
b( 1,10)=  0.0000E+00  0.0000E+00
b( 2,10)=  0.0000E+00  0.0000E+00
b( 3,10)=  0.0000E+00  0.0000E+00
b( 4,10)=  0.0000E+00  0.0000E+00
b( 5,10)=  0.0000E+00  0.0000E+00
b( 6,10)=  0.0000E+00  0.0000E+00
b( 7,10)=  0.0000E+00  0.0000E+00
b( 8,10)=  0.0000E+00  0.0000E+00
b( 9,10)=  0.0000E+00  0.0000E+00
b(10,10)=  1.0000E+00  0.0000E+00
b(11,10)=  0.0000E+00  0.0000E+00
b(12,10)=  0.0000E+00  0.0000E+00
b(13,10)=  0.0000E+00  0.0000E+00
b(14,10)=  0.0000E+00  0.0000E+00
b(15,10)=  0.0000E+00  0.0000E+00
b(16,10)=  0.0000E+00  0.0000E+00
b(17,10)=  0.0000E+00  0.0000E+00
b(18,10)=  0.0000E+00  0.0000E+00
b(19,10)=  0.0000E+00  0.0000E+00
b(20,10)=  0.0000E+00  0.0000E+00
b(21,10)=  0.0000E+00  0.0000E+00
b(22,10)=  0.0000E+00  0.0000E+00
b(23,10)=  0.0000E+00  0.0000E+00
b(24,10)=  0.0000E+00  0.0000E+00
b(25,10)=  0.0000E+00  0.0000E+00
b(26,10)=  0.0000E+00  0.0000E+00
b( 1,11)=  0.0000E+00  0.0000E+00
b( 2,11)=  0.0000E+00  0.0000E+00
b( 3,11)=  0.0000E+00  0.0000E+00
b( 4,11)=  0.0000E+00  0.0000E+00
b( 5,11)=  0.0000E+00  0.0000E+00
b( 6,11)=  0.0000E+00  0.0000E+00
b( 7,11)=  0.0000E+00  0.0000E+00
b( 8,11)=  0.0000E+00  0.0000E+00
b( 9,11)=  0.0000E+00  0.0000E+00
b(10,11)=  0.0000E+00  0.0000E+00
b(11,11)=  1.0000E+00  0.0000E+00
b(12,11)=  0.0000E+00  0.0000E+00
b(13,11)=  0.0000E+00  0.0000E+00
b(14,11)=  0.0000E+00  0.0000E+00
b(15,11)=  0.0000E+00  0.0000E+00
b(16,11)=  0.0000E+00  0.0000E+00
b(17,11)=  0.0000E+00  0.0000E+00
b(18,11)=  0.0000E+00  0.0000E+00
b(19,11)=  0.0000E+00  0.0000E+00
b(20,11)=  0.0000E+00  0.0000E+00
b(21,11)=  0.0000E+00  0.0000E+00
b(22,11)=  0.0000E+00  0.0000E+00
b(23,11)=  0.0000E+00  0.0000E+00
b(24,11)=  0.0000E+00  0.0000E+00
b(25,11)=  0.0000E+00  0.0000E+00
b(26,11)=  0.0000E+00  0.0000E+00
b( 1,12)=  0.0000E+00  0.0000E+00
b( 2,12)=  0.0000E+00  0.0000E+00
b( 3,12)=  0.0000E+00  0.0000E+00
b( 4,12)=  0.0000E+00  0.0000E+00
b( 5,12)=  0.0000E+00  0.0000E+00
b( 6,12)=  0.0000E+00  0.0000E+00
b( 7,12)=  0.0000E+00  0.0000E+00
b( 8,12)=  0.0000E+00  0.0000E+00
b( 9,12)=  0.0000E+00  0.0000E+00
b(10,12)=  0.0000E+00  0.0000E+00
b(11,12)=  0.0000E+00  0.0000E+00
b(12,12)=  1.0000E+00  0.0000E+00
b(13,12)=  0.0000E+00  0.0000E+00
b(14,12)=  0.0000E+00  0.0000E+00
b(15,12)=  0.0000E+00  0.0000E+00
b(16,12)=  0.0000E+00  0.0000E+00
b(17,12)=  0.0000E+00  0.0000E+00
b(18,12)=  0.0000E+00  0.0000E+00
b(19,12)=  0.0000E+00  0.0000E+00
b(20,12)=  0.0000E+00  0.0000E+00
b(21,12)=  0.0000E+00  0.0000E+00
b(22,12)=  0.0000E+00  0.0000E+00
b(23,12)=  0.0000E+00  0.0000E+00
b(24,12)=  0.0000E+00  0.0000E+00
b(25,12)=  0.0000E+00  0.0000E+00
b(26,12)=  0.0000E+00  0.0000E+00
b( 1,13)=  0.0000E+00  0.0000E+00
b( 2,13)=  0.0000E+00  0.0000E+00
b( 3,13)=  0.0000E+00  0.0000E+00
b( 4,13)=  0.0000E+00  0.0000E+00
b( 5,13)=  0.0000E+00  0.0000E+00
b( 6,13)=  0.0000E+00  0.0000E+00
b( 7,13)=  0.0000E+00  0.0000E+00
b( 8,13)=  0.0000E+00  0.0000E+00
b( 9,13)=  0.0000E+00  0.0000E+00
b(10,13)=  0.0000E+00  0.0000E+00
b(11,13)=  0.0000E+00  0.0000E+00
b(12,13)=  0.0000E+00  0.0000E+00
b(13,13)=  1.0000E+00  0.0000E+00
b(14,13)=  0.0000E+00  0.0000E+00
b(15,13)=  0.0000E+00  0.0000E+00
b(16,13)=  0.0000E+00  0.0000E+00
b(17,13)=  0.0000E+00  0.0000E+00
b(18,13)=  0.0000E+00  0.0000E+00
b(19,13)=  0.0000E+00  0.0000E+00
b(20,13)=  0.0000E+00  0.0000E+00
b(21,13)=  0.0000E+00  0.0000E+00
b(22,13)=  0.0000E+00  0.0000E+00
b(23,13)=  0.0000E+00  0.0000E+00
b(24,13)=  0.0000E+00  0.0000E+00
b(25,13)=  0.0000E+00  0.0000E+00
b(26,13)=  0.0000E+00  0.0000E+00
b( 1,14)=  0.0000E+00  0.0000E+00
b( 2,14)=  0.0000E+00  0.0000E+00
b( 3,14)=  0.0000E+00  0.0000E+00
b( 4,14)=  0.0000E+00  0.0000E+00
b( 5,14)=  0.0000E+00  0.0000E+00
b( 6,14)=  0.0000E+00  0.0000E+00
b( 7,14)=  0.0000E+00  0.0000E+00
b( 8,14)=  0.0000E+00  0.0000E+00
b( 9,14)=  0.0000E+00  0.0000E+00
b(10,14)=  0.0000E+00  0.0000E+00
b(11,14)=  0.0000E+00  0.0000E+00
b(12,14)=  0.0000E+00  0.0000E+00
b(13,14)=  0.0000E+00  0.0000E+00
b(14,14)=  1.0000E+00  0.0000E+00
b(15,14)=  0.0000E+00  0.0000E+00
b(16,14)=  0.0000E+00  0.0000E+00
b(17,14)=  0.0000E+00  0.0000E+00
b(18,14)=  0.0000E+00  0.0000E+00
b(19,14)=  0.0000E+00  0.0000E+00
b(20,14)=  0.0000E+00  0.0000E+00
b(21,14)=  0.0000E+00  0.0000E+00
b(22,14)=  0.0000E+00  0.0000E+00
b(23,14)=  0.0000E+00  0.0000E+00
b(24,14)=  0.0000E+00  0.0000E+00
b(25,14)=  0.0000E+00  0.0000E+00
b(26,14)=  0.0000E+00  0.0000E+00
b( 1,15)=  0.0000E+00  0.0000E+00
b( 2,15)=  0.0000E+00  0.0000E+00
b( 3,15)=  0.0000E+00  0.0000E+00
b( 4,15)=  0.0000E+00  0.0000E+00
b( 5,15)=  0.0000E+00  0.0000E+00
b( 6,15)=  0.0000E+00  0.0000E+00
b( 7,15)=  0.0000E+00  0.0000E+00
b( 8,15)=  0.0000E+00  0.0000E+00
b( 9,15)=  0.0000E+00  0.0000E+00
b(10,15)=  0.0000E+00  0.0000E+00
b(11,15)=  0.0000E+00  0.0000E+00
b(12,15)=  0.0000E+00  0.0000E+00
b(13,15)=  0.0000E+00  0.0000E+00
b(14,15)=  0.0000E+00  0.0000E+00
b(15,15)=  1.0000E+00  0.0000E+00
b(16,15)=  0.0000E+00  0.0000E+00
b(17,15)=  0.0000E+00  0.0000E+00
b(18,15)=  0.0000E+00  0.0000E+00
b(19,15)=  0.0000E+00  0.0000E+00
b(20,15)=  0.0000E+00  0.0000E+00
b(21,15)=  0.0000E+00  0.0000E+00
b(22,15)=  0.0000E+00  0.0000E+00
b(23,15)=  0.0000E+00  0.0000E+00
b(24,15)=  0.0000E+00  0.0000E+00
b(25,15)=  0.0000E+00  0.0000E+00
b(26,15)=  0.0000E+00  0.0000E+00
b( 1,16)=  0.0000E+00  0.0000E+00
b( 2,16)=  0.0000E+00  0.0000E+00
b( 3,16)=  0.0000E+00  0.0000E+00
b( 4,16)=  0.0000E+00  0.0000E+00
b( 5,16)=  0.0000E+00  0.0000E+00
b( 6,16)=  0.0000E+00  0.0000E+00
b( 7,16)=  0.0000E+00  0.0000E+00
b( 8,16)=  0.0000E+00  0.0000E+00
b( 9,16)=  0.0000E+00  0.0000E+00
b(10,16)=  0.0000E+00  0.0000E+00
b(11,16)=  0.0000E+00  0.0000E+00
b(12,16)=  0.0000E+00  0.0000E+00
b(13,16)=  0.0000E+00  0.0000E+00
b(14,16)=  0.0000E+00  0.0000E+00
b(15,16)=  0.0000E+00  0.0000E+00
b(16,16)=  1.0000E+00  0.0000E+00
b(17,16)=  0.0000E+00  0.0000E+00
b(18,16)=  0.0000E+00  0.0000E+00
b(19,16)=  0.0000E+00  0.0000E+00
b(20,16)=  0.0000E+00  0.0000E+00
b(21,16)=  0.0000E+00  0.0000E+00
b(22,16)=  0.0000E+00  0.0000E+00
b(23,16)=  0.0000E+00  0.0000E+00
b(24,16)=  0.0000E+00  0.0000E+00
b(25,16)=  0.0000E+00  0.0000E+00
b(26,16)=  0.0000E+00  0.0000E+00
b( 1,17)=  0.0000E+00  0.0000E+00
b( 2,17)=  0.0000E+00  0.0000E+00
b( 3,17)=  0.0000E+00  0.0000E+00
b( 4,17)=  0.0000E+00  0.0000E+00
b( 5,17)=  0.0000E+00  0.0000E+00
b( 6,17)=  0.0000E+00  0.0000E+00
b( 7,17)=  0.0000E+00  0.0000E+00
b( 8,17)=  0.0000E+00  0.0000E+00
b( 9,17)=  0.0000E+00  0.0000E+00
b(10,17)=  0.0000E+00  0.0000E+00
b(11,17)=  0.0000E+00  0.0000E+00
b(12,17)=  0.0000E+00  0.0000E+00
b(13,17)=  0.0000E+00  0.0000E+00
b(14,17)=  0.0000E+00  0.0000E+00
b(15,17)=  0.0000E+00  0.0000E+00
b(16,17)=  0.0000E+00  0.0000E+00
b(17,17)=  1.0000E+00  0.0000E+00
b(18,17)=  0.0000E+00  0.0000E+00
b(19,17)=  0.0000E+00  0.0000E+00
b(20,17)=  0.0000E+00  0.0000E+00
b(21,17)=  0.0000E+00  0.0000E+00
b(22,17)=  0.0000E+00  0.0000E+00
b(23,17)=  0.0000E+00  0.0000E+00
b(24,17)=  0.0000E+00  0.0000E+00
b(25,17)=  0.0000E+00  0.0000E+00
b(26,17)=  0.0000E+00  0.0000E+00
b( 1,18)=  0.0000E+00  0.0000E+00
b( 2,18)=  0.0000E+00  0.0000E+00
b( 3,18)=  0.0000E+00  0.0000E+00
b( 4,18)=  0.0000E+00  0.0000E+00
b( 5,18)=  0.0000E+00  0.0000E+00
b( 6,18)=  0.0000E+00  0.0000E+00
b( 7,18)=  0.0000E+00  0.0000E+00
b( 8,18)=  0.0000E+00  0.0000E+00
b( 9,18)=  0.0000E+00  0.0000E+00
b(10,18)=  0.0000E+00  0.0000E+00
b(11,18)=  0.0000E+00  0.0000E+00
b(12,18)=  0.0000E+00  0.0000E+00
b(13,18)=  0.0000E+00  0.0000E+00
b(14,18)=  0.0000E+00  0.0000E+00
b(15,18)=  0.0000E+00  0.0000E+00
b(16,18)=  0.0000E+00  0.0000E+00
b(17,18)=  0.0000E+00  0.0000E+00
b(18,18)=  1.0000E+00  0.0000E+00
b(19,18)=  0.0000E+00  0.0000E+00
b(20,18)=  0.0000E+00  0.0000E+00
b(21,18)=  0.0000E+00  0.0000E+00
b(22,18)=  0.0000E+00  0.0000E+00
b(23,18)=  0.0000E+00  0.0000E+00
b(24,18)=  0.0000E+00  0.0000E+00
b(25,18)=  0.0000E+00  0.0000E+00
b(26,18)=  0.0000E+00  0.0000E+00
b( 1,19)=  0.0000E+00  0.0000E+00
b( 2,19)=  0.0000E+00  0.0000E+00
b( 3,19)=  0.0000E+00  0.0000E+00
b( 4,19)=  0.0000E+00  0.0000E+00
b( 5,19)=  0.0000E+00  0.0000E+00
b( 6,19)=  0.0000E+00  0.0000E+00
b( 7,19)=  0.0000E+00  0.0000E+00
b( 8,19)=  0.0000E+00  0.0000E+00
b( 9,19)=  0.0000E+00  0.0000E+00
b(10,19)=  0.0000E+00  0.0000E+00
b(11,19)=  0.0000E+00  0.0000E+00
b(12,19)=  0.0000E+00  0.0000E+00
b(13,19)=  0.0000E+00  0.0000E+00
b(14,19)=  0.0000E+00  0.0000E+00
b(15,19)=  0.0000E+00  0.0000E+00
b(16,19)=  0.0000E+00  0.0000E+00
b(17,19)=  0.0000E+00  0.0000E+00
b(18,19)=  0.0000E+00  0.0000E+00
b(19,19)=  1.0000E+00  0.0000E+00
b(20,19)=  0.0000E+00  0.0000E+00
b(21,19)=  0.0000E+00  0.0000E+00
b(22,19)=  0.0000E+00  0.0000E+00
b(23,19)=  0.0000E+00  0.0000E+00
b(24,19)=  0.0000E+00  0.0000E+00
b(25,19)=  0.0000E+00  0.0000E+00
b(26,19)=  0.0000E+00  0.0000E+00
b( 1,20)=  0.0000E+00  0.0000E+00
b( 2,20)=  0.0000E+00  0.0000E+00
b( 3,20)=  0.0000E+00  0.0000E+00
b( 4,20)=  0.0000E+00  0.0000E+00
b( 5,20)=  0.0000E+00  0.0000E+00
b( 6,20)=  0.0000E+00  0.0000E+00
b( 7,20)=  0.0000E+00  0.0000E+00
b( 8,20)=  0.0000E+00  0.0000E+00
b( 9,20)=  0.0000E+00  0.0000E+00
b(10,20)=  0.0000E+00  0.0000E+00
b(11,20)=  0.0000E+00  0.0000E+00
b(12,20)=  0.0000E+00  0.0000E+00
b(13,20)=  0.0000E+00  0.0000E+00
b(14,20)=  0.0000E+00  0.0000E+00
b(15,20)=  0.0000E+00  0.0000E+00
b(16,20)=  0.0000E+00  0.0000E+00
b(17,20)=  0.0000E+00  0.0000E+00
b(18,20)=  0.0000E+00  0.0000E+00
b(19,20)=  0.0000E+00  0.0000E+00
b(20,20)=  1.0000E+00  0.0000E+00
b(21,20)=  0.0000E+00  0.0000E+00
b(22,20)=  0.0000E+00  0.0000E+00
b(23,20)=  0.0000E+00  0.0000E+00
b(24,20)=  0.0000E+00  0.0000E+00
b(25,20)=  0.0000E+00  0.0000E+00
b(26,20)=  0.0000E+00  0.0000E+00
b( 1,21)=  0.0000E+00  0.0000E+00
b( 2,21)=  0.0000E+00  0.0000E+00
b( 3,21)=  0.0000E+00  0.0000E+00
b( 4,21)=  0.0000E+00  0.0000E+00
b( 5,21)=  0.0000E+00  0.0000E+00
b( 6,21)=  0.0000E+00  0.0000E+00
b( 7,21)=  0.0000E+00  0.0000E+00
b( 8,21)=  0.0000E+00  0.0000E+00
b( 9,21)=  0.0000E+00  0.0000E+00
b(10,21)=  0.0000E+00  0.0000E+00
b(11,21)=  0.0000E+00  0.0000E+00
b(12,21)=  0.0000E+00  0.0000E+00
b(13,21)=  0.0000E+00  0.0000E+00
b(14,21)=  0.0000E+00  0.0000E+00
b(15,21)=  0.0000E+00  0.0000E+00
b(16,21)=  0.0000E+00  0.0000E+00
b(17,21)=  0.0000E+00  0.0000E+00
b(18,21)=  0.0000E+00  0.0000E+00
b(19,21)=  0.0000E+00  0.0000E+00
b(20,21)=  0.0000E+00  0.0000E+00
b(21,21)=  1.0000E+00  0.0000E+00
b(22,21)=  0.0000E+00  0.0000E+00
b(23,21)=  0.0000E+00  0.0000E+00
b(24,21)=  0.0000E+00  0.0000E+00
b(25,21)=  0.0000E+00  0.0000E+00
b(26,21)=  0.0000E+00  0.0000E+00
b( 1,22)=  0.0000E+00  0.0000E+00
b( 2,22)=  0.0000E+00  0.0000E+00
b( 3,22)=  0.0000E+00  0.0000E+00
b( 4,22)=  0.0000E+00  0.0000E+00
b( 5,22)=  0.0000E+00  0.0000E+00
b( 6,22)=  0.0000E+00  0.0000E+00
b( 7,22)=  0.0000E+00  0.0000E+00
b( 8,22)=  0.0000E+00  0.0000E+00
b( 9,22)=  0.0000E+00  0.0000E+00
b(10,22)=  0.0000E+00  0.0000E+00
b(11,22)=  0.0000E+00  0.0000E+00
b(12,22)=  0.0000E+00  0.0000E+00
b(13,22)=  0.0000E+00  0.0000E+00
b(14,22)=  0.0000E+00  0.0000E+00
b(15,22)=  0.0000E+00  0.0000E+00
b(16,22)=  0.0000E+00  0.0000E+00
b(17,22)=  0.0000E+00  0.0000E+00
b(18,22)=  0.0000E+00  0.0000E+00
b(19,22)=  0.0000E+00  0.0000E+00
b(20,22)=  0.0000E+00  0.0000E+00
b(21,22)=  0.0000E+00  0.0000E+00
b(22,22)=  1.0000E+00  0.0000E+00
b(23,22)=  0.0000E+00  0.0000E+00
b(24,22)=  0.0000E+00  0.0000E+00
b(25,22)=  0.0000E+00  0.0000E+00
b(26,22)=  0.0000E+00  0.0000E+00
b( 1,23)=  0.0000E+00  0.0000E+00
b( 2,23)=  0.0000E+00  0.0000E+00
b( 3,23)=  0.0000E+00  0.0000E+00
b( 4,23)=  0.0000E+00  0.0000E+00
b( 5,23)=  0.0000E+00  0.0000E+00
b( 6,23)=  0.0000E+00  0.0000E+00
b( 7,23)=  0.0000E+00  0.0000E+00
b( 8,23)=  0.0000E+00  0.0000E+00
b( 9,23)=  0.0000E+00  0.0000E+00
b(10,23)=  0.0000E+00  0.0000E+00
b(11,23)=  0.0000E+00  0.0000E+00
b(12,23)=  0.0000E+00  0.0000E+00
b(13,23)=  0.0000E+00  0.0000E+00
b(14,23)=  0.0000E+00  0.0000E+00
b(15,23)=  0.0000E+00  0.0000E+00
b(16,23)=  0.0000E+00  0.0000E+00
b(17,23)=  0.0000E+00  0.0000E+00
b(18,23)=  0.0000E+00  0.0000E+00
b(19,23)=  0.0000E+00  0.0000E+00
b(20,23)=  0.0000E+00  0.0000E+00
b(21,23)=  0.0000E+00  0.0000E+00
b(22,23)=  0.0000E+00  0.0000E+00
b(23,23)=  1.0000E+00  0.0000E+00
b(24,23)=  0.0000E+00  0.0000E+00
b(25,23)=  0.0000E+00  0.0000E+00
b(26,23)=  0.0000E+00  0.0000E+00
b( 1,24)=  0.0000E+00  0.0000E+00
b( 2,24)=  0.0000E+00  0.0000E+00
b( 3,24)=  0.0000E+00  0.0000E+00
b( 4,24)=  0.0000E+00  0.0000E+00
b( 5,24)=  0.0000E+00  0.0000E+00
b( 6,24)=  0.0000E+00  0.0000E+00
b( 7,24)=  0.0000E+00  0.0000E+00
b( 8,24)=  0.0000E+00  0.0000E+00
b( 9,24)=  0.0000E+00  0.0000E+00
b(10,24)=  0.0000E+00  0.0000E+00
b(11,24)=  0.0000E+00  0.0000E+00
b(12,24)=  0.0000E+00  0.0000E+00
b(13,24)=  0.0000E+00  0.0000E+00
b(14,24)=  0.0000E+00  0.0000E+00
b(15,24)=  0.0000E+00  0.0000E+00
b(16,24)=  0.0000E+00  0.0000E+00
b(17,24)=  0.0000E+00  0.0000E+00
b(18,24)=  0.0000E+00  0.0000E+00
b(19,24)=  0.0000E+00  0.0000E+00
b(20,24)=  0.0000E+00  0.0000E+00
b(21,24)=  0.0000E+00  0.0000E+00
b(22,24)=  0.0000E+00  0.0000E+00
b(23,24)=  0.0000E+00  0.0000E+00
b(24,24)=  1.0000E+00  0.0000E+00
b(25,24)=  0.0000E+00  0.0000E+00
b(26,24)=  0.0000E+00  0.0000E+00
b( 1,25)=  0.0000E+00  0.0000E+00
b( 2,25)=  0.0000E+00  0.0000E+00
b( 3,25)=  0.0000E+00  0.0000E+00
b( 4,25)=  0.0000E+00  0.0000E+00
b( 5,25)=  0.0000E+00  0.0000E+00
b( 6,25)=  0.0000E+00  0.0000E+00
b( 7,25)=  0.0000E+00  0.0000E+00
b( 8,25)=  0.0000E+00  0.0000E+00
b( 9,25)=  0.0000E+00  0.0000E+00
b(10,25)=  0.0000E+00  0.0000E+00
b(11,25)=  0.0000E+00  0.0000E+00
b(12,25)=  0.0000E+00  0.0000E+00
b(13,25)=  0.0000E+00  0.0000E+00
b(14,25)=  0.0000E+00  0.0000E+00
b(15,25)=  0.0000E+00  0.0000E+00
b(16,25)=  0.0000E+00  0.0000E+00
b(17,25)=  0.0000E+00  0.0000E+00
b(18,25)=  0.0000E+00  0.0000E+00
b(19,25)=  0.0000E+00  0.0000E+00
b(20,25)=  0.0000E+00  0.0000E+00
b(21,25)=  0.0000E+00  0.0000E+00
b(22,25)=  0.0000E+00  0.0000E+00
b(23,25)=  0.0000E+00  0.0000E+00
b(24,25)=  0.0000E+00  0.0000E+00
b(25,25)=  1.0000E+00  0.0000E+00
b(26,25)=  0.0000E+00  0.0000E+00
b( 1,26)=  0.0000E+00  0.0000E+00
b( 2,26)=  0.0000E+00  0.0000E+00
b( 3,26)=  0.0000E+00  0.0000E+00
b( 4,26)=  0.0000E+00  0.0000E+00
b( 5,26)=  0.0000E+00  0.0000E+00
b( 6,26)=  0.0000E+00  0.0000E+00
b( 7,26)=  0.0000E+00  0.0000E+00
b( 8,26)=  0.0000E+00  0.0000E+00
b( 9,26)=  0.0000E+00  0.0000E+00
b(10,26)=  0.0000E+00  0.0000E+00
b(11,26)=  0.0000E+00  0.0000E+00
b(12,26)=  0.0000E+00  0.0000E+00
b(13,26)=  0.0000E+00  0.0000E+00
b(14,26)=  0.0000E+00  0.0000E+00
b(15,26)=  0.0000E+00  0.0000E+00
b(16,26)=  0.0000E+00  0.0000E+00
b(17,26)=  0.0000E+00  0.0000E+00
b(18,26)=  0.0000E+00  0.0000E+00
b(19,26)=  0.0000E+00  0.0000E+00
b(20,26)=  0.0000E+00  0.0000E+00
b(21,26)=  0.0000E+00  0.0000E+00
b(22,26)=  0.0000E+00  0.0000E+00
b(23,26)=  0.0000E+00  0.0000E+00
b(24,26)=  0.0000E+00  0.0000E+00
b(25,26)=  0.0000E+00  0.0000E+00
b(26,26)=  1.0000E+00  0.0000E+00

Code: Select all
results with mkl zhegv:
eigen( 1)=  8.6053E-02
eigen( 2)=  1.0131E-01
eigen( 3)=  1.3076E-01
eigen( 4)=  1.4243E-01
eigen( 5)=  1.4720E-01
eigen( 6)=  1.6026E-01
eigen( 7)=  1.6603E-01
eigen( 8)=  1.8213E-01
eigen( 9)=  1.8959E-01
eigen(10)=  2.0850E-01
eigen(11)=  2.6334E-01
eigen(12)=  2.6669E-01
eigen(13)=  2.8979E-01
eigen(14)=  1.1081E+00
eigen(15)=  1.1444E+00
eigen(16)=  1.1819E+00
eigen(17)=  1.3607E+00
eigen(18)=  1.3775E+00
eigen(19)=  1.4243E+00
eigen(20)=  1.4811E+00
eigen(21)=  1.5196E+00
eigen(22)=  1.5665E+00
eigen(23)=  1.6077E+00
eigen(24)=  1.6578E+00
eigen(25)=  1.9861E+00
eigen(26)=  2.4869E+00

Code: Select all
results with magma_zhegvd:
eigen( 1)= -1.4830E+03
eigen( 2)=  0.0000E+00
eigen( 3)=  0.0000E+00
eigen( 4)=  0.0000E+00
eigen( 5)=  0.0000E+00
eigen( 6)=  0.0000E+00
eigen( 7)=  0.0000E+00
eigen( 8)=  0.0000E+00
eigen( 9)=  0.0000E+00
eigen(10)=  0.0000E+00
eigen(11)=  0.0000E+00
eigen(12)=  0.0000E+00
eigen(13)=  0.0000E+00
eigen(14)=  0.0000E+00
eigen(15)=  0.0000E+00
eigen(16)=  0.0000E+00
eigen(17)=  0.0000E+00
eigen(18)=  0.0000E+00
eigen(19)=  0.0000E+00
eigen(20)=  0.0000E+00
eigen(21)=  0.0000E+00
eigen(22)=  0.0000E+00
eigen(23)=  0.0000E+00
eigen(24)=  0.0000E+00
eigen(25)=  0.0000E+00
eigen(26)=  1.4830E+03

Also, the fortan binding for the generalized eigen problem are missing, i join the one i've made for double and double complex version
magma_df77.cpp
(22.52 KiB) Downloaded 51 times
Attachments
magma_dfortran_f90.txt
the .f90 extension was refused by the forum
(17.7 KiB) Downloaded 48 times
magma_zf77.cpp
(23.4 KiB) Downloaded 51 times
fdahm
 
Posts: 8
Joined: Tue Feb 15, 2011 10:04 am

Re: Error in Zhegvd

Postby fdahm » Fri Sep 02, 2011 4:53 am

last file (3 attachments was the max)
Attachments
magma_zfortran_f90.txt
the .f90 extension was refused by the forum
(17.26 KiB) Downloaded 50 times
fdahm
 
Posts: 8
Joined: Tue Feb 15, 2011 10:04 am


Return to User discussion

Who is online

Users browsing this forum: Bing [Bot], Google [Bot], Yahoo [Bot] and 1 guest

cron