Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Annotation of /branches/vis12/src/include/Diderot/diderot.h
ViewVC logotype

Annotation of /branches/vis12/src/include/Diderot/diderot.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1796 - (view) (download) (as text)

1 : jhr 1115 /*! \file diderot.h
2 :     *
3 :     * \author John Reppy
4 :     *
5 :     * This is the interface to the Diderot runtime for the C target. For now,
6 :     * we are targetting single-precision computations.
7 :     */
8 :    
9 :     /*
10 :     * COPYRIGHT (c) 2010 The Diderot Project (http://diderot-language.cs.uchicago.edu)
11 :     * All rights reserved.
12 :     */
13 :    
14 :     #ifndef _DIDEROT_H_
15 :     #define _DIDEROT_H_
16 :    
17 : jhr 1232 #include "Diderot/config.h"
18 :    
19 : jhr 1301 // include OpenCL headers when the target is OpenCL
20 :     #if defined(DIDEROT_TARGET_CL)
21 :     # ifdef HAVE_CL_CL_H
22 :     # include <CL/cl.h>
23 :     # elif defined(HAVE_OPENCL_CL_H)
24 :     # include <OpenCL/cl.h>
25 :     # else
26 :     # error no cl.h
27 :     # endif
28 :     #endif
29 :    
30 : jhr 1115 /* gcc has a bug that breaks code that uses 8-byte vectors (e.g., vec2f_t), but clang
31 :     * handles the code correctly.
32 :     */
33 :     #if defined(__clang__)
34 :     # define VEC2_OK
35 :     #endif
36 :    
37 :     #ifdef NDEBUG
38 : jhr 1295 #define STATIC_INLINE static inline
39 : jhr 1115 #else
40 : jhr 1295 #define STATIC_INLINE static
41 : jhr 1115 #endif
42 :    
43 : jhr 1640 #include <assert.h>
44 : jhr 1115 #include <stdint.h>
45 :     #include <stdbool.h>
46 :     #include <stdlib.h>
47 :     #include <math.h>
48 :     #include <stdio.h> // for printing output
49 : jhr 1640 #include <string.h>
50 :     #include "teem/air.h"
51 : jhr 1716 #include "teem/nrrd.h"
52 : jhr 1115
53 :     #include "types.h"
54 : jhr 1640 #include "world.h"
55 : jhr 1115 #include "strands.h"
56 : jhr 1301 #include "options.h"
57 : jhr 1716 //#include "output.h"
58 : jhr 1115
59 :     /* load image data from Nrrd files */
60 :     extern Status_t Diderot_LoadImage1D (Diderot_string_t name, Diderot_image1D_t **img);
61 :     extern Status_t Diderot_LoadImage2D (Diderot_string_t name, Diderot_image2D_t **img);
62 :     extern Status_t Diderot_LoadImage3D (Diderot_string_t name, Diderot_image3D_t **img);
63 :    
64 : jhr 1380 #ifdef DIDEROT_TARGET_CL
65 :     #include "shadow-types.h"
66 :     #endif
67 : jhr 1232
68 : jhr 1295 //! Summary information about the CPU configuration.
69 : jhr 1232 typedef struct {
70 : jhr 1295 int numHWNodes; //!< \brief number of (possibly multicore) processors
71 :     int numHWCores; //!< \brief total number of (possibly
72 :     //! mulithreaded) cores
73 :     int numHWThreads; //!< \brief total number of hardware threads
74 :     int numCoresPerNode; //!< \brief number of cores per thread
75 :     int numThdsPerCore; //!< \brief number of threads per core
76 : jhr 1232 } CPUInfo_t;
77 :    
78 : jhr 1295 //! \brief function to get information about the CPU configuration
79 : jhr 1232 bool GetNumCPUs (CPUInfo_t *info);
80 :    
81 : jhr 1640 //! Checked memory allocation
82 :     STATIC_INLINE void *CheckedAlloc (size_t szb)
83 :     {
84 :     void *p = malloc(szb);
85 :     if (p == 0) {
86 :     fprintf (stderr, "fatal error: unable to allocate %d bytes of memory\n", (int)szb);
87 :     exit (1);
88 :     }
89 :     return p;
90 :     }
91 : jhr 1232
92 : jhr 1640 #define NEW(ty) (ty *)CheckedAlloc(sizeof(ty))
93 :     #define NEWVEC(ty,n) (ty *)CheckedAlloc(sizeof(ty) * (n))
94 :     #define NEWSTR(s) strcpy((char *)CheckedAlloc(strlen(s)+1), s)
95 : jhr 1718 #define FREE(p) do { if (p != 0) free(p); } while(0)
96 : jhr 1301
97 : jhr 1691 /* dynamic sequences */
98 :     #include "dyn-seq.h"
99 : jhr 1640
100 : jhr 1691
101 : jhr 1115 /********** scalar math functions **********/
102 :    
103 : jhr 1640 STATIC_INLINE Diderot_real_t min (Diderot_real_t a, Diderot_real_t b)
104 : jhr 1115 {
105 : jhr 1295 return (a < b) ? a : b;
106 : jhr 1115 }
107 :    
108 : jhr 1640 STATIC_INLINE Diderot_real_t max (Diderot_real_t a, Diderot_real_t b)
109 : jhr 1115 {
110 : jhr 1295 return (a < b) ? b : a;
111 : jhr 1115 }
112 :    
113 : jhr 1640 STATIC_INLINE Diderot_real_t clamp (Diderot_real_t lo, Diderot_real_t hi, Diderot_real_t x)
114 : jhr 1295 {
115 : jhr 1640 return min(max(lo, x), hi);
116 : jhr 1295 }
117 :    
118 : jhr 1640 STATIC_INLINE Diderot_real_t lerp (Diderot_real_t a, Diderot_real_t b, Diderot_real_t t)
119 : jhr 1115 {
120 :     return a + t*(b - a);
121 :     }
122 :    
123 : jhr 1640 #if defined(DIDEROT_SINGLE_PRECISION)
124 :     # define EPSILON 1.0E-12
125 :     # define ABS(x) fabsf(x)
126 :     # define ATAN2(x,y) atan2f(x,y)
127 :     # define CBRT(x) cbrtf(x)
128 :     # define COS(x) cosf(x)
129 :     # define FLOOR(x) floorf(x)
130 :     # define SIN(x) sinf(x)
131 :     # define SQRT(x) sqrtf(x)
132 :     # define TRUNC(x) truncf(x)
133 :     #elif defined(DIDEROT_DOUBLE_PRECISION)
134 :     # define EPSILON 1.0E-12
135 :     # define ABS(x) fabs(x)
136 :     # define ATAN2(x,y) atan2(x,y)
137 :     # define CBRT(x) cbrt(x)
138 :     # define COS(x) cos(x)
139 :     # define FLOOR(x) floor(x)
140 :     # define SIN(x) sin(x)
141 :     # define SQRT(x) sqrt(x)
142 :     # define TRUNC(x) trunc(x)
143 :     #else
144 :     # error floating-point precision unknown
145 :     #endif
146 :    
147 : jhr 1115 /********** vector math functions **********/
148 :    
149 :     #include "inline-vec2.h"
150 :     #include "inline-vec3.h"
151 :     #include "inline-vec4.h"
152 :     #include "inline-matrix.h"
153 : jhr 1796 #include "inline-tensor3.h"
154 : jhr 1115
155 :     /********** other Diderot support functions **********/
156 :    
157 :     #include "inline-image.h"
158 :    
159 : jhr 1640 /********** other math functions **********/
160 :     int Diderot_evals2x2 (
161 :     Diderot_real_t eval[2],
162 :     const Diderot_real_t _M00, const Diderot_real_t _M01,
163 :     const Diderot_real_t _M11);
164 :     int Diderot_evecs2x2 (
165 :     Diderot_real_t eval[2], Diderot_vec2_t evec[2],
166 :     const Diderot_real_t _M00, const Diderot_real_t _M01,
167 :     const Diderot_real_t _M11);
168 :     int Diderot_evals3x3 (
169 :     Diderot_real_t eval[3],
170 :     const Diderot_real_t _M00, const Diderot_real_t _M01, const Diderot_real_t _M02,
171 :     const Diderot_real_t _M11, const Diderot_real_t _M12,
172 :     const Diderot_real_t _M22);
173 :     int Diderot_evecs3x3 (
174 :     Diderot_real_t eval[3], Diderot_vec3_t evecs[3],
175 :     const Diderot_real_t _M00, const Diderot_real_t _M01, const Diderot_real_t _M02,
176 :     const Diderot_real_t _M11, const Diderot_real_t _M12,
177 :     const Diderot_real_t _M22);
178 :    
179 : jhr 1115 #endif /* !_DIDEROT_H_ */

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0