mpi_cholesky.h 2.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293
  1. /* StarPU --- Runtime system for heterogeneous multicore architectures.
  2. *
  3. * Copyright (C) 2009, 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. #ifndef __DW_CHOLESKY_H__
  18. #define __DW_CHOLESKY_H__
  19. #include <string.h>
  20. #include <math.h>
  21. #include <sys/time.h>
  22. #ifdef STARPU_USE_CUDA
  23. #include <cuda.h>
  24. #include <cuda_runtime.h>
  25. #include <cublas.h>
  26. #endif
  27. #include <common/blas.h>
  28. #include <starpu.h>
  29. #define BLOCKSIZE (size/nblocks)
  30. static unsigned size = 4*1024;
  31. static unsigned nblocks = 16;
  32. static unsigned nbigblocks = 8;
  33. static unsigned noprio = 0;
  34. static unsigned display = 0;
  35. void chol_cpu_codelet_update_u11(void **, void *);
  36. void chol_cpu_codelet_update_u21(void **, void *);
  37. void chol_cpu_codelet_update_u22(void **, void *);
  38. #ifdef STARPU_USE_CUDA
  39. void chol_cublas_codelet_update_u11(void *descr[], void *_args);
  40. void chol_cublas_codelet_update_u21(void *descr[], void *_args);
  41. void chol_cublas_codelet_update_u22(void *descr[], void *_args);
  42. #endif
  43. static void __attribute__((unused)) parse_args(int argc, char **argv)
  44. {
  45. int i;
  46. for (i = 1; i < argc; i++)
  47. {
  48. if (strcmp(argv[i], "-size") == 0)
  49. {
  50. char *argptr;
  51. size = strtol(argv[++i], &argptr, 10);
  52. }
  53. if (strcmp(argv[i], "-nblocks") == 0)
  54. {
  55. char *argptr;
  56. nblocks = strtol(argv[++i], &argptr, 10);
  57. }
  58. if (strcmp(argv[i], "-nbigblocks") == 0)
  59. {
  60. char *argptr;
  61. nbigblocks = strtol(argv[++i], &argptr, 10);
  62. }
  63. if (strcmp(argv[i], "-no-prio") == 0)
  64. {
  65. noprio = 1;
  66. }
  67. if (strcmp(argv[i], "-display") == 0)
  68. {
  69. display = 1;
  70. }
  71. if (strcmp(argv[i], "-h") == 0)
  72. {
  73. printf("usage : %s [-display] [-size size] [-nblocks nblocks]\n", argv[0]);
  74. }
  75. }
  76. if (nblocks > size) nblocks = size;
  77. }
  78. #endif // __DW_CHOLESKY_H__