Skip to content

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

GPU Calc: implement NAN argument handling in PPMT 

AMLOEXT-119 FIX

Change-Id: I43dcd08ff846a4008f13c9d02a69478898363131
Signed-off-by: 's avatarhaochen <haochen@multicorewareinc.com>
Signed-off-by: 's avatarI-Jui (Ray) Sung <ray@multicorewareinc.com>
üst 7acf9f27
Show whitespace changes
Inline Side-by-side
Showing with 123 additions and 70 deletions
sc/source/core/opencl/op_financial.cxx
Dosyayı görüntüle @ acd9cdf4
...@@ -2069,89 +2069,142 @@ void OpPMT::GenSlidingWindowFunction(std::stringstream &ss, ...@@ -2069,89 +2069,142 @@ void OpPMT::GenSlidingWindowFunction(std::stringstream &ss,
funs.insert(GetZw); funs.insert(GetZw);
} }
void OpPPMT::GenSlidingWindowFunction(std::stringstream &ss, void OpPPMT::GenSlidingWindowFunction(std::stringstream &ss,
const std::string sSymName, SubArguments &vSubArguments) const std::string sSymName, SubArguments &vSubArguments)
{ {
ArgVector argVector; ArgVector argVector;
ss << "\ndouble " << sSymName; ss << "\ndouble " << sSymName;
ss << "_"<< BinFuncName() <<"("; ss << "_"<< BinFuncName() <<"(";
for (unsigned i = 0; i < vSubArguments.size(); i++) for (unsigned i = 0; i < vSubArguments.size(); i++)
{ {
if (i) if (i)
ss << ","; ss << ", ";
vSubArguments[i]->GenSlidingWindowDecl(ss); vSubArguments[i]->GenSlidingWindowDecl(ss);
argVector.push_back(vSubArguments[i]->GenSlidingWindowDeclRef());
} }
ss << ") {\n\t"; ss << ") {\n ";
ss << "int gid0 = get_global_id(0);\n\t"; ss << "double tmp = 0;\n ";
ss<<"double tFv=0,tType=0;\n\t"; ss << "int gid0 = get_global_id(0);\n ";
if(vSubArguments.size()==5) ss<<"double tmp0=0;\n ";
ss<<"tFv="<<vSubArguments[4]->GenSlidingWindowDeclRef()<<";\n\t"; ss<<"double tmp1=0;\n ";
else if(vSubArguments.size()==6) ss<<"double tmp2=0;\n ";
ss<<"double tmp3=0;\n ";
ss<<"double tmp4=0,tmp5=0;\n ";
unsigned i = vSubArguments.size();
size_t nItems = 0;
ss <<"\n ";
//while (i-- > 1)
for (size_t i = 0; i < vSubArguments.size(); i++)
{ {
ss<<"tType="<<vSubArguments[5]->GenSlidingWindowDeclRef(); FormulaToken *pCur = vSubArguments[i]->GetFormulaToken();
ss<<";\n\t"; assert(pCur);
ss<<"tFv="<<vSubArguments[4]->GenSlidingWindowDeclRef(); if (pCur->GetType() == formula::svDoubleVectorRef)
ss<<";\n\t"; {
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
} }
ss << "double tmp = 0.0;\n\t"; else {
ss <<"double pmt ;\n\t"; #ifdef ISNAN
ss<<"if("<<vSubArguments[0]->GenSlidingWindowDeclRef()<<"==0.0)\n\t"; ss << "0; i < "<< nCurWindowSize << "; i++){\n ";
ss<<"\treturn ("<<vSubArguments[3]->GenSlidingWindowDeclRef(); #else
ss<<"+tFv)/"; ss << "0; i < "<< nCurWindowSize << "; i++)\n ";
ss<<vSubArguments[2]->GenSlidingWindowDeclRef()<<";\n\t"; #endif
ss <<"double temp1=0;\n\t"; }
ss <<"double abl = pow(1.0+"; nItems += nCurWindowSize;
ss <<vSubArguments[0]->GenSlidingWindowDeclRef(); }
ss<<","; else if (pCur->GetType() == formula::svSingleVectorRef)
ss <<vSubArguments[2]->GenSlidingWindowDeclRef(); {
ss<<");\n\t"; #ifdef ISNAN
ss <<"temp1-=tFv"; const formula::SingleVectorRefToken* pSVR =
ss<<";\n\t"; dynamic_cast< const formula::SingleVectorRefToken* >(pCur);
ss<<"temp1-="; ss << "if (gid0 < " << pSVR->GetArrayLength() << "){\n ";
ss<<vSubArguments[3]->GenSlidingWindowDeclRef(); #else
ss<<"*abl;\n\t"; nItems += 1;
ss <<"pmt = temp1/(1.0+"; #endif
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); }
ss<<"*tType"; else if (pCur->GetType() == formula::svDouble)
ss<<") / ( (abl-1.0)/"; {
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); #ifdef ISNAN
ss<<");\n\t"; ss << "{\n ";
ss <<"double temp = pow( 1+"; #endif
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); nItems += 1;
ss<<" ,"; }
ss<<vSubArguments[1]->GenSlidingWindowDeclRef(); else
ss<<"-2);\n\t"; {
ss <<"double re;\n\t"; #ifdef ISNAN
ss <<"if (tType"; #endif
ss<<" == 0)\n\t"; nItems += 1;
ss <<"{\n\t"; }
ss <<" re = -"; #ifdef ISNAN
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode())
ss<<" * "; {
ss<<vSubArguments[3]->GenSlidingWindowDeclRef(); ss << "if (isNan(";
ss<<"*pow(1+"; ss << vSubArguments[i]->GenSlidingWindowDeclRef();
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); ss << "))\n ";
ss<<" ,"; ss << "tmp"<<i<<"= 0;\n ";
vSubArguments[1]->GenSlidingWindowDeclRef(); ss << "else\n ";
ss<<"-1) - pmt * (pow(1+"; ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); ss << ";\n }\n ";
ss<<" ,"; }
ss<<vSubArguments[1]->GenSlidingWindowDeclRef(); else
ss<<"-1)-1);\n\t"; {
ss <<"}\n\t"; ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
ss <<"else\n\t"; ss <<";\n ";
ss <<"{\n\t"; }
ss <<" re = -"; #else
ss<<vSubArguments[0]->GenSlidingWindowDeclRef(); ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
ss<<" * ("; ss <<";\n ";
ss<<vSubArguments[3]->GenSlidingWindowDeclRef(); #endif
ss<<" + pmt)* temp - pmt * (temp-1);\n\t"; }
ss <<"}\n\t"; ss <<"double pmt ;\n ";
ss <<"tmp = pmt - re;\n\t"; ss<<"if(tmp0==0.0)\n ";
ss<<" return -(tmp3+tmp4)/";
ss<<"tmp2;\n double temp1=0;\n ";
ss <<"double abl = pow(1.0+tmp0,tmp2);\n ";
ss <<"temp1-=tmp4";
ss<<";\n ";
ss<<"temp1-=tmp3*abl;\n ";
ss <<"pmt = temp1/(1.0+tmp0*tmp5";
ss<<") / ( (abl-1.0)/tmp0);\n ";
ss <<"double temp = pow( 1+tmp0,tmp1-2);\n ";
ss <<"double re;\n ";
ss <<"if(tmp1==1.0)\n {\n ";
ss <<"if(tmp5>0.0)\n ";
ss <<" re=0.0;\n ";
ss << "else\n re=-tmp3;\n }\n ";
ss <<"else\n {\n ";
ss <<"if(tmp5>0.0)\n ";
ss <<" re=GetZw(tmp0, tmp1-2.0, pmt, tmp3, 1.0) - pmt;\n ";
ss <<"else\n ";
ss <<" re=GetZw(tmp0, tmp1-1.0, pmt, tmp3, 0.0);\n ";
ss <<"}\n ";
ss <<"re = re * tmp0;\n ";
ss <<"tmp = pmt - re;\n ";
ss << "return tmp;\n"; ss << "return tmp;\n";
ss << "}"; ss << "}";
} }
void OpCoupdaybs::BinInlineFun(std::set<std::string>& decls, void OpCoupdaybs::BinInlineFun(std::set<std::string>& decls,
std::set<std::string>& funs) std::set<std::string>& funs)
{ {
......
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