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])
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|