15 using namespace shogun;
71 float64_t result = compute_recursive1(avec, bvec, alen);
90 float64_t result = compute_recursive2(avec, bvec, alen);
163 for (int32_t j=0; j < len+1; j++)
166 for (int32_t k=1; k < d+1; k++)
173 for (int32_t j=k; j < len+1; j++)
174 DP[k*offs+j]=
DP[k*offs+j-1]+avec[j-1]*bvec[j-1]*
DP[(k-1)*offs+j-1];
189 for (int32_t i=0; i < len; i++)
192 for (int32_t k=1; k < d+1; k++)
195 for (int32_t i=0; i < len; i++)
203 for (int32_t k=1; k < d+1; k++)
206 for (int32_t s=1; s < k+1; s++)
212 sum += sign*
KD[k-s]*
KS[s];