1ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams/* 2ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * Copyright (C) 2012 The Android Open Source Project 3ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * 4ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * Licensed under the Apache License, Version 2.0 (the "License"); 5ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * you may not use this file except in compliance with the License. 6ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * You may obtain a copy of the License at 7ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * 8ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * http://www.apache.org/licenses/LICENSE-2.0 9ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * 10ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * Unless required by applicable law or agreed to in writing, software 11ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * distributed under the License is distributed on an "AS IS" BASIS, 12ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 13ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * See the License for the specific language governing permissions and 14ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams * limitations under the License. 15ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams */ 16ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 17ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 18ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#include "rsdCore.h" 19ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#include "rsdIntrinsics.h" 20ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#include "rsdAllocation.h" 21ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 22ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#include "rsdIntrinsicInlines.h" 23ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 24ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsusing namespace android; 25ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsusing namespace android::renderscript; 26ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 27ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsstruct ConvolveParams { 28ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams float fp[28]; 29ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams short ip[28]; 30ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams ObjectBaseRef<Allocation> alloc; 31ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams}; 32ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 33ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsstatic void Convolve5x5_Bind(const Context *dc, const Script *script, 34ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams void * intrinsicData, uint32_t slot, Allocation *data) { 35ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams ConvolveParams *cp = (ConvolveParams *)intrinsicData; 36ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams rsAssert(slot == 1); 37ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams cp->alloc.set(data); 38ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams} 39ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 40ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsstatic void Convolve5x5_SetVar(const Context *dc, const Script *script, void * intrinsicData, 41ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t slot, void *data, size_t dataLength) { 42ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams ConvolveParams *cp = (ConvolveParams *)intrinsicData; 43ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 44ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams rsAssert(slot == 0); 45ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams memcpy (cp->fp, data, dataLength); 46ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams for(int ct=0; ct < 25; ct++) { 47ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams cp->ip[ct] = (short)(cp->fp[ct] * 255.f + 0.5f); 48ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams } 49ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams} 50ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 51ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 52ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsstatic void One(const RsForEachStubParamStruct *p, uint32_t x, uchar4 *out, 53ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py0, const uchar4 *py1, const uchar4 *py2, const uchar4 *py3, const uchar4 *py4, 54ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const float* coeff) { 55ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 56ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x0 = rsMax((int32_t)x-2, 0); 57ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x1 = rsMax((int32_t)x-1, 0); 58ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x2 = x; 59ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x3 = rsMin((int32_t)x+1, (int32_t)(p->dimX-1)); 60ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x4 = rsMin((int32_t)x+2, (int32_t)(p->dimX-1)); 61ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 62ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams float4 px = convert_float4(py0[x0]) * coeff[0] + 63ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py0[x1]) * coeff[1] + 64ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py0[x2]) * coeff[2] + 65ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py0[x3]) * coeff[3] + 66ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py0[x4]) * coeff[4] + 67ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 68ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py1[x0]) * coeff[5] + 69ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py1[x1]) * coeff[6] + 70ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py1[x2]) * coeff[7] + 71ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py1[x3]) * coeff[8] + 72ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py1[x4]) * coeff[9] + 73ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 74ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py2[x0]) * coeff[10] + 75ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py2[x1]) * coeff[11] + 76ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py2[x2]) * coeff[12] + 77ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py2[x3]) * coeff[13] + 78ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py2[x4]) * coeff[14] + 79ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 80ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py3[x0]) * coeff[15] + 81ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py3[x1]) * coeff[16] + 82ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py3[x2]) * coeff[17] + 83ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py3[x3]) * coeff[18] + 84ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py3[x4]) * coeff[19] + 85ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 86ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py4[x0]) * coeff[20] + 87ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py4[x1]) * coeff[21] + 88ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py4[x2]) * coeff[22] + 89ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py4[x3]) * coeff[23] + 90ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams convert_float4(py4[x4]) * coeff[24]; 91ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 92ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams px = clamp(px, 0.f, 255.f); 93ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uchar4 o = {(uchar)px.x, (uchar)px.y, (uchar)px.z, (uchar)px.w}; 94ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams *out = o; 95ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams} 96ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 97ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsextern "C" void rsdIntrinsicConvolve5x5_K(void *dst, const void *y0, const void *y1, 98ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const void *y2, const void *y3, const void *y4, 99ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const short *coef, uint32_t count); 100ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 101ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsstatic void Convolve5x5_uchar4(const RsForEachStubParamStruct *p, 102ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t xstart, uint32_t xend, 103ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t instep, uint32_t outstep) { 104ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams ConvolveParams *cp = (ConvolveParams *)p->usr; 105ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams DrvAllocation *din = (DrvAllocation *)cp->alloc->mHal.drv; 106ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar *pin = (const uchar *)din->lod[0].mallocPtr; 107ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 108ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t y0 = rsMax((int32_t)p->y-2, 0); 109ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t y1 = rsMax((int32_t)p->y-1, 0); 110ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t y2 = p->y; 111ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t y3 = rsMin((int32_t)p->y+1, (int32_t)(p->dimY-1)); 112ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t y4 = rsMin((int32_t)p->y+2, (int32_t)(p->dimY-1)); 113ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 114ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py0 = (const uchar4 *)(pin + din->lod[0].stride * y0); 115ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py1 = (const uchar4 *)(pin + din->lod[0].stride * y1); 116ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py2 = (const uchar4 *)(pin + din->lod[0].stride * y2); 117ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py3 = (const uchar4 *)(pin + din->lod[0].stride * y3); 118ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams const uchar4 *py4 = (const uchar4 *)(pin + din->lod[0].stride * y4); 119ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 120ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uchar4 *out = (uchar4 *)p->out; 121ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x1 = xstart; 122ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t x2 = xend; 123ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 124ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams while((x1 < x2) && (x1 < 2)) { 125ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams One(p, x1, out, py0, py1, py2, py3, py4, cp->fp); 126ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams out++; 127ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams x1++; 128ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams } 129ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 130ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#if defined(ARCH_ARM_HAVE_NEON) 131ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams if((x1 + 3) < x2) { 132ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams uint32_t len = (x2 - x1 - 3) >> 1; 133ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams rsdIntrinsicConvolve5x5_K(out, py0, py1, py2, py3, py4, cp->ip, len); 134ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams out += len << 1; 135ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams x1 += len << 1; 136ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams } 137ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams#endif 138ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 139ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams while(x1 < x2) { 140ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams One(p, x1, out, py0, py1, py2, py3, py4, cp->fp); 141ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams out++; 142ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams x1++; 143ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams } 144ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams} 145ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 146ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Samsvoid * rsdIntrinsic_InitConvolve5x5(const android::renderscript::Context *dc, 147ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams android::renderscript::Script *script, 148ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams RsdIntriniscFuncs_t *funcs) { 149ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 150ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams script->mHal.info.exportedVariableCount = 2; 151ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams funcs->bind = Convolve5x5_Bind; 152ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams funcs->setVar = Convolve5x5_SetVar; 153ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams funcs->root = Convolve5x5_uchar4; 154ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 155ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams ConvolveParams *cp = (ConvolveParams *)calloc(1, sizeof(ConvolveParams)); 156ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams for(int ct=0; ct < 25; ct++) { 157ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams cp->fp[ct] = 1.f / 25.f; 158ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams cp->ip[ct] = (short)(cp->fp[ct] * 255.f + 0.5f); 159ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams } 160ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams return cp; 161ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams} 162ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 163ca29b8caf56fa4866752f9cea4ec02b2a271dceeJason Sams 164