From d72763ce095f0b879e5baa1713a934225ccb39f7 Mon Sep 17 00:00:00 2001 From: haochen Date: Tue, 13 May 2014 10:27:30 +0800 Subject: GPU Calc: Fix bug For arraybound of slidingwindows parameters in Intercept formula Change-Id: I6c96073e859ef7550033e9a46e6cafb11e8bfe0a --- sc/source/core/opencl/op_statistical.cxx | 267 +++++++++++++++++++++---------- 1 file changed, 179 insertions(+), 88 deletions(-) diff --git a/sc/source/core/opencl/op_statistical.cxx b/sc/source/core/opencl/op_statistical.cxx index 6dfe3444b619..c64588a4fa99 100644 --- a/sc/source/core/opencl/op_statistical.cxx +++ b/sc/source/core/opencl/op_statistical.cxx @@ -5850,121 +5850,212 @@ void OpKurt:: GenSlidingWindowFunction(std::stringstream &ss, void OpIntercept::GenSlidingWindowFunction(std::stringstream &ss, const std::string &sSymName, SubArguments &vSubArguments) { - FormulaToken *pCur = vSubArguments[0]->GetFormulaToken(); - assert(pCur); - const formula::DoubleVectorRefToken* pCurDVR = - static_cast(pCur); - size_t nCurWindowSize = pCurDVR->GetRefRowSize(); - FormulaToken *pCur1 = vSubArguments[1]->GetFormulaToken(); - assert(pCur1); - const formula::DoubleVectorRefToken* pCurDVR1 = - static_cast(pCur1); - size_t nCurWindowSize1 = pCurDVR1->GetRefRowSize(); ss << "\ndouble " << sSymName; - ss << "_"<< BinFuncName() <<"( "; + ss << "_" << BinFuncName() << "("; for (unsigned i = 0; i < vSubArguments.size(); i++) { if (i) ss << ","; vSubArguments[i]->GenSlidingWindowDecl(ss); } - ss << ") {\n"; + ss << "){\n"; ss << " int gid0 = get_global_id(0);\n"; ss << " double fSumX = 0.0;\n"; ss << " double fSumY = 0.0;\n"; + ss << " double fMeanX = 0.0;\n"; + ss << " double fMeanY = 0.0;\n"; ss << " double fSumDeltaXDeltaY = 0.0;\n"; ss << " double fSumSqrDeltaX = 0.0;\n"; - ss << " int length="<GetFormulaToken(); - assert(pCur); - if (pCur->GetType() == formula::svDoubleVectorRef) - { + ss << " return NAN;\n"; + ss << "}\n"; + return ; + } + FormulaToken *pCur = vSubArguments[1]->GetFormulaToken(); + FormulaToken *pCur1 = vSubArguments[0]->GetFormulaToken(); + assert(pCur); + assert(pCur1); + if (pCur->GetType() == formula::svDoubleVectorRef&& + pCur1->GetType() == formula::svDoubleVectorRef) + { + const formula::DoubleVectorRefToken* pDVR = + static_cast(pCur); + const formula::DoubleVectorRefToken* pDVR1 = + static_cast(pCur1); + size_t nCurWindowSize = pDVR->GetRefRowSize(); + size_t nCurWindowSize1 = pDVR1->GetRefRowSize(); + size_t arrayLength = pDVR->GetArrayLength()< + pDVR1->GetArrayLength() ? pDVR->GetArrayLength(): + pDVR1->GetArrayLength(); + if(nCurWindowSize != nCurWindowSize1) + { + ss << " return NAN;\n"; + ss << "}\n"; + return ; } - else if (pCur->GetType() == formula::svSingleVectorRef) + ss << " for (int i = "; + if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed()) + &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed())) { #ifdef ISNAN - ss << "return HUGE_VAL"; + ss << "gid0; i < " << arrayLength; + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "gid0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; #endif } - else if (pCur->GetType() == formula::svDouble) + else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())) { #ifdef ISNAN - ss << "return HUGE_VAL"; + ss << "0; i < " << arrayLength ; + ss << " && i < gid0+" << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < gid0+" << nCurWindowSize << "; i++)\n"; + ss << " {\n"; #endif } - } - ss << " for (int i = 0; i <" << nCurWindowSize << "; i++)\n"; - ss << " {\n"; - ss << " double arg0 = "; - ss << vSubArguments[0]->GenSlidingWindowDeclRef(); - ss << ";\n"; - ss << " double arg1 = "; - ss << vSubArguments[1]->GenSlidingWindowDeclRef(); - ss << ";\n"; -#ifdef ISNAN - ss << " if(isNan(arg0)||((gid0+i)>="; - ss << pCurDVR->GetArrayLength(); - ss << "))\n"; - ss << " {\n"; - ss << " length--;\n"; - ss << " continue;\n"; - ss << " }\n"; + else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())) + { +#ifdef ISNAN + ss << "0; i + gid0 < " << arrayLength; + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; #endif -#ifdef ISNAN - ss << " if(isNan(arg1)||((gid0+i)>="; - ss << pCurDVR1->GetArrayLength(); - ss << "))\n"; - ss << " {\n"; - ss << " length--;\n"; - ss << " continue;\n"; - ss << " }\n"; + } + else if ((pDVR->IsStartFixed() && pDVR->IsEndFixed()) + &&(pDVR1->IsStartFixed() && pDVR1->IsEndFixed())) + { +#ifdef ISNAN + ss << "0; i < " << arrayLength << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << arrayLength << "; i++)\n"; + ss << " {\n"; #endif - ss << " fSumY+=arg0;\n"; - ss << " fSumX+=arg1;\n"; - ss << " }\n"; - ss <<" double fMeanX = fSumX / length;\n"; - ss <<" double fMeanY = fSumY / length;\n"; - ss << " for (int i = 0; i <" << nCurWindowSize << "; i++)\n"; - ss << " {\n"; - ss << " double arg0 = "; - ss << vSubArguments[0]->GenSlidingWindowDeclRef(); - ss << ";\n"; - ss << " double arg1 = "; - ss << vSubArguments[1]->GenSlidingWindowDeclRef(); - ss << ";\n"; -#ifdef ISNAN - ss << " if(isNan(arg0)||((gid0+i)>="; - ss << pCurDVR->GetArrayLength(); - ss << "))\n"; - ss << " {\n"; - ss << " continue;\n"; - ss << " }\n"; + } + else + { +#ifdef ISNAN + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; #endif -#ifdef ISNAN - ss << " if(isNan(arg1)||((gid0+i)>="; - ss <GetArrayLength(); - ss <<"))\n"; - ss << " {\n"; - ss << " continue;\n"; - ss << " }\n"; + ss << " break;\n"; + ss << " }"; + ss << " return NAN;\n"; + ss << "}\n"; + return ; + } + + ss << " argX = "; + ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n"; + ss << " argY = "; + ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n"; +#ifdef ISNAN + ss << " if (isNan(argX) || isNan(argY))\n"; + ss << " continue;\n"; #endif - ss << " fSumDeltaXDeltaY+=(arg1 - fMeanX) * (arg0 - fMeanY);"; - ss << ";\n"; - ss << " fSumSqrDeltaX += pow(arg1 - fMeanX, 2);\n"; - ss << " }\n"; - ss << " tmp = fMeanY - fSumDeltaXDeltaY / fSumSqrDeltaX"; - ss << "* fMeanX;\n"; - ss << " return tmp;\n"; - ss << "}"; + ss << " fSumX += argX;\n"; + ss << " fSumY += argY;\n"; + ss << " fCount += 1.0;\n"; + ss << " }\n"; + + ss << " if (fCount < 1.0)\n"; + ss << " return NAN;\n"; + ss << " else\n"; + ss << " {\n"; + ss << " fMeanX = fSumX * pow(fCount,-1.0);\n"; + ss << " fMeanY = fSumY * pow(fCount,-1.0);\n"; + + ss << " for (int i = "; + if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed()) + &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed())) + { +#ifdef ISNAN + ss << "gid0; i < " << arrayLength; + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "gid0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())) + { +#ifdef ISNAN + ss << "0; i < " << arrayLength ; + ss << " && i < gid0+" << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < gid0+" << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())) + { +#ifdef ISNAN + ss << "0; i + gid0 < " << arrayLength; + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else + { +#ifdef ISNAN + ss << "0; i < " << arrayLength << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << arrayLength << "; i++)\n"; + ss << " {\n"; +#endif + } + + ss << " argX = "; + ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n"; + ss << " argY = "; + ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n"; +#ifdef ISNAN + ss << " if (isNan(argX) || isNan(argY))\n"; + ss << " continue;\n"; +#endif + ss << " fSumDeltaXDeltaY += (argX-fMeanX)*(argY-fMeanY);\n"; + ss << " fSumSqrDeltaX += (argX-fMeanX) * (argX-fMeanX);\n"; + ss << " }\n"; + ss << " if(fSumSqrDeltaX == 0.0)\n"; + ss << " return NAN;\n"; + ss << " else\n"; + ss << " {\n"; + ss << " return fMeanY -"; + ss << " (fSumDeltaXDeltaY*pow(fSumSqrDeltaX,-1.0))*fMeanX;\n"; + ss << " }\n"; + ss << " }\n"; + ss << "}\n"; + } + else + { + ss << " return NAN;\n"; + ss << "}\n"; + } } void OpLogInv:: GenSlidingWindowFunction(std::stringstream &ss, const std::string &sSymName, SubArguments &vSubArguments) -- cgit