sc/source/core/opencl/formulagroupcl.cxx | 16
sc/source/core/opencl/op_statistical.cxx | 1399 ++++---------------------------
sc/source/core/opencl/op_statistical.hxx | 158 ++-
sc/source/core/opencl/opbase.cxx | 7
sc/source/core/opencl/opbase.hxx | 6
5 files changed, 326 insertions(+), 1260 deletions(-)
New commits:
commit 66ad69f5b300060f174fdc6c6e106dbba2b0140a
Author: Luboš Luňák <l.lu...@collabora.com>
AuthorDate: Mon Sep 19 17:51:02 2022 +0200
Commit: Luboš Luňák <l.lu...@collabora.com>
CommitDate: Tue Sep 20 08:46:10 2022 +0200
correct handling of 'nested' in opencl vectors
The 'nested' arguments means the value is used inside a loop, which
already implicitly does range checking.
Change-Id: I71ccbded2f9d623f58d11f11316aab28385f1614
Reviewed-on: https://gerrit.libreoffice.org/c/core/+/140193
Tested-by: Jenkins
Reviewed-by: Luboš Luňák <l.lu...@collabora.com>
diff --git a/sc/source/core/opencl/formulagroupcl.cxx
b/sc/source/core/opencl/formulagroupcl.cxx
index 564ee9cc6393..b4e35e932ba5 100644
--- a/sc/source/core/opencl/formulagroupcl.cxx
+++ b/sc/source/core/opencl/formulagroupcl.cxx
@@ -915,8 +915,8 @@ public:
virtual std::string GenSlidingWindowDeclRef( bool nested ) const override
{
outputstream ss;
- ss << "(!isnan(" << VectorRef::GenSlidingWindowDeclRef();
- ss << ")?" << VectorRef::GenSlidingWindowDeclRef();
+ ss << "(!isnan(" << VectorRef::GenSlidingWindowDeclRef(nested);
+ ss << ")?" << VectorRef::GenSlidingWindowDeclRef(nested);
ss << ":" << mStringArgument.GenSlidingWindowDeclRef(nested);
ss << ")";
return ss.str();
@@ -977,18 +977,18 @@ std::string
DynamicKernelSlidingArgument<Base>::GenSlidingWindowDeclRef( bool ne
outputstream ss;
if (!bIsStartFixed && !bIsEndFixed)
{
- if (nested)
+ if (!nested)
ss << "((i+gid0) <" << nArrayLength << "?";
ss << Base::GetName() << "[i + gid0]";
- if (nested)
+ if (!nested)
ss << ":NAN)";
}
else
{
- if (nested)
+ if (!nested)
ss << "(i <" << nArrayLength << "?";
ss << Base::GetName() << "[i]";
- if (nested)
+ if (!nested)
ss << ":NAN)";
}
return ss.str();
@@ -1146,8 +1146,8 @@ public:
virtual std::string GenSlidingWindowDeclRef( bool nested ) const override
{
outputstream ss;
- ss << "(!isnan(" << mDoubleArgument.GenSlidingWindowDeclRef();
- ss << ")?" << mDoubleArgument.GenSlidingWindowDeclRef();
+ ss << "(!isnan(" << mDoubleArgument.GenSlidingWindowDeclRef(nested);
+ ss << ")?" << mDoubleArgument.GenSlidingWindowDeclRef(nested);
ss << ":" << mStringArgument.GenSlidingWindowDeclRef(nested);
ss << ")";
return ss.str();
commit 4c2b8db7814dead51607935c50c1d0160333d616
Author: Luboš Luňák <l.lu...@collabora.com>
AuthorDate: Mon Sep 19 17:50:58 2022 +0200
Commit: Luboš Luňák <l.lu...@collabora.com>
CommitDate: Tue Sep 20 08:45:48 2022 +0200
reduce opencl copy&paste in op_statistical.cxx (svDoubleVectorRef pair)
Change-Id: I92daab29290c9d72f7531b25950c7f52084cc0bd
Reviewed-on: https://gerrit.libreoffice.org/c/core/+/140192
Tested-by: Jenkins
Reviewed-by: Luboš Luňák <l.lu...@collabora.com>
diff --git a/sc/source/core/opencl/op_statistical.cxx
b/sc/source/core/opencl/op_statistical.cxx
index 29e150de10cd..f4bf77715cde 100644
--- a/sc/source/core/opencl/op_statistical.cxx
+++ b/sc/source/core/opencl/op_statistical.cxx
@@ -586,293 +586,6 @@ void OpStDevP::GenSlidingWindowFunction(outputstream &ss,
ss << "}\n";
}
-void OpSlope::GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments)
-{
- CHECK_PARAMETER_COUNT(2,2);
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0 = get_global_id(0);\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fMeanX = 0.0;\n";
- ss << " double fMeanY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fSumSqrDeltaX = 0.0;\n";
- ss << " double fCount = 0.0;\n";
- ss << " double argX = 0.0;\n";
- ss << " double argY = 0.0;\n";
- FormulaToken *pCur = vSubArguments[1]->GetFormulaToken();
- FormulaToken *pCur1 = vSubArguments[0]->GetFormulaToken();
- assert(pCur);
- assert(pCur1);
- if (pCur->GetType() == formula::svDoubleVectorRef&&
- pCur1->GetType() == formula::svDoubleVectorRef)
- {
- const formula::DoubleVectorRefToken* pDVR =
- static_cast<const formula::DoubleVectorRefToken *>(pCur);
- const formula::DoubleVectorRefToken* pDVR1 =
- static_cast<const formula::DoubleVectorRefToken *>(pCur1);
-
- size_t nCurWindowSize = pDVR->GetRefRowSize();
- size_t nCurWindowSize1 = pDVR1->GetRefRowSize();
- size_t arrayLength = pDVR->GetArrayLength()<
- pDVR1->GetArrayLength() ? pDVR->GetArrayLength():
- pDVR1->GetArrayLength();
- if(nCurWindowSize != nCurWindowSize1)
- throw Unhandled(__FILE__, __LINE__);
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength ;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- throw Unhandled(__FILE__, __LINE__);
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX) || isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumX += argX;\n";
- ss << " fSumY += argY;\n";
- ss << " fCount += 1.0;\n";
- ss << " }\n";
-
- ss << " if (fCount < 1.0)\n";
- ss << " return CreateDoubleError(NoValue);\n";
- ss << " else\n";
- ss << " {\n";
- ss << " fMeanX = fSumX / fCount;\n";
- ss << " fMeanY = fSumY / fCount;\n";
-
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength ;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX) || isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumDeltaXDeltaY += (argX-fMeanX)*(argY-fMeanY);\n";
- ss << " fSumSqrDeltaX += (argX-fMeanX) * (argX-fMeanX);\n";
- ss << " }\n";
- ss << " if(fSumSqrDeltaX == 0.0)\n";
- ss << " return CreateDoubleError(DivisionByZero);\n";
- ss << " else\n";
- ss << " {\n";
- ss << " return fSumDeltaXDeltaY / fSumSqrDeltaX;\n";
- ss << " }\n";
- ss << " }\n";
- ss << "}\n";
- }
- else
- {
- throw Unhandled(__FILE__, __LINE__);
- }
-}
-void OpSTEYX::GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments)
-{
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0 = get_global_id(0);\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fMeanX = 0.0;\n";
- ss << " double fMeanY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fSumSqrDeltaX = 0.0;\n";
- ss << " double fSumSqrDeltaY = 0.0;\n";
- ss << " double fCount = 0.0;\n";
- ss << " double argX = 0.0;\n";
- ss << " double argY = 0.0;\n";
- FormulaToken *pCur = vSubArguments[1]->GetFormulaToken();
- FormulaToken *pCur1 = vSubArguments[0]->GetFormulaToken();
- assert(pCur);
- assert(pCur1);
- if (pCur->GetType() == formula::svDoubleVectorRef&&
- pCur1->GetType() == formula::svDoubleVectorRef)
- {
- const formula::DoubleVectorRefToken* pDVR =
- static_cast<const formula::DoubleVectorRefToken *>(pCur);
- const formula::DoubleVectorRefToken* pDVR1 =
- static_cast<const formula::DoubleVectorRefToken *>(pCur1);
- size_t nCurWindowSize = pDVR->GetRefRowSize();
- size_t nCurWindowSize1 = pDVR1->GetRefRowSize();
- size_t arrayLength = pDVR->GetArrayLength()<
- pDVR1->GetArrayLength() ? pDVR->GetArrayLength():
- pDVR1->GetArrayLength();
- if(nCurWindowSize != nCurWindowSize1)
- {
- ss << " return DBL_MAX;\n";
- ss << "}\n";
- return ;
- }
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- ss << "0; i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- ss << " break;\n";
- ss << " }";
- ss << " return DBL_MAX;\n";
- ss << "}\n";
- return ;
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX) || isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumX += argX;\n";
- ss << " fSumY += argY;\n";
- ss << " fCount += 1.0;\n";
- ss << " }\n";
-
- ss << " if (fCount < 3.0)\n";
- ss << " return DBL_MAX;\n";
- ss << " else\n";
- ss << " {\n";
- ss << " fMeanX = fSumX / fCount;\n";
- ss << " fMeanY = fSumY / fCount;\n";
-
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength ;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())&&
- (!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX)||isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumDeltaXDeltaY +=(argX-fMeanX)*(argY-fMeanY);\n";
- ss << " fSumSqrDeltaX += (argX-fMeanX)*(argX-fMeanX);\n";
- ss << " fSumSqrDeltaY += (argY-fMeanY)*(argY-fMeanY);\n";
- ss << " }\n";
- ss << " if(fSumSqrDeltaX == 0.0)\n";
- ss << " return DBL_MAX;\n";
- ss << " else\n";
- ss << " {\n";
- ss << " return sqrt((fSumSqrDeltaY - fSumDeltaXDeltaY * \n";
- ss << " fSumDeltaXDeltaY / fSumSqrDeltaX)\n";
- ss << " /(fCount - 2.0));\n";
- ss << " }\n";
- ss << " }\n";
- ss << "}\n";
- }
- else
- {
- ss << " return DBL_MAX;\n";
- ss << "}\n";
- }
-}
void OpFisher::GenSlidingWindowFunction( outputstream &ss,
const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -914,196 +627,6 @@ void OpGamma::GenSlidingWindowFunction(
ss << "}\n";
}
-void OpCorrel::GenSlidingWindowFunction(
- outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
-{
- if( vSubArguments.size() !=2 ||vSubArguments[0]->GetFormulaToken()
- ->GetType() != formula::svDoubleVectorRef||vSubArguments[1]
- ->GetFormulaToken()->GetType() != formula::svDoubleVectorRef )
- ///only support DoubleVector in OpCorrelfor GPU calculating.
- throw Unhandled(__FILE__, __LINE__);
- const formula::DoubleVectorRefToken* pCurDVRX =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[0]->GetFormulaToken());
- const formula::DoubleVectorRefToken* pCurDVRY =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[1]->GetFormulaToken());
- if( pCurDVRX->GetRefRowSize() != pCurDVRY->GetRefRowSize() )
- throw Unhandled(__FILE__, __LINE__);
-
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n\t";
- ss << "double vSum = 0.0;\n\t";
- ss << "double vXSum = 0.0;\n\t";
- ss << "double vYSum = 0.0;\n\t";
- ss << "double vXMean = 0.0;\n\t";
- ss << "double vYMean = 0.0;\n\t";
-
- ss << "int gid0 = get_global_id(0);\n\t";
- ss << "double arg0 = 0.0;\n\t";
- ss << "double arg1 = 0.0;\n\t";
- ss << "int cnt = 0;\n\t";
-
- size_t nCurWindowSizeX = pCurDVRY->GetRefRowSize();
-
- ss << "for (int i = ";
- if (!pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "gid0; i < " << nCurWindowSizeX << "; i++) {\n\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t";
- ss << "arg0 = 0.0;\n\t\t\t";
- ss << "arg1 = 0.0;\n\t\t\t";
- ss << "--cnt;\n\t\t";
- ss << "}\n\t\t";
- ss << "++cnt;\n\t\t";
- ss << "vXSum += arg0;\n\t\t";
- ss << "vYSum += arg1;\n\t";
- ss << "}\n\t";
- } else if (pCurDVRX->IsStartFixed() && !pCurDVRX->IsEndFixed()) {
- ss << "0; i < gid0 + " << nCurWindowSizeX << "; i++) {\n\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t";
- ss << "arg0 = 0.0;\n\t\t\t";
- ss << "arg1 = 0.0;\n\t\t\t";
- ss << "--cnt;\n\t\t";
- ss << "}\n\t\t";
- ss << "++cnt;\n\t\t";
- ss << "vXSum += arg0;\n\t\t";
- ss << "vYSum += arg1;\n\t";
- ss << "}\n\t";
- }
- else if (pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "0; i < " << nCurWindowSizeX << "; i++) {\n\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t";
- ss << "arg0 = 0.0;\n\t\t\t";
- ss << "arg1 = 0.0;\n\t\t\t";
- ss << "--cnt;\n\t\t";
- ss << "}\n\t\t";
- ss << "++cnt;\n\t\t";
- ss << "vXSum += arg0;\n\t\t";
- ss << "vYSum += arg1;\n\t";
- ss << "}\n\t";
- } else {
- ss << "0; i < " << nCurWindowSizeX << "; i++) {\n\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i + gid0 >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i + gid0 >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t";
- ss << "arg0 = 0.0;\n\t\t\t";
- ss << "arg1 = 0.0;\n\t\t\t";
- ss << "--cnt;\n\t\t";
- ss << "}\n\t\t";
- ss << "++cnt;\n\t\t";
- ss << "vXSum += arg0;\n\t\t";
- ss << "vYSum += arg1;\n\t";
- ss << "}\n\t";
- }
-
- ss << "if(cnt < 1) {\n\t\t";
- ss << "return DBL_MIN;\n\t";
- ss << "}\n\t";
- ss << "else {\n\t\t";
- ss << "vXMean = vXSum/cnt;\n\t\t";
- ss << "vYMean = vYSum/cnt;\n\t\t";
- ss << "vXSum = 0.0;\n\t\t";
- ss << "vYSum = 0.0;\n\t\t";
-
- ss << "for (int i = ";
- if (!pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "gid0; i < " << nCurWindowSizeX << "; i++) {\n\t\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t\t";
- ss << "arg0 = vXMean;\n\t\t\t\t";
- ss << "arg1 = vYMean;\n\t\t\t";
- ss << "}\n\t\t\t";
- ss << "vXSum += pow(arg0 - vXMean, 2);\n\t\t\t";
- ss << "vYSum += pow(arg1 - vYMean, 2);\n\t\t\t";
- ss << "vSum += (arg0 - vXMean)*(arg1 - vYMean);\n\t\t";
- ss << "}\n\t\t";
- } else if (pCurDVRX->IsStartFixed() && !pCurDVRX->IsEndFixed()) {
- ss << "0; i < gid0 + " << nCurWindowSizeX << "; i++) {\n\t\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t\t";
- ss << "arg0 = vXMean;\n\t\t\t\t";
- ss << "arg1 = vYMean;\n\t\t\t";
- ss << "}\n\t\t\t";
- ss << "vXSum += pow(arg0 - vXMean, 2);\n\t\t\t";
- ss << "vYSum += pow(arg1 - vYMean, 2);\n\t\t\t";
- ss << "vSum += (arg0 - vXMean)*(arg1 - vYMean);\n\t\t";
- ss << "}\n\t\t";
- } else if (pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "0; i < " << nCurWindowSizeX << "; i++) {\n\t\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t\t";
- ss << "arg0 = vXMean;\n\t\t\t\t";
- ss << "arg1 = vYMean;\n\t\t\t";
- ss << "}\n\t\t\t";
- ss << "vXSum += pow(arg0 - vXMean, 2);\n\t\t\t";
- ss << "vYSum += pow(arg1 - vYMean, 2);\n\t\t\t";
- ss << "vSum += (arg0 - vXMean)*(arg1 - vYMean);\n\t\t";
- ss << "}\n\t\t";
- } else {
- ss << "0; i < " << nCurWindowSizeX << "; i++) {\n\t\t\t";
- ss << "arg0 = " << vSubArguments[0]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "arg1 = " << vSubArguments[1]
- ->GenSlidingWindowDeclRef(true) << ";\n\t\t\t";
- ss << "if(isnan(arg0) || isnan(arg1) || (i + gid0 >= ";
- ss << pCurDVRX->GetArrayLength() << ") || (i + gid0 >=";
- ss << pCurDVRY->GetArrayLength() << ")) {\n\t\t\t\t";
- ss << "arg0 = vXMean;\n\t\t\t\t";
- ss << "arg1 = vYMean;\n\t\t\t";
- ss << "}\n\t\t\t";
- ss << "vXSum += ((arg0 - vXMean)*(arg0 - vXMean));\n\t\t\t";
- ss << "vYSum += ((arg1 - vYMean)*(arg1 - vYMean));\n\t\t\t";
- ss << "vSum += (arg0 - vXMean)*(arg1 - vYMean);\n\t\t";
- ss << "}\n\t\t";
- }
-
- ss << "if(vXSum == 0.0 || vYSum == 0.0) {\n\t\t\t";
- ss << "return NAN;\n\t\t";
- ss << "}\n\t\t";
- ss << "else {\n\t\t\t";
- ss << "return vSum/pow(vXSum*vYSum, 0.5);\n\t\t";
- ss << "}\n\t";
- ss << "}\n";
- ss << "}";
-}
-
void OpNegbinomdist::GenSlidingWindowFunction(
outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -1125,95 +648,6 @@ void OpNegbinomdist::GenSlidingWindowFunction(
ss << "}\n";
}
-void OpPearson::GenSlidingWindowFunction(
- outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
-{
- if( vSubArguments.size() !=2 ||vSubArguments[0]->GetFormulaToken()
- ->GetType() != formula::svDoubleVectorRef||vSubArguments[1]
- ->GetFormulaToken()->GetType() != formula::svDoubleVectorRef )
- ///only support DoubleVector in OpPearson for GPU calculating.
- throw Unhandled(__FILE__, __LINE__);
- const formula::DoubleVectorRefToken* pDVR =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[0]->GetFormulaToken());
- const formula::DoubleVectorRefToken* pCurDVRY =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[1]->GetFormulaToken());
- if( pDVR->GetRefRowSize() != pCurDVRY->GetRefRowSize() )
- throw Unhandled(__FILE__, __LINE__);
-
- size_t nCurWindowSize = pDVR->GetRefRowSize();
-
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0=get_global_id(0);\n";
- ss << " double fCount = 0.0;\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fInx;\n";
- ss << " double fIny;\n";
- ss << "for (int i = ";
- if (!pDVR->IsStartFixed() && pDVR->IsEndFixed()) {
- ss << "gid0; i < " << pDVR->GetArrayLength();
- ss << " && i < " << nCurWindowSize << "; i++){\n";
- } else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed()) {
- ss << "0; i < " << pDVR->GetArrayLength();
- ss << " && i < gid0+"<< nCurWindowSize << "; i++){\n";
- } else if (!pDVR->IsStartFixed() && !pDVR->IsEndFixed()){
- ss << "0; i + gid0 < " << pDVR->GetArrayLength();
- ss << " && i < "<< nCurWindowSize << "; i++){\n";
- }
- else {
- ss << "0; i < "<< nCurWindowSize << "; i++){\n";
- }
- ss << " fInx = "<<vSubArguments[0]->GenSlidingWindowDeclRef(true);
- ss << ";\n";
- ss << " fIny = "<<vSubArguments[1]->GenSlidingWindowDeclRef(true);
- ss << " ;\n";
- ss << " if(!isnan(fInx)&&!isnan(fIny)){\n";
- ss << " fSumX += fInx;\n";
- ss << " fSumY += fIny;\n";
- ss << " fCount = fCount + 1;\n";
- ss << " }\n";
- ss << " }\n";
- ss << " if(fCount < 1)\n";
- ss << " return CreateDoubleError(NoValue);\n";
- ss << " double fMeanX = fSumX / fCount;\n";
- ss << " double fMeanY = fSumY / fCount;\n";
- ss << " fSumX = 0.0;\n";
- ss << " fSumY = 0.0;\n";
- ss << "for (int i = ";
- if (!pDVR->IsStartFixed() && pDVR->IsEndFixed()) {
- ss << "gid0; i < " << pDVR->GetArrayLength();
- ss << " && i < " << nCurWindowSize << "; i++){\n";
- } else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed()) {
- ss << "0; i < " << pDVR->GetArrayLength();
- ss << " && i < gid0+"<< nCurWindowSize << "; i++){\n";
- } else if (!pDVR->IsStartFixed() && !pDVR->IsEndFixed()){
- ss << "0; i + gid0 < " << pDVR->GetArrayLength();
- ss << " && i < "<< nCurWindowSize << "; i++){\n";
- }
- else {
- ss << "0; i < "<< nCurWindowSize << "; i++){\n";
- }
- ss << " fInx =
"<<vSubArguments[0]->GenSlidingWindowDeclRef(true);
- ss << " ;\n";
- ss << " fIny =
"<<vSubArguments[1]->GenSlidingWindowDeclRef(true);
- ss << " ;\n";
- ss << " if(!isnan(fInx)&&!isnan(fIny)){\n";
- ss << " fSumDeltaXDeltaY += (fInx - fMeanX) * (fIny -
fMeanY);\n";
- ss << " fSumX += (fInx - fMeanX) * (fInx - fMeanX);\n";
- ss << " fSumY += (fIny - fMeanY) * (fIny - fMeanY);\n";
- ss << " }\n";
- ss << " }\n";
- ss << " if (fSumX == 0 || fSumY == 0)\n";
- ss << " return CreateDoubleError(DivisionByZero);\n";
- ss << " double tmp = ( fSumDeltaXDeltaY / sqrt( fSumX * fSumY));\n";
- ss << " return tmp;\n";
- ss << "}\n";
-}
-
void OpGammaLn::GenSlidingWindowFunction(
outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -1519,87 +953,6 @@ void OpCritBinom::GenSlidingWindowFunction(outputstream&
ss,
ss << "}";
}
-void OpRsq::GenSlidingWindowFunction(
- outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
-{
- if( vSubArguments.size() !=2 ||vSubArguments[0]->GetFormulaToken()
- ->GetType() != formula::svDoubleVectorRef||vSubArguments[1]
- ->GetFormulaToken()->GetType() != formula::svDoubleVectorRef )
- ///only support DoubleVector in OpRsq for GPU calculating.
- throw Unhandled(__FILE__, __LINE__);
- const formula::DoubleVectorRefToken* pCurDVR1 =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[0]->GetFormulaToken());
- const formula::DoubleVectorRefToken* pCurDVR2 =
- static_cast<const formula::DoubleVectorRefToken *>(
- vSubArguments[1]->GetFormulaToken());
- if( pCurDVR1->GetRefRowSize() != pCurDVR2->GetRefRowSize() )
- throw Unhandled(__FILE__, __LINE__);
-
- size_t nCurWindowSize = pCurDVR1->GetRefRowSize();
-
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0=get_global_id(0);\n";
- ss << " double fCount = 0.0;\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fInx;\n";
- ss << " double fIny;\n";
- ss << " double tmp0,tmp1;\n";
-
- ss <<"\n";
-
- ss << " for(int i=0; i<"<<nCurWindowSize<<"; i++)\n";
- ss << " {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef(true);
- ss << "))\n";
- ss << " fInx = 0;\n";
- ss << " else\n";
- ss << " fInx = "<<vSubArguments[0]->GenSlidingWindowDeclRef();
- ss << ";\n";
- ss << " if(isnan(";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef(true);
- ss << "))\n";
- ss << " fIny = 0;\n";
- ss << " else\n";
- ss << " fIny = "<<vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << " ;\n";
- ss << " fSumX += fInx;\n";
- ss << " fSumY += fIny;\n";
- ss << " fCount = fCount + 1;\n";
- ss << " }\n";
- ss << " double fMeanX = fSumX / fCount;\n";
- ss << " double fMeanY = fSumY / fCount;\n";
- ss << " fSumX = 0.0;\n";
- ss << " fSumY = 0.0;\n";
- ss << " for(int i=0; i<"<<nCurWindowSize<<"; i++)\n";
- ss << " {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef(true);
- ss << "))\n";
- ss << " fInx = 0;\n";
- ss << " else\n";
- ss << " fInx = "<<vSubArguments[0]->GenSlidingWindowDeclRef();
- ss << ";\n";
- ss << " if(isnan(";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << "))\n";
- ss << " fIny = 0;\n";
- ss << " else\n";
- ss << " fIny = "<<vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << " ;\n";
- ss << " fSumDeltaXDeltaY += (fInx - fMeanX) * (fIny - fMeanY);\n";
- ss << " fSumX += pow(fInx - fMeanX,2);\n";
- ss << " fSumY += pow(fIny - fMeanY,2);\n";
- ss << " }\n";
- ss << " double tmp = pow( fSumDeltaXDeltaY,2) / (fSumX * fSumY);\n";
- ss << " return tmp ;\n";
- ss << "}\n";
-}
-
void OpChiInv::BinInlineFun(std::set<std::string>& decls,
std::set<std::string>& funs)
{
@@ -2173,160 +1526,6 @@ void OpKurt:: GenSlidingWindowFunction(outputstream &ss,
ss << "}";
}
-void OpIntercept::GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments)
-{
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0 = get_global_id(0);\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fMeanX = 0.0;\n";
- ss << " double fMeanY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fSumSqrDeltaX = 0.0;\n";
- ss << " double fCount = 0.0;\n";
- ss << " double argX = 0.0;\n";
- ss << " double argY = 0.0;\n";
- if(vSubArguments.size() != 2)
- {
- ss << " return NAN;\n";
- ss << "}\n";
- return ;
- }
- FormulaToken *pCur = vSubArguments[1]->GetFormulaToken();
- FormulaToken *pCur1 = vSubArguments[0]->GetFormulaToken();
- assert(pCur);
- assert(pCur1);
- if (pCur->GetType() == formula::svDoubleVectorRef&&
- pCur1->GetType() == formula::svDoubleVectorRef)
- {
- const formula::DoubleVectorRefToken* pDVR =
- static_cast<const formula::DoubleVectorRefToken *>(pCur);
- const formula::DoubleVectorRefToken* pDVR1 =
- static_cast<const formula::DoubleVectorRefToken *>(pCur1);
-
- size_t nCurWindowSize = pDVR->GetRefRowSize();
- size_t nCurWindowSize1 = pDVR1->GetRefRowSize();
- size_t arrayLength = pDVR->GetArrayLength()<
- pDVR1->GetArrayLength() ? pDVR->GetArrayLength():
- pDVR1->GetArrayLength();
- if(nCurWindowSize != nCurWindowSize1)
- {
- ss << " return NAN;\n";
- ss << "}\n";
- return ;
- }
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength ;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- ss << "0; i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- ss << " break;\n";
- ss << " }";
- ss << " return NAN;\n";
- ss << "}\n";
- return ;
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX) || isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumX += argX;\n";
- ss << " fSumY += argY;\n";
- ss << " fCount += 1.0;\n";
- ss << " }\n";
-
- ss << " if (fCount < 1.0)\n";
- ss << " return NAN;\n";
- ss << " else\n";
- ss << " {\n";
- ss << " fMeanX = fSumX / fCount;\n";
- ss << " fMeanY = fSumY / fCount;\n";
-
- ss << " for (int i = ";
- if ((!pDVR->IsStartFixed() && pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && pDVR1->IsEndFixed()))
- {
- ss << "gid0; i < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i < " << arrayLength ;
- ss << " && i < gid0+" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else if ((!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
- &&(!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed()))
- {
- ss << "0; i + gid0 < " << arrayLength;
- ss << " && i < " << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- }
- else
- {
- ss << "0; i < " << arrayLength << "; i++)\n";
- ss << " {\n";
- }
-
- ss << " argX = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " argY = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if (isnan(argX) || isnan(argY))\n";
- ss << " continue;\n";
- ss << " fSumDeltaXDeltaY += (argX-fMeanX)*(argY-fMeanY);\n";
- ss << " fSumSqrDeltaX += (argX-fMeanX) * (argX-fMeanX);\n";
- ss << " }\n";
- ss << " if(fSumSqrDeltaX == 0.0)\n";
- ss << " return NAN;\n";
- ss << " else\n";
- ss << " {\n";
- ss << " return fMeanY -";
- ss << " (fSumDeltaXDeltaY/fSumSqrDeltaX)*fMeanX;\n";
- ss << " }\n";
- ss << " }\n";
- ss << "}\n";
- }
- else
- {
- ss << " return NAN;\n";
- ss << "}\n";
- }
-}
void OpLogInv:: GenSlidingWindowFunction(outputstream &ss,
const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -2502,111 +1701,6 @@ void OpLogInv:: GenSlidingWindowFunction(outputstream
&ss,
ss << "}\n";
}
-void OpForecast::GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments)
-{
- FormulaToken *pCur0 = vSubArguments[0]->GetFormulaToken();
- assert(pCur0);
- const formula::SingleVectorRefToken*pCurDVR0= static_cast<const
- formula::SingleVectorRefToken *>(pCur0);
- FormulaToken *pCur1 = vSubArguments[1]->GetFormulaToken();
- assert(pCur1);
- const formula::DoubleVectorRefToken* pCurDVR1 =
- static_cast<const formula::DoubleVectorRefToken *>(pCur1);
- size_t nCurWindowSize = pCurDVR1->GetRefRowSize();
- FormulaToken *pCur2 = vSubArguments[2]->GetFormulaToken();
- assert(pCur2);
- const formula::DoubleVectorRefToken* pCurDVR2 =
- static_cast<const formula::DoubleVectorRefToken *>(pCur2);
- size_t nCurWindowSize1 = pCurDVR2->GetRefRowSize();
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0 = get_global_id(0);\n";
- ss << " double fSumX = 0.0;\n";
- ss << " double fSumY = 0.0;\n";
- ss << " double fSumDeltaXDeltaY = 0.0;\n";
- ss << " double fSumSqrDeltaX = 0.0;\n";
- if(pCur0->GetType()== formula::svDouble ||
- pCur0->GetType() == formula::svSingleVectorRef)
- {
- ss << " double arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef();
- ss << ";\n";
- }
- else
- throw Unhandled(__FILE__, __LINE__);
- if(pCur1->GetType() != formula::svDoubleVectorRef ||
- pCur2->GetType() != formula::svDoubleVectorRef)
- throw Unhandled(__FILE__, __LINE__);
- else
- {
- ss<< " if(isnan(arg0)||(gid0>=";
- ss<<pCurDVR0->GetArrayLength();
- ss<<"))\n";
- ss<<" arg0 = 0;\n";
- ss << " int length="<<nCurWindowSize;
- ss << ";\n";
- ss << " int length1= "<<nCurWindowSize1;
- ss << ";\n";
- ss << " if(length!=length1)\n";
- ss << " return 0;\n";
- ss << " double tmp = 0;\n";
- ss << " for (int i = 0; i <" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- ss << " double arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef(true);
- ss << ";\n";
- ss << " double arg2 = ";
- ss << vSubArguments[2]->GenSlidingWindowDeclRef(true);
- ss << ";\n";
- ss << " if(isnan(arg1)||((gid0+i)>=";
- ss << pCurDVR1->GetArrayLength();
- ss << "))\n";
- ss << " {\n";
- ss << " length--;\n";
- ss << " continue;\n";
- ss << " }\n";
- ss << " if(isnan(arg2)||((gid0+i)>=";
- ss << pCurDVR2->GetArrayLength();
- ss << "))\n";
- ss << " {\n";
- ss << " length--;\n";
- ss << " continue;\n";
- ss << " }\n";
- ss << " fSumY+=arg1;\n";
- ss << " fSumX+=arg2;\n";
- ss << " }\n";
- ss << " double fMeanX = fSumX / length;\n";
- ss << " double fMeanY = fSumY / length;\n";
- ss << " for (int i = 0; i <" << nCurWindowSize << "; i++)\n";
- ss << " {\n";
- ss << " double arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef(true);
- ss << ";\n";
- ss << " double arg2 = ";
- ss << vSubArguments[2]->GenSlidingWindowDeclRef(true);
- ss << ";\n";
- ss << " if(isnan(arg1)||((gid0+i)>=";
- ss <<pCurDVR1->GetArrayLength();
- ss <<"))\n";
- ss <<" {\n";
- ss <<" continue;\n";
- ss <<" }\n";
- ss << " if(isnan(arg2)||((gid0+i)>=";
- ss <<pCurDVR2->GetArrayLength();
- ss <<"))\n";
- ss <<" {\n";
- ss <<" continue;\n";
- ss <<" }\n";
- ss <<" fSumDeltaXDeltaY+=(arg2 - fMeanX) * (arg1 - fMeanY);\n";
- ss <<" fSumSqrDeltaX+=pow(arg2 - fMeanX, 2);\n";
- ss <<" }\n";
- ss <<" tmp =fMeanY + fSumDeltaXDeltaY / fSumSqrDeltaX *";
- ss <<" (arg0 - fMeanX);\n";
- ss <<" return tmp;\n";
- ss << "}";
- }
-}
void OpLogNormDist::GenSlidingWindowFunction(outputstream &ss,
const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -3381,283 +2475,7 @@ void OpPoisson::GenSlidingWindowFunction(
ss << " }\n";
ss << "}\n";
}
-void OpCovar::GenSlidingWindowFunction(outputstream& ss,
- const std::string &sSymName, SubArguments& vSubArguments)
-{
- CHECK_PARAMETER_COUNT(2,2);
- GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
- ss << "{\n";
- ss << " int gid0 = get_global_id(0);\n";
- ss << " double vSum = 0.0;\n";
- ss << " double vSum0 = 0.0;\n";
- ss << " double vSum1 = 0.0;\n";
- ss << " double vMean0 = 0.0;\n";
- ss << " double vMean1 = 0.0;\n";
- ss << " double arg0 = 0.0;\n";
- ss << " double arg1 = 0.0;\n";
- FormulaToken* pCurX = vSubArguments[0]->GetFormulaToken();
- FormulaToken* pCurY = vSubArguments[1]->GetFormulaToken();
- if ((pCurX->GetType() == formula::svDoubleVectorRef)&&
- (pCurY->GetType() == formula::svDoubleVectorRef)){
- ss << " int cnt = 0;\n";
- const formula::DoubleVectorRefToken* pCurDVRX =
- static_cast<const formula::DoubleVectorRefToken* >(pCurX);
- const formula::DoubleVectorRefToken* pCurDVRY =
- static_cast<const formula::DoubleVectorRefToken* >(pCurY);
- size_t nCurWindowSizeX = pCurDVRX->GetRefRowSize();
- size_t nCurWindowSizeY = pCurDVRY->GetRefRowSize();
- if(nCurWindowSizeX == nCurWindowSizeY)
- {
- ss << " for( ";
- if (!pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "int i = gid0; i < " << nCurWindowSizeX;
- ss << " && i < " << pCurDVRX->GetArrayLength() << " && i < ";
- ss << pCurDVRY->GetArrayLength() << "; i++){\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") ||";
- ss << " isnan("<< vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")) {\n";
- ss << " arg0 = 0.0;\n";
- ss << " arg1 = 0.0;\n";
- ss << " --cnt;\n";
- ss << " }\n";
- ss << "else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << "}\n";
- ss << " ++cnt;\n";
- ss << " vSum0 += arg0;\n";
- ss << " vSum1 += arg1;\n";
- ss << " }\n";
- }
- else if (pCurDVRX->IsStartFixed() && !pCurDVRX->IsEndFixed()) {
- ss << "int i = 0; i < gid0 + " << nCurWindowSizeX << " && ";
- ss << " i < " << pCurDVRX->GetArrayLength() << " && ";
- ss << " i < " << pCurDVRY->GetArrayLength() << "; i++) {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") ||";
- ss << " isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")) {\n";
- ss << " arg0 = 0.0;\n";
- ss << " arg1 = 0.0;\n";
- ss << " --cnt;\n";
- ss << " }\n";
- ss << "else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";}\n";
- ss << " ++cnt;\n";
- ss << " vSum0 += arg0;\n";
- ss << " vSum1 += arg1;\n";
- ss << " }\n";
- }
- else if (pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "int i = 0; i < " << nCurWindowSizeX << " && i < ";
- ss << pCurDVRX->GetArrayLength() << " && i < ";
- ss << pCurDVRY->GetArrayLength() << "; i++) {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") ||";
- ss << " isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")) {\n";
- ss << " arg0 = 0.0;\n";
- ss << " arg1 = 0.0;\n";
- ss << " --cnt;\n";
- ss << " }\n";
- ss << "else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";}\n";
- ss << " ++cnt;\n";
- ss << " vSum0 += arg0;\n";
- ss << " vSum1 += arg1;\n";
- ss << " }\n";
- }
- else {
- ss << "int i = 0; i < " << nCurWindowSizeX << " && ";
- ss << " i + gid0 < " << pCurDVRX->GetArrayLength();
- ss << " && i + gid0 < " << pCurDVRY->GetArrayLength();
- ss << "; i++) {\n";
- ss << "if ((isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ")) || ";
- ss << "(isnan("<< vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << "))) {\n";
- ss << " arg0 = 0.0;\n";
- ss << " arg1 = 0.0;\n";
- ss << " --cnt;\n";
- ss << " }\n";
- ss << " else {\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " }\n";
- ss << " ++cnt;\n";
- ss << " vSum0 += arg0;\n";
- ss << " vSum1 += arg1;\n";
- ss << " }\n";
- }
- ss << " if(cnt < 1) {\n";
- ss << " return CreateDoubleError(NoValue);\n";
- ss << " }\n";
- ss << " else {\n";
- ss << " vMean0 = vSum0 / cnt;\n";
- ss << " vMean1 = vSum1 / cnt;\n";
- ss << " for(";
- if (!pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "int i = gid0; i < " << nCurWindowSizeX;
- ss << " && i < " << pCurDVRX->GetArrayLength() << " && i < ";
- ss << pCurDVRY->GetArrayLength() << "; i++){\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") ||";
- ss << " isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")){\n";
- ss << " arg0 = vMean0;\n";
- ss << " arg1 = vMean1;\n";
- ss << " }\n";
- ss << " else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";}\n";
- ss << " vSum += (arg0 - vMean0) * (arg1 - vMean1);\n";
- ss << " }\n";
- }
- else if (pCurDVRX->IsStartFixed() && !pCurDVRX->IsEndFixed()) {
- ss << "int i = 0; i < gid0 + " << nCurWindowSizeX << " && ";
- ss << " i < " << pCurDVRX->GetArrayLength() << " && ";
- ss << " i < " << pCurDVRY->GetArrayLength() << "; i++) {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") || ";
- ss << "isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")) {\n";
- ss << " arg0 = vMean0;\n";
- ss << " arg1 = vMean1;\n";
- ss << " }\n";
- ss << "else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";}\n";
- ss << " vSum += (arg0 - vMean0) * (arg1 - vMean1);\n";
- ss << " }\n";
- }
- else if (pCurDVRX->IsStartFixed() && pCurDVRX->IsEndFixed()) {
- ss << "int i = 0; i < " << nCurWindowSizeX << " && i < ";
- ss << pCurDVRX->GetArrayLength() << " && i < ";
- ss << pCurDVRY->GetArrayLength() << "; i++) {\n";
- ss << " if(isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ") || ";
- ss << "isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << ")) {\n";
- ss << " arg0 = vMean0;\n";
- ss << " arg1 = vMean1;\n";
- ss << " }\n";
- ss << "else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";}\n";
- ss << " vSum += (arg0 - vMean0) * (arg1 - vMean1);\n";
- ss << " }\n";
- }
- else {
- ss << "int i = 0; i < " << nCurWindowSizeX << " && ";
- ss << " i + gid0 < " << pCurDVRX->GetArrayLength();
- ss << " && i + gid0 < " << pCurDVRX->GetArrayLength();
- ss << "; i++) {\n";
- ss << "if((isnan(";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ")) || ";
- ss << "(isnan(" << vSubArguments[1]->GenSlidingWindowDeclRef();
- ss << "))) {\n";
- ss << " arg0 = vMean0;\n";
- ss << " arg1 = vMean1;\n";
- ss << " }\n";
- ss << " else{\n";
- ss << " arg0 = ";
- ss << vSubArguments[0]->GenSlidingWindowDeclRef() << ";\n";
- ss << " arg1 = ";
- ss << vSubArguments[1]->GenSlidingWindowDeclRef() << ";\n";
- ss << " }\n";
- ss << " vSum += (arg0 - vMean0) * (arg1 - vMean1);\n";
- ss << " }\n";
- }
- ss << " return vSum / cnt;\n";
- ss << " }\n";
- ss << "}";
- }
- else
- {
- throw Unhandled(__FILE__, __LINE__);
- }
- }
- else {
- ss << " int cnt0 = 0,cnt1 = 0;\n";
- for (size_t i = 0; i < vSubArguments.size(); i++)
- {
- FormulaToken* pCur = vSubArguments[i]->GetFormulaToken();
- if (pCur->GetType() == formula::svSingleVectorRef){
- const formula::SingleVectorRefToken* pTVR =
- static_cast< const formula::SingleVectorRefToken* >(pCur);
- ss << " if(isnan(";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef();
- ss << ") || gid0 >= " << pTVR->GetArrayLength() << ")\n";
- ss << " arg" << i << " = 0;\n else\n";
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- ss << " cnt" << i << "++;\n";
- ss << " vSum" << i << " += arg" << i << ";\n";
- }
- else if (pCur->GetType() == formula::svDouble){
- ss << " if(isnan ( ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << "))\n";
- ss << " arg" << i << " = 0;\n else\n";
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- ss << " cnt" << i << "++;\n";
- ss << " vSum" << i << " += arg" << i << ";\n";
- }
- else {
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- ss << " cnt" << i << "++;\n";
- ss << " vSum" << i << " += arg" << i << ";\n";
- }
- }
- ss << " vMean0 = vSum0 / cnt0;\n";
- ss << " vMean1 = vSum0 / cnt1;\n";
- for(size_t i = 0; i < vSubArguments.size(); i++ ) {
- FormulaToken* pCur = vSubArguments[i]->GetFormulaToken();
- if (pCur->GetType() == formula::svSingleVectorRef) {
- const formula::SingleVectorRefToken* pTVR =
- static_cast< const formula::SingleVectorRefToken* >(pCur);
- ss << " if(isnan(";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef();
- ss << ") || gid0 >= " << pTVR->GetArrayLength() << ")\n";
- ss << " arg" << i << " = vMean" << i << ";\n";
- ss << " else\n";
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- }
- else if (pCur->GetType() == formula::svDouble) {
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- ss << " if(isnan(arg" << i << "))\n";
- ss << " arg" << i << " = vMean" << i << ";\n";
- }
- else {
- ss << " arg" << i << " = ";
- ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
- }
- }
- ss << " vSum += (arg0 - vMean0) * ( arg1 - vMean1 );\n";
- ss << " return vSum / cnt0;\n";
- ss << "}";
- }
-}
+
void OpBetaDist::GenSlidingWindowFunction(outputstream &ss,
const std::string &sSymName, SubArguments &vSubArguments)
{
@@ -3866,6 +2684,221 @@ void OpAveDev:: GenSlidingWindowFunction(outputstream
&ss,
ss << "}";
}
+// Block of functions OpCovar-OpRsq that are rather similar.
+
+void OpCovar::GenSlidingWindowFunction(outputstream& ss,
+ const std::string &sSymName, SubArguments& vSubArguments)
+{
+ CHECK_PARAMETER_COUNT( 2, 2 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 0 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 1 );
+ GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
+ ss << "{\n";
+ ss << " int gid0 = get_global_id(0);\n";
+ ss << " double fSumX = 0.0;\n";
+ ss << " double fSumY = 0.0;\n";
+ ss << " double fMeanX = 0.0;\n";
+ ss << " double fMeanY = 0.0;\n";
+ ss << " double fSumDeltaXDeltaY = 0.0;\n";
+ ss << " double fCount = 0.0;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if (!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ " fSumX += arg1;\n"
+ " fSumY += arg2;\n"
+ " fCount += 1.0;\n"
+ " }\n"
+ );
+ ss << " if( fCount < 1 )\n";
+ ss << " return CreateDoubleError(NoValue);\n";
+ ss << " fMeanX = fSumX / fCount;\n";
+ ss << " fMeanY = fSumY / fCount;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if(!isnan(arg1) && !isnan(arg2))\n"
+ " fSumDeltaXDeltaY +=(arg1-fMeanX)*(arg2-fMeanY);\n"
+ );
+ ss << " return fSumDeltaXDeltaY / fCount;\n";
+ ss << "}\n";
+}
+
+void OpForecast::GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments)
+{
+ CHECK_PARAMETER_COUNT( 3, 3 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 1 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 2 );
+ GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
+ ss << "{\n";
+ ss << " int gid0 = get_global_id(0);\n";
+ ss << " double fSumX = 0.0;\n";
+ ss << " double fSumY = 0.0;\n";
+ ss << " double fMeanX = 0.0;\n";
+ ss << " double fMeanY = 0.0;\n";
+ ss << " double fSumDeltaXDeltaY = 0.0;\n";
+ ss << " double fSumSqrDeltaX = 0.0;\n";
+ ss << " double fCount = 0.0;\n";
+ GenerateArg( "arg0", 0, vSubArguments, ss );
+ GenerateRangeArgPair( 1, 2, vSubArguments, ss,
+ " if (!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ // note that arg1 -> Y, arg2 -> X
+ " fSumX += arg2;\n"
+ " fSumY += arg1;\n"
+ " fCount += 1.0;\n"
+ " }\n"
+ );
+ ss << " if( fCount < 1 )\n";
+ ss << " return CreateDoubleError(NoValue);\n";
+ ss << " fMeanX = fSumX / fCount;\n";
+ ss << " fMeanY = fSumY / fCount;\n";
+ GenerateRangeArgPair( 1, 2, vSubArguments, ss,
+ " if(!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ " fSumDeltaXDeltaY +=(arg2-fMeanX)*(arg1-fMeanY);\n"
+ " fSumSqrDeltaX += (arg2-fMeanX)*(arg2-fMeanX);\n"
+ " }\n"
+ );
+ ss << " if(fSumSqrDeltaX == 0.0)\n";
+ ss << " return CreateDoubleError(DivisionByZero);\n";
+ ss << " return fMeanY + fSumDeltaXDeltaY / fSumSqrDeltaX * (arg0 -
fMeanX);\n";
+ ss << "}\n";
+}
+
+void OpInterceptSlopeBase::GenerateCode( outputstream &ss, const std::string
&sSymName,
+ SubArguments &vSubArguments, const char* finalComputeCode )
+{
+ CHECK_PARAMETER_COUNT( 2, 2 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 0 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 1 );
+ GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
+ ss << "{\n";
+ ss << " int gid0 = get_global_id(0);\n";
+ ss << " double fSumX = 0.0;\n";
+ ss << " double fSumY = 0.0;\n";
+ ss << " double fMeanX = 0.0;\n";
+ ss << " double fMeanY = 0.0;\n";
+ ss << " double fSumDeltaXDeltaY = 0.0;\n";
+ ss << " double fSumSqrDeltaX = 0.0;\n";
+ ss << " double fCount = 0.0;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if (!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ // note that arg1 -> Y, arg2 -> X
+ " fSumX += arg2;\n"
+ " fSumY += arg1;\n"
+ " fCount += 1.0;\n"
+ " }\n"
+ );
+ ss << " if( fCount < 1 )\n";
+ ss << " return CreateDoubleError(NoValue);\n";
+ ss << " fMeanX = fSumX / fCount;\n";
+ ss << " fMeanY = fSumY / fCount;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if(!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ " fSumDeltaXDeltaY +=(arg2-fMeanX)*(arg1-fMeanY);\n"
+ " fSumSqrDeltaX += (arg2-fMeanX)*(arg2-fMeanX);\n"
+ " }\n"
+ );
+ ss << finalComputeCode;
+ ss << "}\n";
+}
+
+void OpIntercept::GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments)
+{
+ GenerateCode( ss, sSymName, vSubArguments,
+ " if(fSumSqrDeltaX == 0.0)\n"
+ " return CreateDoubleError(DivisionByZero);\n"
+ " return fMeanY - (fSumDeltaXDeltaY/fSumSqrDeltaX)*fMeanX;\n"
+ );
+}
+
+void OpSlope::GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments)
+{
+ GenerateCode( ss, sSymName, vSubArguments,
+ " if(fSumSqrDeltaX == 0.0)\n"
+ " return CreateDoubleError(DivisionByZero);\n"
+ " return fSumDeltaXDeltaY / fSumSqrDeltaX;\n"
+ );
+}
+
+void OpPearsonCovarBase::GenerateCode( outputstream &ss, const std::string
&sSymName,
+ SubArguments &vSubArguments, double minimalCountValue, const char*
finalComputeCode )
+{
+ CHECK_PARAMETER_COUNT( 2, 2 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 0 );
+ CHECK_PARAMETER_DOUBLEVECTORREF( 1 );
+ GenerateFunctionDeclaration( sSymName, vSubArguments, ss );
+ ss << "{\n";
+ ss << " int gid0 = get_global_id(0);\n";
+ ss << " double fSumX = 0.0;\n";
+ ss << " double fSumY = 0.0;\n";
+ ss << " double fMeanX = 0.0;\n";
+ ss << " double fMeanY = 0.0;\n";
+ ss << " double fSumDeltaXDeltaY = 0.0;\n";
+ ss << " double fSumSqrDeltaX = 0.0;\n";
+ ss << " double fSumSqrDeltaY = 0.0;\n";
+ ss << " double fCount = 0.0;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if (!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ // note that arg1 -> Y, arg2 -> X
+ " fSumX += arg2;\n"
+ " fSumY += arg1;\n"
+ " fCount += 1.0;\n"
+ " }\n"
+ );
+ ss << " if( fCount < " << minimalCountValue <<" )\n";
+ ss << " return CreateDoubleError(NoValue);\n";
+ ss << " fMeanX = fSumX / fCount;\n";
+ ss << " fMeanY = fSumY / fCount;\n";
+ GenerateRangeArgPair( 0, 1, vSubArguments, ss,
+ " if(!isnan(arg1) && !isnan(arg2))\n"
+ " {\n"
+ " fSumDeltaXDeltaY +=(arg2-fMeanX)*(arg1-fMeanY);\n"
+ " fSumSqrDeltaX += (arg2-fMeanX)*(arg2-fMeanX);\n"
+ " fSumSqrDeltaY += (arg1-fMeanY)*(arg1-fMeanY);\n"
+ " }\n"
+ );
+ ss << finalComputeCode;
+ ss << "}\n";
+}
+
+void OpPearson::GenSlidingWindowFunction(
+ outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
+{
+ GenerateCode( ss, sSymName, vSubArguments, 1,
+ " if (fSumSqrDeltaX == 0 || fSumSqrDeltaY == 0)\n"
+ " return CreateDoubleError(DivisionByZero);\n"
+ " return ( fSumDeltaXDeltaY / sqrt( fSumSqrDeltaX *
fSumSqrDeltaY));\n"
+ );
+}
+
+void OpSTEYX::GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments)
+{
+ GenerateCode( ss, sSymName, vSubArguments, 3,
+ " if(fSumSqrDeltaX == 0.0)\n"
+ " return CreateDoubleError(DivisionByZero);\n"
+ " return sqrt((fSumSqrDeltaY - fSumDeltaXDeltaY * \n"
+ " fSumDeltaXDeltaY / fSumSqrDeltaX)\n"
+ " /(fCount - 2.0));\n"
+ );
+}
+
+void OpRsq::GenSlidingWindowFunction(
+ outputstream &ss, const std::string &sSymName, SubArguments &vSubArguments)
+{
+ // This is basically pow(OpPearson,2)
+ GenerateCode( ss, sSymName, vSubArguments, 1,
+ " if (fSumSqrDeltaX == 0 || fSumSqrDeltaY == 0)\n"
+ " return CreateDoubleError(DivisionByZero);\n"
+ " return ( fSumDeltaXDeltaY * fSumDeltaXDeltaY / (fSumSqrDeltaX *
fSumSqrDeltaY));\n"
+ );
+}
+
namespace {
enum MixDoubleString
diff --git a/sc/source/core/opencl/op_statistical.hxx
b/sc/source/core/opencl/op_statistical.hxx
index a433acdbf56d..64089b54f524 100644
--- a/sc/source/core/opencl/op_statistical.hxx
+++ b/sc/source/core/opencl/op_statistical.hxx
@@ -44,13 +44,6 @@ public:
virtual std::string BinFuncName() const override { return "VarP"; }
virtual bool canHandleMultiVector() const override { return true; }
};
-class OpSTEYX: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "STEYX"; }
-};
class OpZTest: public Normal
{
public:
@@ -92,13 +85,6 @@ public:
virtual std::string BinFuncName() const override { return "Skewp"; }
virtual bool canHandleMultiVector() const override { return true; }
};
-class OpSlope: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "Slope"; }
-};
class OpWeibull: public Normal
{
public:
@@ -162,14 +148,6 @@ public:
virtual std::string BinFuncName() const override { return "Gamma"; }
};
-class OpCorrel: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "Correl"; }
-};
-
class OpNegbinomdist: public Normal
{
public:
@@ -178,14 +156,6 @@ public:
virtual std::string BinFuncName() const override { return
"OpNegbinomdist"; }
};
-class OpPearson: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "OpPearson"; }
-};
-
class OpGammaLn: public Normal
{
public:
@@ -222,13 +192,6 @@ public:
virtual bool canHandleMultiVector() const override { return true; }
};
-class OpRsq: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "OpRsq"; }
-};
class OpNormdist:public Normal{
public:
virtual void GenSlidingWindowFunction(outputstream &ss,
@@ -273,13 +236,6 @@ public:
virtual std::string BinFuncName() const override { return "Kurt"; }
virtual bool canHandleMultiVector() const override { return true; }
};
-class OpCovar: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream& ss,
- const std::string &sSymName, SubArguments& vSubArguments) override;
- virtual std::string BinFuncName() const override { return "Covar"; }
-};
class OpPermut:public Normal{
public:
@@ -306,13 +262,6 @@ public:
virtual std::string BinFuncName() const override { return "Confidence"; }
};
-class OpIntercept: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "Intercept"; }
-};
class OpLogInv: public Normal
{
public:
@@ -332,13 +281,6 @@ public:
virtual std::string BinFuncName() const override { return "CritBinom"; }
};
-class OpForecast: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "Forecast"; }
-};
class OpLogNormDist: public Normal
{
public:
@@ -466,6 +408,98 @@ class OpBetainv:public Normal{
virtual void BinInlineFun(std::set<std::string>& ,std::set<std::string>&)
override;
virtual std::string BinFuncName() const override { return "OpBetainv"; }
};
+
+class OpAveDev: public Normal
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "AveDev"; }
+ virtual bool canHandleMultiVector() const override { return true; }
+};
+
+class OpCovar: public Normal
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream& ss,
+ const std::string &sSymName, SubArguments& vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "Covar"; }
+};
+
+class OpForecast: public Normal
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "Forecast"; }
+};
+
+class OpInterceptSlopeBase: public Normal
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override
= 0;
+protected:
+ void GenerateCode( outputstream& ss, const std::string &sSymName,
+ SubArguments &vSubArguments, const char* finalComputeCode );
+};
+
+class OpIntercept: public OpInterceptSlopeBase
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "Intercept"; }
+};
+
+class OpSlope: public OpInterceptSlopeBase
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "Slope"; }
+};
+
+class OpPearsonCovarBase: public Normal
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override
= 0;
+protected:
+ void GenerateCode( outputstream& ss, const std::string &sSymName,
+ SubArguments &vSubArguments, double minimalCountValue, const char*
finalComputeCode );
+};
+
+class OpPearson: public OpPearsonCovarBase
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "OpPearson"; }
+};
+
+class OpCorrel: public OpPearson // they are identical
+{
+public:
+ virtual std::string BinFuncName() const override { return "Correl"; }
+};
+
+class OpSTEYX: public OpPearsonCovarBase
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "STEYX"; }
+};
+
+class OpRsq: public OpPearsonCovarBase
+{
+public:
+ virtual void GenSlidingWindowFunction(outputstream &ss,
+ const std::string &sSymName, SubArguments &vSubArguments) override;
+ virtual std::string BinFuncName() const override { return "OpRsq"; }
+};
+
class OpMinA: public Normal
{
public:
@@ -544,14 +578,6 @@ public:
virtual bool takeString() const override { return true; }
virtual bool takeNumeric() const override { return true; }
};
-class OpAveDev: public Normal
-{
-public:
- virtual void GenSlidingWindowFunction(outputstream &ss,
- const std::string &sSymName, SubArguments &vSubArguments) override;
- virtual std::string BinFuncName() const override { return "AveDev"; }
- virtual bool canHandleMultiVector() const override { return true; }
-};
}
diff --git a/sc/source/core/opencl/opbase.cxx b/sc/source/core/opencl/opbase.cxx
index c73e50b3c593..f66dede91c59 100644
--- a/sc/source/core/opencl/opbase.cxx
+++ b/sc/source/core/opencl/opbase.cxx
@@ -347,22 +347,23 @@ void SlidingFunctionBase::GenerateRangeArgPair( int arg1,
int arg2, SubArguments
throw Unhandled( __FILE__, __LINE__ );
}
+ size_t arrayLength = std::min( pDVR1->GetArrayLength(),
pDVR2->GetArrayLength());
ss << " for (int i = ";
if (!pDVR1->IsStartFixed() && pDVR1->IsEndFixed())
{
- ss << "gid0; i < " << pDVR1->GetArrayLength();
+ ss << "gid0; i < " << arrayLength;
ss << " && i < " << nCurWindowSize1 << "; i++)\n";
ss << " {\n";
}
else if (pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())
{
- ss << "0; i < " << pDVR1->GetArrayLength();
+ ss << "0; i < " << arrayLength;
ss << " && i < gid0+"<< nCurWindowSize1 << "; i++)\n";
ss << " {\n";
}
else if (!pDVR1->IsStartFixed() && !pDVR1->IsEndFixed())
{
- ss << "0; i + gid0 < " << pDVR1->GetArrayLength();
+ ss << "0; i + gid0 < " << arrayLength;
ss << " && i < " << nCurWindowSize1 << "; i++)\n";
ss << " {\n";
}
diff --git a/sc/source/core/opencl/opbase.hxx b/sc/source/core/opencl/opbase.hxx
index c9ab020d6689..9ad6e7480cf0 100644
--- a/sc/source/core/opencl/opbase.hxx
+++ b/sc/source/core/opencl/opbase.hxx
@@ -92,6 +92,12 @@ public:
if( count < ( min )) \
throw InvalidParameterCount( count, __FILE__, __LINE__ ); \
} while( false )
+#define CHECK_PARAMETER_DOUBLEVECTORREF(arg) \
+ do { \
+ FormulaToken *token = vSubArguments[arg]->GetFormulaToken(); \
+ if (token == nullptr || token->GetType() !=
formula::svDoubleVectorRef) \
+ throw Unhandled(__FILE__, __LINE__); \
+ } while( false )
typedef std::shared_ptr<FormulaTreeNode> FormulaTreeNodeRef;
