GPU Calc: implemented LOGNORMDIST

AMLOEXT-112 FIX

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

sc/source/core/opencl/formulagroupcl.cxx

@@ -1184,6 +1184,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments(
                 mvSubArguments.push_back(SoPHelper(ts,
                          ft->Children[i], new OpForecast));
                 break;
+            case ocLogNormDist:
+                mvSubArguments.push_back(SoPHelper(ts,
+                         ft->Children[i], new OpLogNormDist));
+                break;
             case ocExternal:
                 if ( !(pChild->GetExternal().compareTo(OUString(
                     "com.sun.star.sheet.addin.Analysis.getEffect"))))

sc/source/core/opencl/op_statistical.cxx

@@ -2395,6 +2395,82 @@ void OpForecast::GenSlidingWindowFunction(std::stringstream &ss,
     ss <<"    return tmp;\n";
     ss << "}";
 }
+void OpLogNormDist::GenSlidingWindowFunction(std::stringstream &ss,
+            const std::string sSymName, SubArguments &vSubArguments)
+{
+    FormulaToken *tmpCur0 = vSubArguments[0]->GetFormulaToken();
+    const formula::SingleVectorRefToken*tmpCurDVR0= dynamic_cast<const
+        formula::SingleVectorRefToken *>(tmpCur0);
+    FormulaToken *tmpCur1 = vSubArguments[1]->GetFormulaToken();
+    const formula::SingleVectorRefToken*tmpCurDVR1= dynamic_cast<const
+        formula::SingleVectorRefToken *>(tmpCur1);
+    FormulaToken *tmpCur2 = vSubArguments[2]->GetFormulaToken();
+    const formula::SingleVectorRefToken*tmpCurDVR2= dynamic_cast<const
+        formula::SingleVectorRefToken *>(tmpCur2);
+    FormulaToken *tmpCur3 = vSubArguments[3]->GetFormulaToken();
+    const formula::SingleVectorRefToken*tmpCurDVR3= dynamic_cast<const
+        formula::SingleVectorRefToken *>(tmpCur3);
+    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 << "    double arg0 = ";
+    ss << vSubArguments[0]->GenSlidingWindowDeclRef();
+    ss << ";\n";
+    ss << "    double arg1 = ";
+    ss << vSubArguments[1]->GenSlidingWindowDeclRef();
+    ss << ";\n";
+    ss << "    double arg2 = ";
+    ss << vSubArguments[2]->GenSlidingWindowDeclRef();
+    ss << ";\n";
+    ss << "    double arg3 = ";
+    ss << vSubArguments[3]->GenSlidingWindowDeclRef();
+    ss << ";\n";
+    ss << "    double tmp;\n";
+#ifdef ISNAN
+    ss << "    if(isNan(arg0)||(gid0>=";
+    ss << tmpCurDVR0->GetArrayLength();
+    ss << "))\n";
+    ss << "        arg0 = 0;\n";
+#endif
+#ifdef ISNAN
+    ss << "    if(isNan(arg1)||(gid0>=";
+    ss << tmpCurDVR1->GetArrayLength();
+    ss << "))\n";
+    ss << "        arg1 = 0;\n";
+#endif
+#ifdef ISNAN
+    ss << "    if(isNan(arg2)||(gid0>=";
+    ss << tmpCurDVR2->GetArrayLength();
+    ss << "))\n";
+    ss << "        arg2 = 0;\n";
+#endif
+#ifdef ISNAN
+    ss << "    if(isNan(arg3)||(gid0>=";
+    ss << tmpCurDVR3->GetArrayLength();
+    ss << "))\n";
+    ss << "        arg3 = 0;\n";
+#endif
+    ss << "    double temp = (log(arg0)-arg1)/arg2;\n";
+    ss << "    if(arg3)\n";
+    ss << "    {\n";
+    ss << "        if(arg0<=0)\n";
+    ss << "            tmp = 0.0;\n";
+    ss << "        else\n";
+    ss << "            tmp = 0.5 * erfc(-temp * 0.7071067811865475);\n";
+    ss << "    }\n";
+    ss << "    else\n";
+    ss << "        tmp = (0.39894228040143268 * exp(-(temp * temp)";
+    ss << " / 2.0))/arg2/arg0;\n";
+    ss << "    return tmp;\n";
+    ss << "}\n";
+}
 
 }}
 

sc/source/core/opencl/op_statistical.hxx

@@ -218,6 +218,13 @@ public:
             const std::string sSymName, SubArguments &vSubArguments);
     virtual std::string BinFuncName(void) const { return "Forecast"; }
 };
+class OpLogNormDist: public Normal
+{
+public:
+    virtual void GenSlidingWindowFunction(std::stringstream &ss,
+            const std::string sSymName, SubArguments &vSubArguments);
+    virtual std::string BinFuncName(void) const { return "LogNormdist"; }
+};
 }}
 
 #endif