toolbox cleanup and migration from svn

based on commit 13321 from svn

.gitmodules

+[submodule "internal"]
+	path = internal
+	url = https://git.rwth-aachen.de/ita/toolbox-intern.git
+	branch = master

Debug/test_ita_Super.m

+function test_ita_Super()
+%Test itaSuper
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+a = itaSuper([100 1],1); %Create 100x100 itaSupers with 1 data-elem
+b = itaSuper([100 1],1);
+
+tic
+for idx = 1:numel(a)
+    a(idx).data = 1;
+end
+t1 = toc;
+
+tic
+for idx = 1:numel(b)
+    b(idx).data = 1;
+end
+t2 = toc;
+
+disp(t2/t1);

Debug/test_ita_Value.m

+function test_ita_Value()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+a = itaValue('-0.20 kg');
+b = itaValue(2);
+c = itaValue(21,'N');
+d = itaValue();
+s = ita_generate('noise',1,44100,14);
+
+%% stupid arithmetic
+disp('aritmetic test')
+x = a + b;
+x = 2 + b;
+x = b + 2;
+
+x = a - b;
+x = 2 - b;
+x = b - 2;
+
+x = a * b;
+x = 2 * b;
+x = b * 2;
+
+x = a .* b;
+x = 2 .* b;
+x = b .* 2;
+
+x = a / b;
+x = 2 / b;
+x = b / 2;
+
+x = a ./ b;
+x = 2 ./ b;
+x = b ./ 2;
+
+x = sqrt(b);
+x = power(b,2);
+x = power(2,b);
+
+x = log(c);
+x = log10(c);
+x = exp(c);
+
+%% test with itaAudio object
+disp('test with audio objects')
+res = a * s;
+res = a .* s;
+res = a / s;
+res = a ./ s;

Debug/test_ita_all.m

+function test_ita_all(varargin)
+%ITA_TEST_ALL - Test for the ita toolbox
+%  runs all funktions called test_ita_* to find any error within the toolbox
+%  you dont need to change anything in here, all test_ita_ routines will be found automatically
+%
+%  Call: test_ita_all(Options)
+%
+%       Options:
+%           'kernel' [false] - Check only kernel functions, no applications
+%
+%   See also  ita_sqrt, test, ita_sum, ita_audio2struct, test, ita_channelnames_to_numbers.
+%
+%   Reference page in Help browser
+%        <a href="matlab:doc test_ita_all">doc test_ita_all</a>
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+% Author: Roman Scharrer -- Email: rsc@akustik.rwth-aachen.de
+% Created:  15-Jan-2009
+
+%% get the path information
+
+sArgs = struct('kernel',false);
+sArgs = ita_parse_arguments(sArgs,varargin);
+
+
+filename_workspace = 'workspace_tmp.mat';
+test_path = which('test_ita_all.m'); %Find the test directory
+ind = strfind(test_path,'test_ita_all.m');
+test_path(ind-1:end) = []; %Remove filename from path
+
+% save all variables in a temp file
+save([test_path filesep filename_workspace]); 
+% clear everythin including the class definition
+% ccx; 
+clear classes;
+
+% as we cleared all variables we have to look again for the path
+filename_workspace = 'workspace_tmp.mat';
+test_path = which('test_ita_all.m'); %Find the test directory
+ind = strfind(test_path,'test_ita_all.m');
+test_path(ind-1:end) = []; %Remove filename from path
+load([test_path filesep filename_workspace]); % restore arguments
+
+
+tic;
+disp('********************************************************************')
+disp('********************* TESTING ITA TOOLBOX **************************')
+disp('********************************************************************')
+
+%old_path = cd; %Save old working directory
+%cd(test_path); %Change directory to the test directory
+
+%test_files = dir('*test_*.m'); %Find all files with *test_*
+if sArgs.kernel
+    test_files = rdir([ita_toolbox_path filesep 'Debug' filesep '**' filesep 'test_ita_*.m']); % Recursive find of all files with 'test_ita_*.m'
+else
+    test_files = rdir([ita_toolbox_path filesep '**' filesep 'test_ita_*.m']); % Recursive find of all files with 'test_ita_*.m'
+end
+
+errors = {};
+tic;
+for test_number = 1:length(test_files)
+    functionname = strtok(fliplr(strtok(fliplr(test_files(test_number).name),filesep)),'.'); %Find the function name in the filename (everything bevor the first '.')
+
+    if ~strcmp(functionname,'test_ita_all') %Exclude test_ita_all
+        disp(['Testing: ' functionname '...'])
+        try
+            feval(functionname); %Run the test function
+        catch ME
+            ME_test = MException('test_ita_all:error', [functionname ' returned the above error(s)']);
+            ME = addCause(ME,ME_test);
+            errors(end+1) = {ME}; %Get all errors first, we can deal with them later
+        end
+    end
+
+end
+
+toc;
+% cd(old_path); %Restore old working directory
+% restore old workspace variables
+load([test_path filesep filename_workspace]);
+delete([test_path filesep filename_workspace]);
+
+disp('********************************************************************')
+if ~isempty(errors)
+    fprintf(2,'The following errors occured:\n');
+    for i = 1:length(errors)
+        ita_verbose_info('********************************************************************',0);
+        ita_verbose_info(getReport(errors{i}),0);
+        ita_verbose_info('********************************************************************',0);
+    end
+else
+    disp('Congratulations! No Errors occurred!');
+    disp(['Test took ' num2str(toc,3) ' sec']);
+end
+
+commandwindow()
+
+
+
+end
\ No newline at end of file

