On Thu, 17 Sep 2009, Roy Stogner wrote: > On Thu, 17 Sep 2009, Roy Stogner wrote: > >> We use Parallel::datatype<> to wrap a bunch of MPI calls into a much >> nicer interface. Unfortunately, it looks like MPI_DOUBLE_COMPLEX is a >> Fortran-only datatype, and MPI::DOUBLE_COMPLEX requires a C++ MPI2 >> standard whereas we've been trying to stay C MPI1 compatible, so we >> have to write special template definition for Parallel operations >> on std::complex data. If allgather code ever tried to use a >> non-special template definition, this is the error you'd see. >> >> Ben added a new Parallel::allgather() argument a while back, but he >> did it properly on both the regular and complex versions - I'm not >> sure what the problem here could be. I'll keep looking. > > Okay, I've found the problem.
One of the two problems, anyway. Problem two: for some reason the DistributedVector call to Parallel::allgather(std::vector<T>&) with T=std::complex<double> isn't getting resolved to the Parallel::allgather(std::vector<std::complex<T> >&) specialization. I'm not sure how to fix this one... --- Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry® Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9-12, 2009. Register now! http://p.sf.net/sfu/devconf _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
