function [score, warped_points] = eval_wavelet_multi_point(knot_displacements, init_knots, unwarped_image_vec, unwarped_points_vec, image_set, points_set, n_modes, start_points_vec, warp_params, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, wavelet_evaluation_method, filter, filter_level)
% OPTIMISE_WAVELET_MULTI_POINT: 
%
%    GENERAL
%
%
% 
%    INPUT/S
%
%      -X:
%          X
%           
%    OUTPUT/S
%
%      -X:
%           X
%
%    PENDING WORK
%      
%      
%
%    KNOWN BUG/S
%
%      
%
%    COMMENT/S
%
%      
%
%    RELATED FUNCTION/S
%
%      
%
%    ABOUT
%
%      -Created:     February 9th, 2004
%      -Last update: February 2004
%      -Revision:    0.1.0
%      -Author:      R. S. Schestowitz, University of Manchester
% ==============================================================

% scale it all down to try and get it to move more!
% init_knots = init_knots*0.01;
% start_points_vec = start_points_vec*0.01;

warped_points = nrr_trans_1d(start_points_vec, init_knots, init_knots + knot_displacements, warp_params, []);
warped_image = interp1(unwarped_points_vec,unwarped_image_vec,warped_points, 'linear',0);
image_set(:,size(image_set,2)+1) = warped_image;
points_set(:,size(points_set,2)+1) = warped_points;

wavelet_evaluation_method = 'Default';

compressed_data = compute_wavelet(image_set, points_set, filter, filter_level);
score = measure_wavelet(compressed_data, wavelet_evaluation_method, warped_image);

% score = 0; % to be modified w.r.t. eval_groupwise_single_point