pingpong.c 2.1 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788
  1. /* StarPU --- Runtime system for heterogeneous multicore architectures.
  2. *
  3. * Copyright (C) 2009, 2010, 2014-2015 Université de Bordeaux
  4. * Copyright (C) 2010, 2011, 2012, 2013, 2015, 2016 CNRS
  5. *
  6. * StarPU is free software; you can redistribute it and/or modify
  7. * it under the terms of the GNU Lesser General Public License as published by
  8. * the Free Software Foundation; either version 2.1 of the License, or (at
  9. * your option) any later version.
  10. *
  11. * StarPU is distributed in the hope that it will be useful, but
  12. * WITHOUT ANY WARRANTY; without even the implied warranty of
  13. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
  14. *
  15. * See the GNU Lesser General Public License in COPYING.LGPL for more details.
  16. */
  17. #include <starpu_mpi.h>
  18. #include "helper.h"
  19. #ifdef STARPU_QUICK_CHECK
  20. # define NITER 16
  21. #else
  22. # define NITER 2048
  23. #endif
  24. #define SIZE 16
  25. float *tab;
  26. starpu_data_handle_t tab_handle;
  27. int main(int argc, char **argv)
  28. {
  29. int ret, rank, size;
  30. MPI_Init(&argc, &argv);
  31. starpu_mpi_comm_rank(MPI_COMM_WORLD, &rank);
  32. starpu_mpi_comm_size(MPI_COMM_WORLD, &size);
  33. ret = starpu_init(NULL);
  34. STARPU_CHECK_RETURN_VALUE(ret, "starpu_init");
  35. ret = starpu_mpi_init(NULL, NULL, 0);
  36. STARPU_CHECK_RETURN_VALUE(ret, "starpu_mpi_init");
  37. if (size%2 != 0)
  38. {
  39. if (rank == 0)
  40. FPRINTF(stderr, "We need a even number of processes.\n");
  41. starpu_mpi_shutdown();
  42. starpu_shutdown();
  43. MPI_Finalize();
  44. return STARPU_TEST_SKIPPED;
  45. }
  46. tab = calloc(SIZE, sizeof(float));
  47. starpu_vector_data_register(&tab_handle, STARPU_MAIN_RAM, (uintptr_t)tab, SIZE, sizeof(float));
  48. int nloops = NITER;
  49. int loop;
  50. int other_rank = rank%2 == 0 ? rank+1 : rank-1;
  51. for (loop = 0; loop < nloops; loop++)
  52. {
  53. if ((loop % 2) == (rank%2))
  54. {
  55. //FPRINTF_MPI(stderr, "Sending to %d\n", other_rank);
  56. starpu_mpi_send(tab_handle, other_rank, loop, MPI_COMM_WORLD);
  57. }
  58. else
  59. {
  60. MPI_Status status;
  61. //FPRINTF_MPI(stderr, "Receiving from %d\n", other_rank);
  62. starpu_mpi_recv(tab_handle, other_rank, loop, MPI_COMM_WORLD, &status);
  63. }
  64. }
  65. starpu_data_unregister(tab_handle);
  66. free(tab);
  67. starpu_mpi_shutdown();
  68. starpu_shutdown();
  69. MPI_Finalize();
  70. return 0;
  71. }