diff options
author | mingli <mingli@multicorewareinc.com> | 2013-11-07 14:34:23 +0800 |
---|---|---|
committer | I-Jui (Ray) Sung <ray@multicorewareinc.com> | 2013-11-13 17:44:13 -0600 |
commit | a14622c88eda570e029eeb81baadda37d4417196 (patch) | |
tree | 4d0328341ae6b90dc621a73f4067bf24a7d14d7d /sc/source/core/opencl/op_statistical.cxx | |
parent | 848b91fbdf668d0de0511992d8ca22911d63c0b6 (diff) |
GPU Calc: implemented GAMMADIST
AMLOEXT-157 FIX
Change-Id: I806b6522e7481af87c03fbd1b1d500eb5cd1b9bf
Signed-off-by: haochen <haochen@multicorewareinc.com>
Signed-off-by: I-Jui (Ray) Sung <ray@multicorewareinc.com>
Diffstat (limited to 'sc/source/core/opencl/op_statistical.cxx')
-rw-r--r-- | sc/source/core/opencl/op_statistical.cxx | 81 |
1 files changed, 81 insertions, 0 deletions
diff --git a/sc/source/core/opencl/op_statistical.cxx b/sc/source/core/opencl/op_statistical.cxx index a58eee8a91f3..8e39c363b5f7 100644 --- a/sc/source/core/opencl/op_statistical.cxx +++ b/sc/source/core/opencl/op_statistical.cxx @@ -2471,7 +2471,88 @@ void OpLogNormDist::GenSlidingWindowFunction(std::stringstream &ss, ss << " return tmp;\n"; ss << "}\n"; } +void OpGammaDist::BinInlineFun(std::set<std::string>& decls, + std::set<std::string>& funs) +{ + decls.insert(fBigInvDecl);decls.insert(fLogDblMaxDecl); + decls.insert(fHalfMachEpsDecl);decls.insert(fMaxGammaArgumentDecl); + decls.insert(GetGammaSeriesDecl);decls.insert(GetGammaContFractionDecl); + decls.insert(GetLowRegIGammaDecl);decls.insert(GetGammaDistDecl); + decls.insert(GetGammaDistPDFDecl); + funs.insert(GetGammaSeries);funs.insert(GetGammaContFraction); + funs.insert(GetLowRegIGamma);funs.insert(GetGammaDist); + funs.insert(GetGammaDistPDF); +} +void OpGammaDist::GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments) +{ + FormulaToken *tmpCur0 = vSubArguments[0]->GetFormulaToken(); + const formula::SingleVectorRefToken*tmpCurDVR0= dynamic_cast<const +formula::SingleVectorRefToken *>(tmpCur0); + FormulaToken *tmpCur1 = vSubArguments[1]->GetFormulaToken(); + const formula::SingleVectorRefToken*tmpCurDVR1= dynamic_cast<const +formula::SingleVectorRefToken *>(tmpCur1); + FormulaToken *tmpCur2 = vSubArguments[2]->GetFormulaToken(); + const formula::SingleVectorRefToken*tmpCurDVR2= dynamic_cast<const +formula::SingleVectorRefToken *>(tmpCur2); + FormulaToken *tmpCur3 = vSubArguments[3]->GetFormulaToken(); + const formula::SingleVectorRefToken*tmpCurDVR3= dynamic_cast<const +formula::SingleVectorRefToken *>(tmpCur3); + 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 arg0 = "; + ss << vSubArguments[0]->GenSlidingWindowDeclRef(); + ss << ";\n"; + ss << " double arg1 = "; + ss << vSubArguments[1]->GenSlidingWindowDeclRef(); + ss << ";\n"; + ss << " double arg2 = "; + ss << vSubArguments[2]->GenSlidingWindowDeclRef(); + ss << ";\n"; + ss << " double arg3 = "; + ss << vSubArguments[3]->GenSlidingWindowDeclRef(); + ss << ";\n"; + ss << " double tmp;\n"; +#ifdef ISNAN + ss << " if(isNan(arg0)||(gid0>="; + ss << tmpCurDVR0->GetArrayLength(); + ss << "))\n"; + ss << " arg0 = 0;\n"; +#endif +#ifdef ISNAN + ss << " if(isNan(arg1)||(gid0>="; + ss << tmpCurDVR1->GetArrayLength(); + ss << "))\n"; + ss << " arg1 = 0;\n"; +#endif +#ifdef ISNAN + ss << " if(isNan(arg2)||(gid0>="; + ss << tmpCurDVR2->GetArrayLength(); + ss << "))\n"; + ss << " arg2 = 0;\n"; +#endif +#ifdef ISNAN + ss << " if(isNan(arg3)||(gid0>="; + ss << tmpCurDVR3->GetArrayLength(); + ss << "))\n"; + ss << " arg3 = 0;\n"; +#endif + ss << " if (arg3)\n"; + ss << " tmp=GetGammaDist( arg0, arg1, arg2);\n"; + ss << " else\n"; + ss << " tmp=GetGammaDistPDF( arg0, arg1, arg2);\n"; + ss << " return tmp;\n"; + ss << "}\n"; +} }} /* vim:set shiftwidth=4 softtabstop=4 expandtab: */ |