diff options
Diffstat (limited to 'sc')
| -rw-r--r-- | sc/source/core/opencl/op_financial.cxx | 137 |
1 files changed, 49 insertions, 88 deletions
diff --git a/sc/source/core/opencl/op_financial.cxx b/sc/source/core/opencl/op_financial.cxx index 216e7e906cee..71ea85837c1b 100644 --- a/sc/source/core/opencl/op_financial.cxx +++ b/sc/source/core/opencl/op_financial.cxx @@ -3402,129 +3402,90 @@ void OpPPMT::GenSlidingWindowFunction(std::stringstream &ss, ss << ", "; vSubArguments[i]->GenSlidingWindowDecl(ss); } - ss << ") {\n "; - ss << "double tmp = 0;\n "; - ss << "int gid0 = get_global_id(0);\n "; - ss<<"double tmp0=0;\n "; - ss<<"double tmp1=0;\n "; - ss<<"double tmp2=0;\n "; - ss<<"double tmp3=0;\n "; - ss<<"double tmp4=0,tmp5=0;\n "; - size_t nItems = 0; + ss<<") {\n"; + ss<<" double tmp = 0;\n"; + ss<<" int gid0 = get_global_id(0);\n"; + ss<<" double arg=0;\n"; + ss<<" double tmp0=0;\n"; + ss<<" double tmp1=0;\n"; + ss<<" double tmp2=0;\n"; + ss<<" double tmp3=0;\n"; + ss<<" double tmp4=0,tmp5=0;\n"; ss <<"\n "; //while (i-- > 1) for (size_t i = 0; i < vSubArguments.size(); i++) { FormulaToken *pCur = vSubArguments[i]->GetFormulaToken(); assert(pCur); - 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 "; -#else - ss << "gid0; i < "<< nCurWindowSize << "; i++)\n "; -#endif - } else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed()) { -#ifdef ISNAN - ss << "0; i < " << pDVR->GetArrayLength(); - ss << " && i < gid0+"<< nCurWindowSize << "; i++){\n "; -#else - ss << "0; i < gid0+"<< nCurWindowSize << "; i++)\n "; -#endif - } else if (!pDVR->IsStartFixed() && !pDVR->IsEndFixed()){ -#ifdef ISNAN - ss << "0; i + gid0 < " << pDVR->GetArrayLength(); - ss << " && i < "<< nCurWindowSize << "; i++){\n "; -#else - ss << "0; i < "<< nCurWindowSize << "; i++)\n "; -#endif - } - else { -#ifdef ISNAN - ss << "0; i < "<< nCurWindowSize << "; i++){\n "; -#else - ss << "0; i < "<< nCurWindowSize << "; i++)\n "; -#endif - } - nItems += nCurWindowSize; - } - else if (pCur->GetType() == formula::svSingleVectorRef) + if (pCur->GetType() == formula::svSingleVectorRef) { #ifdef ISNAN const formula::SingleVectorRefToken* pSVR = dynamic_cast< const formula::SingleVectorRefToken* >(pCur); - ss << "if (gid0 < " << pSVR->GetArrayLength() << "){\n "; + ss << "if (gid0 < " << pSVR->GetArrayLength() << "){\n"; #else - nItems += 1; #endif } else if (pCur->GetType() == formula::svDouble) { #ifdef ISNAN - ss << "{\n "; + ss << "{\n"; #endif - nItems += 1; } else { #ifdef ISNAN #endif - nItems += 1; } #ifdef ISNAN if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode()) { - ss << "if (isNan("; + ss << " arg="; ss << vSubArguments[i]->GenSlidingWindowDeclRef(); - ss << "))\n "; - ss << "tmp"<<i<<"= 0;\n "; - ss << "else\n "; - ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef(); - ss << ";\n }\n "; + ss << ";\n"; + ss << " if (isNan(arg))\n"; + ss << " tmp"<<i<<"= 0;\n"; + ss << " else\n"; + ss << " tmp"<<i<<"=arg;\n"; + ss << " }\n"; } else { - ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef(); - ss <<";\n "; + ss<<" tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef(); + ss<<";\n"; } #else - ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef(); - ss <<";\n "; + ss<<"tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef(); + ss<<";\n"; #endif } - ss <<"double pmt ;\n "; - ss<<"if(tmp0==0.0)\n "; - ss<<" return -(tmp3+tmp4)/"; - ss<<"tmp2;\n double temp1=0;\n "; - ss <<"double abl = pow(1.0+tmp0,tmp2);\n "; - ss <<"temp1-=tmp4"; - ss<<";\n "; - ss<<"temp1-=tmp3*abl;\n "; - ss <<"pmt = temp1/(1.0+tmp0*tmp5"; - ss<<") / ( (abl-1.0)/tmp0);\n "; - ss <<"double temp = pow( 1+tmp0,tmp1-2);\n "; - ss <<"double re;\n "; - ss <<"if(tmp1==1.0)\n {\n "; - ss <<"if(tmp5>0.0)\n "; - ss <<" re=0.0;\n "; - ss << "else\n re=-tmp3;\n }\n "; - ss <<"else\n {\n "; - ss <<"if(tmp5>0.0)\n "; - ss <<" re=GetZw(tmp0, tmp1-2.0, pmt, tmp3, 1.0) - pmt;\n "; - ss <<"else\n "; - ss <<" re=GetZw(tmp0, tmp1-1.0, pmt, tmp3, 0.0);\n "; - ss <<"}\n "; - ss <<"re = re * tmp0;\n "; - ss <<"tmp = pmt - re;\n "; - ss << "return tmp;\n"; - ss << "}"; + ss<<" double pmt=0 ;\n"; + ss<<" if(tmp0==0.0)\n"; + ss<<" return -(tmp3+tmp4)/tmp2;\n"; + ss<<" pmt=pmt-tmp4-tmp3*pow(1.0+tmp0,tmp2);\n"; + ss<<" pmt=pmt*pow(( (1.0+tmp0*tmp5)* "; + ss<<"( (pow(1.0+tmp0,tmp2)-1.0)/tmp0)),-1);\n"; + ss<<" double temp = pow( 1+tmp0,tmp1-2);\n"; + ss<<" double re;\n"; + ss<<" if(tmp1==1.0){\n"; + ss<<" if(tmp5>0.0)\n"; + ss<<" re=0.0;\n"; + ss<<" else\n"; + ss<<" re=-tmp3;\n"; + ss<<" }\n"; + ss<<" else\n"; + ss<<" {\n"; + ss<<" if(tmp5>0.0)\n "; + ss<<" re=GetZw(tmp0, tmp1-2.0, pmt, tmp3, 1.0) - pmt;\n"; + ss<<" else\n"; + ss<<" re=GetZw(tmp0, tmp1-1.0, pmt, tmp3, 0.0);\n"; + ss<<" }\n "; + ss<<" re = re * tmp0;\n"; + ss<<" tmp = pmt - re;\n"; + ss<<" return tmp;\n"; + ss<<"}"; } + void OpCoupdaybs::BinInlineFun(std::set<std::string>& decls, std::set<std::string>& funs) { |