R_dissimilarity (t_float * const X_,
const int N_,
const int dim_,
t_float * const members_,
const unsigned char method,
const unsigned char metric,
const t_float p,
bool make_row_repr)
: Xa(X_),
dim(dim_),
members(members_),
postprocessfn(NULL),
postprocessarg(p),
N(N_)
{
switch (method) {
case METHOD_VECTOR_SINGLE:
switch (metric) {
case METRIC_R_EUCLIDEAN:
distfn = &R_dissimilarity::sqeuclidean<false>;
postprocessfn = &cluster_result::sqrt;
break;
case METRIC_R_MAXIMUM:
distfn = &R_dissimilarity::maximum;
break;
case METRIC_R_MANHATTAN:
distfn = &R_dissimilarity::manhattan;
break;
case METRIC_R_CANBERRA:
distfn = &R_dissimilarity::canberra;
break;
case METRIC_R_BINARY:
distfn = &R_dissimilarity::dist_binary;
break;
case METRIC_R_MINKOWSKI:
distfn = &R_dissimilarity::minkowski;
postprocessfn = &cluster_result::power;
break;
default:
throw std::runtime_error(std::string("Invalid method."));
}
break;
case METHOD_VECTOR_WARD:
postprocessfn = &cluster_result::sqrtdouble;
break;
default:
postprocessfn = &cluster_result::sqrt;
}
if (make_row_repr) {
row_repr.init(2*N-1);
for (t_index i=0; i<N; ++i) {
row_repr[i] = i;
}
}
}
void init(const t_index size) {
parent.init(2*size-1, 0);
nextparent = size;
}
