GPU Calc: implemented for DSTDEVP

AMLOEXT-340 FIX

Change-Id: I180595e08858b91fc84c58abfaa1a7b2b05b199b
Signed-off-by: haochen <haochen@multicorewareinc.com>
Signed-off-by: Wei Wei <weiwei@multicorewareinc.com>
Signed-off-by: I-Jui (Ray) Sung <ray@multicorewareinc.com>

3 files changed, 217 insertions, 0 deletions

diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx
index b88ce9722486..f13ee029c085 100644
--- a/sc/source/core/opencl/formulagroupcl.cxx
+++ b/sc/source/core/opencl/formulagroupcl.cxx
@@ -2392,6 +2392,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments(
                 mvSubArguments.push_back(SoPHelper(ts,
                          ft->Children[i], new OpDstdev));
                  break;
+            case ocDBStdDevP:
+                mvSubArguments.push_back(SoPHelper(ts,
+                         ft->Children[i], new OpDstdevp));
+                 break;
             case ocExternal:
                 if ( !(pChild->GetExternal().compareTo(OUString(
                     "com.sun.star.sheet.addin.Analysis.getEffect"))))
diff --git a/sc/source/core/opencl/op_database.cxx b/sc/source/core/opencl/op_database.cxx
index 40fd25b14916..2907a22de41f 100644
--- a/sc/source/core/opencl/op_database.cxx
+++ b/sc/source/core/opencl/op_database.cxx
@@ -790,6 +790,212 @@ void OpDstdev::GenSlidingWindowFunction(std::stringstream &ss,
     ss << "}";
 }
 
+void OpDstdevp::GenSlidingWindowFunction(std::stringstream &ss,
+            const std::string sSymName, SubArguments &vSubArguments)
+{
+    ss << "\ndouble " << sSymName;
+    ss << "_"<< BinFuncName() <<"(";
+    for (unsigned i = 0; i < vSubArguments.size(); i++)
+    {
+        if (i)
+            ss << ",";
+        vSubArguments[i]->GenSlidingWindowDecl(ss);
+    }
+    ss << ")\n    {\n";
+    ss << "    int gid0=get_global_id(0);\n";
+    ss << "    double var = 0;\n";
+    ss << "    double mean = 0;\n";
+    ss << "    double value =0;\n";
+    ss << "    int count = 0;\n";
+    GenTmpVariables(ss,vSubArguments);
+    int dataCol = 0;
+    int dataRow = 0;
+    if(vSubArguments[0]->GetFormulaToken()->GetType() ==
+    formula::svDoubleVectorRef)
+    {
+        formula::FormulaToken *tmpCur = vSubArguments[0]->GetFormulaToken();
+        const formula::DoubleVectorRefToken*pCurDVR1= dynamic_cast<const
+            formula::DoubleVectorRefToken *>(tmpCur);
+        dataCol = pCurDVR1->GetArrays().size();
+        dataRow = pCurDVR1->GetArrayLength();
+    }
+    else
+        throw Unhandled();
+    if(vSubArguments[dataCol]->GetFormulaToken()->GetType() ==
+     formula::svSingleVectorRef)
+     {
+         const formula::SingleVectorRefToken*pTmpDVR1= dynamic_cast<const
+         formula::SingleVectorRefToken *>(vSubArguments[dataCol]->
+         GetFormulaToken());
+         ss << "    tmp"<<dataCol<<"=";
+         ss << vSubArguments[dataCol]->GenSlidingWindowDeclRef()<<";\n";
+         ss << "    if(gid0>="<<pTmpDVR1->GetArrayLength()<<" ||isNan(";
+         ss << "tmp"<<dataCol<<"))\n";
+         ss << "        tmp"<<dataCol<<"=0;\n";
+     }
+    else
+        throw Unhandled();
+    int conditionCol = 0;
+    int conditionRow = 0;
+    if(vSubArguments[dataCol + 1]->GetFormulaToken()->GetType() ==
+    formula::svDoubleVectorRef)
+    {
+        formula::FormulaToken *tmpCur = vSubArguments[dataCol + 1]->
+            GetFormulaToken();
+        const formula::DoubleVectorRefToken*pCurDVR2= dynamic_cast<const
+            formula::DoubleVectorRefToken *>(tmpCur);
+        conditionCol = pCurDVR2->GetArrays().size();
+        conditionRow = pCurDVR2->GetArrayLength();
+    }
+    else{
+        throw Unhandled();
+    }
+    if(dataCol!=conditionCol)
+            throw Unhandled();
+    if(dataCol > 0 && dataRow > 0)
+    {
+        formula::FormulaToken *tmpCur1 = vSubArguments[0]->GetFormulaToken();
+        formula::FormulaToken *tmpCur2 = vSubArguments[dataCol + 1]->
+            GetFormulaToken();
+        const formula::DoubleVectorRefToken*pCurDVR1= dynamic_cast<const
+            formula::DoubleVectorRefToken *>(tmpCur1);
+        const formula::DoubleVectorRefToken*pCurDVR2= dynamic_cast<const
+            formula::DoubleVectorRefToken *>(tmpCur2);
+
+        if(pCurDVR1->IsStartFixed() && pCurDVR1->IsEndFixed() &&
+            pCurDVR2->IsStartFixed() && pCurDVR2->IsEndFixed())
+        {
+            ss << "    int i,j,p;\n";
+            ss << "    bool flag;\n";
+
+            ss << "    for(p = 1;p < " << dataRow << ";++p)\n";
+            ss << "    {\n";
+            ss << "        i = p;\n";
+            for(int i = 0; i < dataCol; ++i){
+                if(vSubArguments[i]->GetFormulaToken()->GetType() !=
+                        formula::svDoubleVectorRef)
+                    throw Unhandled();
+                formula::FormulaToken *tmpCur_tmp =
+                    vSubArguments[i]->GetFormulaToken();
+                const formula::DoubleVectorRefToken*pCurDVR_tmp=
+                    dynamic_cast<const formula::DoubleVectorRefToken *>
+                    (tmpCur_tmp);
+                ss << "        tmp"<<i<<"=";
+                ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n";
+                ss << "        if(isNan(tmp"<<i<<"))\n";
+                ss <<"            tmp"<<i<<" = 0;\n";
+            }
+            ss << "        flag = false;\n";
+            ss << "        for(j = 1; j < " << conditionRow << ";++j)\n";
+            ss << "        {\n";
+            ss << "            i = j;\n";
+            ss << "            if (flag)\n";
+            ss << "                break;\n";
+            ss << "            else{\n";
+            for(int i = dataCol + 1; i < dataCol + 1 + conditionCol; ++i){
+                if(vSubArguments[i]->GetFormulaToken()->GetType() !=
+                        formula::svDoubleVectorRef)
+                    throw Unhandled();
+                formula::FormulaToken *tmpCur_tmp =
+                    vSubArguments[i]->GetFormulaToken();
+                const formula::DoubleVectorRefToken*pCurDVR_tmp=
+                    dynamic_cast<const formula::DoubleVectorRefToken *>
+                    (tmpCur_tmp);
+                ss << "                tmp"<<i<<"=";
+                ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n";
+                ss << "                if(!isNan(tmp"<<i<<")){\n";
+                ss << "                    if(tmp"<<(i-dataCol-1)<<"!=tmp";
+                ss << i<<"){\n";
+                ss << "                        continue;\n";
+                ss << "                    }\n";
+                ss << "                }\n";
+            }
+            ss << "                flag=true;\n";
+            ss << "            }\n";
+            ss << "        }\n";
+            ss << "        if (flag){\n";
+            ss << "            count++;\n";
+            for(int i = 0; i < dataCol; ++i){
+                ss << "            if(tmp"<<dataCol<<"=="<<(i+1)<<"){\n";
+                ss << "                value=tmp"<<i<<";\n";
+                ss << "            }\n";
+            }
+            ss << "            mean+=value;\n";
+            ss << "        }\n";
+            ss << "    }\n";
+
+            ss << "    if(count<=1)\n";
+            ss << "        return 0;\n";
+
+            ss << "    mean/=count;\n";
+
+            ss << "    for(p = 1;p < " << dataRow << ";++p)\n";
+            ss << "    {\n";
+            ss << "        i = p;\n";
+            for(int i = 0; i < dataCol; ++i){
+                if(vSubArguments[i]->GetFormulaToken()->GetType() !=
+                        formula::svDoubleVectorRef)
+                    throw Unhandled();
+                formula::FormulaToken *tmpCur_tmp =
+                    vSubArguments[i]->GetFormulaToken();
+                const formula::DoubleVectorRefToken*pCurDVR_tmp=
+                    dynamic_cast<const formula::DoubleVectorRefToken *>
+                    (tmpCur_tmp);
+                ss << "        tmp"<<i<<"=";
+                ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n";
+                ss << "        if(isNan(tmp"<<i<<"))\n";
+                ss <<"            tmp"<<i<<" = 0;\n";
+            }
+            ss << "        flag = false;\n";
+            ss << "        for(j = 1; j < " << conditionRow << ";++j)\n";
+            ss << "        {\n";
+            ss << "            i = j;\n";
+            ss << "            if (flag)\n";
+            ss << "                break;\n";
+            ss << "            else{\n";
+            for(int i = dataCol + 1; i < dataCol + 1 + conditionCol; ++i){
+                if(vSubArguments[i]->GetFormulaToken()->GetType() !=
+                        formula::svDoubleVectorRef)
+                    throw Unhandled();
+                formula::FormulaToken *tmpCur_tmp =
+                    vSubArguments[i]->GetFormulaToken();
+                const formula::DoubleVectorRefToken*pCurDVR_tmp=
+                    dynamic_cast<const formula::DoubleVectorRefToken *>
+                    (tmpCur_tmp);
+                ss << "                tmp"<<i<<"=";
+                ss << vSubArguments[i]->GenSlidingWindowDeclRef()<<";\n";
+                ss << "                if(!isNan(tmp"<<i<<")){\n";
+                ss << "                    if(tmp"<<(i-dataCol-1)<<"!=tmp";
+                ss << i<<"){\n";
+                ss << "                        continue;\n";
+                ss << "                    }\n";
+                ss << "                }\n";
+            }
+            ss << "                flag=true;\n";
+            ss << "            }\n";
+            ss << "        }\n";
+            ss << "        if (flag){\n";
+            for(int i = 0; i < dataCol; ++i){
+                ss << "            if(tmp"<<dataCol<<"=="<<(i+1)<<"){\n";
+                ss << "                value=tmp"<<i<<";\n";
+                ss << "            }\n";
+            }
+            ss << "            var+=pow(mean-value,2);\n";
+            ss << "        }\n";
+            ss << "    }\n";
+
+            ss << "    var = sqrt( var/count );\n";
+        }
+        else
+            ss << "var = -1;\n";
+    }
+    else
+        ss << "var = -1;\n";
+    ss << "    return var;\n";
+    ss << "}";
+}
+
+
 
 
 }}
diff --git a/sc/source/core/opencl/op_database.hxx b/sc/source/core/opencl/op_database.hxx
index 01a1bf2ed805..3630657abade 100644
--- a/sc/source/core/opencl/op_database.hxx
+++ b/sc/source/core/opencl/op_database.hxx
@@ -55,6 +55,13 @@ public:
     virtual std::string BinFuncName(void) const { return "Dstdev"; }
 };
 
+class OpDstdevp: public CheckVariables
+{
+public:
+    virtual void GenSlidingWindowFunction(std::stringstream &ss,
+             const std::string sSymName, SubArguments &vSubArguments);
+    virtual std::string BinFuncName(void) const { return "Dstdevp"; }
+};
 
 
 }}