Debug/test_ita_analysis.m

+function test_ita_analysis()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+%% get two signals and a tp
+sr = 44100;
+fftdeg = 15;
+a       = ita_generate('sine',1,1000,sr,fftdeg);
+b       = ita_generate('noise',1,sr,fftdeg);
+c       = ita_merge(ita_fft(a),b);
+tp      = ita_fft(ita_merge(a,a,a,a));
+tp.signalType = 'energy';
+
+
+%% used to test the transformation at the beginning of functions
+timeData = b;
+freqData = ita_fft(b);
+
+
+%% fft transformation test - energy-signal
+timeData = ita_generate('impulse', 1, sr, fftdeg);
+freqData = ita_fft(timeData);
+timeData_b = ita_ifft(freqData);
+
+if sum(find(timeData.timeData-timeData_b.timeData))
+   error('FFT or IFFT ERROR with energy-signal'); 
+end
+
+%% fft transformation test - power-signal
+timeData = ita_generate('sine', 1, 1000, sr,fftdeg);
+freqData = ita_fft(timeData);
+timeData_b = ita_ifft(freqData);
+if 0.0001 < sum(abs(timeData.timeData.^2-timeData_b.timeData.^2))
+   error('FFT or IFFT ERROR with power-signal');
+end
+
+%% test functions
+c = ita_spk2frequencybands(b);
+
+
+end
\ No newline at end of file

Debug/test_ita_arithmetic.m

+function test_ita_arithmetic()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+%% get two signals and a tp
+sr = 44100;
+fftdeg = 15;
+a       = ita_generate('sine',1,1000,sr,fftdeg);
+b       = ita_generate('noise',1,sr,fftdeg);
+c       = merge(a,b);
+tp      = fft(merge(a,a,a));
+tp.signalType = 'energy';
+imp = ita_generate('impulse',1,sr,fftdeg-5);
+
+%% used to test the transformation at the beginning of functions
+a_spk    = fft(a);
+timeData = b;
+freqData = fft(b);
+
+%% test functions
+d = ita_add(timeData,timeData);
+d = ita_add(freqData,freqData);
+d = ita_add(freqData,timeData);
+d = ita_add(timeData,freqData);
+
+d = ita_conj(freqData);
+
+d = ita_divide_spk(a_spk,tp);
+
+d = ita_invert_spk(timeData);
+d = ita_invert_spk(freqData);
+
+d = ita_multiply_dat(a,timeData);
+
+d = ita_multiply_spk(a_spk,tp);
+if ~(tp == ita_divide_spk(ita_multiply_spk(tp,freqData),freqData)) 
+    error('ITA_DIVIDE_SPK/ITA_MULTIPLY_SPK:Oh Lord. Function seems not to work properly.')
+end
+
+d = ita_negate(a);
+d = ita_negate(a_spk);
+if ~(b == ita_negate(ita_negate(b))) 
+    error('ITA_NEGATE:Oh Lord. Function seems not to work properly.')
+end
+
+d = ita_negTozero(freqData);
+
+d = ita_posTozero(freqData);
+
+d = timeData.^2;
+d = freqData^2;
+
+d = timeData.^(0.5);
+d = freqData^(0.5);
+
+d = ita_subtract(timeData,timeData);
+d = ita_subtract(freqData,freqData);
+if ~(timeData == ita_subtract(ita_add(timeData,freqData),freqData)) 
+    error('ITA_SUBTRACT/ITA_ADD:Oh Lord. Function seems not to work properly.')
+end
+
+
+d = ita_sum(freqData);
+d = ita_sum(timeData);
+
+
+d = ita_amplify(timeData,2);
+d = ita_amplify(timeData,'2');
+d = ita_amplify(timeData,'2dB');
+
+
+%overloaded stuff
+d = a .* tp;
+d = a * 2;
+d = a .* 2;
+d = 2 * a;
+d = 2 .* a;
+
+d = a ./ b;
+d = a / 2;
+d = a ./ 2;
+d = 2 / a;
+d = 2 ./ a;
+
+d = a + a;
+d = a - a;
+d = 1 + a;
+d = 1 - a;
+d = a + 1;
+d = a - 1;
+
+d = a^2;
+d = sqrt(a_spk);
+
+d = abs(freqData);
+d = abs(timeData);
+
+d = sum(freqData);
+d = sum(timeData);
+
+d = mean(freqData);
+d = mean(timeData);
+
+end
\ No newline at end of file

