subgraph_repeat.c 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148
  1. /* StarPU --- Runtime system for heterogeneous multicore architectures.
  2. *
  3. * Copyright (C) 2010 Université de Bordeaux 1
  4. * Copyright (C) 2010, 2011 Centre National de la Recherche Scientifique
  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 <sys/time.h>
  18. #include <starpu.h>
  19. #include <pthread.h>
  20. #include "../common/helper.h"
  21. static unsigned niter = 16384;
  22. /*
  23. *
  24. * /-->B--\
  25. * | |
  26. * -----> A D---\--->
  27. * ^ | | |
  28. * | \-->C--/ |
  29. * | |
  30. * \--------------/
  31. *
  32. * - {B, C} depend on A
  33. * - D depends on {B, C}
  34. * - A, B, C and D are resubmitted at the end of the loop (or not)
  35. */
  36. static struct starpu_task taskA, taskB, taskC, taskD;
  37. static unsigned loop_cnt = 0;
  38. static unsigned check_cnt = 0;
  39. static pthread_cond_t cond = PTHREAD_COND_INITIALIZER;
  40. static pthread_mutex_t mutex = PTHREAD_MUTEX_INITIALIZER;
  41. static void dummy_func(void *descr[] __attribute__ ((unused)), void *arg __attribute__ ((unused)))
  42. {
  43. STARPU_ATOMIC_ADD(&check_cnt, 1);
  44. }
  45. static starpu_codelet dummy_codelet =
  46. {
  47. .where = STARPU_CPU|STARPU_CUDA|STARPU_OPENCL,
  48. .cpu_func = dummy_func,
  49. .cuda_func = dummy_func,
  50. .opencl_func = dummy_func,
  51. .model = NULL,
  52. .nbuffers = 0
  53. };
  54. static void callback_task_D(void *arg __attribute__((unused)))
  55. {
  56. loop_cnt++;
  57. if (loop_cnt == niter)
  58. {
  59. /* We are done */
  60. pthread_mutex_lock(&mutex);
  61. pthread_cond_signal(&cond);
  62. pthread_mutex_unlock(&mutex);
  63. }
  64. else {
  65. /* Let's go for another iteration */
  66. starpu_task_submit(&taskA);
  67. starpu_task_submit(&taskB);
  68. starpu_task_submit(&taskC);
  69. starpu_task_submit(&taskD);
  70. }
  71. }
  72. int main(int argc, char **argv)
  73. {
  74. // unsigned i;
  75. // double timing;
  76. // struct timeval start;
  77. // struct timeval end;
  78. int ret;
  79. ret = starpu_init(NULL);
  80. STARPU_CHECK_RETURN_VALUE(ret, "starpu_init");
  81. starpu_task_init(&taskA);
  82. taskA.cl = &dummy_codelet;
  83. taskA.cl_arg = &taskA;
  84. starpu_task_init(&taskB);
  85. taskB.cl = &dummy_codelet;
  86. taskB.cl_arg = &taskB;
  87. starpu_task_init(&taskC);
  88. taskC.cl = &dummy_codelet;
  89. taskC.cl_arg = &taskC;
  90. starpu_task_init(&taskD);
  91. taskD.cl = &dummy_codelet;
  92. taskD.cl_arg = &taskD;
  93. taskD.callback_func = callback_task_D;
  94. struct starpu_task *depsBC_array[1] = {&taskA};
  95. starpu_task_declare_deps_array(&taskB, 1, depsBC_array);
  96. starpu_task_declare_deps_array(&taskC, 1, depsBC_array);
  97. struct starpu_task *depsD_array[2] = {&taskB, &taskC};
  98. starpu_task_declare_deps_array(&taskD, 2, depsD_array);
  99. ret = starpu_task_submit(&taskA); if (ret == -ENODEV) goto enodev; STARPU_CHECK_RETURN_VALUE(ret, "starpu_task_submit");
  100. ret = starpu_task_submit(&taskB); if (ret == -ENODEV) goto enodev; STARPU_CHECK_RETURN_VALUE(ret, "starpu_task_submit");
  101. ret = starpu_task_submit(&taskC); if (ret == -ENODEV) goto enodev; STARPU_CHECK_RETURN_VALUE(ret, "starpu_task_submit");
  102. ret = starpu_task_submit(&taskD); if (ret == -ENODEV) goto enodev; STARPU_CHECK_RETURN_VALUE(ret, "starpu_task_submit");
  103. /* Wait for the termination of all loops */
  104. pthread_mutex_lock(&mutex);
  105. if (loop_cnt < niter)
  106. pthread_cond_wait(&cond, &mutex);
  107. pthread_mutex_unlock(&mutex);
  108. STARPU_ASSERT(check_cnt == (4*loop_cnt));
  109. /* Cleanup the statically allocated tasks */
  110. starpu_task_deinit(&taskA);
  111. starpu_task_deinit(&taskB);
  112. starpu_task_deinit(&taskC);
  113. starpu_task_deinit(&taskD);
  114. starpu_shutdown();
  115. return 0;
  116. enodev:
  117. fprintf(stderr, "WARNING: No one can execute this task\n");
  118. /* yes, we do not perform the computation but we did detect that no one
  119. * could perform the kernel, so this is not an error from StarPU */
  120. starpu_shutdown();
  121. return 77;
  122. }