final int N = szData - 1; // vector/matrix dimensions
final double c = Math.sqrt(2.0/N); // normalization constant
// Create matrix kernel and compute element values
DoubleSquareMatrix K = new DoubleSquareMatrix( N + 1 );
for (m=0; m<=N; m++)
for (n=m; n<=N; n++) {
k = c*Math.sin(Math.PI*m*n/N);
K.setElement(m, n, k);
K.setElement(n, m, k);
}
this.szData = N + 1;
this.matKer = K;
}