J'utilise python + numpy + scipy pour faire du filtrage par convolution sur un tableau de nombres complexes.
field = np.zeros((field_size, field_size), dtype=complex)
...
field = scipy.signal.convolve(field, kernel, 'same')Ainsi, lorsque je veux utiliser un tableau complexe dans numpy, tout ce que je dois faire est de passer le paramètre dtype=complex. Pour mes recherches, j'ai besoin d'implémenter deux autres types de nombres complexes : dual (i*i=0) et double (i*i=1). Ce n'est pas un gros problème - je prends simplement le code source python pour les nombres complexes et je change la fonction de multiplication. Le problème : comment faire un tableau numpy de ces types numériques exotiques ?
