Jaehoon Seol asked: >1. Can anyone tell me why the ouput of the array content a in testGSolver > is different from that of the array content a in dgeco. No, but your program certainly appears to have a mistake in it. Anyway, I don't see any Fortran 77 here, both program units are Fortran 90. >2. I also appreciate if you can tell me what difference there is between > double precision a(lda,1) > and > double precision a(lda,n) The first one declares an array with 1 column of LDA elements. The second one declares an array with N columns of LDA elements. > call dgeco(at,nrow,nrow,pivot_info,cond_num,extras) ^^ This variable does not exist - presumably you meant "a". (implicit none ought to have warned you about this - indeed the NAGWare compiler gives an error here). >=========== output result ==================== > dgeco : lda = 2 > dgeco : n = 2 > dgeco : shape of a = 2 2 > dgeco : 11.0 10.0 > dgeco : 10.0 10.0 > 1.0 3.0 > 2.0 5.0 > pivot information: 0 0 > condition number : 0.0 The output I get, after changing "at" to "a", is: dgeco : lda = 2 dgeco : n = 2 dgeco : shape of a =: 2 2 dgeco : 11.0000000000000000 10.0000000000000000 dgeco : 10.0000000000000000 10.0000000000000000 11.0000000000000000 10.0000000000000000 10.0000000000000000 10.0000000000000000 pivot information: 1 134281576 condition number : 5.1528479173536188E-270 ...I don't see any problem with this output. Cheers, -- ...........................Malcolm Cohen, NAG Ltd., Oxford, U.K. ([log in to unmask]) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%