<p>In <a href="https://github.com/mlpack/mlpack/pull/742#discussion_r72753244">src/mlpack/methods/preprocess/preprocess_describe_main.cpp</a>:</p>
<pre style='color:#555'>> +* Calculates the sum of deviations to the Nth Power
> +*
> +* @param input Vector that captures a dimension of a dataset
> +* @param rowMean Mean of the given vector.
> +* @return sum of nth power deviations
> +*/
> +double SumNthPowerDeviations(const arma::rowvec& input,
> + const double& rowMean,
> + const size_t Nth) // Degree of Power
> +{
> + double sum = 0;
> + for (size_t i = 0; i < input.n_elem; ++i)
> + {
> + sum += pow(input(i) - rowMean, Nth);
> + }
> + return sum;
</pre>
<p>This one could be vectorize, which performance is better need to measure, but it should make the codes shorter and easier to read.</p>
<pre><code>return arma::sum(arma::pow(input - rowMean,
static_cast<double>(Nth)));
</code></pre>
<p style="font-size:small;-webkit-text-size-adjust:none;color:#666;">—<br />You are receiving this because you are subscribed to this thread.<br />Reply to this email directly, <a href="https://github.com/mlpack/mlpack/pull/742/files/27ac82e876c3d68637d9e72d28bf0dadbc229ace#r72753244">view it on GitHub</a>, or <a href="https://github.com/notifications/unsubscribe-auth/AJ4bFJZqy_irKi1050lEFsm02VImjVz6ks5qaa-vgaJpZM4JWax8">mute the thread</a>.<img alt="" height="1" src="https://github.com/notifications/beacon/AJ4bFBd03S0_C8T8BaDiqKLdWXkM2_58ks5qaa-vgaJpZM4JWax8.gif" width="1" /></p>
<div itemscope itemtype="http://schema.org/EmailMessage">
<div itemprop="action" itemscope itemtype="http://schema.org/ViewAction">
<link itemprop="url" href="https://github.com/mlpack/mlpack/pull/742/files/27ac82e876c3d68637d9e72d28bf0dadbc229ace#r72753244"></link>
<meta itemprop="name" content="View Pull Request"></meta>
</div>
<meta itemprop="description" content="View this Pull Request on GitHub"></meta>
</div>