Intel® Fortran Compiler
Build applications that can scale for the future with optimized code designed for Intel® Xeon® and compatible processors.

array question

liatchin
Beginner
1,109 Views

I have a two-dimensional matrix where one of its elements depends on subject I

I create a three-dimensional matrix A where the first subscript is I.

DIMENSION A(N_SUB, N_PAR, N_PAR+1)

where N_SUB = number of subjects and N_PAR = 7

I am trying to multiply the transpose of A to a constant matrix B (7 by 7).

I would like to use the routine dmxtyf to obtain the product. This routine is supposed to multiply two two-dimensional matrices.

How can I convert A to a two-dimensional matrix so I can use the routine dmxtyf?

I tried

DIMENSION MT(N_PAR, N_PAR)

DIMENSION E(N_PAR+1,N_PAR+1)

DO 10 I = 1, N_SUB

DO 9 J = 1, N_PAR

DO 8 K =1, N_PAR+1

MT(J,K) = A(I,J,K)

8 CONTINUE

9 CONTINUE

CALLDMXTYF(7,8,MT,7,7,7,SMATINV,7,8,7,E,8)

DO 7 J=1,N_PAR+1

WRITE(20,130)(E(J,K),K=1,N_PAR+1)

7 CONTINUE

10 CONTINUE

But it gives me all zeros!

SOMEONE CAN HELP ME PLEASE.

0 Kudos
0 Replies
Reply