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
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
/*
* -- High Performance Computing Linpack Benchmark (HPL)
* HPL - 2.0 - September 10, 2008
* Antoine P. Petitet
* University of Tennessee, Knoxville
* Innovative Computing Laboratory
* (C) Copyright 2000-2008 All Rights Reserved
*
* -- Copyright notice and Licensing terms:
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* are met:
*
* 1. Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
*
* 2. Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions, and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
*
* 3. All advertising materials mentioning features or use of this
* software must display the following acknowledgement:
* This product includes software developed at the University of
* Tennessee, Knoxville, Innovative Computing Laboratory.
*
* 4. The name of the University, the name of the Laboratory, or the
* names of its contributors may not be used to endorse or promote
* products derived from this software without specific written
* permission.
*
* -- Disclaimer:
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
* ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE UNIVERSITY
* OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
* DATA OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
#ifndef HPL_MISC_H
#define HPL_MISC_H
/*
* ---------------------------------------------------------------------
* Include files
* ---------------------------------------------------------------------
*/
#ifdef __STDC__
#define HPL_STDC_HEADERS
#endif
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#ifdef HPL_STDC_HEADERS
#include <stdarg.h>
#define STDC_ARGS(p) p
#else
#include <varargs.h>
#define STDC_ARGS(p) ()
#endif
#ifdef HPL_CALL_VSIPL
#include <vsip.h>
#endif
/*
* ---------------------------------------------------------------------
* #define macro constants
* ---------------------------------------------------------------------
*/
#define HPL_rone 1.0
#define HPL_rtwo 2.0
#define HPL_rzero 0.0
/*
* ---------------------------------------------------------------------
* #define macros definitions
* ---------------------------------------------------------------------
*/
#define Mabs( a_ ) ( ( (a_) < 0 ) ? -(a_) : (a_) )
#define Mmin( a_, b_ ) ( ( (a_) < (b_) ) ? (a_) : (b_) )
#define Mmax( a_, b_ ) ( ( (a_) > (b_) ) ? (a_) : (b_) )
#define Mfloor(a,b) (((a)>0) ? (((a)/(b))) : (-(((-(a))+(b)-1)/(b))))
#define Mceil(a,b) ( ( (a)+(b)-1 ) / (b) )
#define Miceil(a,b) (((a)>0) ? ((((a)+(b)-1)/(b))) : (-((-(a))/(b))))
#define Mupcase(C) (((C)>96 && (C)<123) ? (C) & 0xDF : (C))
#define Mlowcase(C) (((C)>64 && (C)< 91) ? (C) | 32 : (C))
/*
* Mptr returns a pointer to a_( i_, j_ ) for readability reasons and
* also less silly errors ...
*/
#define Mptr( a_, i_, j_, lda_ ) \
( (a_) + (size_t)(i_) + (size_t)(j_)*(size_t)(lda_) )
/*
* Align pointer
*/
#define HPL_PTR( ptr_, al_ ) \
( ( ( (size_t)(ptr_)+(al_)-1 ) / (al_) ) * (al_) )
#endif
/*
* End of hpl_misc.h
*/
|
