cholesky.h 1.9 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879
  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 <semaphore.h>
  20. #include <string.h>
  21. #include <math.h>
  22. #include <sys/time.h>
  23. #ifdef STARPU_USE_CUDA
  24. #include <cuda.h>
  25. #include <cuda_runtime.h>
  26. #include <cublas.h>
  27. #endif
  28. #include <common/blas.h>
  29. #include <starpu.h>
  30. #define NMAXBLOCKS 32
  31. #define BLOCKSIZE (size/nblocks)
  32. static unsigned size = 4*1024;
  33. static unsigned nblocks = 16;
  34. static unsigned nbigblocks = 8;
  35. static unsigned noprio = 0;
  36. static unsigned display = 0;
  37. static void __attribute__((unused)) parse_args(int argc, char **argv)
  38. {
  39. int i;
  40. for (i = 1; i < argc; i++) {
  41. if (strcmp(argv[i], "-size") == 0) {
  42. char *argptr;
  43. size = strtol(argv[++i], &argptr, 10);
  44. }
  45. if (strcmp(argv[i], "-nblocks") == 0) {
  46. char *argptr;
  47. nblocks = strtol(argv[++i], &argptr, 10);
  48. }
  49. if (strcmp(argv[i], "-nbigblocks") == 0) {
  50. char *argptr;
  51. nbigblocks = strtol(argv[++i], &argptr, 10);
  52. }
  53. if (strcmp(argv[i], "-no-prio") == 0) {
  54. noprio = 1;
  55. }
  56. if (strcmp(argv[i], "-display") == 0) {
  57. display = 1;
  58. }
  59. if (strcmp(argv[i], "-h") == 0) {
  60. printf("usage : %s [-display] [-size size] [-nblocks nblocks]\n", argv[0]);
  61. }
  62. }
  63. if (nblocks > size) nblocks = size;
  64. }
  65. #endif // __DW_CHOLESKY_H__