mpi_cholesky.h 2.3 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586
  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. if (strcmp(argv[i], "-size") == 0) {
  48. char *argptr;
  49. size = strtol(argv[++i], &argptr, 10);
  50. }
  51. if (strcmp(argv[i], "-nblocks") == 0) {
  52. char *argptr;
  53. nblocks = strtol(argv[++i], &argptr, 10);
  54. }
  55. if (strcmp(argv[i], "-nbigblocks") == 0) {
  56. char *argptr;
  57. nbigblocks = strtol(argv[++i], &argptr, 10);
  58. }
  59. if (strcmp(argv[i], "-no-prio") == 0) {
  60. noprio = 1;
  61. }
  62. if (strcmp(argv[i], "-display") == 0) {
  63. display = 1;
  64. }
  65. if (strcmp(argv[i], "-h") == 0) {
  66. printf("usage : %s [-display] [-size size] [-nblocks nblocks]\n", argv[0]);
  67. }
  68. }
  69. if (nblocks > size) nblocks = size;
  70. }
  71. #endif // __DW_CHOLESKY_H__