test_535_eigen.c running, check 535_double_f.out 535_2.dat 4 4 A real, A imag, B real, B ima A matrix ( 1, 1)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 1, 2)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 1, 3)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 1, 4)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 2, 1)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 2, 2)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 2, 3)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 2, 4)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 3, 1)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 3, 2)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 3, 3)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 3, 4)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 4, 1)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 4, 2)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 4, 3)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 4, 4)=( 2.00000000000000E+00,-2.00000000000000E+00) 4 B matrix ( 1, 1)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 1, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 2, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 3, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=( 1.00000000000000E+00, 0.00000000000000E+00) cxhes results A matrix ( 1, 1)=( 2.00000000000000E+00,-2.00000000000000E+00) ( 1, 2)=( 3.46410161513775E+00,-3.46410161513775E+00) ( 1, 3)=(-4.44089209850063E-16, 4.44089209850063E-16) ( 1, 4)=(-1.35525271560688E-17, 1.35525271560688E-17) ( 2, 1)=( 4.89897948556636E+00, 0.00000000000000E+00) ( 2, 2)=( 8.48528137423857E+00, 0.00000000000000E+00) ( 2, 3)=(-1.08779196440841E-15, 0.00000000000000E+00) ( 2, 4)=(-1.82579645846559E-16, 0.00000000000000E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 5.43895982204207E-16, 0.00000000000000E+00) ( 3, 3)=(-6.97261119368420E-32, 3.39198063679146E-66) ( 3, 4)=(-2.12787206838508E-33,-3.74187059412907E-50) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=(-0.00000000000000E+00, 0.00000000000000E+00) cxhes results B matrix ( 1, 1)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 1, 2)=(-0.00000000000000E+00, 0.00000000000000E+00) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 7.07106781186548E-01, 7.07106781186548E-01) ( 2, 3)=(-9.06493303673679E-17,-9.06493303673679E-17) ( 2, 4)=( 2.12774676350280E-17, 2.12774676350280E-17) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=(-7.07106781186547E-01,-7.07106781186547E-01) ( 3, 4)=(-6.40987562127854E-17,-6.40987562127854E-17) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=( 1.00000000000000E+00, 0.00000000000000E+00) cxhes results Z matrix ( 1, 1)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 1, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 5.77350269189626E-01, 0.00000000000000E+00) ( 2, 3)=(-7.40148683083438E-17, 0.00000000000000E+00) ( 2, 4)=( 8.16496580927726E-01, 0.00000000000000E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 5.77350269189626E-01, 0.00000000000000E+00) ( 3, 3)=(-7.07106781186548E-01, 0.00000000000000E+00) ( 3, 4)=(-4.08248290463863E-01, 0.00000000000000E+00) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 5.77350269189626E-01, 0.00000000000000E+00) ( 4, 3)=( 7.07106781186547E-01, 0.00000000000000E+00) ( 4, 4)=(-4.08248290463863E-01, 0.00000000000000E+00) cxval results A matrix ( 1, 1)=( 2.71947991102104E-16,-2.71947991102104E-16) ( 1, 2)=( 1.57009245868377E-16,-1.57009245868377E-16) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=(-5.28148200397360E-17, 5.28148200397360E-17) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 8.00000000000000E+00,-8.00000000000000E+00) ( 2, 3)=( 8.88178419700125E-16,-8.88178419700125E-16) ( 2, 4)=( 1.18583006013567E-16,-1.18583006013567E-16) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=( 4.93038065763132E-32,-4.93038065763132E-32) ( 3, 4)=( 1.50463276905253E-33,-1.50463276905253E-33) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) cxval results B matrix ( 1, 1)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 1, 2)=(-7.85046229341888E-17, 0.00000000000000E+00) ( 1, 3)=(-6.40987562127855E-17, 0.00000000000000E+00) ( 1, 4)=( 1.50454416512056E-17, 0.00000000000000E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 2, 3)=( 1.11022302462516E-16, 0.00000000000000E+00) ( 2, 4)=(-2.60594693622011E-17, 0.00000000000000E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 3, 4)=( 9.06493303673679E-17, 0.00000000000000E+00) ( 4, 1)=( 3.14018491736755E-16, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=( 1.00000000000000E+00, 0.00000000000000E+00) cxval results Z matrix ( 1, 1)=(-8.66025403784439E-01, 0.00000000000000E+00) ( 1, 2)=(-5.00000000000000E-01,-0.00000000000000E+00) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 1)=( 2.88675134594813E-01, 0.00000000000000E+00) ( 2, 2)=(-5.00000000000000E-01,-0.00000000000000E+00) ( 2, 3)=(-7.40148683083438E-17, 0.00000000000000E+00) ( 2, 4)=( 8.16496580927726E-01, 0.00000000000000E+00) ( 3, 1)=( 2.88675134594813E-01, 0.00000000000000E+00) ( 3, 2)=(-5.00000000000000E-01,-0.00000000000000E+00) ( 3, 3)=(-7.07106781186548E-01, 0.00000000000000E+00) ( 3, 4)=(-4.08248290463863E-01, 0.00000000000000E+00) ( 4, 1)=( 2.88675134594813E-01, 0.00000000000000E+00) ( 4, 2)=(-5.00000000000000E-01,-0.00000000000000E+00) ( 4, 3)=( 7.07106781186547E-01, 0.00000000000000E+00) ( 4, 4)=(-4.08248290463863E-01, 0.00000000000000E+00) cxval error= 0 ALFR(I) ALFI(I) BETA(I) 1 2.71947991102104E-16 -2.71947991102104E-16 1.00000000000000E+00 2 8.00000000000000E+00 -8.00000000000000E+00 1.00000000000000E+00 3 4.93038065763132E-32 -4.93038065763132E-32 1.00000000000000E+00 4 0.00000000000000E+00 0.00000000000000E+00 1.00000000000000E+00 cxvec final eigenvctors ( 1, 1)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 1, 2)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 1, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 1, 4)=(-1.92754842380641E-01, 0.00000000000000E+00) ( 2, 1)=( 3.33333333333333E-01, 0.00000000000000E+00) ( 2, 2)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 2, 3)=(-8.71576399210524E-33, 0.00000000000000E+00) ( 2, 4)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 3, 1)=( 3.33333333333333E-01, 0.00000000000000E+00) ( 3, 2)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 3, 4)=(-3.78891686667455E-01, 0.00000000000000E+00) ( 4, 1)=( 3.33333333333333E-01, 0.00000000000000E+00) ( 4, 2)=(-1.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 1.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=(-4.28353470951904E-01, 0.00000000000000E+00) cxcheck resdul= 1.17756934401283E-16 eigdet check eigenvalues 2 = SINGULARITY 1 = SINGULARITY 3 = SINGULARITY 3 = SINGULARITY evalvec check eigenvectors SUM( 1)= 3.55271367880050E-15 SUM( 2)= 1.42108547152020E-14 SUM( 3)= 5.32907051820075E-15 SUM( 4)= 0.00000000000000E+00 TOTAL EIGENVECTOR ERROR = 2.30926389122033E-14 test_535_eigen.adb finished