Debug/test_ita_backwardcompatibility.m

+function test_ita_backwardcompatibility
+% %Should test if we are still able to read and work with old MF files 
+% 
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+a = ita_read(which('AD22-1.SPK'));
+% 
+% close(f1);
+% close(f2);
+end
\ No newline at end of file

Debug/test_ita_class.m

+function test_ita_class()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+%% RSC - test ita_struct2audio and ita_audio2struct, test if conversion is transparent
+noise = ita_generate('noise',1,44100,3); 
+
+timedata = noise.timeData;
+time = noise.time;
+freqdata = noise.freqData;
+freq = noise.freq;
+
+freqVector = noise.freqVector;
+samplingRate = noise.samplingRate;
+signalType = noise.signalType;
+fftDegree = noise.fftDegree;
+trackLength = noise.trackLength;
+
+clear timedata time freqdata freq freqVector samplingRate signalType fftDegree trackLength
+end
\ No newline at end of file

Debug/test_ita_convolve.m

+function test_ita_convolve
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+% Testinstanzen
+I = [6 6
+    8 6
+    12 6
+    ];
+
+k = size(I,1);
+
+for i=1:k
+   m = I(i,1);
+   n = I(i,2);
+
+   A = ita_generate('noise',1,44100,m);
+   B = ita_generate('impulse',1,44100,n);
+
+   C1 = ita_convolve(A,B,'overlap_add',false);
+   C1 = ita_ifft(C1);
+   
+   C2 = ita_convolve(A,B,'overlap_add',true);
+   C2 = ita_ifft(C2);
+
+   if ~(C1 == C2)
+      error('test_convolve_overlap_add: Overlap-add and linear convolution dont return same results') 
+   end
+   
+end
+end
\ No newline at end of file

Debug/test_ita_coordinates.m

+function test_ita_coordinates()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+
+a = ita_generate('noise',1,44100,15);
+a = ita_split(a,[1 1 1 1]);
+a.channelCoordinates = itaCoordinates([1 5 9; 2 6 10; 3 7 11; 4 8 12],'cart');
+%a.channelCoordinates = build_search_database(a.channelCoordinates);
+
+[ind, dist] = findnearest(a.channelCoordinates,[2 2 2; 3 3 3],'cart',1); %Check matlab search (few elements)
+if ~ind == [2; 3] %#ok<BDSCA>
+    error('findnearest@itaCoordinates does not work');
+end
+
+[ind, dist] = findnearest([a.channelCoordinates a.channelCoordinates a.channelCoordinates],[2 2 2; 3 3 3],'cart',1); % Check mex-search (lots of elements)
+
+
+
+coords = a.channelCoordinates;
+tmp = coords.makeSph;
+coords_new = tmp.makeCart;
+coords_new.cart - coords.cart
+
+
+tmp = coords.makePol;
+coords_new = tmp.makeCart;
+coords_new.cart - coords.cart
+
+coords = coords.makeSph;
+tmp = coords.makePol;
+coords_new = tmp.makeSph;
+coords_new.sph - coords.sph
+
+end
\ No newline at end of file

Debug/test_ita_dsp.m

+function test_ita_dsp()
+%get two signals and a tp
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved. 
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder. 
+% </ITA-Toolbox>
+
+sr = 44100;
+fftdeg = 10;
+a       = ita_generate('flatnoise',1,sr,fftdeg);
+a       = ita_merge(a,a,a);
+
+b = ita_mpb_filter(a,[20,7000]);
+b = ita_mpb_filter(a,[20,30000]);
+b = ita_mpb_filter(a,[0,7000]);
+b = ita_mpb_filter(a,[20 0]);
+a = ita_time_shift(a,20,'samples');
+a = ita_time_window(a,[0.001 0.002],'time');
+a = ita_minimumphase(a);
+a = ita_zerophase(a);
+a = ita_integrate(a);
+% a = ita_apparentmass2impedance(a);
+
+%test stfft and inverse
+% if ~(ita_ifft(a) == ita_istfft(ita_stfft(a,'blocksize',256)))
+%     disp(['Difference is ' num2str(max(max(res.dat)))]);
+%     error(['ita_stfft or ita_istfft are not working']);
+%     
+% end
+end
\ No newline at end of file

