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
|
/*
* Copyright 2011-2013 Blender Foundation
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#include "stdosl.h"
float safe_divide(float a, float b)
{
return (b != 0.0) ? a / b : 0.0;
}
float safe_modulo(float a, float b)
{
return (b != 0.0) ? fmod(a, b) : 0.0;
}
float safe_sqrt(float a)
{
return (a > 0.0) ? sqrt(a) : 0.0;
}
float safe_log(float a, float b)
{
return (a > 0.0 && b > 0.0) ? log(a) / log(b) : 0.0;
}
shader node_math(string type = "add",
float A = 0.5,
float B = 0.5,
output float Result = 0.0)
{
/* OSL asin, acos, pow check for values that could give rise to nan. */
if (type == "add")
Result = A + B;
else if (type == "subtract")
Result = A - B;
else if (type == "multiply")
Result = A * B;
else if (type == "divide")
Result = safe_divide(A, B);
else if (type == "power")
Result = pow(A, B);
else if (type == "logarithm")
Result = safe_log(A, B);
else if (type == "sqrt")
Result = safe_sqrt(A);
else if (type == "absolute")
Result = fabs(A);
else if (type == "minimum")
Result = min(A, B);
else if (type == "maximum")
Result = max(A, B);
else if (type == "less_than")
Result = A < B;
else if (type == "greater_than")
Result = A > B;
else if (type == "round")
Result = floor(A + 0.5);
else if (type == "floor")
Result = floor(A);
else if (type == "ceil")
Result = ceil(A);
else if (type == "fraction")
Result = A - floor(A);
else if (type == "modulo")
Result = safe_modulo(A, B);
else if (type == "sine")
Result = sin(A);
else if (type == "cosine")
Result = cos(A);
else if (type == "tangent")
Result = tan(A);
else if (type == "arcsine")
Result = asin(A);
else if (type == "arccosine")
Result = acos(A);
else if (type == "arctangent")
Result = atan(A);
else if (type == "arctan2")
Result = atan2(A, B);
else
Result = 0.0;
}
|