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
110
111
112
113
114
|
/**************************************************************************
*
* Copyright 2007 Tungsten Graphics, Inc., Cedar Park, Texas.
* All Rights Reserved.
*
* Permission is hereby granted, free of charge, to any person obtaining a
* copy of this software and associated documentation files (the
* "Software"), to deal in the Software without restriction, including
* without limitation the rights to use, copy, modify, merge, publish,
* distribute, sub license, and/or sell copies of the Software, and to
* permit persons to whom the Software is furnished to do so, subject to
* the following conditions:
*
* The above copyright notice and this permission notice (including the
* next paragraph) shall be included in all copies or substantial portions
* of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
* MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NON-INFRINGEMENT.
* IN NO EVENT SHALL TUNGSTEN GRAPHICS AND/OR ITS SUPPLIERS BE LIABLE FOR
* ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
* TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
* SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
*
**************************************************************************/
/*
* Authors:
* Zack Rusin zack@tungstengraphics.com
*/
typedef __attribute__(( ext_vector_type(4) )) float float4;
extern float powf(float a, float b);
inline float approx(float a, float b)
{
if (b < -128.0f) b = -128.0f;
if (b > 128.0f) b = 128.0f;
if (a < 0) a = 0;
return powf(a, b);
}
inline float4 lit(float4 tmp)
{
float4 result;
result.x = 1.0;
result.w = 1.0;
if (tmp.x > 0) {
result.y = tmp.x;
result.z = approx(tmp.y, tmp.w);
} else {
result.y = 0;
result.z = 0;
}
return result;
}
inline float4 cmp(float4 tmp0, float4 tmp1, float4 tmp2)
{
float4 result;
result.x = (tmp0.x < 0.0) ? tmp1.x : tmp2.x;
result.y = (tmp0.y < 0.0) ? tmp1.y : tmp2.y;
result.z = (tmp0.z < 0.0) ? tmp1.z : tmp2.z;
result.w = (tmp0.w < 0.0) ? tmp1.w : tmp2.w;
return result;
}
extern float cosf(float val);
extern float sinf(float val);
inline float4 vcos(float4 val)
{
float4 result;
printf("VEC IN is %f %f %f %f\n", val.x, val.y, val.z, val.w);
result.x = cosf(val.x);
result.y = cosf(val.x);
result.z = cosf(val.x);
result.w = cosf(val.x);
printf("VEC OUT is %f %f %f %f\n", result.x, result.y, result.z, result.w);
return result;
}
inline float4 scs(float4 val)
{
float4 result;
float tmp = val.x;
result.x = cosf(tmp);
result.y = sinf(tmp);
return result;
}
inline float4 vsin(float4 val)
{
float4 result;
float tmp = val.x;
float res = sinf(tmp);
result.x = res;
result.y = res;
result.z = res;
result.w = res;
return result;
}
inline int kil(float4 val)
{
if (val.x < 0 || val.y < 0 || val.z < 0 || val.w < 0)
return 1;
else
return 0;
}
|
