GPU Calc: implemented FINV

AMLOEXT-159 FIX

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

3 files changed, 215 insertions, 0 deletions

diff --git a/sc/source/core/opencl/formulagroupcl.cxx b/sc/source/core/opencl/formulagroupcl.cxx
index 3f7af7267e7e..00b9124db047 100644
--- a/sc/source/core/opencl/formulagroupcl.cxx
+++ b/sc/source/core/opencl/formulagroupcl.cxx
@@ -1241,6 +1241,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments(
                 mvSubArguments.push_back(SoPHelper(ts,
                          ft->Children[i], new OpFloor));
                 break;
+            case ocFInv:
+                mvSubArguments.push_back(SoPHelper(ts,
+                         ft->Children[i], new OpFInv));
+                break;
             case ocExternal:
                 if ( !(pChild->GetExternal().compareTo(OUString(
                     "com.sun.star.sheet.addin.Analysis.getEffect"))))
diff --git a/sc/source/core/opencl/op_statistical.cxx b/sc/source/core/opencl/op_statistical.cxx
index c3ffb5e0352f..3d8aa2e600bf 100644
--- a/sc/source/core/opencl/op_statistical.cxx
+++ b/sc/source/core/opencl/op_statistical.cxx
@@ -3264,6 +3264,208 @@ void OpGammaInv::GenSlidingWindowFunction(std::stringstream &ss,
     "    }\n"
     "}\n";
 }
+void OpFInv::BinInlineFun(std::set<std::string>& decls,
+    std::set<std::string>& funs)
+{
+    decls.insert(fMachEpsDecl);decls.insert(fMaxGammaArgumentDecl);
+    decls.insert(lcl_getLanczosSumDecl);decls.insert(GetBetaDecl);
+    decls.insert(GetLogBetaDecl);decls.insert(GetBetaDistPDFDecl);
+    decls.insert(lcl_GetBetaHelperContFracDecl);decls.insert(GetFInvValueDecl);
+    funs.insert(lcl_getLanczosSum);funs.insert(GetBeta);
+    funs.insert(GetLogBeta);funs.insert(GetBetaDistPDF);
+    funs.insert(lcl_GetBetaHelperContFrac);funs.insert(GetFInvValue);
+}
+
+void OpFInv::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";
+    ss << "    int gid0=get_global_id(0);\n";
+    ss << "    double tmp;\n";
+    ss << "    double arg0,arg1,arg2;\n";
+    size_t i = vSubArguments.size();
+    size_t nItems = 0;
+    for (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)
+        {
+#ifdef  ISNAN
+            const formula::SingleVectorRefToken* pSVR =
+                dynamic_cast< const formula::SingleVectorRefToken* >(pCur);
+            ss << "    if (gid0 < " << pSVR->GetArrayLength() << ")\n";
+            ss << "    {\n";
+            ss << "        if (isNan(";
+            ss << vSubArguments[i]->GenSlidingWindowDeclRef();
+            ss << "))\n";
+            ss << "            arg"<<i<<"= 0;\n";
+            ss << "        else\n";
+            ss << "            arg"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
+            ss << ";\n";
+            ss << "    }\n";
+            ss << "    else\n";
+            ss << "        arg"<<i<<"= 0;\n";
+#endif
+        }
+        else if (pCur->GetType() == formula::svDouble)
+        {
+#ifdef  ISNAN
+            ss << "    if (isNan(";
+            ss << vSubArguments[i]->GenSlidingWindowDeclRef();
+            ss << "))\n";
+            ss << "        arg"<<i<<"= 0;\n";
+            ss << "    else\n";
+            ss << "        arg"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
+            ss << ";\n";
+#endif
+        }
+    }
+    ss << "    double fF2=floor(arg2);\n"
+    "    double fF1=floor(arg1);\n"
+    "    bool bConvError;\n"
+    "    double fAx=fF1*0.5;\n"
+    "    double fBx=fF1;\n"
+    "    bConvError = false;\n"
+    "    const double fYEps = 1.0E-307;\n"
+    "    const double fXEps = 2.22045e-016;\n"
+    "    double fAy = arg0-GetFInvValue(fF1,fF2,fAx);\n"
+    "    double fBy = arg0-GetFInvValue(fF1,fF2,fBx);\n"
+    "    double fTemp;\n"
+    "    unsigned short nCount;\n"
+    "    for (nCount = 0; nCount < 1000 && !((fAy < 0.0 && fBy > 0.0)"
+    " || (fAy > 0.0 && fBy < 0.0)); nCount++)\n"
+    "    {\n"
+    "        if (fabs(fAy) <= fabs(fBy))\n"
+    "        {\n"
+    "            fTemp = fAx;\n"
+    "            fAx += 2.0 * (fAx - fBx);\n"
+    "            if (fAx < 0.0)\n"
+    "                fAx = 0.0;\n"
+    "            fBx = fTemp;\n"
+    "            fBy = fAy;\n"
+    "            fAy = arg0-GetFInvValue(fF1,fF2,fAx);\n"
+    "        }\n"
+    "        else\n"
+    "        {\n"
+    "            fTemp = fBx;\n"
+    "            fBx += 2.0 * (fBx - fAx);\n"
+    "            fAx = fTemp;\n"
+    "            fAy = fBy;\n"
+    "            fBy = arg0-GetFInvValue(fF1,fF2,fBx);\n"
+    "        }\n"
+    "    }\n"
+    "    if (fAy == 0.0)\n"
+    "    {\n"
+    "        tmp = fAx;\n"
+    "        return tmp;\n"
+    "    }\n"
+    "    if (fBy == 0.0)\n"
+    "    {\n"
+    "        tmp = fBx;\n"
+    "        return tmp;\n"
+    "    }\n"
+    "    if (!((fAy < 0.0 && fBy > 0.0) || (fAy > 0.0 && fBy < 0.0)))\n"
+    "    {\n"
+    "        bConvError = true;\n"
+    "        tmp = 0.0;\n"
+    "        return tmp;\n"
+    "    }\n"
+    "    double fPx = fAx;\n"
+    "    double fPy = fAy;\n"
+    "    double fQx = fBx;\n"
+    "    double fQy = fBy;\n"
+    "    double fRx = fAx;\n"
+    "    double fRy = fAy;\n"
+    "    double fSx = 0.5 * (fAx + fBx);\n"
+    "    bool bHasToInterpolate = true;\n"
+    "    nCount = 0;\n"
+    "    while ( nCount < 500 && fabs(fRy) > fYEps &&"
+    "(fBx-fAx) > fmax( fabs(fAx), fabs(fBx)) * fXEps )\n"
+    "    {\n"
+    "        if (bHasToInterpolate)\n"
+    "        {\n"
+    "            if (fPy!=fQy && fQy!=fRy && fRy!=fPy)\n"
+    "            {\n"
+    "                fSx = fPx * fRy * fQy / (fRy-fPy) / (fQy-fPy)+"
+    " fRx * fQy * fPy"
+    "/ (fQy-fRy) / (fPy-fRy)+ fQx * fPy * fRy / (fPy-fQy) / (fRy-fQy);\n"
+    "                bHasToInterpolate = (fAx < fSx) && (fSx < fBx);\n"
+    "            }\n"
+    "            else\n"
+    "                bHasToInterpolate = false;\n"
+    "        }\n"
+    "        if(!bHasToInterpolate)\n"
+    "        {\n"
+    "            fSx = 0.5 * (fAx + fBx);\n"
+    "            fPx = fAx; fPy = fAy;\n"
+    "            fQx = fBx; fQy = fBy;\n"
+    "            bHasToInterpolate = true;\n"
+    "        }\n"
+    "        fPx = fQx; fQx = fRx; fRx = fSx;\n"
+    "        fPy = fQy; fQy = fRy;\n"
+    "        fRy = arg0-GetFInvValue(fF1,fF2,fSx);\n"
+    "        if ((fAy < 0.0 && fRy > 0.0) || (fAy > 0.0 && fRy < 0.0))\n"
+    "        {\n"
+    "            fBx = fRx; fBy = fRy;\n"
+    "        }\n"
+    "        else\n"
+    "        {\n"
+    "            fAx = fRx; fAy = fRy;\n"
+    "        }\n"
+    "        bHasToInterpolate = bHasToInterpolate && (fabs(fRy)"
+    " * 2.0 <= fabs(fQy));\n"
+    "        ++nCount;\n"
+    "    }\n"
+    "    tmp = fRx;\n"
+    "    return tmp;"
+    "}";
+}
 
 
 }}
diff --git a/sc/source/core/opencl/op_statistical.hxx b/sc/source/core/opencl/op_statistical.hxx
index 153d1e926e25..ca8e343fbfbb 100644
--- a/sc/source/core/opencl/op_statistical.hxx
+++ b/sc/source/core/opencl/op_statistical.hxx
@@ -279,6 +279,15 @@ public:
 );
     virtual std::string BinFuncName(void) const { return "GammaInv"; }
 };
+class OpFInv: public Normal
+{
+public:
+    virtual void GenSlidingWindowFunction(std::stringstream &ss,
+            const std::string sSymName, SubArguments &vSubArguments);
+    void BinInlineFun(std::set<std::string>& decls,std::set<std::string>& funs
+);
+    virtual std::string BinFuncName(void) const { return "FInv"; }
+};
 
 }}