search

pchanial / pyoperators

Operators and solvers for high-performance computing.
http://pchanial.github.com/pyoperators
Other
14 stars 9 forks source link

cg is not compatible with complex dtype'd operators #11

Open ghisvail opened 11 years ago

ghisvail commented 11 years ago

In the current implementation of CG, there is an explicit type check on the dtype to exclude complex dtype'd computation. What forces such restriction ?

pchanial commented 11 years ago

Hi, I haven't used cg for complex dtype so I'd rather have something not implemented than wrong :-) Feel free to modify the small bits and have a separate _norm2_complex function for complex dtype.