blob: 373699604621b8d7106d2bcda24b11de703f980d [file] [log] [blame]
joshualittac977922014-07-22 09:52:11 -07001/*
2 * Copyright 2014 Google Inc.
3 *
4 * Use of this source code is governed by a BSD-style license that can be
5 * found in the LICENSE file.
6 */
Mike Kleinc0bd9f92019-04-23 12:05:21 -05007#include "src/gpu/effects/GrMatrixConvolutionEffect.h"
Robert Phillips296b1cc2017-03-15 10:42:12 -04008
Adlai Holler00ddb002020-05-11 19:25:52 -04009#include "include/private/SkHalf.h"
10#include "src/gpu/GrBitmapTextureMaker.h"
11#include "src/gpu/GrContextPriv.h"
12#include "src/gpu/GrProxyProvider.h"
13#include "src/gpu/GrRecordingContextPriv.h"
Greg Daniel456f9b52020-03-05 19:14:18 +000014#include "src/gpu/GrTexture.h"
Greg Danielf91aeb22019-06-18 09:58:02 -040015#include "src/gpu/GrTextureProxy.h"
Brian Salomon694ec492020-04-14 13:39:31 -040016#include "src/gpu/effects/GrTextureEffect.h"
Mike Kleinc0bd9f92019-04-23 12:05:21 -050017#include "src/gpu/glsl/GrGLSLFragmentProcessor.h"
18#include "src/gpu/glsl/GrGLSLFragmentShaderBuilder.h"
19#include "src/gpu/glsl/GrGLSLProgramDataManager.h"
20#include "src/gpu/glsl/GrGLSLUniformHandler.h"
joshualittac977922014-07-22 09:52:11 -070021
egdaniel64c47282015-11-13 06:54:19 -080022class GrGLMatrixConvolutionEffect : public GrGLSLFragmentProcessor {
joshualittac977922014-07-22 09:52:11 -070023public:
robertphillips9cdb9922016-02-03 12:25:40 -080024 void emitCode(EmitArgs&) override;
joshualittac977922014-07-22 09:52:11 -070025
Brian Salomon94efbf52016-11-29 13:43:05 -050026 static inline void GenKey(const GrProcessor&, const GrShaderCaps&, GrProcessorKeyBuilder*);
joshualittac977922014-07-22 09:52:11 -070027
wangyixb1daa862015-08-18 11:29:31 -070028protected:
Brian Salomonab015ef2017-04-04 10:15:51 -040029 void onSetData(const GrGLSLProgramDataManager&, const GrFragmentProcessor&) override;
joshualittac977922014-07-22 09:52:11 -070030
31private:
egdaniel018fb622015-10-28 07:26:40 -070032 typedef GrGLSLProgramDataManager::UniformHandle UniformHandle;
joshualittac977922014-07-22 09:52:11 -070033
Adlai Holler00ddb002020-05-11 19:25:52 -040034 void emitKernelBlock(EmitArgs&, SkIPoint);
35
joshualitt5ae5fc52014-07-29 12:59:27 -070036 UniformHandle fKernelUni;
joshualitt5ae5fc52014-07-29 12:59:27 -070037 UniformHandle fKernelOffsetUni;
38 UniformHandle fGainUni;
39 UniformHandle fBiasUni;
Adlai Holler00ddb002020-05-11 19:25:52 -040040 UniformHandle fKernelBiasUni;
joshualittac977922014-07-22 09:52:11 -070041
egdaniel64c47282015-11-13 06:54:19 -080042 typedef GrGLSLFragmentProcessor INHERITED;
joshualittac977922014-07-22 09:52:11 -070043};
44
Adlai Holler00ddb002020-05-11 19:25:52 -040045GrMatrixConvolutionEffect::KernelWrapper GrMatrixConvolutionEffect::KernelWrapper::Make(
46 GrRecordingContext* context, SkISize size, const GrCaps& caps, const SkScalar* values) {
47 if (nullptr == context || nullptr == values || size.isEmpty()) {
48 return {};
49 }
50 const int length = size.area();
51 // Small kernel -> just fill the array.
52 KernelWrapper result(size);
53 if (length <= kMaxUniformSize) {
54 for (int i = 0; i < length; i++) {
55 result.fArray[i] = SkScalarToFloat(values[i]);
56 }
57 return result;
58 }
59
60 ScalableSampler& scalableSampler = result.fScalableSampler;
61 bool useA16 =
62 context->defaultBackendFormat(kA16_float_SkColorType, GrRenderable::kNo).isValid();
63 SkScalar min = values[0];
64 if (!useA16) {
65 // Determine min and max values to figure out inner gain & bias.
66 SkScalar max = values[0];
67 for (int i = 1; i < length; i++) {
68 if (values[i] < min) {
69 min = values[i];
70 }
71 if (values[i] > max) {
72 max = values[i];
73 }
74 }
75 // Treat near-0 gain (i.e. box blur) as 1, and let the kernelBias
76 // move everything up to the final value.
77 const SkScalar computedGain = max - min;
78 scalableSampler.fGain =
79 SkScalarNearlyZero(computedGain) ? 1.0f : SkScalarToFloat(computedGain);
80 // Inner bias is pre-inner-gain so we divide that out.
81 scalableSampler.fBias = SkScalarToFloat(min) / scalableSampler.fGain;
82 }
83
84 // TODO: Enable kernel caching and check perf.
85 static constexpr bool kCacheKernelTexture = false;
86
87 GrUniqueKey key;
88 if (kCacheKernelTexture) {
89 static const GrUniqueKey::Domain kDomain = GrUniqueKey::GenerateDomain();
90 GrUniqueKey::Builder builder(&key, kDomain, length, "Matrix Convolution Kernel");
91 // Texture cache key is the exact content of the kernel.
92 static_assert(sizeof(float) == 4);
93 for (int i = 0; i < length; i++) {
94 builder[i] = *(const uint32_t*)&values[i];
95 }
96 builder.finish();
97 }
98
99 // Find or create a texture.
100 GrProxyProvider* proxyProvider = context->priv().proxyProvider();
101 GrSurfaceProxyView view;
102 SkColorType colorType = useA16 ? kA16_float_SkColorType : kAlpha_8_SkColorType;
103 sk_sp<GrTextureProxy> cachedKernel;
104 if (kCacheKernelTexture && (cachedKernel = proxyProvider->findOrCreateProxyByUniqueKey(key))) {
105 GrSwizzle swizzle =
106 context->priv().caps()->getReadSwizzle(cachedKernel->backendFormat(),
107 SkColorTypeToGrColorType(colorType));
108 view = {std::move(cachedKernel), kTopLeft_GrSurfaceOrigin, swizzle};
109 } else {
110 SkBitmap bm;
111 auto info = SkImageInfo::Make({(int)GrNextPow2(length), 1}, colorType,
112 kPremul_SkAlphaType, nullptr);
113 if (!bm.tryAllocPixels(info)) {
114 return {};
115 }
116 for (int i = 0; i < length; i++) {
117 if (useA16) {
118 *bm.getAddr16(i, 0) = SkFloatToHalf(values[i]);
119 } else {
120 *bm.getAddr8(i, 0) =
121 SkScalarRoundToInt((values[i] - min) / scalableSampler.fGain * 255);
122 }
123 }
124 bm.setImmutable();
125 GrBitmapTextureMaker maker(context, bm, GrImageTexGenPolicy::kNew_Uncached_Budgeted);
126 view = maker.view(GrMipMapped::kNo);
127 if (!view) {
128 return {};
129 }
130 if (kCacheKernelTexture) {
131 proxyProvider->assignUniqueKeyToProxy(key, view.asTextureProxy());
132 }
133 }
134 scalableSampler.fSampler = { std::move(view) };
135 return result;
136}
137
138bool GrMatrixConvolutionEffect::KernelWrapper::operator==(const KernelWrapper& k) const {
139 if (fSize != k.fSize) {
140 return false;
141 } else if (this->isSampled()) {
142 return fScalableSampler == k.fScalableSampler;
143 } else {
144 return std::equal(fArray.begin(), fArray.begin() + fSize.area(), k.fArray.begin());
145 }
146}
147
148bool GrMatrixConvolutionEffect::KernelWrapper::ScalableSampler::operator==(
149 const ScalableSampler& k) const {
150 return fSampler == k.fSampler && fGain == k.fGain && fBias == k.fBias;
151}
152
153// For sampled kernels, emit a for loop that does all the kernel accumulation.
154// For uniform kernels, emit a single iteration. Function is called repeatedly in a for loop.
155// loc is ignored for sampled kernels.
156void GrGLMatrixConvolutionEffect::emitKernelBlock(EmitArgs& args, SkIPoint loc) {
157 const GrMatrixConvolutionEffect& mce = args.fFp.cast<GrMatrixConvolutionEffect>();
158 GrGLSLFPFragmentBuilder* fragBuilder = args.fFragBuilder;
159 GrGLSLUniformHandler* uniformHandler = args.fUniformHandler;
160 int kernelWidth = mce.kernelSize().width();
161 int kernelHeight = mce.kernelSize().height();
162 int kernelArea = kernelWidth * kernelHeight;
163
164 if (mce.kernelIsSampled()) {
165 fragBuilder->codeAppendf("half2 kernelCoord = half2(0, 0);");
166 fragBuilder->codeAppendf("for (int i = 0; i < %d; ++i)", (int)kernelArea);
167 }
168
169 GrGLSLShaderBuilder::ShaderBlock block(fragBuilder);
170
171 fragBuilder->codeAppend("half k;");
172 fragBuilder->codeAppend("half2 sourceOffset;");
173 if (mce.kernelIsSampled()) {
174 const char* kernelBias = uniformHandler->getUniformCStr(fKernelBiasUni);
175 fragBuilder->codeAppend("k = ");
176 fragBuilder->appendTextureLookup(args.fTexSamplers[0], "kernelCoord");
177 fragBuilder->codeAppendf(".w + %s;", kernelBias);
178 fragBuilder->codeAppendf("sourceOffset.y = floor(i / %d);", kernelWidth);
179 fragBuilder->codeAppendf("sourceOffset.x = i - sourceOffset.y * %d;", kernelWidth);
180 float kernelStride = 1.0f / (float)GrNextPow2(kernelArea);
181 fragBuilder->codeAppendf("kernelCoord.x += %f;", kernelStride);
182 } else {
183 fragBuilder->codeAppendf("sourceOffset = half2(%d, %d);", loc.x(), loc.y());
184 int offset = loc.y() * kernelWidth + loc.x();
185 static constexpr const char kVecSuffix[][4] = { ".x", ".y", ".z", ".w" };
186 const char* kernel = uniformHandler->getUniformCStr(fKernelUni);
187 fragBuilder->codeAppendf("k = %s[%d]%s;", kernel, offset / 4,
188 kVecSuffix[offset & 0x3]);
189 }
190
191 auto sample = this->invokeChild(0, args, "coord + sourceOffset");
192 fragBuilder->codeAppendf("half4 c = %s;", sample.c_str());
193 if (!mce.convolveAlpha()) {
194 fragBuilder->codeAppend("c.rgb /= max(c.a, 0.0001);");
195 fragBuilder->codeAppend("c.rgb = saturate(c.rgb);");
196 }
197 fragBuilder->codeAppend("sum += c * k;");
198}
199
wangyix7c157a92015-07-22 15:08:53 -0700200void GrGLMatrixConvolutionEffect::emitCode(EmitArgs& args) {
robertphillipsbf536af2016-02-04 06:11:53 -0800201 const GrMatrixConvolutionEffect& mce = args.fFp.cast<GrMatrixConvolutionEffect>();
robertphillipsbf536af2016-02-04 06:11:53 -0800202
Adlai Holler00ddb002020-05-11 19:25:52 -0400203 int kernelWidth = mce.kernelSize().width();
204 int kernelHeight = mce.kernelSize().height();
robertphillipsbf536af2016-02-04 06:11:53 -0800205
Adlai Holler00ddb002020-05-11 19:25:52 -0400206 int arrayCount = (kernelWidth * kernelHeight + 3) / 4;
207 SkASSERT(4 * arrayCount >= kernelWidth * kernelHeight);
jvanverth78d6eb02016-03-02 13:21:16 -0800208
egdaniel7ea439b2015-12-03 09:20:44 -0800209 GrGLSLUniformHandler* uniformHandler = args.fUniformHandler;
Adlai Holler00ddb002020-05-11 19:25:52 -0400210 if (mce.kernelIsSampled()) {
211 fKernelBiasUni = uniformHandler->addUniform(&mce, kFragment_GrShaderFlag,
212 kHalf_GrSLType, "KernelBias");
213 } else {
214 fKernelUni = uniformHandler->addUniformArray(&mce, kFragment_GrShaderFlag,
215 kHalf4_GrSLType, "Kernel", arrayCount);
216 }
Ethan Nicholas16464c32020-04-06 13:53:05 -0400217 fKernelOffsetUni = uniformHandler->addUniform(&mce, kFragment_GrShaderFlag, kHalf2_GrSLType,
egdaniel7ea439b2015-12-03 09:20:44 -0800218 "KernelOffset");
Brian Salomon41e377d2020-05-05 01:27:57 +0000219 fGainUni = uniformHandler->addUniform(&mce, kFragment_GrShaderFlag, kHalf_GrSLType, "Gain");
220 fBiasUni = uniformHandler->addUniform(&mce, kFragment_GrShaderFlag, kHalf_GrSLType, "Bias");
joshualittac977922014-07-22 09:52:11 -0700221
egdaniel7ea439b2015-12-03 09:20:44 -0800222 const char* kernelOffset = uniformHandler->getUniformCStr(fKernelOffsetUni);
egdaniel7ea439b2015-12-03 09:20:44 -0800223 const char* gain = uniformHandler->getUniformCStr(fGainUni);
224 const char* bias = uniformHandler->getUniformCStr(fBiasUni);
joshualittac977922014-07-22 09:52:11 -0700225
cdalton85285412016-02-18 12:37:07 -0800226 GrGLSLFPFragmentBuilder* fragBuilder = args.fFragBuilder;
Ethan Nicholas58430122020-04-14 09:54:02 -0400227 SkString coords2D = fragBuilder->ensureCoords2D(args.fTransformedCoords[0].fVaryingPoint,
228 mce.sampleMatrix());
Brian Salomon41e377d2020-05-05 01:27:57 +0000229 fragBuilder->codeAppend("half4 sum = half4(0, 0, 0, 0);");
Brian Salomon694ec492020-04-14 13:39:31 -0400230 fragBuilder->codeAppendf("float2 coord = %s - %s;", coords2D.c_str(), kernelOffset);
joshualitt5ae5fc52014-07-29 12:59:27 -0700231
Adlai Holler00ddb002020-05-11 19:25:52 -0400232 if (mce.kernelIsSampled()) {
233 this->emitKernelBlock(args, {});
234 } else {
235 for (int x = 0; x < kernelWidth; ++x) {
236 for (int y = 0; y < kernelHeight; ++y) {
237 this->emitKernelBlock(args, SkIPoint::Make(x, y));
joshualitt5ae5fc52014-07-29 12:59:27 -0700238 }
joshualitt5ae5fc52014-07-29 12:59:27 -0700239 }
joshualittac977922014-07-22 09:52:11 -0700240 }
Adlai Holler00ddb002020-05-11 19:25:52 -0400241
robertphillipsbf536af2016-02-04 06:11:53 -0800242 if (mce.convolveAlpha()) {
Brian Salomon41e377d2020-05-05 01:27:57 +0000243 fragBuilder->codeAppendf("%s = sum * %s + %s;", args.fOutputColor, gain, bias);
Ethan Nicholas12fb9cf2018-08-03 16:16:57 -0400244 fragBuilder->codeAppendf("%s.a = saturate(%s.a);", args.fOutputColor, args.fOutputColor);
egdaniel4ca2e602015-11-18 08:01:26 -0800245 fragBuilder->codeAppendf("%s.rgb = clamp(%s.rgb, 0.0, %s.a);",
246 args.fOutputColor, args.fOutputColor, args.fOutputColor);
joshualittac977922014-07-22 09:52:11 -0700247 } else {
Brian Salomon694ec492020-04-14 13:39:31 -0400248 auto sample = this->invokeChild(0, args, coords2D.c_str());
Adlai Holler00ddb002020-05-11 19:25:52 -0400249 fragBuilder->codeAppendf("half4 c = %s;", sample.c_str());
Brian Salomon41e377d2020-05-05 01:27:57 +0000250 fragBuilder->codeAppendf("%s.a = c.a;", args.fOutputColor);
251 fragBuilder->codeAppendf("%s.rgb = saturate(sum.rgb * %s + %s);", args.fOutputColor, gain, bias);
egdaniel4ca2e602015-11-18 08:01:26 -0800252 fragBuilder->codeAppendf("%s.rgb *= %s.a;", args.fOutputColor, args.fOutputColor);
joshualittac977922014-07-22 09:52:11 -0700253 }
Ethan Nicholas2983f402017-05-08 09:36:08 -0400254 fragBuilder->codeAppendf("%s *= %s;\n", args.fOutputColor, args.fInputColor);
joshualittac977922014-07-22 09:52:11 -0700255}
256
joshualittb0a8a372014-09-23 09:50:21 -0700257void GrGLMatrixConvolutionEffect::GenKey(const GrProcessor& processor,
Brian Salomon94efbf52016-11-29 13:43:05 -0500258 const GrShaderCaps&, GrProcessorKeyBuilder* b) {
joshualittb0a8a372014-09-23 09:50:21 -0700259 const GrMatrixConvolutionEffect& m = processor.cast<GrMatrixConvolutionEffect>();
joshualitt5ae5fc52014-07-29 12:59:27 -0700260 SkASSERT(m.kernelSize().width() <= 0x7FFF && m.kernelSize().height() <= 0xFFFF);
261 uint32_t key = m.kernelSize().width() << 16 | m.kernelSize().height();
caryclark952538e2016-02-26 05:01:42 -0800262 key |= m.convolveAlpha() ? 1U << 31 : 0;
joshualittac977922014-07-22 09:52:11 -0700263 b->add32(key);
264}
265
egdaniel018fb622015-10-28 07:26:40 -0700266void GrGLMatrixConvolutionEffect::onSetData(const GrGLSLProgramDataManager& pdman,
Brian Salomonab015ef2017-04-04 10:15:51 -0400267 const GrFragmentProcessor& processor) {
joshualittb0a8a372014-09-23 09:50:21 -0700268 const GrMatrixConvolutionEffect& conv = processor.cast<GrMatrixConvolutionEffect>();
Brian Salomon694ec492020-04-14 13:39:31 -0400269 pdman.set2fv(fKernelOffsetUni, 1, conv.kernelOffset().ptr());
Adlai Holler00ddb002020-05-11 19:25:52 -0400270 float totalGain = conv.gain();
271 if (conv.kernelIsSampled()) {
272 totalGain *= conv.kernelSampleGain();
273 pdman.set1f(fKernelBiasUni, conv.kernelSampleBias());
274 } else {
275 int kernelCount = conv.kernelSize().area();
276 int arrayCount = (kernelCount + 3) / 4;
277 SkASSERT(4 * arrayCount >= kernelCount);
278 pdman.set4fv(fKernelUni, arrayCount, conv.kernel());
279 }
kkinnunen7510b222014-07-30 00:04:16 -0700280 pdman.set1f(fBiasUni, conv.bias());
Adlai Holler00ddb002020-05-11 19:25:52 -0400281 pdman.set1f(fGainUni, totalGain);
joshualittac977922014-07-22 09:52:11 -0700282}
283
Brian Salomon694ec492020-04-14 13:39:31 -0400284GrMatrixConvolutionEffect::GrMatrixConvolutionEffect(std::unique_ptr<GrFragmentProcessor> child,
Adlai Holler00ddb002020-05-11 19:25:52 -0400285 KernelWrapper kernel,
Robert Phillips40fd7c92017-01-30 08:06:27 -0500286 SkScalar gain,
287 SkScalar bias,
288 const SkIPoint& kernelOffset,
Robert Phillips40fd7c92017-01-30 08:06:27 -0500289 bool convolveAlpha)
Brian Salomon6cd51b52017-07-26 19:07:15 -0400290 // To advertise either the modulation or opaqueness optimizations we'd have to examine the
291 // parameters.
Ethan Nicholasabff9562017-10-09 10:54:08 -0400292 : INHERITED(kGrMatrixConvolutionEffect_ClassID, kNone_OptimizationFlags)
Adlai Holler00ddb002020-05-11 19:25:52 -0400293 , fKernel(std::move(kernel))
Brian Salomon6cd51b52017-07-26 19:07:15 -0400294 , fGain(SkScalarToFloat(gain))
295 , fBias(SkScalarToFloat(bias) / 255.0f)
296 , fConvolveAlpha(convolveAlpha) {
Brian Salomon694ec492020-04-14 13:39:31 -0400297 child->setSampledWithExplicitCoords();
298 this->registerChildProcessor(std::move(child));
Adlai Holler00ddb002020-05-11 19:25:52 -0400299 if (fKernel.isSampled()) {
300 this->setTextureSamplerCnt(1);
Robert Phillips40fd7c92017-01-30 08:06:27 -0500301 }
Brian Salomon694ec492020-04-14 13:39:31 -0400302 fKernelOffset = {static_cast<float>(kernelOffset.x()),
303 static_cast<float>(kernelOffset.y())};
304 this->addCoordTransform(&fCoordTransform);
Robert Phillips40fd7c92017-01-30 08:06:27 -0500305}
306
Brian Salomon3f6f9652017-07-28 07:34:05 -0400307GrMatrixConvolutionEffect::GrMatrixConvolutionEffect(const GrMatrixConvolutionEffect& that)
Ethan Nicholasabff9562017-10-09 10:54:08 -0400308 : INHERITED(kGrMatrixConvolutionEffect_ClassID, kNone_OptimizationFlags)
Adlai Holler00ddb002020-05-11 19:25:52 -0400309 , fKernel(that.fKernel)
Brian Salomon3f6f9652017-07-28 07:34:05 -0400310 , fGain(that.fGain)
311 , fBias(that.fBias)
Brian Salomon694ec492020-04-14 13:39:31 -0400312 , fKernelOffset(that.fKernelOffset)
Brian Salomon3f6f9652017-07-28 07:34:05 -0400313 , fConvolveAlpha(that.fConvolveAlpha) {
Brian Salomon694ec492020-04-14 13:39:31 -0400314 auto child = that.childProcessor(0).clone();
315 child->setSampledWithExplicitCoords();
316 this->registerChildProcessor(std::move(child));
Adlai Holler00ddb002020-05-11 19:25:52 -0400317 if (fKernel.isSampled()) {
318 this->setTextureSamplerCnt(1);
319 }
Brian Salomon3f6f9652017-07-28 07:34:05 -0400320 this->addCoordTransform(&fCoordTransform);
Brian Salomon3f6f9652017-07-28 07:34:05 -0400321}
322
Brian Salomonaff329b2017-08-11 09:40:37 -0400323std::unique_ptr<GrFragmentProcessor> GrMatrixConvolutionEffect::clone() const {
324 return std::unique_ptr<GrFragmentProcessor>(new GrMatrixConvolutionEffect(*this));
Brian Salomon3f6f9652017-07-28 07:34:05 -0400325}
326
Brian Salomon94efbf52016-11-29 13:43:05 -0500327void GrMatrixConvolutionEffect::onGetGLSLProcessorKey(const GrShaderCaps& caps,
egdaniel57d3b032015-11-13 11:57:27 -0800328 GrProcessorKeyBuilder* b) const {
joshualitteb2a6762014-12-04 11:35:33 -0800329 GrGLMatrixConvolutionEffect::GenKey(*this, caps, b);
330}
331
egdaniel57d3b032015-11-13 11:57:27 -0800332GrGLSLFragmentProcessor* GrMatrixConvolutionEffect::onCreateGLSLInstance() const {
robertphillipsbf536af2016-02-04 06:11:53 -0800333 return new GrGLMatrixConvolutionEffect;
joshualittac977922014-07-22 09:52:11 -0700334}
335
bsalomon0e08fc12014-10-15 08:19:04 -0700336bool GrMatrixConvolutionEffect::onIsEqual(const GrFragmentProcessor& sBase) const {
joshualitt49586be2014-09-16 08:21:41 -0700337 const GrMatrixConvolutionEffect& s = sBase.cast<GrMatrixConvolutionEffect>();
Adlai Holler00ddb002020-05-11 19:25:52 -0400338 return fKernel == s.fKernel &&
joshualittac977922014-07-22 09:52:11 -0700339 fGain == s.gain() &&
340 fBias == s.bias() &&
Brian Salomon694ec492020-04-14 13:39:31 -0400341 fKernelOffset == s.kernelOffset() &&
342 fConvolveAlpha == s.convolveAlpha();
joshualittac977922014-07-22 09:52:11 -0700343}
344
Adlai Holler00ddb002020-05-11 19:25:52 -0400345const GrFragmentProcessor::TextureSampler& GrMatrixConvolutionEffect::onTextureSampler(
346 int index) const {
347 return IthTextureSampler(index, fKernel.scalableSampler().fSampler);
348}
349
Robert Phillips4e962c62018-06-20 07:46:19 -0400350static void fill_in_1D_gaussian_kernel_with_stride(float* kernel, int size, int stride,
351 float twoSigmaSqrd) {
352 SkASSERT(!SkScalarNearlyZero(twoSigmaSqrd, SK_ScalarNearlyZero));
353
354 const float sigmaDenom = 1.0f / twoSigmaSqrd;
355 const int radius = size / 2;
356
357 float sum = 0.0f;
358 for (int i = 0; i < size; ++i) {
359 float term = static_cast<float>(i - radius);
360 // Note that the constant term (1/(sqrt(2*pi*sigma^2)) of the Gaussian
361 // is dropped here, since we renormalize the kernel below.
362 kernel[i * stride] = sk_float_exp(-term * term * sigmaDenom);
363 sum += kernel[i * stride];
364 }
365 // Normalize the kernel
366 float scale = 1.0f / sum;
367 for (int i = 0; i < size; ++i) {
368 kernel[i * stride] *= scale;
369 }
370}
371
Robert Phillips40fd7c92017-01-30 08:06:27 -0500372static void fill_in_2D_gaussian_kernel(float* kernel, int width, int height,
373 SkScalar sigmaX, SkScalar sigmaY) {
Greg Daniel4eda8d92018-04-03 14:03:15 -0400374 const float twoSigmaSqrdX = 2.0f * SkScalarToFloat(SkScalarSquare(sigmaX));
375 const float twoSigmaSqrdY = 2.0f * SkScalarToFloat(SkScalarSquare(sigmaY));
376
Robert Phillips4e962c62018-06-20 07:46:19 -0400377 // TODO: in all of these degenerate cases we're uploading (and using) a whole lot of zeros.
Greg Daniel4eda8d92018-04-03 14:03:15 -0400378 if (SkScalarNearlyZero(twoSigmaSqrdX, SK_ScalarNearlyZero) ||
379 SkScalarNearlyZero(twoSigmaSqrdY, SK_ScalarNearlyZero)) {
Robert Phillips4e962c62018-06-20 07:46:19 -0400380 // In this case the 2D Gaussian degenerates to a 1D Gaussian (in X or Y) or a point
381 SkASSERT(3 == width || 3 == height);
Brian Salomon694ec492020-04-14 13:39:31 -0400382 std::fill_n(kernel, width*height, 0);
Robert Phillips4e962c62018-06-20 07:46:19 -0400383
384 if (SkScalarNearlyZero(twoSigmaSqrdX, SK_ScalarNearlyZero) &&
385 SkScalarNearlyZero(twoSigmaSqrdY, SK_ScalarNearlyZero)) {
386 // A point
387 SkASSERT(3 == width && 3 == height);
388 kernel[4] = 1.0f;
389 } else if (SkScalarNearlyZero(twoSigmaSqrdX, SK_ScalarNearlyZero)) {
390 // A 1D Gaussian in Y
391 SkASSERT(3 == width);
392 // Down the middle column of the kernel with a stride of width
393 fill_in_1D_gaussian_kernel_with_stride(&kernel[1], height, width, twoSigmaSqrdY);
394 } else {
395 // A 1D Gaussian in X
396 SkASSERT(SkScalarNearlyZero(twoSigmaSqrdY, SK_ScalarNearlyZero));
397 SkASSERT(3 == height);
398 // Down the middle row of the kernel with a stride of 1
399 fill_in_1D_gaussian_kernel_with_stride(&kernel[width], width, 1, twoSigmaSqrdX);
Greg Daniel4eda8d92018-04-03 14:03:15 -0400400 }
401 return;
402 }
403
404 const float sigmaXDenom = 1.0f / twoSigmaSqrdX;
405 const float sigmaYDenom = 1.0f / twoSigmaSqrdY;
Robert Phillips40fd7c92017-01-30 08:06:27 -0500406 const int xRadius = width / 2;
407 const int yRadius = height / 2;
408
joshualitt5acfea72014-08-11 13:55:34 -0700409 float sum = 0.0f;
joshualitt5acfea72014-08-11 13:55:34 -0700410 for (int x = 0; x < width; x++) {
Robert Phillips40fd7c92017-01-30 08:06:27 -0500411 float xTerm = static_cast<float>(x - xRadius);
412 xTerm = xTerm * xTerm * sigmaXDenom;
413 for (int y = 0; y < height; y++) {
414 float yTerm = static_cast<float>(y - yRadius);
415 float xyTerm = sk_float_exp(-(xTerm + yTerm * yTerm * sigmaYDenom));
416 // Note that the constant term (1/(sqrt(2*pi*sigma^2)) of the Gaussian
417 // is dropped here, since we renormalize the kernel below.
418 kernel[y * width + x] = xyTerm;
419 sum += xyTerm;
420 }
joshualitt5acfea72014-08-11 13:55:34 -0700421 }
422 // Normalize the kernel
423 float scale = 1.0f / sum;
424 for (int i = 0; i < width * height; ++i) {
425 kernel[i] *= scale;
426 }
Robert Phillips40fd7c92017-01-30 08:06:27 -0500427}
428
Adlai Holler00ddb002020-05-11 19:25:52 -0400429std::unique_ptr<GrFragmentProcessor> GrMatrixConvolutionEffect::Make(GrRecordingContext* context,
430 GrSurfaceProxyView srcView,
Brian Salomon694ec492020-04-14 13:39:31 -0400431 const SkIRect& srcBounds,
432 const SkISize& kernelSize,
433 const SkScalar* kernel,
434 SkScalar gain,
435 SkScalar bias,
436 const SkIPoint& kernelOffset,
437 GrSamplerState::WrapMode wm,
438 bool convolveAlpha,
439 const GrCaps& caps) {
Adlai Holler00ddb002020-05-11 19:25:52 -0400440 auto kw = KernelWrapper::Make(context, kernelSize, caps, kernel);
441 if (!kw.isValid()) {
442 return nullptr;
443 }
Brian Salomon694ec492020-04-14 13:39:31 -0400444 GrSamplerState sampler(wm, GrSamplerState::Filter::kNearest);
445 auto child = GrTextureEffect::MakeSubset(std::move(srcView), kPremul_SkAlphaType, SkMatrix::I(),
446 sampler, SkRect::Make(srcBounds), caps);
447 return std::unique_ptr<GrFragmentProcessor>(new GrMatrixConvolutionEffect(
Adlai Holler00ddb002020-05-11 19:25:52 -0400448 std::move(child), std::move(kw), gain, bias, kernelOffset, convolveAlpha));
Brian Salomon694ec492020-04-14 13:39:31 -0400449}
450
Brian Salomonaff329b2017-08-11 09:40:37 -0400451std::unique_ptr<GrFragmentProcessor> GrMatrixConvolutionEffect::MakeGaussian(
Adlai Holler00ddb002020-05-11 19:25:52 -0400452 GrRecordingContext* context,
Greg Daniel5c082492020-01-29 15:06:49 -0500453 GrSurfaceProxyView srcView,
Robert Phillips5140f9a2018-05-11 16:11:45 -0400454 const SkIRect& srcBounds,
Brian Salomonaff329b2017-08-11 09:40:37 -0400455 const SkISize& kernelSize,
456 SkScalar gain,
457 SkScalar bias,
458 const SkIPoint& kernelOffset,
Brian Salomon694ec492020-04-14 13:39:31 -0400459 GrSamplerState::WrapMode wm,
Brian Salomonaff329b2017-08-11 09:40:37 -0400460 bool convolveAlpha,
461 SkScalar sigmaX,
Brian Salomon694ec492020-04-14 13:39:31 -0400462 SkScalar sigmaY,
463 const GrCaps& caps) {
Adlai Holler00ddb002020-05-11 19:25:52 -0400464 SkAutoSTMalloc<32, float> kernel(kernelSize.area());
465 fill_in_2D_gaussian_kernel(kernel.get(), kernelSize.width(), kernelSize.height(),
466 sigmaX, sigmaY);
467 return Make(context, std::move(srcView), srcBounds, kernelSize, kernel.get(),
468 gain, bias, kernelOffset, wm, convolveAlpha, caps);
Robert Phillips40fd7c92017-01-30 08:06:27 -0500469}
470
joshualittb0a8a372014-09-23 09:50:21 -0700471GR_DEFINE_FRAGMENT_PROCESSOR_TEST(GrMatrixConvolutionEffect);
joshualittac977922014-07-22 09:52:11 -0700472
Hal Canary6f6961e2017-01-31 13:50:44 -0500473#if GR_TEST_UTILS
Brian Salomonaff329b2017-08-11 09:40:37 -0400474std::unique_ptr<GrFragmentProcessor> GrMatrixConvolutionEffect::TestCreate(GrProcessorTestData* d) {
Greg Daniel026a60c2020-02-12 10:53:51 -0500475 auto [view, ct, at] = d->randomView();
Robert Phillipsdbc8eeb2017-02-21 10:04:31 -0500476
Adlai Holler00ddb002020-05-11 19:25:52 -0400477 static constexpr size_t kMaxTestKernelSize = 2 * kMaxUniformSize;
478 int width = d->fRandom->nextRangeU(1, kMaxTestKernelSize);
479 int height = d->fRandom->nextRangeU(1, kMaxTestKernelSize / width);
joshualittac977922014-07-22 09:52:11 -0700480 SkISize kernelSize = SkISize::Make(width, height);
Ben Wagner7ecc5962016-11-02 17:07:33 -0400481 std::unique_ptr<SkScalar[]> kernel(new SkScalar[width * height]);
joshualittac977922014-07-22 09:52:11 -0700482 for (int i = 0; i < width * height; i++) {
joshualitt0067ff52015-07-08 14:26:19 -0700483 kernel.get()[i] = d->fRandom->nextSScalar1();
joshualittac977922014-07-22 09:52:11 -0700484 }
joshualitt0067ff52015-07-08 14:26:19 -0700485 SkScalar gain = d->fRandom->nextSScalar1();
486 SkScalar bias = d->fRandom->nextSScalar1();
Greg Daniel026a60c2020-02-12 10:53:51 -0500487
488 uint32_t kernalOffsetX = d->fRandom->nextRangeU(0, kernelSize.width());
489 uint32_t kernalOffsetY = d->fRandom->nextRangeU(0, kernelSize.height());
490 SkIPoint kernelOffset = SkIPoint::Make(kernalOffsetX, kernalOffsetY);
491
492 uint32_t boundsX = d->fRandom->nextRangeU(0, view.width());
493 uint32_t boundsY = d->fRandom->nextRangeU(0, view.height());
494 uint32_t boundsW = d->fRandom->nextRangeU(0, view.width());
495 uint32_t boundsH = d->fRandom->nextRangeU(0, view.height());
496 SkIRect bounds = SkIRect::MakeXYWH(boundsX, boundsY, boundsW, boundsH);
497
Brian Salomon694ec492020-04-14 13:39:31 -0400498 auto wm = static_cast<GrSamplerState::WrapMode>(
499 d->fRandom->nextULessThan(GrSamplerState::kWrapModeCount));
joshualitt0067ff52015-07-08 14:26:19 -0700500 bool convolveAlpha = d->fRandom->nextBool();
Adlai Holler00ddb002020-05-11 19:25:52 -0400501 return GrMatrixConvolutionEffect::Make(d->context()->priv().asRecordingContext(),
502 std::move(view),
bungeman06ca8ec2016-06-09 08:01:03 -0700503 bounds,
504 kernelSize,
505 kernel.get(),
506 gain,
507 bias,
508 kernelOffset,
Brian Salomon694ec492020-04-14 13:39:31 -0400509 wm,
510 convolveAlpha,
511 *d->caps());
joshualittac977922014-07-22 09:52:11 -0700512}
Hal Canary6f6961e2017-01-31 13:50:44 -0500513#endif