function o = horzcat(varargin)

% Overloads the horzcat method for dates objects.
%
% INPUTS
% - varargin [dates]
%
% OUTPUTS
% - o [dates] object containing dates defined in varargin{:}

% Copyright © 2013-2023 Dynare Team
%
% This code 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 dates submodule 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/>.

if ~all(cellfun(@isdates, varargin))
    error('dates:horzcat:ArgCheck', 'All input arguments must be dates objects.')
end

o = copy(varargin{1});
for i=2:nargin
    p = varargin{i};
    if isequal(o.freq, p.freq) || isempty(o.freq) || isempty(p.freq)
        if isempty(o)
            o = copy(p);
            continue
        end
        if ~isempty(p)
            o.time = [o.time; p.time];
        end
    else
        error('dates:horzcat', 'All input arguments must have the same frequency.')
    end
end

return % --*-- Unit tests --*--

%@test:1
% Define some dates
B1 = '1953Q4';
B2 = '1950Q2';
B3 = '1950Q1';
B4 = '1945Q3';
B5 = '2009Q2';

% Define expected results.
e.time = [1945*4+3; 1950*4+1; 1950*4+2; 1953*4+4; 2009*4+2];
e.freq = 4;

% Call the tested routine.
d = dates(B4,B3,B2,B1);
d = [d, dates(B5)];

% Check the results.
t(1) = dassert(d.time,e.time);
t(2) = dassert(d.freq,e.freq);
t(3) = size(e.time,1)==d.ndat();
T = all(t);
%@eof:1

%@test:2
% Define some dates
B1 = '1953Q4';
B2 = '1950Q2';
B3 = '1950Q1';
B4 = '1945Q3';
B5 = '2009Q2';

% Define expected results.
e.time = [1945*4+3; 1950*4+1; 1950*4+2; 1953*4+4; 2009*4+2];
e.freq = 4;

% Call the tested routine.
d = dates(B4,B3,B2);
d = [d, dates(B1), dates(B5)];

% Check the results.
t(1) = dassert(d.time,e.time);
t(2) = dassert(d.freq,e.freq);
t(3) = size(e.time,1)==d.ndat();
T = all(t);
%@eof:2

%@test:3
% Define some dates
B1 = '1953Q4';
B2 = '1950Q2';
B3 = '1950Q1';
B4 = '1945Q3';
B5 = '2009Q2';

% Define expected results.
e.time = [1945*4+3; 1950*4+1; 1950*4+2; 1953*4+4; 2009*4+2];
e.freq = 4;

% Call the tested routine.
d = dates(B4,B3,B2);
d = [d, dates(B1,B5)];

% Check the results.
t(1) = dassert(d.time,e.time);
t(2) = dassert(d.freq,e.freq);
t(3) = size(e.time,1)==d.ndat();
T = all(t);
%@eof:3

%@test:4
% Define some dates
B1 = '1953Q4';
B2 = '1950Q2';
B3 = '1950Q1';
B4 = '1945Q3';
B5 = '2009Q2';

% Define expected results.
e.time = [1945*4+3; 1950*4+1; 1950*4+2; 1953*4+4; 2009*4+2];
e.freq = 4;

% Call the tested routine.
d = dates(B4,B3,B2);
d = [d, [dates(B1), dates(B5)]];

% Check the results.
t(1) = dassert(d.time,e.time);
t(2) = dassert(d.freq,e.freq);
t(3) = size(e.time,1)==d.ndat();
T = all(t);
%@eof:4

%@test:5
try
    e = [dates('1938M11') dates()];
    t(1) = true;
catch
    t(1) = false;
end

if t(1)
    t(2) = e==dates('1938M11');
end

T = all(t);
%@eof:5

%@test:6
try
    e = [dates() dates('1938M11')];
    t(1) = true;
catch
    t(1) = false;
end

if t(1)
    t(2) = e==dates('1938M11');
end

T = all(t);
%@eof:6