Debug/test_ita_edit.m

+function varargout = test_ita_edit()
+
+% <ITA-Toolbox>
+% This file is part of the ITA-Toolbox. Some rights reserved.
+% You can find the license for this m-file in the license.txt file in the ITA-Toolbox folder.
+% </ITA-Toolbox>
+
+%% get three signals and a tp
+sr = 44100;
+fftdeg = 10;
+a       = ita_generate('sine',1,1000,sr,fftdeg);
+b       = ita_generate('noise',1,sr,fftdeg);
+c       = ita_merge(a,b);
+d       = ita_generate_sweep('mode','exp','freqRange',[2 22000],'samplingRate',44100,'fftDegree',15);
+tp      = ita_fft(ita_merge(a,a,a,a));
+tp.signalType = 'energy';
+tparray(1) = tp;
+tparray(2) = tp;
+
+%% used to test the transformation at the beginning of functions
+timeData = b;
+freqData = ita_fft(b);
+
+%% test functions
+%x = ita_import(a.dat,44100,'t');
+
+e = ita_abs(a);
+
+e = ita_amplify(c,6);
+e = ita_amplify(c,'6');
+e = ita_amplify(c,'6dB');
+if ~(a == ita_amplify(ita_amplify(a,'3dB'),'-3dB')) 
+    error('ITA_AMPLIFY:Oh Lord. Function seems not to work properly.')
+end
+
+e = ita_append(timeData,freqData);
+e = ita_append(a,a);
+if ~([a.dat a.dat] == e.dat) 
+    error('ITA_APPEND:Oh Lord. Function seems not to work properly.')
+end
+e = ita_convolve(a,b); 
+% 
+% e = ita_CTC_filter(tp);
+
+e = ita_envelope(c);
+
+e = ita_extend_dat(c);
+e = ita_extend_dat(c,10000);             %indirect call of extract
+e = ita_extend_dat(c,100000);
+e = ita_extend_dat(c,10000,'symmetric'); %indirect call of extract
+e = ita_extend_dat(c,100000,'symmetric');
+e = ita_extend_dat(a,d);
+
+e = ita_extract_dat(c,50);
+e = ita_extract_dat(c,5,'symmetric');
+e = ita_extract_dat(c,5);
+e = ita_extract_dat(c,4,'firstsample',16);
+e = ita_extract_dat(c,6,'forcesamples');
+if e.nSamples ~= 6
+   error('forcesamples does not work') 
+end
+e = ita_extract_dat(c,4,'random');
+if ~(a == ita_extract_dat(ita_extend_dat(a)))
+    error('ITA_EXTRACT_DAT/ITA_EXTEND_DAT:Oh Lord. Function seems not to work properly.')
+end
+
+e = ita_imag(a);
+
+e = ita_make_filter([100 1000],sr,20);
+e = ita_make_filter([100 1000],sr,fftdeg);
+e = ita_make_filter([100 1000],c);
+
+e = ita_mean(tp);
+tp.channelNames{1} = 'test';
+e = ita_mean(ita_merge(tp,tp),'same_channelnames_only');
+if ~(a == ita_mean(ita_merge(a,a,a,a)))
+    error('ITA_MEAN:Oh Lord. Function seems not to work properly.')
+end
+e = ita_mean(tparray);
+
+
+e = ita_median(tp);
+if ~(a == ita_median(ita_merge(a,a,a,a)))
+    ita_verbose_info('ITA_MEDIAN:Oh Lord. Function seems not to work properly.',0)
+end
+
+e = ita_merge(ita_fft(a),a);
+e = ita_merge(ita_fft(a),b);
+e = ita_merge(a,d);
+
+e = ita_minimumphase(c);
+
+e = ita_mpb_filter(c,[0 2000]);
+e = ita_mpb_filter(c,[0 2000],'order',2);
+e = ita_mpb_filter(c,[0 2000],'order',2,'class',0);
+e = ita_mpb_filter(c,[0 2000],'zerophase');
+e = ita_mpb_filter(c,[0 2000],'minimumphase'); 
+ 
+e = ita_mpb_filter(a,'octaves',3);
+e = ita_mpb_filter(a,'oct',3);
+
+e = ita_mpb_filter(c,'c-weight');
+e = ita_mpb_filter(c,'a-weight');
+
+e = ita_mpb_filter(c,[2000 0]);
+
+% e = ita_multiple_time_windows(a,'blocksize',1024);
+% e = ita_multiple_time_windows(ita_merge(a,b),'blocksize',1024,'overlap',0.1);
+% e = ita_multiple_time_windows(ita_merge(a,b),'blocksize',1024,'overlap',0.1,'channels',[2]);