1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
|
/* StarPU --- Runtime system for heterogeneous multicore architectures.
*
* Copyright (C) 2009, 2010 Université de Bordeaux 1
* Copyright (C) 2010, 2011 Centre National de la Recherche Scientifique
*
* StarPU is free software; you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation; either version 2.1 of the License, or (at
* your option) any later version.
*
* StarPU is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
*
* See the GNU Lesser General Public License in COPYING.LGPL for more details.
*/
#ifndef __DW_CHOLESKY_H__
#define __DW_CHOLESKY_H__
#include <string.h>
#include <math.h>
#include <sys/time.h>
#ifdef STARPU_USE_CUDA
#include <cuda.h>
#include <cuda_runtime.h>
#include <cublas.h>
#endif
#include <common/blas.h>
#include <starpu.h>
#define BLOCKSIZE (size/nblocks)
static unsigned size = 4*1024;
static unsigned nblocks = 16;
static unsigned nbigblocks = 8;
static unsigned noprio = 0;
static unsigned display = 0;
void chol_cpu_codelet_update_u11(void **, void *);
void chol_cpu_codelet_update_u21(void **, void *);
void chol_cpu_codelet_update_u22(void **, void *);
#ifdef STARPU_USE_CUDA
void chol_cublas_codelet_update_u11(void *descr[], void *_args);
void chol_cublas_codelet_update_u21(void *descr[], void *_args);
void chol_cublas_codelet_update_u22(void *descr[], void *_args);
#endif
static void __attribute__((unused)) parse_args(int argc, char **argv)
{
int i;
for (i = 1; i < argc; i++)
{
if (strcmp(argv[i], "-size") == 0)
{
char *argptr;
size = strtol(argv[++i], &argptr, 10);
}
if (strcmp(argv[i], "-nblocks") == 0)
{
char *argptr;
nblocks = strtol(argv[++i], &argptr, 10);
}
if (strcmp(argv[i], "-nbigblocks") == 0)
{
char *argptr;
nbigblocks = strtol(argv[++i], &argptr, 10);
}
if (strcmp(argv[i], "-no-prio") == 0)
{
noprio = 1;
}
if (strcmp(argv[i], "-display") == 0)
{
display = 1;
}
if (strcmp(argv[i], "-h") == 0)
{
printf("usage : %s [-display] [-size size] [-nblocks nblocks]\n", argv[0]);
}
}
if (nblocks > size) nblocks = size;
}
#endif // __DW_CHOLESKY_H__
|
