Ice Sheet System Model  4.18
Code documentation
Macros | Typedefs | Functions | Variables
issmmpi.h File Reference
#include <cstddef>
#include <cassert>
#include "../../shared/Numerics/types.h"

Go to the source code of this file.

Macros

#define ISSM_MPI_CHAR   1
 
#define ISSM_MPI_DOUBLE   2
 
#define ISSM_MPI_PDOUBLE   3
 
#define ISSM_MPI_INT   4
 
#define ISSM_MPI_LONG_LONG_INT   5
 
#define ISSM_MPI_MAX   1
 
#define ISSM_MPI_MIN   2
 
#define ISSM_MPI_PROD   3
 
#define ISSM_MPI_SUM   4
 
#define ISSM_MPI_COMM_WORLD   1
 
#define ISSM_MPI_STATUS_IGNORE   &ourIssmMPIStatusIgnore
 
#define ISSM_MPI_ANY_TAG   2
 
#define ISSM_MPI_ANY_SOURCE   3
 

Typedefs

typedef int ISSM_MPI_Comm
 
typedef int ISSM_MPI_Datatype
 
typedef int ISSM_MPI_Op
 
typedef int ISSM_MPI_Status
 

Functions

template<class T >
ISSM_MPI_Datatype TypeToMPIType ()
 
template<>
ISSM_MPI_Datatype TypeToMPIType< IssmDouble > ()
 
template<>
ISSM_MPI_Datatype TypeToMPIType< int > ()
 
template<>
ISSM_MPI_Datatype TypeToMPIType< char > ()
 
