EMD distance - example of use

Rubner Of C Language implementation
We try to use Rubner Open C The language code (example1.c),
Compile dependencies emd.c and emd.h.
Type of characteristic quantity feature_t stay emd.h Is defined as follows ：
typedef struct { int X,Y,Z; } feature_t;
The specific implementation code is shown in emd.c.

The answer to the above example is as follows ：

# include <stdio.h>
# include <math.h>
# include "emd.h"

/*  Euclid distance  */
float dist(feature_t *F1, feature_t *F2) {
    
    int dX = F1->X - F2->X;
    int dY = F1->Y - F2->Y;
    int dZ = F1->Z - F2->Z;
    return sqrt(dXdX + dY*dY + dZ*dZ);
}

int main() {
    

    /*  Distribution P The eigenvector of the  */
    feature_t f1[4] = {
     {
    100,40,22}, {
    211,20,2}, {
    32,190,150}, {
    2,100,100} };
    /* Distribution Q The eigenvector of the  */
    feature_t f2[3] = {
     {
    0,0,0}, {
    50,100,80}, {
    255,255,255} };
    /* Distribution P The weight of  */
    float w1[5] = {
     0.4, 0.3, 0.2, 0.1 };
    /* Distribution Q The weight of  */
    float w2[3] = {
     0.5, 0.3, 0.2 };
    /* Distribution P The signature of the  */
    signature_t s1 = {
     4, f1, w1 };
    /* Distribution Q The signature of the  */
    signature_t s2 = {
     3, f2, w2};
    /*  Calculation EMD */
    float e;
    e = emd(&s1, &s2, dist, 0, 0);
    printf("emd = %f\n", e); return 0;

}