@MSGID:
<e54a9d2d-7ae1-4e2b-ac8f-ddc6daf8d1f8n@googlegroups.com> ec0cdbba
@REPLY:
<be249762-6f1d-4b0d-8c37-bfd262721bd8n@googlegroups.com> 3a24e025
@REPLYADDR gah4 <gah4@u.washington.edu>
@REPLYTO 2:5075/128 gah4
@CHRS: CP866 2
@RFC: 1 0
@RFC-References:
<f8466a68-55d1-4d8c-a607-5c41c85404d2n@googlegroups.com> <4d061409-7c22-4e8c-a020-6e4476e8f165n@googlegroups.com>
<be249762-6f1d-4b0d-8c37-bfd262721bd8n@googlegroups.com>
@RFC-Message-ID:
<e54a9d2d-7ae1-4e2b-ac8f-ddc6daf8d1f8n@googlegroups.com>
@TZUTC: -0700
@PID: G2/1.0
@TID: FIDOGATE-5.12-ge4e8b94
On Friday, September 1, 2023 at 12:37:16 AM UTC-7, Farzad Tatar wrote:

(snip)

> real:: TEST(2,2) = TRANSPOSE(reshape((/1.,2.,3.,4./), (/2,2/))) 
> real:: TEST2(2) = (1.,1.) 
> integer :: ipiv2(2) 
> integer, parameter :: sizet=2 
> goto 999 
> ... 
> 
> 999 NRHS=1 
> print*, `before solving`, test, test2 
> CALL DGESV(sizet, nrhs, TEST, sizet, IPIV2, TEST2, sizet, INFO) 
> print*, ` after solving:`, test, test2 
> print*, ` ipiv2 vector: ` ipiv2 

(snip)

DGESV needs double precision 3rd and 6th argument.

If you give it single precision, it will try to access, in this case,
four double precision variables in an array big enough for two.

The  -1.7014126E+38 is the most negative single precision value
on many systems.  That is a sign that something is wrong.

If you change to SGESV it should work.




--- G2/1.0
 * Origin: usenet.network (2:5075/128)
SEEN-BY: 5001/100 5005/49 5015/255 5019/40 5020/715
848 1042 4441 12000
SEEN-BY: 5030/49 1081 5058/104 5075/128
@PATH: 5075/128 5020/1042 4441