Fortran element by element array multiplication

[Kwetla]

Homework Statement

Hello,

I'm having a problem in multiplying two vectors together in a specific way in Fortran. I can do it in Matlab, but can't work out how to do it in Fortran.

The problem is that i want to multiply two vectors together, but only each element by it's corresponding element in the other vector. In Matlab you simply use A.*B which works perfectly.

Homework Equations

For example if A = [1 2 3]'

and B = [4 5 6]'

Then A.*B = [4 10 18]

The Attempt at a Solution

Any attempt i make in Fortran tries to multiply all the elements together (usually by dot product etc), but i just want the first element multiplied by the first element, and then the second mulitplied by the second etc.

Sorry about the strange formatting, i tried to crowbar this question into the template...

 

[Mark44]

Homework Statement

Hello,

I'm having a problem in multiplying two vectors together in a specific way in Fortran. I can do it in Matlab, but can't work out how to do it in Fortran.

The problem is that i want to multiply two vectors together, but only each element by it's corresponding element in the other vector. In Matlab you simply use A.*B which works perfectly.

Homework Equations

For example if A = [1 2 3]'

and B = [4 5 6]'

Then A.*B = [4 10 18]

The Attempt at a Solution

Any attempt i make in Fortran tries to multiply all the elements together (usually by dot product etc), but i just want the first element multiplied by the first element, and then the second mulitplied by the second etc.

Sorry about the strange formatting, i tried to crowbar this question into the template...

Assuming that all variables and the three arrays are declared and that A and B are initialized, this should work.

N = 3
DO I = 1, N
  C(I) = A(I) * B(I)
END DO

 

[Kwetla]

Ah yes, that would work. Although I'm surprised there isn't a function for it.

Thanks a lot though!

 

[LawrenceC]

Vector times vector can also be a scalar, 32.

 

[hotvette]

If you are using at least Fortran 90, I believe the following will work:

C = A*B

or

C(1:N) = A(1:N)*B(1:N)

The 1st version assumes A,B,C have the same dimension. The 2nd version is safer.