int PowerPolyByDouble(Polynom* _pPoly, Double* _pDouble, InternalType** _pOut)
{
bool bComplex1 = _pPoly->isComplex();
bool bScalar1 = _pPoly->isScalar();
double* bImg = _pDouble->getImg();
bool bNumericallyComplex1 = _pDouble->isNumericallyComplex();
if(!bNumericallyComplex1)
{
return 2;
}
if (_pDouble->isEmpty())
{
//p ** []
*_pOut = Double::Empty();
return 0;
}
if (bScalar1)
{
//p ^ x or p ^ X
int iRank = 0;
int* piRank = new int[_pDouble->getSize()];
_pPoly->getRank(&iRank);
for (int i = 0 ; i < _pDouble->getSize() ; i++)
{
int iInputRank = (int)_pDouble->get(i);
if (iInputRank < 0)
{
//call overload
_pOut = NULL;
delete[] piRank;
return 0;
}
piRank[i] = iRank * iInputRank;
}
Polynom* pOut = new Polynom(_pPoly->getVariableName(), _pDouble->getRows(), _pDouble->getCols(), piRank);
delete[] piRank;
pOut->setComplex(bComplex1);
for (int i = 0 ; i < _pDouble->getSize() ; i++)
{
SinglePoly* pCoeffOut = pOut->get(i);
int iCurrentRank = 0;
int iLoop = (int)_pDouble->get(i);
//initialize Out to 1
pCoeffOut->set(0, 1);
//get a copy of p
Polynom* pP = _pPoly->clone()->getAs<Polynom>();
pP->setComplex(_pPoly->isComplex());
while (iLoop)
{
SinglePoly* ps = pP->get()[0];
if (iLoop % 2)
{
int iRank = pP->getMaxRank();
if (bComplex1)
{
C2F(wpmul1)(pCoeffOut->get(), pCoeffOut->getImg(), &iCurrentRank, ps->get(), ps->getImg(), &iRank, pCoeffOut->get(), pCoeffOut->getImg());
}
else
{
C2F(dpmul1)(pCoeffOut->get(), &iCurrentRank, ps->get(), &iRank, pCoeffOut->get());
}
iCurrentRank += iRank;
}
iLoop /= 2;
if (iLoop)
{
//p = p * p
Polynom* pTemp = NULL;
MultiplyPolyByPoly(pP, pP, &pTemp);
pP->killMe();
pP = pTemp;
}
}
pP->killMe();
}
*_pOut = pOut;
}
return 0;
}
