Skip to content

C
core
Kaydet (Commit) 31cc96c5 authored tarafından shiming zhang's avatar shiming zhang Kaydeden (comit) I-Jui (Ray) Sung
Dosyalara gözat
Seçenekler

GPU Calc: implemented VARP 

AMLOEXT-183 FIX

Change-Id: I65a26ecdac93e741a6fcaf40a9f135ac1758d875
Signed-off-by: 's avatarhaochen <haochen@multicorewareinc.com>
Signed-off-by: 's avatarI-Jui (Ray) Sung <ray@multicorewareinc.com>
üst a30a12b1
Hide whitespace changes
Inline Side-by-side
Showing with 244 additions and 0 deletions
sc/source/core/opencl/formulagroupcl.cxx
Dosyayı görüntüle @ 31cc96c5
...@@ -1345,6 +1345,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments( ...@@ -1345,6 +1345,10 @@ DynamicKernelSoPArguments::DynamicKernelSoPArguments(
mvSubArguments.push_back(SoPHelper(ts, mvSubArguments.push_back(SoPHelper(ts,
ft->Children[i],new OpBinomdist)); ft->Children[i],new OpBinomdist));
break; break;
case ocVarP:
mvSubArguments.push_back(SoPHelper(ts,
ft->Children[i], new OpVarP));
break;
case ocExternal: case ocExternal:
if ( !(pChild->GetExternal().compareTo(OUString( if ( !(pChild->GetExternal().compareTo(OUString(
"com.sun.star.sheet.addin.Analysis.getEffect")))) "com.sun.star.sheet.addin.Analysis.getEffect"))))
......
sc/source/core/opencl/op_statistical.cxx
Dosyayı görüntüle @ 31cc96c5
...@@ -254,6 +254,239 @@ void OpVar::GenSlidingWindowFunction(std::stringstream &ss, ...@@ -254,6 +254,239 @@ void OpVar::GenSlidingWindowFunction(std::stringstream &ss,
ss << " return vSum / (fCount - 1.0);\n"; ss << " return vSum / (fCount - 1.0);\n";
ss << "}\n"; ss << "}\n";
} }
void OpVarP::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 fSum=0.0;\n";
ss << " double fMean=0.0;\n";
ss << " double vSum=0.0;\n";
ss << " double fCount=0.0;\n";
ss << " double arg=0.0;\n";
unsigned i = vSubArguments.size();
while (i--)
{
FormulaToken *pCur = vSubArguments[i]->GetFormulaToken();
assert(pCur);
if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode())
{
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";
ss << " {\n";
#else
ss << "gid0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed())
{
#ifdef ISNAN
ss << "0; i < " << pDVR->GetArrayLength();
ss << " && i < gid0+"<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < gid0+"<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else if (!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
{
#ifdef ISNAN
ss << "0; i + gid0 < " << pDVR->GetArrayLength();
ss << " && i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else
{
#ifdef ISNAN
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
ss << " arg = ";
ss << vSubArguments[i]->GenSlidingWindowDeclRef();
ss << ";\n";
#ifdef ISNAN
ss << " if (isNan(arg))\n";
ss << " continue;\n";
#endif
ss << " fSum += arg;\n";
ss << " fCount += 1;\n";
ss << " }\n";
}
else if (pCur->GetType() == formula::svSingleVectorRef)
{
const formula::SingleVectorRefToken* pSVR =
dynamic_cast< const formula::SingleVectorRefToken*>(pCur);
assert(pSVR);
#ifdef ISNAN
ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n";
ss << " {\n";
ss << " if (!isNan(";
ss << vSubArguments[i]->GenSlidingWindowDeclRef() << "))\n";
ss << " {\n";
#endif
ss << " arg=";
ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
ss << " fSum += arg;\n";
ss << " fCount += 1;\n";
#ifdef ISNAN
ss << " }\n";
ss << " }\n";
#endif
}
else
{
ss << " arg=" << pCur->GetDouble() << ";\n";
ss << " fSum += arg;\n";
ss << " fCount += 1;\n";
}
}
else
{
ss << " arg = ";
ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
ss << " fSum += arg;\n";
ss << " fCount += 1;\n";
}
if (i == 0)
{
ss << " fMean=fSum/" << "fCount" << ";\n";
}
}
i = vSubArguments.size();
while (i--)
{
FormulaToken *pCur = vSubArguments[i]->GetFormulaToken();
assert(pCur);
if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode())
{
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";
ss << " {\n";
#else
ss << "gid0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else if (pDVR->IsStartFixed() && !pDVR->IsEndFixed())
{
#ifdef ISNAN
ss << "0; i < " << pDVR->GetArrayLength();
ss << " && i < gid0+"<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < gid0+"<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else if (!pDVR->IsStartFixed() && !pDVR->IsEndFixed())
{
#ifdef ISNAN
ss << "0; i + gid0 < " << pDVR->GetArrayLength();
ss << " && i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
else
{
#ifdef ISNAN
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#else
ss << "0; i < "<< nCurWindowSize << "; i++)\n";
ss << " {\n";
#endif
}
ss << " arg = ";
ss << vSubArguments[i]->GenSlidingWindowDeclRef();
ss << ";\n";
#ifdef ISNAN
ss << " if (isNan(arg))\n";
ss << " continue;\n";
#endif
ss << " vSum += (arg - fMean) * (arg - fMean);\n";
ss << " }\n";
}
else if (pCur->GetType() == formula::svSingleVectorRef)
{
const formula::SingleVectorRefToken* pSVR =
dynamic_cast< const formula::SingleVectorRefToken*>(pCur);
assert(pSVR);
#ifdef ISNAN
ss << " if (gid0 < " << pSVR->GetArrayLength() << ")\n";
ss << " {\n";
ss << " if (!isNan(";
ss << vSubArguments[i]->GenSlidingWindowDeclRef();
ss << "))\n";
ss << " {\n";
#endif
ss << " arg=";
ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
ss << " vSum += (arg - fMean) * (arg - fMean);\n";
#ifdef ISNAN
ss << " }\n";
ss << " }\n";
#endif
}
else
{
ss << " arg=" << pCur->GetDouble() << ";\n";
ss << " vSum += (arg - fMean) * (arg - fMean);\n";
}
}
else
{
ss << " arg = ";
ss << vSubArguments[i]->GenSlidingWindowDeclRef() << ";\n";
ss << " vSum += (arg - fMean) * (arg - fMean);\n";
}
}
ss << " if (fCount == 0.0)\n";
ss << " return DBL_MAX;\n";
ss << " else\n";
ss << " return vSum / fCount;\n";
ss << "}\n";
}
void OpExponDist::GenSlidingWindowFunction(std::stringstream &ss, void OpExponDist::GenSlidingWindowFunction(std::stringstream &ss,
const std::string sSymName, SubArguments &vSubArguments) const std::string sSymName, SubArguments &vSubArguments)
{ {
......
sc/source/core/opencl/op_statistical.hxx
Dosyayı görüntüle @ 31cc96c5
...@@ -35,6 +35,13 @@ public: ...@@ -35,6 +35,13 @@ public:
const std::string sSymName, SubArguments &vSubArguments); const std::string sSymName, SubArguments &vSubArguments);
virtual std::string BinFuncName(void) const { return "Var"; } virtual std::string BinFuncName(void) const { return "Var"; }
}; };
class OpVarP: public Normal
{
public:
virtual void GenSlidingWindowFunction(std::stringstream &ss,
const std::string sSymName, SubArguments &vSubArguments);
virtual std::string BinFuncName(void) const { return "VarP"; }
};
class OpSkewp: public Normal class OpSkewp: public Normal
{ {
public: public:
......
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