pfft icon indicating copy to clipboard operation
pfft copied to clipboard

Distributed Matrix Transposes

Open croor opened this issue 7 years ago • 1 comments

Is there a PFFT way to transpose pencil-decomposed 3D array globally?

Is using a pfft_plan_many_*_skipped function, and skipping all three transforms a viable option?

croor avatar May 02 '18 12:05 croor

I've been looking into this feature as well. However, this (skip_trafo) is probably different from the kind of transpose you really wanted. I think it would only do the kind of transpose due to TRANSPOSED_IN and TRANSPOSED_OUT; and even for that I suspect it won't work for the half-plane hermitian complex field.

A proper solution will probably need one to at least provide the desired new axes ordering, like this

https://docs.scipy.org/doc/numpy/reference/generated/numpy.transpose.html

and it shall still accept flags like TRANSPOSED_IN and TRANSPOSED_OUT flags, which really is mostly concerned about the data distribution scheme more than a transpose. I believe the global transp routines in pfft can do this, but I don't think I know how to implement this.

rainwoodman avatar May 25 '18 23:05 rainwoodman