24 using namespace shogun;
56 return "DiffusionMaps";
76 bool use_arpack =
true;
78 bool use_arpack =
false;
90 p_vector[i] += kernel_matrix.
matrix[j*N+i];
108 p_vector[i] += kernel_matrix.
matrix[j*N+i];
117 kernel_matrix.
matrix[i*N+j] /= p_vector[i]*p_vector[j];
130 arpack_dsxupd(kernel_matrix.
matrix,NULL,
false,N,m_target_dim,
"LA",
false,1,
false,
true,0.0,0.0,s_values,kernel_matrix.
matrix,info);
135 new_feature_matrix[j*m_target_dim+i] = kernel_matrix[j*m_target_dim+i];
140 SG_WARNING(
"LAPACK does not provide efficient routines to construct embedding (this may take time). Consider installing ARPACK.");
141 wrap_dgesvd(
'O',
'N',N,N,kernel_matrix.
matrix,N,s_values,NULL,1,NULL,1,&info);
145 new_feature_matrix[j*m_target_dim+i] =
146 kernel_matrix[(m_target_dim-i-1)*N+j];
150 SG_ERROR(
"Eigenproblem solving failed with %d code", info);