diff options
author | Lukas Tönne <lukas.toenne@gmail.com> | 2016-04-20 17:45:29 +0300 |
---|---|---|
committer | Lukas Tönne <lukas.toenne@gmail.com> | 2016-04-20 17:45:29 +0300 |
commit | aae0598aa233d43d0ffc81a070b2ce20d6871c1d (patch) | |
tree | 378100f45cf1db972c671a6886496a5604f7ab29 /intern/cycles/util/util_math_cdf.h | |
parent | 95d7d3c2a6aa91f1dcb33ff0cdfc751e2886d8cd (diff) | |
parent | d7e4f920fd93a4ae5679e0eb6b228a349ab3b082 (diff) |
Merge branch 'master' into temp_depsgraph_split_ubereval
Diffstat (limited to 'intern/cycles/util/util_math_cdf.h')
-rw-r--r-- | intern/cycles/util/util_math_cdf.h | 78 |
1 files changed, 78 insertions, 0 deletions
diff --git a/intern/cycles/util/util_math_cdf.h b/intern/cycles/util/util_math_cdf.h new file mode 100644 index 00000000000..47dfb68ba44 --- /dev/null +++ b/intern/cycles/util/util_math_cdf.h @@ -0,0 +1,78 @@ +/* + * Copyright 2011-2015 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. + */ + +#ifndef __UTIL_MATH_CDF_H__ +#define __UTIL_MATH_CDF_H__ + +#include "util_algorithm.h" +#include "util_math.h" +#include "util_vector.h" + +CCL_NAMESPACE_BEGIN + +/* Evaluate CDF of a given functor with given range and resolution. */ +template <typename Functor> +void util_cdf_evaluate(const int resolution, + const float from, + const float to, + Functor functor, + vector<float> &cdf) +{ + const int cdf_count = resolution + 1; + const float range = to - from; + cdf.resize(cdf_count); + cdf[0] = 0.0f; + /* Actual CDF evaluation. */ + for(int i = 0; i < resolution; ++i) { + float x = from + range * (float)i / (resolution - 1); + float y = functor(x); + cdf[i + 1] = cdf[i] + fabsf(y); + } + /* Normalize the CDF. */ + for(int i = 0; i <= resolution; i++) { + cdf[i] /= cdf[resolution]; + } +} + +/* Invert pre-calculated CDF function. */ +void util_cdf_invert(const int resolution, + const float from, + const float to, + const vector<float> &cdf, + const bool make_symmetric, + vector<float> &inv_cdf); + +/* Evaluate inverted CDF of a given functor with given range and resolution. */ +template <typename Functor> +void util_cdf_inverted(const int resolution, + const float from, + const float to, + Functor functor, + const bool make_symmetric, + vector<float> &inv_cdf) +{ + vector<float> cdf; + /* There is no much smartness going around lower resolution for the CDF table, + * this just to match the old code from pixel filter so it all stays exactly + * the same and no regression tests are failed. + */ + util_cdf_evaluate(resolution - 1, from, to, functor, cdf); + util_cdf_invert(resolution, from, to, cdf, make_symmetric, inv_cdf); +} + +CCL_NAMESPACE_END + +#endif /* __UTIL_MATH_H_CDF__ */ |