GPU Calc: Optimized PPMT

AMLOEXT-277

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

sc/source/core/opencl/op_financial.cxx

@@ -3402,129 +3402,90 @@ void OpPPMT::GenSlidingWindowFunction(std::stringstream &ss,
             ss << ", ";
         vSubArguments[i]->GenSlidingWindowDecl(ss);
     }
-    ss << ") {\n    ";
-    ss << "double tmp = 0;\n    ";
-    ss << "int gid0 = get_global_id(0);\n    ";
-    ss<<"double tmp0=0;\n    ";
-    ss<<"double tmp1=0;\n    ";
-    ss<<"double tmp2=0;\n    ";
-    ss<<"double tmp3=0;\n    ";
-    ss<<"double tmp4=0,tmp5=0;\n    ";
-    size_t nItems = 0;
+    ss<<") {\n";
+    ss<<"    double tmp = 0;\n";
+    ss<<"    int gid0 = get_global_id(0);\n";
+    ss<<"    double arg=0;\n";
+    ss<<"    double tmp0=0;\n";
+    ss<<"    double tmp1=0;\n";
+    ss<<"    double tmp2=0;\n";
+    ss<<"    double tmp3=0;\n";
+    ss<<"    double tmp4=0,tmp5=0;\n";
     ss <<"\n    ";
     //while (i-- > 1)
     for (size_t 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)
+        if (pCur->GetType() == formula::svSingleVectorRef)
         {
 #ifdef  ISNAN
             const formula::SingleVectorRefToken* pSVR =
                 dynamic_cast< const formula::SingleVectorRefToken* >(pCur);
-            ss << "if (gid0 < " << pSVR->GetArrayLength() << "){\n        ";
+            ss << "if (gid0 < " << pSVR->GetArrayLength() << "){\n";
 #else
-            nItems += 1;
 #endif
         }
         else if (pCur->GetType() == formula::svDouble)
         {
 #ifdef  ISNAN
-            ss << "{\n        ";
+            ss << "{\n";
 #endif
-            nItems += 1;
         }
         else
         {
 #ifdef  ISNAN
 #endif
-            nItems += 1;
         }
 #ifdef  ISNAN
         if(ocPush==vSubArguments[i]->GetFormulaToken()->GetOpCode())
         {
-            ss << "if (isNan(";
+            ss << "        arg=";
             ss << vSubArguments[i]->GenSlidingWindowDeclRef();
-            ss << "))\n            ";
-            ss << "tmp"<<i<<"= 0;\n        ";
-            ss << "else\n            ";
-            ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
-            ss << ";\n    }\n    ";
+            ss << ";\n";
+            ss << "        if (isNan(arg))\n";
+            ss << "            tmp"<<i<<"= 0;\n";
+            ss << "        else\n";
+            ss << "            tmp"<<i<<"=arg;\n";
+            ss << "    }\n";
         }
         else
         {
-            ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
-            ss <<";\n    ";
+            ss<<"    tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
+            ss<<";\n";
         }
 #else
-        ss << "tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
-        ss <<";\n    ";
+        ss<<"tmp"<<i<<"="<<vSubArguments[i]->GenSlidingWindowDeclRef();
+        ss<<";\n";
 #endif
     }
-    ss <<"double pmt ;\n    ";
-    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 << "}";
+    ss<<"    double pmt=0 ;\n";
+    ss<<"    if(tmp0==0.0)\n";
+    ss<<"        return -(tmp3+tmp4)/tmp2;\n";
+    ss<<"    pmt=pmt-tmp4-tmp3*pow(1.0+tmp0,tmp2);\n";
+    ss<<"    pmt=pmt*pow(( (1.0+tmp0*tmp5)* ";
+    ss<<"( (pow(1.0+tmp0,tmp2)-1.0)/tmp0)),-1);\n";
+    ss<<"    double temp = pow( 1+tmp0,tmp1-2);\n";
+    ss<<"    double re;\n";
+    ss<<"    if(tmp1==1.0){\n";
+    ss<<"        if(tmp5>0.0)\n";
+    ss<<"            re=0.0;\n";
+    ss<<"        else\n";
+    ss<<"            re=-tmp3;\n";
+    ss<<"    }\n";
+    ss<<"    else\n";
+    ss<<"    {\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<<"}";
 }
+
 void OpCoupdaybs::BinInlineFun(std::set<std::string>& decls,
     std::set<std::string>& funs)
 {