bug fix and readjustment of scaling factor

gui/src/de/mpg/mpimet/yacgui/utils/Default.java

@@ -51,7 +51,7 @@ public class Default {
 	public static final String INTERPOLATION_ATTR_WEIGHTED = Dom.Value.DISTANCE_WEIGHTED;
 	public static final String INTERPOLATION_ATTR_RBF_N = "9";
 	public static final String INTERPOLATION_ATTR_RBF_SCALE = "14.87973";
-	public static final String INTERPOLATION_ATTR_GAUSS_SCALE = "0.000002";
+	public static final String INTERPOLATION_ATTR_GAUSS_SCALE = "0.1";
 	public static final String INTERPOLATION_ATTR_RBF_KERNEL = Dom.Value.GAUSS_KERNEL;
 	
 	public static final String INTERPOLATION_ATTR_CONSERVATIVE_ORDER = "1";

src/interpolation_method_nnn.c

@@ -800,7 +800,7 @@ generate_interpolation_plan_nnn(struct intermediate_method_data * data) {
                double sum_weights = 0.0;
 
                for (unsigned k = 0; k < num_src_per_tgt; ++k) {
-                  weights[i*num_src_per_tgt+k] = exp ( -0.5 * distances[k] / gauss_scale * sum_d );
+                  weights[i*num_src_per_tgt+k] = exp ( -1.0 * distances[k] / ( gauss_scale * sum_d ) );
                   sum_weights += weights[i*num_src_per_tgt+k];
                }
                for (unsigned k = 0; k < num_src_per_tgt; ++k)

src/interpolation_method_nnn.h

@@ -54,7 +54,7 @@ enum yac_nnn_reduction_type {
 };
 
 #define YAC_DEFAULT_RBF_SCALE (1.487973e+01)
-#define YAC_DEFAULT_GAUSS_SCALE (2.0e-6)
+#define YAC_DEFAULT_GAUSS_SCALE (0.1)
 
 /**
  * constructor for a interpolation method of type interpolation_method_nnn