Complex arrays support

[seadra]

Apologies if there's already a way to do this, but it there a way remove the assumption that arrays are real numbers?

When doing

@variables u[1:2,1:2]
u'*u

the resulting expression doesn't use complex conjugation, and similarly, conj(u) just gives the matrix itself.

I tried @variables u[1:2,1:2]::Complex but it gives the error

LoadError: AssertionError: @variables expects a tuple of expressions or an expression of a tuple (`@variables x y z(t) v[1:3] w[1:2,1:4]` or `@variables x y z(t) v[1:3] w[1:2,1:4] k=1.0`)
in expression starting at In[42]:1

Similarly, the following

@syms u[1:2,1:2]::Complex

gives

MethodError: Cannot `convert` an object of type Expr to an object of type Symbol
Closest candidates are:
  convert(::Type{T}, ::T) where T at essentials.jl:205
  Symbol(::Any...) at strings/basic.jl:229

[leonandonayre]

This could help you

@variables a[1:2,1:2] b[1:2,1:2]
u = a + im*b

[shashi]

Sorry there's no way to do this properly right now, although @leonandonayre 's trick might give you some success. It will currently not give you an array of eltype <: Number.

[seadra]

Ah, sure, I see.

The workaround is fine on its own, but in my case it's not compatible with the rest of code which doesn't use separate parameters for real and imaginary parts, so I'll put it on hold until this works.

[shashi]

This is in #301

Try

@variables (u::Complex)[1:2,1:2]

on the master branch!