complex arithmetic in C

[Cleve Ashcraft]

i've got to do some complex analysis in C, and this is the one time
i look back kindly to all the fortran programming which i've done.

i've seen examples of creating a complex struct and writing
procedures to operate on complex numbers, but those always struck
me as lessons on structures rather then a good implementation.

i've got a bit of algorithm development to do, for which ease of
coding is most important. later i'll optimize the computations.
the kernels are CAXPY's and CDOT's, and their sparse indexed versions.

dealing with struct's may be easy, but will it later have to be
hardcoded into two arrays for efficiency?

any suggestions/experience to pass on? thanks.

cleve ashcraft
ashcraft at cs.yale.edu