[Fortran] Quick question about complex exponentials

[Matterwave]

Hey, so I just have a quick question. I am trying to set a complex variable (in an array) as $e^{i\alpha_1}$ and the line I used in my code looks like this:

  hmajphasemix(2,2)=(cos(alpha1),sin(alpha1))

But the compiler is telling me that it "expects a right parenthesis" at this line. I'm kind of confused cus I'm counting the right number of parenthesis. Is this the correct way to set this variable or have I messed something up?

 

[gsal]

Only literal constants can be assigned that way:

cnumber = (1.2, 2.3)

When using variables, as you are doing, then you need:

hmajphasemix(2,2) = cmplx(cos(alpha1),sin(alpha1))

 

[Matterwave]

Ah, so I can only use (X,Y) when X and Y are just numbers, but if X and Y are some functions of some variable then I need to use cmplx(X,Y) is that right?

 

[gsal]

Yeap

 

[gsal]

Not even functions, even if they just variables... Without cmplx they need to be literal contacts

 

[gsal]

Contacts->constants

 

[Matterwave]

Ah, thanks! What if X and Y are parameters? Would the simpler one work then?

 

[gsal]

Don't think so...anyway, you can do your own testing now that you know what's going on.

 

[Mark44]

anyway, you can do your own testing now that you know what's going on.

This^

 

[Matterwave]

Alrighty. Thanks again. :)

 

[serp777]

One question: why would you ever still program in fortran lol?

 

[Mark44]

One question: why would you ever still program in fortran lol?

Apparently there are good reasons to do so. There is a ton of existing Fortran code for scientific and engineering applications. Also, Fortran code lends itself more easily to highly parallelized applications than do C and C++ programs, or so I understand. One of the drawbacks in C/C++ code is that data can be accessed indirectly (via pointers), which can be problematic where a block of data has to be accessed by multiple threads. The relatively new restrict keyword in C (C99) is an effort to restrict the access to data by only a single pointer.

 

[meBigGuy]

The last time I programmed in FORTRAN (an FM multipath simulation) it was because of the ease of working with complex numbers.
(and I happened to have a graphing package)
I'd probably do it in MATLAB now and it would take forever to run.

Maybe there are some good complex math packages for C or C++ now.