diff options
author | shiming zhang <shiming@multicorewareinc.com> | 2013-12-16 11:39:14 +0800 |
---|---|---|
committer | I-Jui (Ray) Sung <ray@multicorewareinc.com> | 2013-12-19 17:54:00 -0600 |
commit | f034d1a2f532f43f8abb4b020d6e4e91f19ffc8c (patch) | |
tree | 3f4c7d007af11c7965633ce974e138193e904746 | |
parent | 49235d00f6f19bd622eba3be6f11878ac870ceca (diff) |
GPU Calc: implemented VARPA
AMLOEXT-219 FIX
Change-Id: Id88635589ea23a9b44b027d5037b33febbd82cab
Signed-off-by: haochen <haochen@multicorewareinc.com>
Signed-off-by: Wei Wei <weiwei@multicorewareinc.com>
Signed-off-by: I-Jui (Ray) Sung <ray@multicorewareinc.com>
-rw-r--r-- | sc/source/core/opencl/formulagroupcl.cxx | 4 | ||||
-rw-r--r-- | sc/source/core/opencl/op_statistical.cxx | 421 | ||||
-rw-r--r-- | sc/source/core/opencl/op_statistical.hxx | 9 |
3 files changed, 434 insertions, 0 deletions
diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx index 1f0b3c452069..f2dea68168d0 100644 --- a/sc/source/core/opencl/formulagroupcl.cxx +++ b/sc/source/core/opencl/formulagroupcl.cxx @@ -2565,6 +2565,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments( mvSubArguments.push_back(SoPHelper(ts, ft->Children[i], new OpVarA)); break; + case ocVarPA: + mvSubArguments.push_back(SoPHelper(ts, + ft->Children[i], new OpVarPA)); + break; case ocExternal: if ( !(pChild->GetExternal().compareTo(OUString( "com.sun.star.sheet.addin.Analysis.getEffect")))) diff --git a/sc/source/core/opencl/op_statistical.cxx b/sc/source/core/opencl/op_statistical.cxx index bd8ac70ad023..3068d8a089b2 100644 --- a/sc/source/core/opencl/op_statistical.cxx +++ b/sc/source/core/opencl/op_statistical.cxx @@ -9984,6 +9984,427 @@ void OpVarA::GenSlidingWindowFunction(std::stringstream &ss, ss << " return vSum * pow(fCount - 1.0,-1.0);\n"; ss << "}\n"; } +void OpVarPA::GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments) +{ + int isMixedDV = 0; + int isMixedSV = 0; + ss << "\ndouble " << sSymName; + ss << "_" << BinFuncName() << "("; + for (unsigned i = 0; i < vSubArguments.size(); i++) + { + if (i) + ss << ","; + vSubArguments[i]->GenSlidingWindowDecl(ss); + } + ss << "){\n"; + ss << " int gid0 = get_global_id(0);\n"; + ss << " double fSum = 0.0;\n"; + ss << " double fMean = 0.0;\n"; + ss << " double vSum = 0.0;\n"; + ss << " double fCount = 0.0;\n"; + ss << " double arg = 0.0;\n"; + unsigned i = vSubArguments.size(); + while (i--) + { + FormulaToken *pCur = vSubArguments[i]->GetFormulaToken(); + assert(pCur); + + if(ocPush == vSubArguments[i]->GetFormulaToken()->GetOpCode()) + { + if (pCur->GetType() == formula::svDoubleVectorRef) + { + const formula::DoubleVectorRefToken* pDVR = + dynamic_cast<const formula::DoubleVectorRefToken *>(pCur); + if(pDVR->GetArrays()[0].mpNumericArray + && pDVR->GetArrays()[0].mpStringArray) + isMixedDV = svDoubleVectorRefDoubleString; + else if(pDVR->GetArrays()[0].mpNumericArray) + isMixedDV = svDoubleVectorRefDouble; + else if(pDVR->GetArrays()[0].mpStringArray) + isMixedDV = svDoubleVectorRefString; + else + isMixedDV = svDoubleVectorRefNULL; + + size_t nCurWindowSize = pDVR->GetRefRowSize(); + ss << " for (int i = "; + if (!pDVR->IsStartFixed() && pDVR->IsEndFixed()) + { +#ifdef ISNAN + ss << "gid0; i < " << pDVR->GetArrayLength(); + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "gid0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + { +#ifdef ISNAN + ss << "0; i < " << pDVR->GetArrayLength(); + 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()) + { +#ifdef ISNAN + ss << "0; i + gid0 < " << pDVR->GetArrayLength(); + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else + { +#ifdef ISNAN + ss << "0; i < " << pDVR->GetArrayLength() << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << pDVR->GetArrayLength() << "; i++)\n"; + ss << " {\n"; +#endif + } + + if(isMixedDV == svDoubleVectorRefDoubleString) + { + ss << " arg = "; + ss << vSubArguments[i]->GenDoubleSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " == 0)\n"; + ss << " continue;\n"; + ss << " if(isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " {\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " continue;\n"; + ss << " }\n"; +#endif + ss << " fSum += arg;\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " }\n"; + } + else if(isMixedDV == svDoubleVectorRefDouble) + { + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (isNan(arg))\n"; + ss << " continue;\n"; +#endif + ss << " fSum += arg;\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " }\n"; + } + else if(isMixedDV == svDoubleVectorRefString) + { + ss << " if ("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << " == 0)\n"; + ss << " continue;\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " }\n"; + } + else + { + ss << " continue;\n"; + ss << " }\n"; + } + + } + else if (pCur->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken* pSVR = + dynamic_cast< const formula::SingleVectorRefToken* >(pCur); + assert(pSVR); + if(pSVR->GetArray().mpNumericArray + && pSVR->GetArray().mpStringArray) + isMixedSV = svSingleVectorRefDoubleString; + else if(pSVR->GetArray().mpNumericArray) + isMixedSV = svSingleVectorRefDouble; + else if(pSVR->GetArray().mpStringArray) + isMixedSV = svSingleVectorRefString; + else + isMixedSV = svSingleVectorRefNULL; + + if(isMixedSV == svSingleVectorRefDoubleString) + { +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; +#endif + ss << " arg = "; + ss << vSubArguments[i]->GenDoubleSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (!isNan(arg))\n"; + ss << " {\n"; +#endif + ss << " fSum += arg;\n"; + ss << " fCount = fCount + 1.0;\n"; +#ifdef ISNAN + ss << " }\n"; + ss << " if (isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " }\n"; +#endif + } + else if(isMixedSV == svSingleVectorRefDouble) + { +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; +#endif + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; +#ifdef ISNAN + ss << " if (!isNan(arg))\n"; + ss << " {\n"; +#endif + ss << " fSum += arg;\n"; + ss << " fCount += 1.0;\n"; +#ifdef ISNAN + ss << " }\n"; + ss << " }\n"; +#endif + } + else if(isMixedSV == svSingleVectorRefString) + { + + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " if ("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " fCount = fCount + 1.0;\n"; + ss << " }\n"; + } + else + { + ss << " arg =0.0;\n"; + } + } + else + { + ss << " arg = " << pCur->GetDouble() << ";\n"; + ss << " fSum += arg;\n"; + ss << " fCount = fCount + 1.0;\n"; + } + } + else + { + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; + ss << " fSum += arg;\n"; + ss << " fCount = fCount + 1.0;\n"; + } + if (i == 0) + { + ss << " fMean = fSum * pow(fCount,-1.0);\n"; + } + } + i = vSubArguments.size(); + while (i--) + { + FormulaToken *pCur = vSubArguments[i]->GetFormulaToken(); + assert(pCur); + + if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode()) + { + if (pCur->GetType() == formula::svDoubleVectorRef) + { + const formula::DoubleVectorRefToken* pDVR = + dynamic_cast<const formula::DoubleVectorRefToken *>(pCur); + size_t nCurWindowSize = pDVR->GetRefRowSize(); + ss << " for (int i = "; + if (!pDVR->IsStartFixed() && pDVR->IsEndFixed()) + { +#ifdef ISNAN + ss << "gid0; i < " << pDVR->GetArrayLength(); + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "gid0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed()) + { +#ifdef ISNAN + ss << "0; i < " << pDVR->GetArrayLength(); + 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()) + { +#ifdef ISNAN + ss << "0; i + gid0 < " << pDVR->GetArrayLength(); + ss << " && i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + else + { +#ifdef ISNAN + ss << "0; i < " << pDVR->GetArrayLength() << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < " << pDVR->GetArrayLength() << "; i++)\n"; + ss << " {\n"; +#endif + } + if(isMixedDV == svDoubleVectorRefDoubleString) + { + ss << " arg = "; + ss << vSubArguments[i]->GenDoubleSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " == 0)\n"; + ss << " continue;\n"; + ss << " if(isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " arg = 0.0;\n"; +#endif + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; + ss << " }\n"; + + } + else if(isMixedDV == svDoubleVectorRefDouble) + { + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (isNan(arg))\n"; + ss << " continue;\n"; +#endif + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; + ss << " }\n"; + + } + else if(isMixedDV == svDoubleVectorRefString) + { + ss << " if ("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << " == 0)\n"; + ss << " continue;\n"; + ss << " arg = 0.0;\n"; + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; + ss << " }\n"; + } + else + { + ss << " continue;\n"; + ss << " }\n"; + } + } + else if (pCur->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken* pSVR = + dynamic_cast< const formula::SingleVectorRefToken* >(pCur); + assert(pSVR); + if(isMixedSV == svSingleVectorRefDoubleString) + { +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; +#endif + ss << " arg = "; + ss << vSubArguments[i]->GenDoubleSlidingWindowDeclRef(); + ss << ";\n"; +#ifdef ISNAN + ss << " if (!isNan(arg))\n"; +#endif + ss << " vSum += (arg - fMean)*(arg - fMean);\n"; +#ifdef ISNAN + ss << " if (isNan(arg) && "; + ss << vSubArguments[i]->GenStringSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " {\n"; + ss << " arg = 0.0;\n"; + ss << " vSum += (arg - fMean)*(arg - fMean);\n"; + ss << " }\n"; + ss << " }\n"; +#endif + } + else if(isMixedSV == svSingleVectorRefDouble) + { +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; +#endif + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; +#ifdef ISNAN + ss << " if (!isNan(arg))\n"; + ss << " {\n"; +#endif + ss << " vSum += (arg - fMean)*(arg - fMean);\n"; +#ifdef ISNAN + ss << " }\n"; + ss << " }\n"; +#endif + } + else if(isMixedSV == svSingleVectorRefString) + { + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " if ("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << " != 0)\n"; + ss << " {\n"; + ss << " arg = 0.0;\n"; + ss << " vSum += (arg - fMean)*(arg - fMean);\n"; + ss << " }\n"; + ss << " }\n"; + } + else + { + ss << " arg = 0.0;\n"; + } + } + else + { + ss << " arg = " << pCur->GetDouble() << ";\n"; + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; + } + } + else + { + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; + } + } + ss << " if (fCount == 0.0)\n"; + ss << " return DBL_MAX;\n"; + ss << " else\n"; + ss << " return vSum * pow(fCount,-1.0);\n"; + ss << "}\n"; +} }} /* vim:set shiftwidth=4 softtabstop=4 expandtab: */ diff --git a/sc/source/core/opencl/op_statistical.hxx b/sc/source/core/opencl/op_statistical.hxx index 45cbf6798c33..3e3f1961d6ef 100644 --- a/sc/source/core/opencl/op_statistical.hxx +++ b/sc/source/core/opencl/op_statistical.hxx @@ -496,6 +496,15 @@ public: virtual bool takeString() const { return true; } virtual bool takeNumeric() const { return true; } }; +class OpVarPA: public Normal +{ +public: + virtual void GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments); + virtual std::string BinFuncName(void) const { return "OpVarPA"; } + virtual bool takeString() const { return true; } + virtual bool takeNumeric() const { return true; } +}; class OpAverageA: public Normal { public: |