Skip to content

C
core
Kaydet (Commit) 6c006e1c authored tarafından mingli's avatar mingli Kaydeden (comit) I-Jui (Ray) Sung
Dosyalara gözat
Seçenekler

GPU Calc: implemented FINV 

AMLOEXT-159 FIX

Change-Id: Ib15803d87b2f108bc1a9f77eeab908b41e3eb402
Signed-off-by: 's avatarhaochen <haochen@multicorewareinc.com>
Signed-off-by: 's avatarI-Jui (Ray) Sung <ray@multicorewareinc.com>
üst a8bb0352
Hide whitespace changes
Inline Side-by-side
Showing with 215 additions and 0 deletions
sc/source/core/opencl/formulagroupcl.cxx
Dosyayı görüntüle @ 6c006e1c
......@@ -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"))))
......
sc/source/core/opencl/op_statistical.cxx
Dosyayı görüntüle @ 6c006e1c
......@@ -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;"
"}";
}
}}
......
sc/source/core/opencl/op_statistical.hxx
Dosyayı görüntüle @ 6c006e1c
......@@ -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"; }
};
}}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment