diff options
author | dechuang <dechuang@multicorewareinc.com> | 2013-12-11 16:51:05 +0800 |
---|---|---|
committer | I-Jui (Ray) Sung <ray@multicorewareinc.com> | 2013-12-19 17:53:47 -0600 |
commit | 68226cbd578a316c8190867daf76f9b6f5d435aa (patch) | |
tree | e232e58751b8985dae5804a6eb2a2267ecf4170d /sc/source | |
parent | 68401ee7c2b7c50e26e6ed11e12eb7cd8ca23adb (diff) |
GPU Calc: implemented for DVARP
AMLOEXT-343 FIX
Change-Id: I8a56eeb016992b113d8b161df6256a4630249dfe
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>
Diffstat (limited to 'sc/source')
-rw-r--r-- | sc/source/core/opencl/formulagroupcl.cxx | 4 | ||||
-rw-r--r-- | sc/source/core/opencl/op_database.cxx | 205 | ||||
-rw-r--r-- | sc/source/core/opencl/op_database.hxx | 7 |
3 files changed, 216 insertions, 0 deletions
diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx index 62329a611c84..b1ed9f44bde1 100644 --- a/sc/source/core/opencl/formulagroupcl.cxx +++ b/sc/source/core/opencl/formulagroupcl.cxx @@ -2404,6 +2404,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments( mvSubArguments.push_back(SoPHelper(ts, ft->Children[i], new OpDvar)); break; + case ocDBVarP: + mvSubArguments.push_back(SoPHelper(ts, + ft->Children[i], new OpDvarp)); + break; case ocExternal: if ( !(pChild->GetExternal().compareTo(OUString( "com.sun.star.sheet.addin.Analysis.getEffect")))) diff --git a/sc/source/core/opencl/op_database.cxx b/sc/source/core/opencl/op_database.cxx index 70e7941b74a7..e8e5b66c643d 100644 --- a/sc/source/core/opencl/op_database.cxx +++ b/sc/source/core/opencl/op_database.cxx @@ -1340,6 +1340,211 @@ void OpDvar::GenSlidingWindowFunction(std::stringstream &ss, ss << "}"; } +void OpDvarp::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 {\n"; + ss << " int gid0=get_global_id(0);\n"; + ss << " double var = 0;\n"; + ss << " double mean = 0;\n"; + ss << " double value =0;\n"; + ss << " int count = 0;\n"; + GenTmpVariables(ss,vSubArguments); + int dataCol = 0; + int dataRow = 0; + if(vSubArguments[0]->GetFormulaToken()->GetType() == + formula::svDoubleVectorRef) + { + formula::FormulaToken *tmpCur = vSubArguments[0]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR1= dynamic_cast<const + formula::DoubleVectorRefToken *>(tmpCur); + dataCol = pCurDVR1->GetArrays().size(); + dataRow = pCurDVR1->GetArrayLength(); + } + else + throw Unhandled(); + if(vSubArguments[dataCol]->GetFormulaToken()->GetType() == + formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken*pTmpDVR1= dynamic_cast<const + formula::SingleVectorRefToken *>(vSubArguments[dataCol]-> + GetFormulaToken()); + ss << " tmp"<<dataCol<<"="; + ss << vSubArguments[dataCol]->GenSlidingWindowDeclRef()<<";\n"; + ss << " if(gid0>="<<pTmpDVR1->GetArrayLength()<<" ||isNan("; + ss << "tmp"<<dataCol<<"))\n"; + ss << " tmp"<<dataCol<<"=0;\n"; + } + else + throw Unhandled(); + int conditionCol = 0; + int conditionRow = 0; + if(vSubArguments[dataCol + 1]->GetFormulaToken()->GetType() == + formula::svDoubleVectorRef) + { + formula::FormulaToken *tmpCur = vSubArguments[dataCol + 1]-> + GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR2= dynamic_cast<const + formula::DoubleVectorRefToken *>(tmpCur); + conditionCol = pCurDVR2->GetArrays().size(); + conditionRow = pCurDVR2->GetArrayLength(); + } + else{ + throw Unhandled(); + } + if(dataCol!=conditionCol) + throw Unhandled(); + if(dataCol > 0 && dataRow > 0) + { + formula::FormulaToken *tmpCur1 = vSubArguments[0]->GetFormulaToken(); + formula::FormulaToken *tmpCur2 = vSubArguments[dataCol + 1]-> + GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR1= dynamic_cast<const + formula::DoubleVectorRefToken *>(tmpCur1); + const formula::DoubleVectorRefToken*pCurDVR2= dynamic_cast<const + formula::DoubleVectorRefToken *>(tmpCur2); + + if(pCurDVR1->IsStartFixed() && pCurDVR1->IsEndFixed() && + pCurDVR2->IsStartFixed() && pCurDVR2->IsEndFixed()) + { + ss << " int i,j,p;\n"; + ss << " bool flag;\n"; + + ss << " for(p = 1;p < " << dataRow << ";++p)\n"; + ss << " {\n"; + ss << " i = p;\n"; + for(int i = 0; i < dataCol; ++i){ + if(vSubArguments[i]->GetFormulaToken()->GetType() != + formula::svDoubleVectorRef) + throw Unhandled(); + formula::FormulaToken *tmpCur_tmp = + vSubArguments[i]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR_tmp= + dynamic_cast<const formula::DoubleVectorRefToken *> + (tmpCur_tmp); + ss << " tmp"<<i<<"="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n"; + ss << " if(isNan(tmp"<<i<<"))\n"; + ss <<" tmp"<<i<<" = 0;\n"; + } + ss << " flag = false;\n"; + ss << " for(j = 1; j < " << conditionRow << ";++j)\n"; + ss << " {\n"; + ss << " i = j;\n"; + ss << " if (flag)\n"; + ss << " break;\n"; + ss << " else{\n"; + for(int i = dataCol + 1; i < dataCol + 1 + conditionCol; ++i){ + if(vSubArguments[i]->GetFormulaToken()->GetType() != + formula::svDoubleVectorRef) + throw Unhandled(); + formula::FormulaToken *tmpCur_tmp = + vSubArguments[i]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR_tmp= + dynamic_cast<const formula::DoubleVectorRefToken *> + (tmpCur_tmp); + ss << " tmp"<<i<<"="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n"; + ss << " if(!isNan(tmp"<<i<<")){\n"; + ss << " if(tmp"<<(i-dataCol-1)<<"!=tmp"; + ss << i<<"){\n"; + ss << " continue;\n"; + ss << " }\n"; + ss << " }\n"; + } + ss << " flag=true;\n"; + ss << " }\n"; + ss << " }\n"; + ss << " if (flag){\n"; + ss << " count++;\n"; + for(int i = 0; i < dataCol; ++i){ + ss << " if(tmp"<<dataCol<<"=="<<(i+1)<<"){\n"; + ss << " value=tmp"<<i<<";\n"; + ss << " }\n"; + } + ss << " mean+=value;\n"; + ss << " }\n"; + ss << " }\n"; + + ss << " if(count<=0)\n"; + ss << " return 0;\n"; + + ss << " mean/=count;\n"; + + ss << " for(p = 1;p < " << dataRow << ";++p)\n"; + ss << " {\n"; + ss << " i = p;\n"; + for(int i = 0; i < dataCol; ++i){ + if(vSubArguments[i]->GetFormulaToken()->GetType() != + formula::svDoubleVectorRef) + throw Unhandled(); + formula::FormulaToken *tmpCur_tmp = + vSubArguments[i]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR_tmp= + dynamic_cast<const formula::DoubleVectorRefToken *> + (tmpCur_tmp); + ss << " tmp"<<i<<"="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n"; + ss << " if(isNan(tmp"<<i<<"))\n"; + ss <<" tmp"<<i<<" = 0;\n"; + } + ss << " flag = false;\n"; + ss << " for(j = 1; j < " << conditionRow << ";++j)\n"; + ss << " {\n"; + ss << " i = j;\n"; + ss << " if (flag)\n"; + ss << " break;\n"; + ss << " else{\n"; + for(int i = dataCol + 1; i < dataCol + 1 + conditionCol; ++i){ + if(vSubArguments[i]->GetFormulaToken()->GetType() != + formula::svDoubleVectorRef) + throw Unhandled(); + formula::FormulaToken *tmpCur_tmp = + vSubArguments[i]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR_tmp= + dynamic_cast<const formula::DoubleVectorRefToken *> + (tmpCur_tmp); + ss << " tmp"<<i<<"="; + ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n"; + ss << " if(!isNan(tmp"<<i<<")){\n"; + ss << " if(tmp"<<(i-dataCol-1)<<"!=tmp"; + ss << i<<"){\n"; + ss << " continue;\n"; + ss << " }\n"; + ss << " }\n"; + } + ss << " flag=true;\n"; + ss << " }\n"; + ss << " }\n"; + ss << " if (flag){\n"; + for(int i = 0; i < dataCol; ++i){ + ss << " if(tmp"<<dataCol<<"=="<<(i+1)<<"){\n"; + ss << " value=tmp"<<i<<";\n"; + ss << " }\n"; + } + ss << " var+=pow(mean-value,2);\n"; + ss << " }\n"; + ss << " }\n"; + + ss << " var = var/count;\n"; + } + else + ss << "var = -1;\n"; + } + else + ss << "var = -1;\n"; + ss << " return var;\n"; + ss << "}"; +} + }} diff --git a/sc/source/core/opencl/op_database.hxx b/sc/source/core/opencl/op_database.hxx index df464ac64362..bd668d321862 100644 --- a/sc/source/core/opencl/op_database.hxx +++ b/sc/source/core/opencl/op_database.hxx @@ -79,6 +79,13 @@ public: virtual std::string BinFuncName(void) const { return "Dvar"; } }; +class OpDvarp: public CheckVariables +{ +public: + virtual void GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments); + virtual std::string BinFuncName(void) const { return "Dvarp"; } +}; }} |