From e8a755e782de1a843a4874d067e38522e1a5911b Mon Sep 17 00:00:00 2001 From: shiming zhang Date: Thu, 7 Nov 2013 17:11:47 +0800 Subject: GPU Calc: implemented VAR AMLOEXT-79 FIX Change-Id: I048d02b752b9d995f424ca63ca481fe57ab9a648 Signed-off-by: haochen Signed-off-by: I-Jui (Ray) Sung --- sc/source/core/opencl/formulagroupcl.cxx | 4 + sc/source/core/opencl/op_statistical.cxx | 232 +++++++++++++++++++++++++++++++ sc/source/core/opencl/op_statistical.hxx | 9 +- 3 files changed, 243 insertions(+), 2 deletions(-) diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx index 8337e58fe5c8..2c07ab38462c 100644 --- a/sc/source/core/opencl/formulagroupcl.cxx +++ b/sc/source/core/opencl/formulagroupcl.cxx @@ -1212,6 +1212,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments( mvSubArguments.push_back(SoPHelper(ts, ft->Children[i], new OpFdist)); break; + case ocVar: + mvSubArguments.push_back(SoPHelper(ts, + ft->Children[i], new OpVar)); + 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 3ab6ea93f65b..4c2f272672cb 100644 --- a/sc/source/core/opencl/op_statistical.cxx +++ b/sc/source/core/opencl/op_statistical.cxx @@ -22,6 +22,238 @@ using namespace formula; namespace sc { namespace opencl { +void OpVar::GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments) +{ + 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(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 < "<< nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < "<< nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + 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 += 1;\n"; + ss << " }\n"; + } + else if (pCur->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken* pSVR = + dynamic_cast< const formula::SingleVectorRefToken*>(pCur); + assert(pSVR); +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " if (!isNan("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << "))\n"; + ss << " {\n"; +#endif + ss << " arg="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; + ss << " fSum += arg;\n"; + ss << " fCount += 1;\n"; +#ifdef ISNAN + ss << " }\n"; + ss << " }\n"; +#endif + } + else + { + ss << " arg=" << pCur->GetDouble() << ";\n"; + ss << " fSum += arg;\n"; + ss << " fCount += 1;\n"; + } + } + else + { + ss << " arg = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; + ss << " fSum += arg;\n"; + ss << " fCount += 1;\n"; + } + if (i == 0) + { + ss << " fMean=fSum/" << "fCount" << ";\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(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 < "<< nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#else + ss << "0; i < "<< nCurWindowSize << "; i++)\n"; + ss << " {\n"; +#endif + } + + 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 (pCur->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken* pSVR = + dynamic_cast< const formula::SingleVectorRefToken*>(pCur); + assert(pSVR); +#ifdef ISNAN + ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " if (!isNan("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << "))\n"; + ss << " {\n"; +#endif + ss << " arg="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n"; + ss << " vSum += (arg - fMean) * (arg - fMean);\n"; +#ifdef ISNAN + ss << " }\n"; + ss << " }\n"; +#endif + } + 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 <= 1.0)\n"; + ss << " return DBL_MAX;\n"; + ss << " else\n"; + ss << " return vSum / (fCount - 1.0);\n"; + ss << "}\n"; +} void OpFdist::BinInlineFun(std::set& decls, std::set& funs) { diff --git a/sc/source/core/opencl/op_statistical.hxx b/sc/source/core/opencl/op_statistical.hxx index cbd276a62ed8..33d596ab32b9 100644 --- a/sc/source/core/opencl/op_statistical.hxx +++ b/sc/source/core/opencl/op_statistical.hxx @@ -21,7 +21,13 @@ public: const std::string sSymName, SubArguments &vSubArguments); virtual std::string BinFuncName(void) const { return "Standard"; } }; - +class OpVar: public Normal +{ +public: + virtual void GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments); + virtual std::string BinFuncName(void) const { return "Var"; } +}; class OpWeibull: public Normal { public: @@ -37,7 +43,6 @@ public: virtual void BinInlineFun(std::set& ,std::set&); virtual std::string BinFuncName(void) const { return "Fdist"; } }; - class OpFisher: public Normal { public: -- cgit