template<class T >
int ISSM_MPI_Bcast (T *buffer, int count, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Allgather (void *sendbuf, int sendcount, ISSM_MPI_Datatype sendtype, void *recvbuf, int recvcount, ISSM_MPI_Datatype recvtype, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Allgatherv (void *sendbuf, int sendcount, ISSM_MPI_Datatype sendtype, void *recvbuf, int *recvcounts, int *displs, ISSM_MPI_Datatype recvtype, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Allreduce (void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Barrier (ISSM_MPI_Comm comm)
 
int ISSM_MPI_Bcast (void *buffer, int count, ISSM_MPI_Datatype datatype, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Comm_free (ISSM_MPI_Comm *comm)
 
int ISSM_MPI_Comm_rank (ISSM_MPI_Comm comm, int *rank)
 
int ISSM_MPI_Comm_size (ISSM_MPI_Comm comm, int *size)
 
int ISSM_MPI_Finalize (void)
 
int ISSM_MPI_Gather (void *sendbuf, int sendcnt, ISSM_MPI_Datatype sendtype, void *recvbuf, int recvcnt, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Gatherv (void *sendbuf, int sendcnt, ISSM_MPI_Datatype sendtype, void *recvbuf, int *recvcnts, int *displs, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Init (int *argc, char ***argv)
 
int ISSM_MPI_Recv (void *buf, int count, ISSM_MPI_Datatype datatype, int source, int tag, ISSM_MPI_Comm comm, ISSM_MPI_Status *status)
 
int ISSM_MPI_Reduce (void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Scatter (void *sendbuf, int sendcnt, ISSM_MPI_Datatype sendtype, void *recvbuf, int recvcnt, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Scatterv (void *sendbuf, int *sendcnts, int *displs, ISSM_MPI_Datatype sendtype, void *recvbuf, int recvcnt, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm)
 
int ISSM_MPI_Send (void *buf, int count, ISSM_MPI_Datatype datatype, int dest, int tag, ISSM_MPI_Comm comm)
 
double ISSM_MPI_Wtime (void)
 
int ISSM_MPI_Comm_split (ISSM_MPI_Comm comm, int color, int key, ISSM_MPI_Comm *newcomm)
 
int ISSM_MPI_Intercomm_create (ISSM_MPI_Comm comm, int local_leader, ISSM_MPI_Comm peer_comm, int remote_leader, int tag, ISSM_MPI_Comm *newintercomm)
 
void ISSM_MPI_ContiguousInAdolc (size_t aSize)
 

Variables

ISSM_MPI_Status ourIssmMPIStatusIgnore
 

Macro Definition Documentation

◆ ISSM_MPI_CHAR

#define ISSM_MPI_CHAR   1

Definition at line 124 of file issmmpi.h.

◆ ISSM_MPI_DOUBLE

#define ISSM_MPI_DOUBLE   2

Definition at line 125 of file issmmpi.h.

◆ ISSM_MPI_PDOUBLE

#define ISSM_MPI_PDOUBLE   3

Definition at line 126 of file issmmpi.h.

◆ ISSM_MPI_INT

#define ISSM_MPI_INT   4

Definition at line 127 of file issmmpi.h.

◆ ISSM_MPI_LONG_LONG_INT

#define ISSM_MPI_LONG_LONG_INT   5

Definition at line 128 of file issmmpi.h.

◆ ISSM_MPI_MAX

#define ISSM_MPI_MAX   1

Definition at line 131 of file issmmpi.h.

◆ ISSM_MPI_MIN

#define ISSM_MPI_MIN   2

Definition at line 132 of file issmmpi.h.

◆ ISSM_MPI_PROD

#define ISSM_MPI_PROD   3

Definition at line 133 of file issmmpi.h.

◆ ISSM_MPI_SUM

#define ISSM_MPI_SUM   4

Definition at line 134 of file issmmpi.h.

◆ ISSM_MPI_COMM_WORLD

#define ISSM_MPI_COMM_WORLD   1

Definition at line 137 of file issmmpi.h.

◆ ISSM_MPI_STATUS_IGNORE

#define ISSM_MPI_STATUS_IGNORE   &ourIssmMPIStatusIgnore

Definition at line 139 of file issmmpi.h.

◆ ISSM_MPI_ANY_TAG

#define ISSM_MPI_ANY_TAG   2

Definition at line 140 of file issmmpi.h.

◆ ISSM_MPI_ANY_SOURCE

#define ISSM_MPI_ANY_SOURCE   3

Definition at line 141 of file issmmpi.h.

Typedef Documentation

◆ ISSM_MPI_Comm

typedef int ISSM_MPI_Comm

Definition at line 118 of file issmmpi.h.

◆ ISSM_MPI_Datatype

typedef int ISSM_MPI_Datatype

Definition at line 119 of file issmmpi.h.

◆ ISSM_MPI_Op

typedef int ISSM_MPI_Op

Definition at line 120 of file issmmpi.h.

◆ ISSM_MPI_Status

typedef int ISSM_MPI_Status

Definition at line 121 of file issmmpi.h.

Function Documentation

◆ TypeToMPIType()

template<class T >
ISSM_MPI_Datatype TypeToMPIType ( )

Definition at line 146 of file issmmpi.h.

146 {assert(false);};

◆ TypeToMPIType< IssmDouble >()

template<>
ISSM_MPI_Datatype TypeToMPIType< IssmDouble > ( )
inline

Definition at line 147 of file issmmpi.h.

147 {return ISSM_MPI_DOUBLE;};

◆ TypeToMPIType< int >()

template<>
ISSM_MPI_Datatype TypeToMPIType< int > ( )
inline

Definition at line 151 of file issmmpi.h.

151 {return ISSM_MPI_INT;};

◆ TypeToMPIType< char >()

template<>
ISSM_MPI_Datatype TypeToMPIType< char > ( )
inline

Definition at line 152 of file issmmpi.h.

152 {return ISSM_MPI_CHAR;};

◆ ISSM_MPI_Bcast() [1/2]

template<class T >
int ISSM_MPI_Bcast ( T *  buffer,
int  count,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 154 of file issmmpi.h.

154  { /*{{{*/
155 
156  int rc=0;
157 
158  /*Get MPI type*/
159  ISSM_MPI_Datatype datatype = TypeToMPIType<T>();
160 
161 #ifdef _HAVE_MPI_
162 # ifdef _HAVE_AMPI_
163  rc=AMPI_Bcast(buffer,
164  count,
165  datatype,
166  root,
167  comm);
168 # else
169  rc=MPI_Bcast(buffer,
170  count,
171  datatype,
172  root,
173  comm);
174 # endif
175 #else
176  // nothing to be done here
177 #endif
178  return rc;
179 }/*}}}*/

◆ ISSM_MPI_Allgather()

int ISSM_MPI_Allgather ( void *  sendbuf,
int  sendcount,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int  recvcount,
ISSM_MPI_Datatype  recvtype,
ISSM_MPI_Comm  comm 
)

Definition at line 45 of file issmmpi.cpp.

45  { /*{{{*/
46  int rc=0;
47  assert(sendcount==recvcount || sendtype==recvtype); // we handle only identical representations
48 #ifdef _HAVE_MPI_
49 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
50  rc=AMPI_Allgather(sendbuf,
51  sendcount,
52  sendtype,
53  recvbuf,
54  recvcount,
55  recvtype,
56  comm);
57 # else
58  rc=MPI_Allgather(sendbuf,
59  sendcount,
60  sendtype,
61  recvbuf,
62  recvcount,
63  recvtype,
64  comm);
65 # endif
66 #else
67 # ifdef _HAVE_AD_
68  if (sendtype==ISSM_MPI_DOUBLE) {
69  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
70  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
71  for(int i=0;i<sendcount;++i) activeRecvBuf[i]=activeSendBuf[i];
72  }
73  else
74 # endif
75  memcpy(recvbuf,sendbuf,sizeHelper(sendtype)*sendcount);
76 #endif
77  return rc;
78 } /*}}}*/

◆ ISSM_MPI_Allgatherv()

int ISSM_MPI_Allgatherv ( void *  sendbuf,
int  sendcount,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int *  recvcounts,
int *  displs,
ISSM_MPI_Datatype  recvtype,
ISSM_MPI_Comm  comm 
)

Definition at line 79 of file issmmpi.cpp.

79  { /*{{{*/
80  int rc=0;
81  assert(sendtype==recvtype); // we handle only identical representations
82 #ifdef _HAVE_MPI_
83 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
84  rc=AMPI_Allgatherv(sendbuf,
85  sendcount,
86  sendtype,
87  recvbuf,
88  recvcounts,
89  displs,
90  recvtype,
91  comm);
92 # else
93  rc=MPI_Allgatherv(sendbuf,
94  sendcount,
95  sendtype,
96  recvbuf,
97  recvcounts,
98  displs,
99  recvtype,
100  comm);
101 # endif
102 #else
103  assert(sendcount==recvcounts[0]); // we handle only identical representations
104 # ifdef _HAVE_AD_
105  if (sendtype==ISSM_MPI_DOUBLE) {
106  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
107  IssmDouble* activeRecvBuf=(IssmDouble*)(recvbuf)+displs[0];
108  for(int i=0;i<sendcount;++i) activeRecvBuf[i]=activeSendBuf[i];
109  }
110  else
111 # endif
112  memcpy((char*)recvbuf+(sizeHelper(recvtype)*displs[0]),sendbuf,sizeHelper(sendtype)*sendcount);
113 #endif
114  return rc;
115 }/*}}}*/

◆ ISSM_MPI_Allreduce()

int ISSM_MPI_Allreduce ( void *  sendbuf,
void *  recvbuf,
int  count,
ISSM_MPI_Datatype  datatype,
ISSM_MPI_Op  op,
ISSM_MPI_Comm  comm 
)

Definition at line 116 of file issmmpi.cpp.

116  {/*{{{*/
117 
118  int rc=0;
119 #ifdef _HAVE_MPI_
120 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
121  rc=AMPI_Allreduce(sendbuf,
122  recvbuf,
123  count,
124  datatype,
125  op,
126  comm);
127 # else
128  rc=MPI_Allreduce(sendbuf,
129  recvbuf,
130  count,
131  datatype,
132  op,
133  comm);
134 # endif
135 #else
136 #ifdef _HAVE_AD_
137  if (datatype==ISSM_MPI_DOUBLE) {
138  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
139  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
140  for(int i=0;i<count;++i) activeRecvBuf[i]=activeSendBuf[i];
141  }
142  else
143 # endif
144  memcpy(recvbuf,sendbuf,sizeHelper(datatype)*count);
145 #endif
146  return rc;
147 }/*}}}*/

◆ ISSM_MPI_Barrier()

int ISSM_MPI_Barrier ( ISSM_MPI_Comm  comm)

Definition at line 148 of file issmmpi.cpp.

148  { /*{{{*/
149 
150  int rc=0;
151 #ifdef _HAVE_MPI_
152 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
153  rc=AMPI_Barrier(comm);
154 # else
155  rc=MPI_Barrier(comm);
156 # endif
157 #else
158 // do nothing
159 #endif
160  return rc;
161 }/*}}}*/

◆ ISSM_MPI_Bcast() [2/2]

int ISSM_MPI_Bcast ( void *  buffer,
int  count,
ISSM_MPI_Datatype  datatype,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 162 of file issmmpi.cpp.

162  { /*{{{*/
163 
164  int rc=0;
165 #ifdef _HAVE_MPI_
166 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
167  rc=AMPI_Bcast(buffer,
168  count,
169  datatype,
170  root,
171  comm);
172 # else
173  rc=MPI_Bcast(buffer,
174  count,
175  datatype,
176  root,
177  comm);
178 # endif
179 #else
180 // nothing to be done here
181 #endif
182  return rc;
183 }/*}}}*/

◆ ISSM_MPI_Comm_free()

int ISSM_MPI_Comm_free ( ISSM_MPI_Comm *  comm)

Definition at line 184 of file issmmpi.cpp.

184  { /*{{{*/
185 
186  int rc=0;
187 #ifdef _HAVE_MPI_
188 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
189  assert(0); // to be implemented
190 # else
191  rc=MPI_Comm_free(comm);
192 # endif
193 #else
194 // do nothing
195 #endif
196  return rc;
197 }/*}}}*/

◆ ISSM_MPI_Comm_rank()

int ISSM_MPI_Comm_rank ( ISSM_MPI_Comm  comm,
int *  rank 
)

Definition at line 198 of file issmmpi.cpp.

198  { /*{{{*/
199 
200  int rc=0;
201 #ifdef _HAVE_MPI_
202  rc=MPI_Comm_rank(comm,
203  rank);
204 #else
205  *rank=0;
206 #endif
207  return rc;
208 }/*}}}*/

◆ ISSM_MPI_Comm_size()

int ISSM_MPI_Comm_size ( ISSM_MPI_Comm  comm,
int *  size 
)

Definition at line 209 of file issmmpi.cpp.

209  { /*{{{*/
210 
211  int rc=0;
212 #ifdef _HAVE_MPI_
213  rc=MPI_Comm_size(comm,
214  size);
215 #else
216  *size=1;
217 #endif
218  return rc;
219 }/*}}}*/

◆ ISSM_MPI_Finalize()

int ISSM_MPI_Finalize ( void  )

Definition at line 220 of file issmmpi.cpp.

220  { /*{{{*/
221 
222  int rc=0;
223  #ifdef _HAVE_MPI_
224  #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
225  #if defined(_HAVE_ADJOINTMPI_)
226  rc=AMPI_Finalize();
227  #elif defined(_HAVE_MEDIPACK_)
228  /*Old implementation*/
229  //TOOL::finalize();
230  /*New implementation*/
231  delete mpiTypes;
232  rc=AMPI_Finalize();
233  #else
234  rc=AMPI_Finalize_NT();
235  #endif
236  #else
237  rc=MPI_Finalize();
238  #endif
239  #endif
240  return rc;
241 }/*}}}*/

◆ ISSM_MPI_Gather()

int ISSM_MPI_Gather ( void *  sendbuf,
int  sendcnt,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int  recvcnt,
ISSM_MPI_Datatype  recvtype,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 242 of file issmmpi.cpp.

242  { /*{{{*/
243 
244  int rc=0;
245  assert(sendtype==recvtype && sendcnt==recvcnt); // we handle only identical representations
246 #ifdef _HAVE_MPI_
247 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
248  rc=AMPI_Gather(sendbuf,
249  sendcnt,
250  sendtype,
251  recvbuf,
252  recvcnt,
253  recvtype,
254  root,
255  comm);
256 # else
257  rc=MPI_Gather(sendbuf,
258  sendcnt,
259  sendtype,
260  recvbuf,
261  recvcnt,
262  recvtype,
263  root,
264  comm);
265 # endif
266 #else
267 # ifdef _HAVE_AD_
268  if (sendtype==ISSM_MPI_DOUBLE) {
269  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
270  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
271  for(int i=0;i<sendcnt;++i) activeRecvBuf[i]=activeSendBuf[i];
272  }
273  else
274 # endif
275  memcpy(recvbuf,sendbuf,sizeHelper(sendtype)*sendcnt);
276 #endif
277  return rc;
278 }/*}}}*/

◆ ISSM_MPI_Gatherv()

int ISSM_MPI_Gatherv ( void *  sendbuf,
int  sendcnt,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int *  recvcnts,
int *  displs,
ISSM_MPI_Datatype  recvtype,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 279 of file issmmpi.cpp.

279  {/*{{{*/
280 
281  int rc=0;
282  assert(sendtype==recvtype); // we handle only identical representations
283 #ifdef _HAVE_MPI_
284 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
285  rc=AMPI_Gatherv(sendbuf,
286  sendcnt,
287  sendtype,
288  recvbuf,
289  recvcnts,
290  displs,
291  recvtype,
292  root,
293  comm);
294 # else
295  rc=MPI_Gatherv(sendbuf,
296  sendcnt,
297  sendtype,
298  recvbuf,
299  recvcnts,
300  displs,
301  recvtype,
302  root,
303  comm);
304 # endif
305 #else
306  assert(sendcnt==recvcnts[0]); // we handle only identical representations
307 #ifdef _HAVE_AD_
308  if (sendtype==ISSM_MPI_DOUBLE) {
309  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
310  IssmDouble* activeRecvBuf=(IssmDouble*)(recvbuf)+displs[0];
311  for(int i=0;i<sendcnt;++i) activeRecvBuf[i]=activeSendBuf[i];
312  }
313  else
314 # endif
315  memcpy((char*)recvbuf+(sizeHelper(recvtype)*displs[0]),sendbuf,sizeHelper(sendtype)*sendcnt);
316 #endif
317  return rc;
318 }/*}}}*/

◆ ISSM_MPI_Init()

int ISSM_MPI_Init ( int *  argc,
char ***  argv 
)

Definition at line 319 of file issmmpi.cpp.

319  { /*{{{*/
320 
321  int rc=0;
322  #ifdef _HAVE_MPI_
323  #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
324  #if defined(_HAVE_ADJOINTMPI_)
325  rc=AMPI_Init(argc,argv);
326  #elif defined(_HAVE_MEDIPACK_)
327  rc=AMPI_Init(argc,argv);
328  /*Old implementation of Medipack*/
329  //TOOL::init();
330  /*New*/
331  //MpiTypes* mpiTypes;
332  mpiTypes = new MpiTypes();
333  #else
334  rc=AMPI_Init_NT(argc,argv);
335  #endif
336  #else
337  rc=MPI_Init(argc,argv);
338  #endif
339  #endif
340  return rc;
341 }/*}}}*/

◆ ISSM_MPI_Recv()

int ISSM_MPI_Recv ( void *  buf,
int  count,
ISSM_MPI_Datatype  datatype,
int  source,
int  tag,
ISSM_MPI_Comm  comm,
ISSM_MPI_Status *  status 
)

Definition at line 342 of file issmmpi.cpp.

342  { /*{{{*/
343 
344  int rc=0;
345 #ifdef _HAVE_MPI_
346 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
347  rc=AMPI_Recv(buf,
348  count,
349  datatype,
350  source,
351  tag,
352  #if !defined(_HAVE_ADJOINTMPI_) && !defined(_HAVE_MEDIPACK_)
353  AMPI_FROM_SEND, // as long as there are no other variants
354  #endif
355  comm,
356  status);
357 # else
358  rc=MPI_Recv(buf,
359  count,
360  datatype,
361  source,
362  tag,
363  comm,
364  status);
365 # endif
366 #else
367 // nothing to be done here
368 // as long as nobody tries to do anything with 'status'
369 // we won't do anything to it here either
370 #endif
371  return rc;
372 }/*}}}*/

◆ ISSM_MPI_Reduce()

int ISSM_MPI_Reduce ( void *  sendbuf,
void *  recvbuf,
int  count,
ISSM_MPI_Datatype  datatype,
ISSM_MPI_Op  op,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 373 of file issmmpi.cpp.

373  { /*{{{*/
374 
375  int rc=0;
376 #ifdef _HAVE_MPI_
377 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
378  rc=AMPI_Reduce(sendbuf,
379  recvbuf,
380  count,
381  datatype,
382  op,
383  root,
384  comm);
385 # else
386  rc=MPI_Reduce(sendbuf,
387  recvbuf,
388  count,
389  datatype,
390  op,
391  root,
392  comm);
393 # endif
394 #else
395 # ifdef _HAVE_AD_
396  if (datatype==ISSM_MPI_DOUBLE) {
397  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
398  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
399  for(int i=0;i<count;++i) activeRecvBuf[i]=activeSendBuf[i];
400  }
401  else
402 # endif
403  memcpy(recvbuf,sendbuf,sizeHelper(datatype)*count);
404 #endif
405  return rc;
406 }/*}}}*/

◆ ISSM_MPI_Scatter()

int ISSM_MPI_Scatter ( void *  sendbuf,
int  sendcnt,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int  recvcnt,
ISSM_MPI_Datatype  recvtype,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 407 of file issmmpi.cpp.

407  { /*{{{*/
408 
409  int rc=0;
410  assert(sendtype==recvtype && sendcnt==recvcnt); // we handle only identical representations
411 #ifdef _HAVE_MPI_
412 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
413  rc=AMPI_Scatter(sendbuf,
414  sendcnt,
415  sendtype,
416  recvbuf,
417  recvcnt,
418  recvtype,
419  root,
420  comm);
421 # else
422  rc=MPI_Scatter(sendbuf,
423  sendcnt,
424  sendtype,
425  recvbuf,
426  recvcnt,
427  recvtype,
428  root,
429  comm);
430 # endif
431 #else
432 # ifdef _HAVE_AD_
433  if (sendtype==ISSM_MPI_DOUBLE) {
434  IssmDouble* activeSendBuf=(IssmDouble*)sendbuf;
435  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
436  for(int i=0;i<recvcnt;++i) activeRecvBuf[i]=activeSendBuf[i];
437  }
438  else
439 # endif
440  memcpy(recvbuf,sendbuf,sizeHelper(sendtype)*recvcnt);
441 #endif
442  return rc;
443 }/*}}}*/

◆ ISSM_MPI_Scatterv()

int ISSM_MPI_Scatterv ( void *  sendbuf,
int *  sendcnts,
int *  displs,
ISSM_MPI_Datatype  sendtype,
void *  recvbuf,
int  recvcnt,
ISSM_MPI_Datatype  recvtype,
int  root,
ISSM_MPI_Comm  comm 
)

Definition at line 444 of file issmmpi.cpp.

444  { /*{{{*/
445 
446  int rc=0;
447  assert(sendtype==recvtype); // we handle only identical representations
448 #ifdef _HAVE_MPI_
449 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
450  rc=AMPI_Scatterv(sendbuf,
451  sendcnts,
452  displs,
453  sendtype,
454  recvbuf,
455  recvcnt,
456  recvtype,
457  root,
458  comm);
459 # else
460  rc=MPI_Scatterv(sendbuf,
461  sendcnts,
462  displs,
463  sendtype,
464  recvbuf,
465  recvcnt,
466  recvtype,
467  root,
468  comm);
469 # endif
470 #else
471  assert(sendcnts[0]==recvcnt); // we handle only identical representations
472 # ifdef _HAVE_AD_
473  if (sendtype==ISSM_MPI_DOUBLE) {
474  IssmDouble* activeSendBuf=(IssmDouble*)(sendbuf)+displs[0];
475  IssmDouble* activeRecvBuf=(IssmDouble*)recvbuf;
476  for(int i=0;i<recvcnt;++i) activeRecvBuf[i]=activeSendBuf[i];
477  }
478  else
479 # endif
480  memcpy(recvbuf,(char*)sendbuf+(sizeHelper(sendtype)*displs[0]),sizeHelper(sendtype)*recvcnt);
481 #endif
482  return rc;
483 }/*}}}*/

◆ ISSM_MPI_Send()

int ISSM_MPI_Send ( void *  buf,
int  count,
ISSM_MPI_Datatype  datatype,
int  dest,
int  tag,
ISSM_MPI_Comm  comm 
)

Definition at line 484 of file issmmpi.cpp.

484  { /*{{{*/
485 
486  int rc=0;
487 #ifdef _HAVE_MPI_
488 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
489  rc=AMPI_Send(buf,
490  count,
491  datatype,
492  dest,
493  tag,
494  #if !defined(_HAVE_ADJOINTMPI_) && !defined(_HAVE_MEDIPACK_)
495  AMPI_TO_RECV, // as long as there are no other variants
496  #endif
497  comm);
498 # else
499  rc=MPI_Send(buf,
500  count,
501  datatype,
502  dest,
503  tag,
504  comm);
505 # endif
506 #else
507 // nothing to be done here
508 #endif
509  return rc;
510 }/*}}}*/

◆ ISSM_MPI_Wtime()

double ISSM_MPI_Wtime ( void  )

Definition at line 511 of file issmmpi.cpp.

511  {/*{{{*/
512 
513 #ifdef _HAVE_MPI_
514  return MPI_Wtime();
515 #else
516  assert(0); // to be implemented
517  return 0.0;
518 #endif
519 }/*}}}*/

◆ ISSM_MPI_Comm_split()

int ISSM_MPI_Comm_split ( ISSM_MPI_Comm  comm,
int  color,
int  key,
ISSM_MPI_Comm *  newcomm 
)

Definition at line 528 of file issmmpi.cpp.

528  { /*{{{*/
529 
530 int rc=0;
531 #ifdef _HAVE_MPI_
532 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
533 rc=MPI_Comm_split(comm, color, key, newcomm);
534 #else
535 rc=MPI_Comm_split(comm, color, key, newcomm);
536 #endif
537 #else
538 // nothing to be done here
539 #endif
540  return rc;
541 }/*}}}*/

◆ ISSM_MPI_Intercomm_create()

int ISSM_MPI_Intercomm_create ( ISSM_MPI_Comm  comm,
int  local_leader,
ISSM_MPI_Comm  peer_comm,
int  remote_leader,
int  tag,
ISSM_MPI_Comm *  newintercomm 
)

Definition at line 542 of file issmmpi.cpp.

542  { /*{{{*/
543 
544  int rc=0;
545 #ifdef _HAVE_MPI_
546 #if defined(_HAVE_AMPI_) && !defined(_WRAPPERS_)
547  rc=MPI_Intercomm_create(comm,local_leader,peer_comm,remote_leader,tag,newintercomm);
548 #else
549  rc=MPI_Intercomm_create(comm,local_leader,peer_comm,remote_leader,tag,newintercomm);
550 #endif
551 #else
552  // nothing to be done here
553 #endif
554  return rc;
555 }/*}}}*/

◆ ISSM_MPI_ContiguousInAdolc()

void ISSM_MPI_ContiguousInAdolc ( size_t  aSize)

Definition at line 520 of file issmmpi.cpp.

520  { /*{{{*/
521 
522 #ifdef _HAVE_ADOLC_
523  ensureContiguousLocations(aSize);
524 #else
525  fprintf(stderr, "*** Codipack ISSM_MPI_ContiguousInAdolc()\n");
526 #endif
527 }/*}}}*/

Variable Documentation

◆ ourIssmMPIStatusIgnore

ISSM_MPI_Status ourIssmMPIStatusIgnore
extern

Definition at line 24 of file issmmpi.cpp.

IssmDouble
double IssmDouble
Definition: types.h:37
ISSM_MPI_Datatype
int ISSM_MPI_Datatype
Definition: issmmpi.h:119
ISSM_MPI_DOUBLE
#define ISSM_MPI_DOUBLE
Definition: issmmpi.h:125
ISSM_MPI_INT
#define ISSM_MPI_INT
Definition: issmmpi.h:127
sizeHelper
size_t sizeHelper(ISSM_MPI_Datatype type)
Definition: issmmpi.cpp:25
ISSM_MPI_CHAR
#define ISSM_MPI_CHAR
Definition: issmmpi.h:124
Generated on Thu Jul 2 2020 08:09:25 for Ice Sheet System Model by   doxygen 1.8.19