function [ftest_all, numFeatures] = ...
    featExtrGaborAdapt(imagesCellTest, V, PCANet, bestWavelets, param, numImagesTest, stepPrint)

%one iteration to get size and init
ftest = Gabor_FeaExt(imagesCellTest(1), V, PCANet, bestWavelets, param);
numFeatures = size(ftest,1);
%init feature matrix
ftest_all = sparse(numFeatures, numImagesTest);
ftest_all(:, 1) = ftest;

%start pool
%parto da 2 perchŔ 1 ho giÓ fatto
parfor j = 2 : numel(imagesCellTest)
% for j = 2 : length(vectorIndexTest)
    
    %get id of current worker
%     t = getCurrentTask();
    
    %display progress
    if mod(j, stepPrint) == 0
%         fprintf(1, ['\t\tCore ' num2str(t.ID) ': ' num2str(j) ' / ' num2str(numImagesTest) '\n'])
    end %if mod(i, 100) == 0
    
    %imt = im2double(imread([dirDB filenameTest{j}]));
    %image size must be a power of 2
    %imt = imresize(imt, imageSize);
    
    %get image
    im = imagesCellTest(j);
    
    %PCANet output
    ftest = Gabor_FeaExt(im, V, PCANet, bestWavelets, param);
    
    %whos ftest
    
    %save descriptor
    %w/out wpca
    ftest_all(:, j) = ftest;
    
end %parfor i = 1 : numImagesTest