22 using namespace shogun;
53 return "DiffusionMaps";
82 p_vector[i] += kernel_matrix.
matrix[j*N+i];
87 cblas_dger(CblasColMajor,N,N,1.0,p_vector,1,p_vector,1,p_matrix,N);
99 p_vector[i] += kernel_matrix.
matrix[j*N+i];
103 float64_t ppt = cblas_ddot(N,p_vector,1,p_vector,1);
106 for (i=0; i<N*N; i++)
108 kernel_matrix.
matrix[i] /= ppt;
116 cblas_dgemm(CblasColMajor,CblasTrans,CblasNoTrans,
118 1.0,kernel_matrix.
matrix,N,
126 SG_ERROR(
"DGESVD failed with %d code", info);
140 new_feature_matrix[j*m_target_dim+i] = kernel_matrix.
matrix[(m_target_dim-i-1)*N+j]/kernel_matrix.
matrix[(
m_target_dim)*N+j];