From 64527a0090f440389891d0a29bfe36ce832f8ec5 Mon Sep 17 00:00:00 2001 From: yiming ju Date: Tue, 5 Nov 2013 17:16:36 +0800 Subject: GPU Calc: implemented for XIRR AMLOEXT-137 FIX Change-Id: I108a55037dbbff18848c0166604b316548e1f3c0 Signed-off-by: haochen Signed-off-by: I-Jui (Ray) Sung --- sc/source/core/opencl/formulagroupcl.cxx | 6 ++ sc/source/core/opencl/op_financial.cxx | 120 +++++++++++++++++++++++++++++++ sc/source/core/opencl/op_financial.hxx | 10 +++ 3 files changed, 136 insertions(+) diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx index 1f53072588fb..fbe43c896d3d 100644 --- a/sc/source/core/opencl/formulagroupcl.cxx +++ b/sc/source/core/opencl/formulagroupcl.cxx @@ -1243,6 +1243,12 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments( mvSubArguments.push_back(SoPHelper(ts, ft->Children[i], new OpMDuration)); } + else if ( !(pChild->GetExternal().compareTo(OUString( + "com.sun.star.sheet.addin.Analysis.getXirr")))) + { + mvSubArguments.push_back(SoPHelper(ts, ft->Children[i], + new OpXirr)); + } break; default: throw UnhandledToken(pChild, "unhandled opcode"); diff --git a/sc/source/core/opencl/op_financial.cxx b/sc/source/core/opencl/op_financial.cxx index 0e619332ebb2..58b0a87cca00 100644 --- a/sc/source/core/opencl/op_financial.cxx +++ b/sc/source/core/opencl/op_financial.cxx @@ -4254,6 +4254,126 @@ void OpVDB::GenSlidingWindowFunction( } } + +void OpXirr::GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments) +{ + FormulaToken *tmpCur = vSubArguments[0]->GetFormulaToken(); + const formula::DoubleVectorRefToken*pCurDVR= dynamic_cast(tmpCur); + size_t nCurWindowSize = pCurDVR->GetArrayLength() < + pCurDVR->GetRefRowSize() ? pCurDVR->GetArrayLength(): + pCurDVR->GetRefRowSize() ; + 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 << " int doubleIndex = gid0;\n"; + ss << " int singleIndex = gid0;\n"; + ss << " double result = 0;\n"; + ss << " int i=0;\n"; + if(vSubArguments.size()<2) + { + ss << " result = -DBL_MAX;\n"; + ss << " return result;\n"; + }else + { + GenTmpVariables(ss,vSubArguments); + if(vSubArguments.size() == 2) + { + ss << " double tmp2 = 0.1;\n"; + }else + { + CheckSubArgumentIsNan(ss,vSubArguments,2); + } + ss << " if(tmp2<=-1)\n"; + ss << " result = -DBL_MAX;\n"; + ss << " else\n"; + ss << " {\n"; + ss << " double fMaxEps = 1e-10;\n"; + ss << " int nMaxIter = 50;\n"; + ss << " double fNewRate, fRateEps, fResultValue, fResultValue2;\n"; + ss << " int nIter = 0;\n"; + ss << " int bContLoop;\n"; + ss << " int windowsSize = "; + ss << nCurWindowSize; + ss << ";\n"; + CheckSubArgumentIsNan(ss,vSubArguments,0); + CheckSubArgumentIsNan(ss,vSubArguments,1); + ss << " double D_0 = tmp1;\n"; + ss << " double V_0 = tmp0;\n"; + ss << " double fResultRate = tmp2;\n"; + ss << " double r;\n"; + ss << " double fResult;\n"; + ss << " do\n"; + ss << " {\n"; + ss << " fResultValue = V_0;\n"; + ss << " r = fResultRate + 1;\n"; + ss << " for (i = "; + if (!pCurDVR->IsStartFixed() && pCurDVR->IsEndFixed()) { + ss << "gid0+1; i < "<< nCurWindowSize <<"; i++)\n"; + } else if (pCurDVR->IsStartFixed() && !pCurDVR->IsEndFixed()) { + ss << "1; i < gid0+"<< nCurWindowSize <<"; i++)\n"; + } else { + ss << "1; i < "<< nCurWindowSize <<"; i++)\n"; + } + ss << " {\n"; + if(!pCurDVR->IsStartFixed() && !pCurDVR->IsEndFixed()) + { + ss<< " doubleIndex =i+gid0;\n"; + }else + { + ss<< " doubleIndex =i;\n"; + } + CheckSubArgumentIsNan(ss,vSubArguments,0); + CheckSubArgumentIsNan(ss,vSubArguments,1); + ss << " fResultValue += tmp0/pow(r,(tmp1 - D_0)/365.0);\n"; + ss << " }\n"; + ss << " fResultValue2 = 0;\n"; + + ss << " for (i = "; + if (!pCurDVR->IsStartFixed() && pCurDVR->IsEndFixed()) { + ss << "gid0+1; i < "<< nCurWindowSize <<"; i++)\n"; + } else if (pCurDVR->IsStartFixed() && !pCurDVR->IsEndFixed()) { + ss << "1; i < gid0+"<< nCurWindowSize <<"; i++)\n"; + } else { + ss << "1; i < "<< nCurWindowSize <<"; i++)\n"; + } + ss << " {\n"; + if(!pCurDVR->IsStartFixed() && !pCurDVR->IsEndFixed()) + { + ss<< " doubleIndex =i+gid0;\n"; + }else + { + ss<< " doubleIndex =i;\n"; + } + CheckSubArgumentIsNan(ss,vSubArguments,0); + CheckSubArgumentIsNan(ss,vSubArguments,1); + ss << " double E_i = (tmp1 - D_0)/365.0;\n"; + ss << " fResultValue2 -= E_i * tmp0 / pow(r,E_i + 1.0);\n"; + ss << " }\n"; + ss << " fNewRate = fResultRate - fResultValue / fResultValue2;\n"; + ss << " fRateEps = fabs( fNewRate - fResultRate );\n"; + ss << " fResultRate = fNewRate;\n"; + ss << " bContLoop = (fRateEps > fMaxEps) && (fabs( fResultValue ) > fMaxEps);\n"; + ss << " }\n"; + ss << " while( bContLoop && (++nIter < nMaxIter) );\n"; + ss << " if( bContLoop )\n"; + ss << " result = -DBL_MAX;\n"; + ss << " result = fResultRate;\n"; + ss << " }\n"; + ss << " return result;\n"; + ss << "}"; + } + +} + }} /* vim:set shiftwidth=4 softtabstop=4 expandtab: */ diff --git a/sc/source/core/opencl/op_financial.hxx b/sc/source/core/opencl/op_financial.hxx index c7450bd84aee..68f291590250 100644 --- a/sc/source/core/opencl/op_financial.hxx +++ b/sc/source/core/opencl/op_financial.hxx @@ -481,6 +481,16 @@ public: virtual std::string BinFuncName(void) const { return "VDB"; } virtual void BinInlineFun(std::set& ,std::set& ); }; + +class OpXirr: public CheckVariables +{ +public: + virtual void GenSlidingWindowFunction(std::stringstream &ss, + const std::string sSymName, SubArguments &vSubArguments); + + virtual std::string BinFuncName(void) const { return "Xirr"; } +}; + }} #endif -- cgit