[SOLVED] vectors in fortran90

aihaike · 04-23-2009, 09:07 AM

Dear all,

I'm not an expert of Fortran 90.
I'm wondering how to perform such kind of thing :

Code:

(/l,m,n/) = matmul(mat,(/1.,2.,1./)) (doesn't work)

where mat is a well defined array.

I just want tot avoid the creation of a vector results

Code:

res = matmul(mat,(/1.,2.,1./)) (works, of course)

Any idea?

Eric.

paulsm4 · 04-24-2009, 11:36 AM

Hi -

What's wrong with the good old-fashioned FORTRAN IV technique of declaring a function with three arguments: the two arrays you're multiplying, and a third array to hold the product?

aihaike · 04-25-2009, 06:53 AM

Hi,

you're right, I can do that.
Thanks!,

Eric.