Complex support

[jjjkkkjjj]

To use vdsp, DSPSplitComplex seems to be needed according to documents

To use blas package, DSPComplex will be needed according to this discussion

to achieve to support complex type, using DSPComplex is ideal? I must check the difference between DSPComplex and DSPSplitComplex at first.

DSPComplex is the consecutive float values.

Complex data are stored as ordered pairs of floating-point numbers. Because they are stored as ordered pairs, complex vectors require address strides that are multiples of two.

by document

On the other hand, DSPSplitComplex is stored in different memories

A structure that represents a single-precision complex vector with the real and imaginary parts stored in separate arrays.

by document

So I need to implement the function to connect this difference of memory layout

[jjjkkkjjj]

Basically I store complex data as DSPComplex, and then pass DSPSplitComplex(realptr, imgptr)with 2 times strides to vDSP function

[jjjkkkjjj]

To convert DSPComplex into DSPSplitComplex, use vDSP_ctoz

[jjjkkkjjj]

Hmm, Create another object as MfArrayComplex?

[jjjkkkjjj]

Create MfComplexArray with realdata and imagdata simply!

[jjjkkkjjj]

ToDo

• [x] Basic Arithmetic
• [x] Angle
• [x] Conjugate

[jjjkkkjjj]

MfComplexArray inherent from MfArray and MfComplexData inherit from MfData?

[jjjkkkjjj]

How about withUnsafeComplexPointer?

[jjjkkkjjj]

Complex addition was implemented by ca948f4. But, this implementation will be hard for other modification... So, should i integrate MfComplexArray into MfArray?

[jjjkkkjjj]

ToDo list

• [x] Arithmetic Operation
• [x] Angle, Conjugate and Absolute
• [x] Math (partial: sin,cos,tan,exp,log)
• [x] Basic Subscription Getter
• [x] Basic Subscription Setter
• [x] Boolean Indexing Getter
• [ ] Boolean Indexing Setter
• [x] Fancy Indexing Getter
• [ ] Fancy Indexing Setter