function [report, time] = run_unitary_tests(listoffiles)

% Runs unitary tests defined in a collection of files.
%
% INPUTS
%  - listoffiles [cell], The list of m files (with path) where the unitary tests are written.
%                This cell, with n elements, is the output of get_directory_description routine.
%
% OUTPUTS
%  - report      [cell], Results of the unitary tests (n rows and 5 columns). Each row stores
%                the second output argument of mtest routine (info).
%  - time        [double], Current date and time as date vector.
%
%
% See also get_directory_description, mtest

% Copyright (C) 2013-2017 Dynare Team
%
% This file is part of Dynare.
%
% Dynare is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% Dynare is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with Dynare.  If not, see <https://www.gnu.org/licenses/>.

report = {};

testcoverage = zeros(2,1);

if nargout>2
    list_of_routines_without_unit_tests = {};
end

skipline()

for f=1:length(listoffiles)
    if isempty(strfind(listoffiles{f},'.#'))
        if is_unitary_test_available(listoffiles{f})
            testcoverage(1) = testcoverage(1) + 1;
            [check, info] = mtest(listoffiles{f});
            r0 = sprintf('[%s/%s]',num2str(sum([info{:,3}])),num2str(size(info,1)));
            if check
                r1 = 'PASSED';
            else
                r1 = 'FAILED';
            end
            disp(sprintf('***** Unitary tests in %s %s %s!', listoffiles{f}, r0, r1));
            report = [report; info];
        else
            testcoverage(2) = testcoverage(2) + 1;
            list_of_routines_without_unit_tests(testcoverage(2)) = { listoffiles{f} };
        end
    end
end

atc = 100*testcoverage(1)/sum(testcoverage);

skipline(2)
disp(sprintf('Approximated test coverage is %s%%',num2str(atc)))

if nargout>1
    time = clock;
end
