10 using namespace shogun;
11 using namespace Eigen;
19 double eps,
int itermax)
37 for (
int p = 0; p < m; p++)
39 for (
int q = p+1; q < m; q++)
45 if (fabs(theta) > eps)
78 for (k=0, cm=A; k<K; k++, cm+=M2)
80 diff_on = cm[pp] - cm[qq];
81 sum_off = cm[pq] + cm[qp];
83 G11 += diff_on * diff_on;
84 G22 += sum_off * sum_off;
85 G12 += diff_on * sum_off;
91 return -0.5 * atan2 ( toff , ton+sqrt(ton*ton+toff*toff) );
105 for (k=0, kMN=0; k<K; k++, kMN+=MN)
107 for (cpt=0, ix=p+kMN, iy=q+kMN; cpt<N; cpt++, ix+=M, iy+=M)
121 int k, ix, iy, cpt, kMN;
126 for (k=0, kMN=0; k<K; k++, kMN+=M*N)
128 for (cpt=0, ix=pM+kMN, iy=qM+kMN; cpt<M; cpt++)
132 A[ix++] = c*nx - s*ny;
133 A[iy++] = s*nx + c*ny;
147 for (
int j = 0; j < n; j++, ix+=m, iy+=m)