diff options
Diffstat (limited to 'sc/source/core/opencl/op_addin.cxx')
-rw-r--r-- | sc/source/core/opencl/op_addin.cxx | 236 |
1 files changed, 236 insertions, 0 deletions
diff --git a/sc/source/core/opencl/op_addin.cxx b/sc/source/core/opencl/op_addin.cxx new file mode 100644 index 000000000..8a42de323 --- /dev/null +++ b/sc/source/core/opencl/op_addin.cxx @@ -0,0 +1,236 @@ +/* -*- Mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */ +/* + * This file is part of the LibreOffice project. + * + * This Source Code Form is subject to the terms of the Mozilla Public + * License, v. 2.0. If a copy of the MPL was not distributed with this + * file, You can obtain one at http://mozilla.org/MPL/2.0/. + */ + +#include "op_addin.hxx" + +#include <formula/vectortoken.hxx> +#include <sstream> + +using namespace formula; + +namespace sc::opencl { + +void OpBesselj::GenSlidingWindowFunction(std::stringstream &ss, + const std::string &sSymName, SubArguments &vSubArguments) +{ + CHECK_PARAMETER_COUNT( 2, 2 ); + ss << "\ndouble " << sSymName; + ss << "_" << BinFuncName() << "("; + for (size_t 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 x = 0.0;\n"; + ss << " double N = 0.0;\n"; + FormulaToken *tmpCur0 = vSubArguments[0]->GetFormulaToken(); + assert(tmpCur0); + if(ocPush == vSubArguments[0]->GetFormulaToken()->GetOpCode()) + { + if(tmpCur0->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken*tmpCurSVR0 = + static_cast<const formula::SingleVectorRefToken *>(tmpCur0); + ss << " if (gid0 < " << tmpCurSVR0->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " x = "; + ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n"; + ss << " if (isnan(x))\n"; + ss << " x = 0.0;\n"; + ss << " }\n"; + } + else if(tmpCur0->GetType() == formula::svDouble) + { + ss << " x = " << tmpCur0->GetDouble() << ";\n"; + } + else + { + throw Unhandled(__FILE__, __LINE__); + } + } + else + { + ss << " x = "; + ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n"; + } + + FormulaToken *tmpCur1 = vSubArguments[1]->GetFormulaToken(); + assert(tmpCur1); + if(ocPush == vSubArguments[1]->GetFormulaToken()->GetOpCode()) + { + if(tmpCur1->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken*tmpCurSVR1 = + static_cast<const formula::SingleVectorRefToken *>(tmpCur1); + ss << " if (gid0 < " << tmpCurSVR1->GetArrayLength() << ")\n"; + ss << " {\n"; + ss << " N = "; + ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n"; + ss << " if (isnan(N))\n"; + ss << " N = 0.0;\n"; + ss << " }\n"; + } + else if(tmpCur1->GetType() == formula::svDouble) + { + ss << " N = " << tmpCur1->GetDouble() << ";\n"; + } + else + { + throw Unhandled(__FILE__, __LINE__); + } + } + else + { + ss << " N = "; + ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n"; + } + ss << " double f_PI = 3.1415926535897932385;\n"; + ss << " double f_2_DIV_PI = 2.0 / f_PI;\n"; + ss << " double f_PI_DIV_2 = f_PI / 2.0;\n"; + ss << " double f_PI_DIV_4 = f_PI / 4.0;\n"; + ss << " if( N < 0.0 )\n"; + ss << " return CreateDoubleError(IllegalArgument);\n"; + ss << " if (x == 0.0)\n"; + ss << " return (N == 0.0) ? 1.0 : 0.0;\n"; + ss << " double fSign = ((int)N % 2 == 1 && x < 0.0) ? -1.0 : 1.0;\n"; + ss << " double fX = fabs(x);\n"; + ss << " double fMaxIteration = 9000000.0;\n"; + ss << " double fEstimateIteration = fX * 1.5 + N;\n"; + ss << " bool bAsymptoticPossible = pow(fX,0.4) > N;\n"; + ss << " if (fEstimateIteration > fMaxIteration)\n"; + ss << " {\n"; + ss << " if (bAsymptoticPossible)\n"; + ss << " return fSign * sqrt(f_2_DIV_PI/fX)"; + ss << "* cos(fX-N*f_PI_DIV_2-f_PI_DIV_4);\n"; + ss << " else\n"; + ss << " return CreateDoubleError(NoConvergence);\n"; + ss << " }\n"; + ss << " double epsilon = 1.0e-15;\n"; + ss << " bool bHasfound = false;\n"; + ss << " double k= 0.0;\n"; + ss << " double u ;\n"; + ss << " double m_bar;\n"; + ss << " double g_bar;\n"; + ss << " double g_bar_delta_u;\n"; + ss << " double g = 0.0;\n"; + ss << " double delta_u = 0.0;\n"; + ss << " double f_bar = -1.0;\n"; + ss << " if (N==0)\n"; + ss << " {\n"; + ss << " u = 1.0;\n"; + ss << " g_bar_delta_u = 0.0;\n"; + ss << " g_bar = - 2.0/fX; \n"; + ss << " delta_u = g_bar_delta_u / g_bar;\n"; + ss << " u = u + delta_u ;\n"; + ss << " g = -1.0 / g_bar; \n"; + ss << " f_bar = f_bar * g;\n"; + ss << " k = 2.0;\n"; + ss << " }\n"; + ss << " if (N!=0)\n"; + ss << " {\n"; + ss << " u=0.0;\n"; + ss << " for (k =1.0; k<= N-1; k = k + 1.0)\n"; + ss << " {\n"; + ss << " m_bar=2.0 * fmod(k-1.0, 2.0) * f_bar;\n"; + ss << " g_bar_delta_u = - g * delta_u - m_bar * u;\n"; + ss << " g_bar = m_bar - 2.0*k/fX + g;\n"; + ss << " delta_u = g_bar_delta_u / g_bar;\n"; + ss << " u = u + delta_u;\n"; + ss << " g = -1.0/g_bar;\n"; + ss << " f_bar=f_bar * g;\n"; + ss << " }\n"; + ss << " m_bar=2.0 * fmod(k-1.0, 2.0) * f_bar;\n"; + ss << " g_bar_delta_u = f_bar - g * delta_u - m_bar * u;\n"; + ss << " g_bar = m_bar - 2.0*k/fX + g;\n"; + ss << " delta_u = g_bar_delta_u / g_bar;\n"; + ss << " u = u + delta_u;\n"; + ss << " g = -1.0/g_bar;\n"; + ss << " f_bar = f_bar * g;\n"; + ss << " k = k + 1.0;\n"; + ss << " }\n"; + ss << " do\n"; + ss << " {\n"; + ss << " m_bar = 2.0 * fmod(k-1.0, 2.0) * f_bar;\n"; + ss << " g_bar_delta_u = - g * delta_u - m_bar * u;\n"; + ss << " g_bar = m_bar - 2.0*k/fX + g;\n"; + ss << " delta_u = g_bar_delta_u / g_bar;\n"; + ss << " u = u + delta_u;\n"; + ss << " g = -pow(g_bar,-1.0);\n"; + ss << " f_bar = f_bar * g;\n"; + ss << " bHasfound = (fabs(delta_u)<=fabs(u)*epsilon);\n"; + ss << " k = k + 1.0;\n"; + ss << " }\n"; + ss << " while (!bHasfound && k <= fMaxIteration);\n"; + ss << " if (bHasfound)\n"; + ss << " return u * fSign;\n"; + ss << " else\n"; + ss << " return CreateDoubleError(NoConvergence);\n"; + ss << "}"; +} +void OpGestep::GenSlidingWindowFunction( + std::stringstream &ss,const std::string &sSymName, + SubArguments &vSubArguments) +{ + ss << "\ndouble " << sSymName; + ss << "_"<< BinFuncName() <<"("; + for (size_t i = 0; i < vSubArguments.size(); i++) + { + if (i) + ss << ","; + vSubArguments[i]->GenSlidingWindowDecl(ss); + } + ss << ")\n"; + ss << "{\n"; + ss << " double tmp=0,tmp0 =0,tmp1 = 0;\n"; + ss << " int gid0=get_global_id(0);\n"; + ss <<"\n"; + for (size_t i = 0; i < vSubArguments.size(); i++) + { + FormulaToken *pCur = vSubArguments[i]->GetFormulaToken(); + assert(pCur); + if (pCur->GetType() == formula::svSingleVectorRef) + { + const formula::SingleVectorRefToken& rSVR = + dynamic_cast< const formula::SingleVectorRefToken& >(*pCur); + ss << " if (gid0 < " << rSVR.GetArrayLength() << ")\n"; + ss << " {\n"; + } + else if (pCur->GetType() == formula::svDouble) + { + ss << " {\n"; + } + + if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode()) + { + ss << " if (isnan("; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << "))\n"; + ss << " tmp"<<i<<" = 0;\n"; + ss << " else\n"; + ss << " tmp"<<i<<" = "; + ss << vSubArguments[i]->GenSlidingWindowDeclRef(); + ss << ";\n }\n"; + } + else + { + ss << "tmp"<<i<<" ="<<vSubArguments[i]->GenSlidingWindowDeclRef(); + ss <<";\n"; + } + } + ss << " tmp =tmp0 >= tmp1 ? 1 : 0;\n"; + ss << " return tmp;\n"; + ss << "}\n"; +} + +} + +/* vim:set shiftwidth=4 softtabstop=4 expandtab: */ |