Changeset 17669

Timestamp:

04/08/14 10:29:27 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: removing template that makes AD mpi off crash

File:

• issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h (modified) (1 diff)

issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h

@@ -89 +89 @@
 template <> inline ISSM_MPI_Datatype TypeToMPIType<int>(){return ISSM_MPI_INT;};
 template <> inline ISSM_MPI_Datatype TypeToMPIType<char>(){return ISSM_MPI_CHAR;};
-template <class T> ISSM_MPI_Datatype ISSM_MPI_Reduce(T *sendbuf,T *recvbuf,int count,ISSM_MPI_Op op, int root, ISSM_MPI_Comm comm){ /*{{{*/
-
-        /*Get MPI type*/
-        ISSM_MPI_Datatype datatype = TypeToMPIType<T>();
-
-        int rc=0;
-#ifdef _HAVE_MPI_
-# ifdef _HAVE_AMPI_
-        rc=AMPI_Reduce(sendbuf,
-                                recvbuf,
-                                count,
-                                datatype,
-                                op,
-                                root,
-                                comm);
-# else
-        rc=MPI_Reduce(sendbuf,
-                                recvbuf,
-                                count,
-                                datatype,
-                                op,
-                                root,
-                                comm);
-# endif
-#else
-# ifdef _HAVE_ADOLC_
-        if (datatype==ISSM_MPI_DOUBLE) { 
-                IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
-                IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
-                for(int i=0;i<count;++i) activeRecvBuf[i]=activeSendBuf[i];
-        }
-        else 
-# endif
-         memcpy(recvbuf,sendbuf,sizeHelper(datatype)*count);
-#endif
-        return rc;
-}/*}}}*/
 template <class T> int ISSM_MPI_Bcast(T *buffer, int count,int root, ISSM_MPI_Comm comm){  /*{{{*/