test_535_eigen.c running, check 535_double_f.out 535b.dat 4 0 A real, A imag, B real, B imag A matrix ( 1, 1)=( 1.00000000000000E+00, 4.00000000000000E+00) ( 1, 2)=( 2.00000000000000E+00, 4.00000000000000E+00) ( 1, 3)=( 3.00000000000000E+00, 4.00000000000000E+00) ( 1, 4)=( 4.00000000000000E+00, 4.00000000000000E+00) ( 2, 1)=( 2.00000000000000E+00, 4.00000000000000E+00) ( 2, 2)=( 2.00000000000000E+00, 3.00000000000000E+00) ( 2, 3)=( 3.00000000000000E+00, 3.00000000000000E+00) ( 2, 4)=( 4.00000000000000E+00, 3.00000000000000E+00) ( 3, 1)=( 3.00000000000000E+00, 4.00000000000000E+00) ( 3, 2)=( 3.00000000000000E+00, 3.00000000000000E+00) ( 3, 3)=( 3.00000000000000E+00, 2.00000000000000E+00) ( 3, 4)=( 4.00000000000000E+00, 2.00000000000000E+00) ( 4, 1)=( 4.00000000000000E+00, 4.00000000000000E+00) ( 4, 2)=( 4.00000000000000E+00, 3.00000000000000E+00) ( 4, 3)=( 4.00000000000000E+00, 2.00000000000000E+00) ( 4, 4)=( 4.00000000000000E+00, 1.00000000000000E+00) 4 0 B matrix ( 1, 1)=( 3.00000000000000E+00, 5.00000000000000E+00) ( 1, 2)=( 2.00000000000000E+00, 4.00000000000000E+00) ( 1, 3)=( 3.00000000000000E+00, 4.00000000000000E+00) ( 1, 4)=( 2.00000000000000E+00, 3.00000000000000E+00) ( 2, 1)=( 4.00000000000000E+00, 4.00000000000000E+00) ( 2, 2)=( 4.00000000000000E+00, 5.00000000000000E+00) ( 2, 3)=( 2.00000000000000E+00, 3.00000000000000E+00) ( 2, 4)=( 5.00000000000000E+00, 2.00000000000000E+00) ( 3, 1)=( 5.00000000000000E+00, 3.00000000000000E+00) ( 3, 2)=( 1.00000000000000E+00, 2.00000000000000E+00) ( 3, 3)=( 3.00000000000000E+00, 5.00000000000000E+00) ( 3, 4)=( 3.00000000000000E+00, 4.00000000000000E+00) ( 4, 1)=( 2.00000000000000E+00, 2.00000000000000E+00) ( 4, 2)=( 5.00000000000000E+00, 5.00000000000000E+00) ( 4, 3)=( 5.00000000000000E+00, 4.00000000000000E+00) ( 4, 4)=( 3.00000000000000E+00, 5.00000000000000E+00) cxhes results A matrix ( 1, 1)=( 8.66025403784438E+00, 2.50185116648838E+00) ( 1, 2)=( 1.02048517075613E+01, 3.06484766924932E+00) ( 1, 3)=(-5.53442031808211E+00,-3.87111467981838E-01) ( 1, 4)=( 5.72835331690388E+00, 6.85456973361309E+00) ( 2, 1)=( 3.56941742315756E+00, 0.00000000000000E+00) ( 2, 2)=( 2.43342125839683E+00, 5.72183237055104E-01) ( 2, 3)=(-2.09874895849320E+00, 5.85441109710767E-01) ( 2, 4)=( 2.39243684955681E+00, 7.51040934757044E-01) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 1.55013000069668E+00, 0.00000000000000E+00) ( 3, 3)=(-8.59900003064926E-01,-2.75481567205145E-01) ( 3, 4)=( 6.10761316860971E-01, 7.08544799427637E-01) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 4.40401361558925E-01, 0.00000000000000E+00) ( 4, 4)=(-3.75646427119184E-01,-4.19882505172425E-01) cxhes results B matrix ( 1, 1)=( 1.03923048454133E+01, 0.00000000000000E+00) ( 1, 2)=( 1.05293483205003E+01, 5.91347705313051E+00) ( 1, 3)=(-6.25805577535850E+00,-9.98562894355634E-01) ( 1, 4)=( 4.52354342706592E+00, 6.54217374187679E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 4.20186054369101E+00, 3.86669229775551E+00) ( 2, 3)=(-3.32642492170112E+00, 1.04868287534201E+00) ( 2, 4)=( 9.74427460650081E-01, 3.27370224613930E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=(-1.90898440955069E+00, 1.55346184421301E+00) ( 3, 4)=( 1.79131211755246E+00, 7.16471541885188E-01) ( 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.28769891869547E+00, 3.58672013489079E+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)=( 1.01537258082017E-01, 5.85157899660972E-01) ( 2, 3)=(-2.65509949735913E-01, 3.32119205985510E-01) ( 2, 4)=( 6.82994668568015E-01, 0.00000000000000E+00) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 5.33180493738035E-01,-3.59336400355185E-02) ( 3, 3)=(-7.82979986967023E-01,-1.46463880144728E-01) ( 3, 4)=(-2.81636666175644E-01, 2.44702819477580E-02) ( 4, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 2)=( 6.01421443163285E-01, 0.00000000000000E+00) ( 4, 3)=( 4.07075028255037E-01,-1.37774790146888E-01) ( 4, 4)=( 1.35833329567750E-01, 6.59658686635044E-01) cxval results A matrix ( 1, 1)=( 1.29794445704901E+01,-2.36750869459400E-01) ( 1, 2)=(-5.92776562009671E+00, 1.23550258505706E-01) ( 1, 3)=( 2.39460226660901E+00,-7.78288493520875E+00) ( 1, 4)=(-7.76436807205454E+00,-1.93884868458262E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=(-1.30514105346515E+00, 1.27114472644927E+00) ( 2, 3)=(-1.31414155878309E-01,-5.13077226696182E-01) ( 2, 4)=(-1.14330702677213E+00, 3.07113665080051E-01) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=( 5.46024298976512E-01, 6.28694792886101E-01) ( 3, 4)=( 1.41131789343449E-01, 9.22289018225629E-04) ( 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)=(-5.02851902561303E-01,-7.04665265289702E-02) cxval results B matrix ( 1, 1)=( 1.60662476435062E+01, 0.00000000000000E+00) ( 1, 2)=(-4.05898419395243E+00,-2.49633019064254E+00) ( 1, 3)=( 1.47062087286451E+00,-7.23778203021444E+00) ( 1, 4)=(-6.75481867102203E+00,-3.74688179872881E+00) ( 2, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 2, 2)=( 2.40217058602559E+00, 0.00000000000000E+00) ( 2, 3)=(-2.01327084762040E+00, 6.74498987529713E-01) ( 2, 4)=( 1.39387438461376E+00, 7.35560110174409E-01) ( 3, 1)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 3, 3)=( 3.46379397111464E+00, 0.00000000000000E+00) ( 3, 4)=(-4.33514550341068E-01, 8.18565462261393E-01) ( 4, 1)=( 1.00269717101977E-14, 0.00000000000000E+00) ( 4, 2)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 3)=( 0.00000000000000E+00, 0.00000000000000E+00) ( 4, 4)=( 4.16357265380934E+00, 0.00000000000000E+00) cxval results Z matrix ( 1, 1)=(-5.12471025813368E-01,-1.49425728517140E-01) ( 1, 2)=(-1.64331232884351E-01,-4.33704490547799E-01) ( 1, 3)=(-2.22967754186674E-01, 6.31763655373095E-01) ( 1, 4)=(-1.65704361729520E-01,-1.53763418870315E-01) ( 2, 1)=(-4.65659836507567E-01,-1.48280782311016E-01) ( 2, 2)=(-2.65857395868339E-01, 8.22665105079545E-02) ( 2, 3)=( 8.72989874721428E-02,-2.01559596694678E-01) ( 2, 4)=( 7.94678721361265E-01,-6.29614020981949E-02) ( 3, 1)=(-6.21658763519715E-01,-3.34439098217670E-02) ( 3, 2)=( 6.44572458097741E-01, 1.34003996136518E-01) ( 3, 3)=( 2.48041451308364E-01,-1.80077211959532E-01) ( 3, 4)=(-2.27185294264982E-01, 1.82826803701503E-01) ( 4, 1)=( 1.58290673876538E-02, 2.97305319293995E-01) ( 4, 2)=( 2.21114971773556E-01, 4.74473573159973E-01) ( 4, 3)=(-4.11298976336278E-01, 4.89687112695082E-01) ( 4, 4)=( 2.89137481037440E-01, 3.80501230168746E-01) cxval error= 0 ALFR(I) ALFI(I) BETA(I) 1 1.29794445704901E+01 -2.36750869459400E-01 1.60662476435062E+01 2 -1.30514105346515E+00 1.27114472644927E+00 2.40217058602559E+00 3 5.46024298976512E-01 6.28694792886101E-01 3.46379397111464E+00 4 -5.02851902561303E-01 -7.04665265289702E-02 4.16357265380934E+00 cxvec final eigenvctors ( 1, 1)=(-8.23170302533290E-01,-2.40019076111644E-01) ( 1, 2)=(-5.32895546167177E-01,-8.46181030793758E-01) ( 1, 3)=(-1.59597085685057E-01, 3.76089639051483E-01) ( 1, 4)=(-4.16156033993671E-01, 5.33943988159099E-02) ( 2, 1)=(-7.47978576714956E-01,-2.38179975621264E-01) ( 2, 2)=(-6.61112168696531E-01,-4.43017703852335E-02) ( 2, 3)=( 2.50687236804254E-01,-5.54003217320536E-01) ( 2, 4)=( 9.61914772993826E-01,-7.88011754199423E-02) ( 3, 1)=(-9.98556028854121E-01,-5.37201753449100E-02) ( 3, 2)=( 6.24477063963740E-01, 7.64570919378335E-02) ( 3, 3)=( 6.64279548428045E-01,-7.09970257609573E-01) ( 3, 4)=(-9.96354459756025E-01,-8.53099673208182E-02) ( 4, 1)=( 2.54258631883320E-02, 4.77554627092461E-01) ( 4, 2)=( 2.96931358406698E-01, 9.09871648319794E-01) ( 4, 3)=(-6.86775248977112E-01, 7.26869835247292E-01) ( 4, 4)=( 4.36804119443368E-01, 1.36815217634592E-01) cxcheck resdul= 1.55129584733915E-16 eigdet check eigenvalues ( 5.36860057401588E+03,-3.15347284006587E+04) = DETERMINANT (-1.29722258417982E+06, 5.11824960758573E+06) = DETERMINANT ( 1.58587013218631E+08,-4.34203383332994E+08) = DETERMINANT (-5.50554599905819E+08,-1.64504847996577E+09) = DETERMINANT evalvec check eigenvectors SUM( 1)= 8.18698667521032E+02 SUM( 2)= 2.47666714249396E+01 SUM( 3)= 1.13583074317665E+01 SUM( 4)= 6.77428884751150E+00 TOTAL EIGENVECTOR ERROR = 8.61597935225250E+02 test_535_eigen.adb finished