From f2960bf774113d278eeb43faf9a30c31495e6843 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Ry=C3=BBk=29?=
 <stepan@adjemian.eu>
Date: Sun, 1 Jan 2023 20:16:22 +0100
Subject: [PATCH] Use new unit test interface.

---
 src/@arima/arima.m                            |   6 +-
 src/@dates/append.m                           |   6 +-
 src/@dates/append_.m                          |   6 +-
 src/@dates/char.m                             |   6 +-
 src/@dates/colon.m                            |   8 +-
 src/@dates/copy.m                             |   8 +-
 src/@dates/dates.m                            |   9 +-
 src/@dates/double.m                           |   6 +-
 src/@dates/end.m                              |   6 +-
 src/@dates/eq.m                               |   8 +-
 src/@dates/ge.m                               |   6 +-
 src/@dates/gt.m                               |   6 +-
 src/@dates/horzcat.m                          |   8 +-
 src/@dates/intersect.m                        |   6 +-
 src/@dates/isempty.m                          |   6 +-
 src/@dates/isequal.m                          |   6 +-
 src/@dates/le.m                               |   6 +-
 src/@dates/length.m                           |   6 +-
 src/@dates/lt.m                               |   6 +-
 src/@dates/max.m                              |   6 +-
 src/@dates/min.m                              |   6 +-
 src/@dates/minus.m                            |   6 +-
 src/@dates/mtimes.m                           |   6 +-
 src/@dates/ndat.m                             |   8 +-
 src/@dates/ne.m                               |   6 +-
 src/@dates/plus.m                             |   6 +-
 src/@dates/pop.m                              |   6 +-
 src/@dates/pop_.m                             |   6 +-
 src/@dates/private/comparison_arg_checks.m    |   8 +-
 src/@dates/private/dates4display.m            | 190 +++++++++---------
 src/@dates/remove.m                           |   6 +-
 src/@dates/remove_.m                          |   6 +-
 src/@dates/setdiff.m                          |   6 +-
 src/@dates/sort.m                             |   6 +-
 src/@dates/sort_.m                            |   8 +-
 src/@dates/strings.m                          |   6 +-
 src/@dates/subsasgn.m                         |   8 +-
 src/@dates/subsref.m                          |   8 +-
 src/@dates/uminus.m                           |   8 +-
 src/@dates/union.m                            |  10 +-
 src/@dates/unique.m                           |   8 +-
 src/@dates/unique_.m                          |   6 +-
 src/@dates/uplus.m                            |   8 +-
 src/@dates/vertcat.m                          |   6 +-
 src/@dates/year.m                             |   6 +-
 src/@dseries/abs.m                            |   6 +-
 src/@dseries/abs_.m                           |   8 +-
 src/@dseries/align.m                          |   8 +-
 src/@dseries/align_.m                         |   8 +-
 src/@dseries/backcast.m                       |   8 +-
 src/@dseries/backcast_.m                      |   8 +-
 src/@dseries/baxter_king_filter.m             |   6 +-
 src/@dseries/baxter_king_filter_.m            |   6 +-
 src/@dseries/center.m                         |   6 +-
 src/@dseries/center_.m                        |   6 +-
 src/@dseries/chain.m                          |   8 +-
 src/@dseries/chain_.m                         |   8 +-
 src/@dseries/copy.m                           |   8 +-
 src/@dseries/cumprod.m                        |   6 +-
 src/@dseries/cumprod_.m                       |   8 +-
 src/@dseries/cumsum.m                         |   8 +-
 src/@dseries/cumsum_.m                        |   8 +-
 src/@dseries/detrend.m                        |   8 +-
 src/@dseries/detrend_.m                       |   8 +-
 src/@dseries/dgrowth.m                        |   8 +-
 src/@dseries/dgrowth_.m                       |   8 +-
 src/@dseries/diff.m                           |   8 +-
 src/@dseries/diff_.m                          |   8 +-
 src/@dseries/dseries.m                        |   6 +-
 src/@dseries/eq.m                             |   8 +-
 src/@dseries/exist.m                          |   6 +-
 src/@dseries/exp.m                            |   8 +-
 src/@dseries/exp_.m                           |   8 +-
 src/@dseries/extract.m                        |   8 +-
 src/@dseries/firstdate.m                      |   6 +-
 src/@dseries/firstobservedperiod.m            |   6 +-
 src/@dseries/flip.m                           |   9 +-
 src/@dseries/flip_.m                          |  71 +++----
 src/@dseries/ge.m                             |   6 +-
 src/@dseries/gt.m                             |   6 +-
 src/@dseries/hdiff.m                          |   8 +-
 src/@dseries/hdiff_.m                         |   8 +-
 src/@dseries/hgrowth.m                        |   8 +-
 src/@dseries/hgrowth_.m                       |   6 +-
 src/@dseries/horzcat.m                        |   6 +-
 src/@dseries/hpcycle.m                        |   8 +-
 src/@dseries/hpcycle_.m                       |   8 +-
 src/@dseries/hptrend.m                        |   8 +-
 src/@dseries/hptrend_.m                       |   8 +-
 src/@dseries/insert.m                         |   8 +-
 src/@dseries/isinf.m                          |   6 +-
 src/@dseries/isnan.m                          |   6 +-
 src/@dseries/isreal.m                         |   6 +-
 src/@dseries/lag.m                            |   8 +-
 src/@dseries/lag_.m                           |   8 +-
 src/@dseries/lastobservedperiod.m             |   8 +-
 src/@dseries/le.m                             |   6 +-
 src/@dseries/lead.m                           |   8 +-
 src/@dseries/lead_.m                          |   8 +-
 src/@dseries/length.m                         |   8 +-
 src/@dseries/lineartrend.m                    |  12 +-
 src/@dseries/log.m                            |   8 +-
 src/@dseries/log_.m                           |   8 +-
 src/@dseries/lt.m                             |   6 +-
 src/@dseries/mdiff.m                          |   8 +-
 src/@dseries/mdiff_.m                         |   6 +-
 src/@dseries/mean.m                           |   8 +-
 src/@dseries/merge.m                          |  10 +-
 src/@dseries/mgrowth.m                        |   8 +-
 src/@dseries/mgrowth_.m                       |   6 +-
 src/@dseries/minus.m                          |   6 +-
 src/@dseries/mpower.m                         |   8 +-
 src/@dseries/mrdivide.m                       |   6 +-
 src/@dseries/mtimes.m                         |   6 +-
 src/@dseries/nanmean.m                        |   8 +-
 src/@dseries/nanstd.m                         |   8 +-
 src/@dseries/ne.m                             |   8 +-
 src/@dseries/onesidedhpcycle.m                |   8 +-
 src/@dseries/onesidedhpcycle_.m               |  10 +-
 src/@dseries/onesidedhptrend.m                |   8 +-
 src/@dseries/onesidedhptrend_.m               |  10 +-
 src/@dseries/plus.m                           |  10 +-
 src/@dseries/pop.m                            |   6 +-
 src/@dseries/pop_.m                           |  10 +-
 src/@dseries/projection.m                     |  23 ++-
 src/@dseries/qdiff.m                          |   8 +-
 src/@dseries/qdiff_.m                         |   6 +-
 src/@dseries/qgrowth.m                        |   8 +-
 src/@dseries/qgrowth_.m                       |   6 +-
 src/@dseries/remove.m                         |   6 +-
 src/@dseries/remove_.m                        |   6 +-
 src/@dseries/rename.m                         |   8 +-
 src/@dseries/rename_.m                        |   6 +-
 src/@dseries/round.m                          |   8 +-
 src/@dseries/round_.m                         |   8 +-
 src/@dseries/save.m                           |   8 +-
 src/@dseries/set_names.m                      |   8 +-
 src/@dseries/size.m                           |  10 +-
 src/@dseries/std.m                            |   8 +-
 src/@dseries/subsample.m                      |   6 +-
 src/@dseries/subsasgn.m                       |   8 +-
 src/@dseries/subsref.m                        |   6 +-
 src/@dseries/tag.m                            |  10 +-
 src/@dseries/tex_rename.m                     |   8 +-
 src/@dseries/tex_rename_.m                    |   8 +-
 src/@dseries/uminus.m                         |   6 +-
 src/@dseries/vertcat.m                        |   6 +-
 src/@dseries/ydiff.m                          |   8 +-
 src/@dseries/ydiff_.m                         |   8 +-
 src/@dseries/ygrowth.m                        |   8 +-
 src/@dseries/ygrowth_.m                       |  10 +-
 src/@x13/clean.m                              |   8 +-
 src/@x13/print.m                              |   8 +-
 src/@x13/subsasgn.m                           |   8 +-
 src/@x13/subsref.m                            |   8 +-
 src/@x13/x13.m                                |   9 +-
 src/mdbnomics2dseries/mdbnomics2dseries.m     | 183 ++++++++---------
 src/read/load_csv_file_data.m                 |   6 +-
 src/read/load_m_file_data.m                   |   6 +-
 src/read/load_mat_file_data.m                 |   8 +-
 src/utilities/convert/date2string.m           |   8 +-
 src/utilities/convert/dates2M.m               |   8 +-
 src/utilities/convert/dates2Q.m               |   8 +-
 src/utilities/convert/dates2S.m               |   8 +-
 src/utilities/convert/dates2Y.m               |   8 +-
 src/utilities/convert/dseries2M.m             |   6 +-
 src/utilities/convert/dseries2Q.m             |   8 +-
 src/utilities/convert/dseries2S.m             |   8 +-
 src/utilities/convert/dseries2Y.m             |   8 +-
 src/utilities/convert/freq2string.m           |   6 +-
 src/utilities/convert/string2date.m           |   8 +-
 src/utilities/convert/string2freq.m           |   6 +-
 src/utilities/from/from.m                     |   6 +-
 src/utilities/get/get_series_start_date.m     |   8 +-
 .../insert/insert_column_vector_in_a_matrix.m |   8 +-
 ...t_object_in_a_one_dimensional_cell_array.m |   8 +-
 src/utilities/is/isbiannual.m                 |   8 +-
 src/utilities/is/isdaily.m                    |   8 +-
 src/utilities/is/isdate.m                     |   8 +-
 src/utilities/is/isdates.m                    |   8 +-
 src/utilities/is/isfreq.m                     |   8 +-
 src/utilities/is/isleapyear.m                 |   8 +-
 src/utilities/is/islongyear.m                 |   8 +-
 src/utilities/is/ismonthly.m                  |   8 +-
 src/utilities/is/isquarterly.m                |   8 +-
 src/utilities/is/isstringdate.m               |   8 +-
 src/utilities/is/issubperiod.m                |   8 +-
 src/utilities/is/isvalidday.m                 |   8 +-
 src/utilities/is/isweekly.m                   |   8 +-
 src/utilities/is/isyearly.m                   |   8 +-
 src/utilities/misc/dummy.m                    |   8 +-
 src/utilities/misc/firstdayofyear.m           |   8 +-
 src/utilities/misc/firstweekofyear.m          |  14 +-
 src/utilities/misc/weeksince.m                |  12 +-
 src/utilities/missing/isint/isint.m           |   8 +-
 src/utilities/missing/shiftS/shiftS.m         |   6 +-
 .../op/add_periods_to_array_of_dates.m        |   8 +-
 src/utilities/op/add_periods_to_date.m        |   8 +-
 src/utilities/str/name2tex.m                  |   6 +-
 src/utilities/x13/checkcommandcompatibility.m |   8 +-
 src/utilities/x13/checkoptioncompatibility.m  |   8 +-
 src/utilities/x13/isoption.m                  |   8 +-
 src/utilities/x13/printoption.m               |   8 +-
 src/utilities/x13/printspan.m                 |   8 +-
 src/utilities/x13/printstart.m                |   6 +-
 src/utilities/x13/setdefaultmember.m          |   8 +-
 tests/runalltests.m                           |   4 +-
 207 files changed, 998 insertions(+), 978 deletions(-)

diff --git a/src/@arima/arima.m b/src/@arima/arima.m
index d6dc8da..4254641 100644
--- a/src/@arima/arima.m
+++ b/src/@arima/arima.m
@@ -1,8 +1,8 @@
-classdef arima<handle % --*-- Unitary tests --*--
+classdef arima<handle
 
 % Class for ARIMA models (interface to X13).
 
-% Copyright © 2017-2022 Dynare Team
+% Copyright © 2017-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
@@ -151,7 +151,7 @@ methods
         end
         end % arima
 end % methods
-end % classdef
+end % classdef --*-- Unit tests --*--
 
 %@test:1
 %$ try
diff --git a/src/@dates/append.m b/src/@dates/append.m
index 6e3f412..696a321 100644
--- a/src/@dates/append.m
+++ b/src/@dates/append.m
@@ -1,4 +1,4 @@
-function o = append(o, d) % --*-- Unitary tests --*--
+function o = append(o, d)
 
 % append method for dates class.
 %
@@ -9,7 +9,7 @@ function o = append(o, d) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates]       containing dates defined in o and d.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2022 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
@@ -42,7 +42,7 @@ end
 o = copy(o);
 o.append_(d);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/append_.m b/src/@dates/append_.m
index 398f032..db4e1db 100644
--- a/src/@dates/append_.m
+++ b/src/@dates/append_.m
@@ -1,4 +1,4 @@
-function o = append_(o, d) % --*-- Unitary tests --*--
+function o = append_(o, d)
 
 % append method for dates class (in place modification).
 %
@@ -9,7 +9,7 @@ function o = append_(o, d) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates] dates object containing dates defined in o and d.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2022 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
@@ -41,7 +41,7 @@ end
 
 o.time = [o.time; d.time];
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/char.m b/src/@dates/char.m
index b68b2e1..b6f03f5 100644
--- a/src/@dates/char.m
+++ b/src/@dates/char.m
@@ -1,4 +1,4 @@
-function s = char(o) % --*-- Unitary tests --*--
+function s = char(o)
 
 % Given a one element dates object, returns a string with the formatted date.
 %
@@ -8,7 +8,7 @@ function s = char(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - s  [string]
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-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
@@ -29,7 +29,7 @@ end
 
 s = date2string(o.time, o.freq);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
diff --git a/src/@dates/colon.m b/src/@dates/colon.m
index 0e0035c..e400eee 100644
--- a/src/@dates/colon.m
+++ b/src/@dates/colon.m
@@ -1,4 +1,4 @@
-function q = colon(varargin) % --*-- Unitary tests --*--
+function q = colon(varargin)
 
 % Overloads the colon operator (:). This method can be used to create ranges of dates.
 %
@@ -14,7 +14,7 @@ function q = colon(varargin) % --*-- Unitary tests --*--
 % 1. p must be greater than o if d>0.
 % 2. p and q are dates objects with one element.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -72,7 +72,7 @@ q.freq = o.freq;
 % Set time
 q.time = transpose(o.time:d:p.time);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define two dates
@@ -172,4 +172,4 @@ d = dd(1950, 01, 28):dd(1950, 01, 30);
 t(1) = isequal(d.time(:,1), e.time(:,1));
 t(2) = isequal(d.freq, e.freq);
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/@dates/copy.m b/src/@dates/copy.m
index 59a79ab..2abcca9 100644
--- a/src/@dates/copy.m
+++ b/src/@dates/copy.m
@@ -1,4 +1,4 @@
-function p = copy(o) % --*-- Unitary tests --*--
+function p = copy(o)
 
 % Do a copy of a dates object.
 %
@@ -8,7 +8,7 @@ function p = copy(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - p [dates]
 
-% Copyright © 2014-2020 Dynare Team
+% Copyright © 2014-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
@@ -27,7 +27,7 @@ p = dates();
 p.freq = o.freq;
 p.time = o.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dates/dates.m b/src/@dates/dates.m
index fe98faa..549bb24 100644
--- a/src/@dates/dates.m
+++ b/src/@dates/dates.m
@@ -1,6 +1,6 @@
-classdef dates<handle % --*-- Unitary tests --*--
+classdef dates<handle
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-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
@@ -193,8 +193,7 @@ methods
     n = ndat(o);
     n = length(o);
 end % methods
-end % classdef
-
+end % classdef --*-- Unit tests --*--
 
 %@test:1
 %$ % Define some dates
@@ -378,4 +377,4 @@ end % classdef
 %$ end
 %$
 %$ T = all(t);
-%@eof:10
\ No newline at end of file
+%@eof:10
diff --git a/src/@dates/double.m b/src/@dates/double.m
index 8a96eed..26a4043 100644
--- a/src/@dates/double.m
+++ b/src/@dates/double.m
@@ -1,4 +1,4 @@
-function [m, f]  = double(o) % --*-- Unitary tests --*--
+function [m, f]  = double(o)
 
 % Returns a vector of doubles with the fractional part corresponding
 % to the subperiod. Used for plots and to store dates in a matrix.
@@ -13,7 +13,7 @@ function [m, f]  = double(o) % --*-- Unitary tests --*--
 % REMARKS
 %  Obviously the frequency is lost during the conversion.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -45,7 +45,7 @@ if nargout>1
     f = o.freq;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
diff --git a/src/@dates/end.m b/src/@dates/end.m
index b0913dd..4c7b928 100644
--- a/src/@dates/end.m
+++ b/src/@dates/end.m
@@ -1,4 +1,4 @@
-function lastIndex = end(o, k, n) % --*-- Unitary tests --*--
+function lastIndex = end(o, k, n)
 
 % Overloads end keyword.
 %
@@ -10,7 +10,7 @@ function lastIndex = end(o, k, n) % --*-- Unitary tests --*--
 % OUTPUTS
 %   lastIndex [integer] last dates index
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -28,7 +28,7 @@ function lastIndex = end(o, k, n) % --*-- Unitary tests --*--
 assert(k==1 && n==1, 'dates:end:ArgCheck', 'dates only has one dimension');
 lastIndex = o.ndat();
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
diff --git a/src/@dates/eq.m b/src/@dates/eq.m
index 9d64d88..f8b6362 100644
--- a/src/@dates/eq.m
+++ b/src/@dates/eq.m
@@ -1,4 +1,4 @@
-function l = eq(varargin) % --*-- Unitary tests --*--
+function l = eq(varargin)
 
 % Overloads == operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = eq(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2020 Dynare Team
+% Copyright © 2013-2022 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
@@ -37,7 +37,7 @@ else
     l = bsxfun(@eq,o.time(:,1),p.time(:,1));
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
@@ -134,4 +134,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dates/ge.m b/src/@dates/ge.m
index 96b201e..0c1d7c8 100644
--- a/src/@dates/ge.m
+++ b/src/@dates/ge.m
@@ -1,4 +1,4 @@
-function l = ge(varargin)  % --*-- Unitary tests --*--
+function l = ge(varargin)
 
 % Overloads the >= operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = ge(varargin)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -28,7 +28,7 @@ function l = ge(varargin)  % --*-- Unitary tests --*--
 
 l = o.time>=p.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/gt.m b/src/@dates/gt.m
index 52a6aed..cb73869 100644
--- a/src/@dates/gt.m
+++ b/src/@dates/gt.m
@@ -1,4 +1,4 @@
-function l = gt(varargin)  % --*-- Unitary tests --*--
+function l = gt(varargin)
 
 % Overloads the > operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = gt(varargin)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -28,7 +28,7 @@ function l = gt(varargin)  % --*-- Unitary tests --*--
 
 l = o.time>p.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/horzcat.m b/src/@dates/horzcat.m
index abda27b..0247033 100644
--- a/src/@dates/horzcat.m
+++ b/src/@dates/horzcat.m
@@ -1,4 +1,4 @@
-function o = horzcat(varargin) % --*-- Unitary tests --*--
+function o = horzcat(varargin)
 
 % Overloads the horzcat method for dates objects.
 %
@@ -8,7 +8,7 @@ function o = horzcat(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates] object containing dates defined in varargin{:}
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -43,7 +43,7 @@ for i=2:nargin
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -165,4 +165,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/@dates/intersect.m b/src/@dates/intersect.m
index 61aca8c..cace905 100644
--- a/src/@dates/intersect.m
+++ b/src/@dates/intersect.m
@@ -1,4 +1,4 @@
-function q = intersect(o, p) % --*-- Unitary tests --*--
+function q = intersect(o, p)
 
 % Overloads intersect function for dates objects.
 %
@@ -9,7 +9,7 @@ function q = intersect(o, p) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q [dates] All the common elements in o and p.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -48,7 +48,7 @@ end
 q.freq = o.freq;
 q.time = time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/isempty.m b/src/@dates/isempty.m
index 770d0a4..97e59bf 100644
--- a/src/@dates/isempty.m
+++ b/src/@dates/isempty.m
@@ -1,4 +1,4 @@
-function l = isempty(o) % --*-- Unitary tests --*--
+function l = isempty(o)
 
 % Returns true (1) if and only if o dates object is empty.
 %
@@ -8,7 +8,7 @@ function l = isempty(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -25,7 +25,7 @@ function l = isempty(o) % --*-- Unitary tests --*--
 
 l = isequal(o.ndat(), 0);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Instantiate an empty dates object
diff --git a/src/@dates/isequal.m b/src/@dates/isequal.m
index 6f923a6..9df7829 100644
--- a/src/@dates/isequal.m
+++ b/src/@dates/isequal.m
@@ -1,4 +1,4 @@
-function l = isequal(o, p, fake) % --*-- Unitary tests --*--
+function l = isequal(o, p, fake)
 
 % Overloads isequal function for dates objects. Returns true true iff o and p have the same elements.
 %
@@ -9,7 +9,7 @@ function l = isequal(o, p, fake) % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -30,7 +30,7 @@ end
 
 l = isequal(o.freq, p.freq) && isequal(o.time, p.time);
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 d1 = dates('1938Q1');
diff --git a/src/@dates/le.m b/src/@dates/le.m
index 8368746..2d1f307 100644
--- a/src/@dates/le.m
+++ b/src/@dates/le.m
@@ -1,4 +1,4 @@
-function l = le(varargin)  % --*-- Unitary tests --*--
+function l = le(varargin)
 
 % Overloads the <= operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = le(varargin)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -28,7 +28,7 @@ function l = le(varargin)  % --*-- Unitary tests --*--
 
 l = o.time<=p.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/length.m b/src/@dates/length.m
index 0307cb6..0e2a33f 100644
--- a/src/@dates/length.m
+++ b/src/@dates/length.m
@@ -1,4 +1,4 @@
-function n = length(o) % --*-- Unitary tests --*--
+function n = length(o)
 
 % Returns the number of elements in a dates object.
 %
@@ -8,7 +8,7 @@ function n = length(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - n [integer] Number of elements in o.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -25,7 +25,7 @@ function n = length(o) % --*-- Unitary tests --*--
 
 n = o.ndat();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 d = dates('1938Q1');
diff --git a/src/@dates/lt.m b/src/@dates/lt.m
index 3dc091e..67b07a5 100644
--- a/src/@dates/lt.m
+++ b/src/@dates/lt.m
@@ -1,4 +1,4 @@
-function l = lt(varargin)  % --*-- Unitary tests --*--
+function l = lt(varargin)
 
 % Overloads the < operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = lt(varargin)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -28,7 +28,7 @@ function l = lt(varargin)  % --*-- Unitary tests --*--
 
 l = o.time<p.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/max.m b/src/@dates/max.m
index ecaa3b7..1148bfd 100644
--- a/src/@dates/max.m
+++ b/src/@dates/max.m
@@ -1,4 +1,4 @@
-function q = max(varargin)  % --*-- Unitary tests --*--
+function q = max(varargin)
 
 % Overloads the max function for dates objects.
 %
@@ -8,7 +8,7 @@ function q = max(varargin)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - q        [dates]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -36,7 +36,7 @@ switch nargin
     q = max(horzcat(varargin{:}));
 end
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/min.m b/src/@dates/min.m
index c993a28..70300ab 100644
--- a/src/@dates/min.m
+++ b/src/@dates/min.m
@@ -1,4 +1,4 @@
-function q = min(varargin) % --*-- Unitary tests --*--
+function q = min(varargin)
 
 % Overloads the min function for dates objects.
 %
@@ -8,7 +8,7 @@ function q = min(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q        [dates]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -36,7 +36,7 @@ switch nargin
     q = min(horzcat(varargin{:}));
 end
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/minus.m b/src/@dates/minus.m
index 6f622c2..26fa04a 100644
--- a/src/@dates/minus.m
+++ b/src/@dates/minus.m
@@ -1,4 +1,4 @@
-function q = minus(o,p) % --*-- Unitary tests --*--
+function q = minus(o,p)
 
 % Overloads the minus operator (-). If o and p are dates objects, the method . If
 % one of the inputs is an integer or a vector of integers, the method shifts the dates object by X (the interger input) periods backward.
@@ -18,7 +18,7 @@ function q = minus(o,p) % --*-- Unitary tests --*--
 %    p periods backward.
 % 3. If o is not a dates object, an error is returned.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -68,7 +68,7 @@ else
     error('dates:minus:ArgCheck', 'You should read the manual.')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/mtimes.m b/src/@dates/mtimes.m
index d59bafa..1f8be8d 100644
--- a/src/@dates/mtimes.m
+++ b/src/@dates/mtimes.m
@@ -1,4 +1,4 @@
-function p = mtimes(o,n) % --*-- Unitary tests --*--
+function p = mtimes(o,n)
 
 % Overloads the times operator (*). Returns dates object o replicated n times.
 %
@@ -12,7 +12,7 @@ function p = mtimes(o,n) % --*-- Unitary tests --*--
 % 1. If A = dates('2000Q1'), then B=A*3 is a dates object equal to dates('2000Q1','2000Q1','2000Q1')
 % 2. If A = dates('2003Q1','2009Q2'), then B=A*2 is a dates object equal to dates('2003Q1','2009Q2','2003Q1','2009Q2')
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -34,7 +34,7 @@ end
 p = copy(o);
 p.time = repmat(p.time, n, 1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/ndat.m b/src/@dates/ndat.m
index 83e303a..c791363 100644
--- a/src/@dates/ndat.m
+++ b/src/@dates/ndat.m
@@ -1,4 +1,4 @@
-function s = ndat(o) % --*-- Unitary tests --*--
+function s = ndat(o)
 
 % Given a one element dates object, returns a string with the formatted date.
 %
@@ -8,7 +8,7 @@ function s = ndat(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - s  [integer]
 
-% Copyright © 2015-2021 Dynare Team
+% Copyright © 2015-2022 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
@@ -25,7 +25,7 @@ function s = ndat(o) % --*-- Unitary tests --*--
 
 s = length(o.time);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -79,4 +79,4 @@ if t(1)
     t(2) = dassert(card,10);
 end
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dates/ne.m b/src/@dates/ne.m
index 98e0777..af27f4c 100644
--- a/src/@dates/ne.m
+++ b/src/@dates/ne.m
@@ -1,4 +1,4 @@
-function l = ne(varargin) % --*-- Unitary tests --*--
+function l = ne(varargin)
 
 % Overloads ~= operator for dates objects.
 %
@@ -9,7 +9,7 @@ function l = ne(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - l [logical] column vector of max(n,1) elements (zeros or ones).
 
-% Copyright © 2013-2020 Dynare Team
+% 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
@@ -28,7 +28,7 @@ function l = ne(varargin) % --*-- Unitary tests --*--
 
 l = o.time~=p.time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/plus.m b/src/@dates/plus.m
index d1bc9d4..6ebc307 100644
--- a/src/@dates/plus.m
+++ b/src/@dates/plus.m
@@ -1,4 +1,4 @@
-function q = plus(o, p) % --*-- Unitary tests --*--
+function q = plus(o, p)
 
 % Overloads the plus (+) binary operator.
 %
@@ -14,7 +14,7 @@ function q = plus(o, p) % --*-- Unitary tests --*--
 % 2. If  one of the inputs is an integer or a vector of integers, the method shifts the dates
 %    object by X (the interger input) periods forward.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -65,7 +65,7 @@ else
     error('dates:plus:ArgCheck','Please read the manual.')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/pop.m b/src/@dates/pop.m
index 3d48776..821af24 100644
--- a/src/@dates/pop.m
+++ b/src/@dates/pop.m
@@ -1,4 +1,4 @@
-function o = pop(o, p) % --*-- Unitary tests --*--
+function o = pop(o, p)
 
 % pop method for dates class (removes a date).
 %
@@ -15,7 +15,7 @@ function o = pop(o, p) % --*-- Unitary tests --*--
 %
 % See also remove, setdiff.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -38,7 +38,7 @@ else
     o.pop_();
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/pop_.m b/src/@dates/pop_.m
index 068ce4c..4a359e8 100644
--- a/src/@dates/pop_.m
+++ b/src/@dates/pop_.m
@@ -1,4 +1,4 @@
-function o = pop_(o, p) % --*-- Unitary tests --*--
+function o = pop_(o, p)
 
 % pop method for dates class (in place modification).
 %
@@ -15,7 +15,7 @@ function o = pop_(o, p) % --*-- Unitary tests --*--
 %
 % See also remove, setdiff.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -66,7 +66,7 @@ else
     o.time = o.time(jdx);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/private/comparison_arg_checks.m b/src/@dates/private/comparison_arg_checks.m
index 3a225ae..adc7041 100644
--- a/src/@dates/private/comparison_arg_checks.m
+++ b/src/@dates/private/comparison_arg_checks.m
@@ -1,4 +1,4 @@
-function [o, p] = comparison_arg_checks(varargin) % --*-- Unitary tests --*--
+function [o, p] = comparison_arg_checks(varargin)
 
 % Returns two dates objects or an error if objects to be compared are not compatible.
 %
@@ -9,7 +9,7 @@ function [o, p] = comparison_arg_checks(varargin) % --*-- Unitary tests --*--
 % - o [dates] dates object with n or 1 elements.
 % - p [dates] dates object with n or 1 elements.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-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
@@ -49,7 +49,7 @@ end
 o = varargin{1};
 p = varargin{2};
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 OPATH = pwd();
@@ -166,4 +166,4 @@ end
 
 T = all(t);
 cd(OPATH);
-%@eof:7
\ No newline at end of file
+%@eof:7
diff --git a/src/@dates/private/dates4display.m b/src/@dates/private/dates4display.m
index 1a7ccd5..12df3e8 100644
--- a/src/@dates/private/dates4display.m
+++ b/src/@dates/private/dates4display.m
@@ -1,4 +1,4 @@
-function str = dates4display(o, name, max_number_of_elements) % --*-- Unitary tests --*--
+function str = dates4display(o, name, max_number_of_elements)
 
 % Converts a list object to a string.
 %
@@ -10,7 +10,7 @@ function str = dates4display(o, name, max_number_of_elements) % --*-- Unitary te
 % OUTPUTS
 % - str  [string] Representation of the dates object as a string.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-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
@@ -47,106 +47,108 @@ else
 end
 str = sprintf('%s%s>\n', str, date2string(o.time(o.length()), o.freq));
 
+return % --*-- Unit tests --*--
+
 %@test:1
-%$ OPATH = pwd();
-%$ DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
-%$ cd([DSERIES_PATH '/@dates/private']);
-%$
-%$ try
-%$     toto = dates();
-%$     str = dates4display(toto, 'toto', 5);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     expected_str_1 = sprintf('toto is an empty dates object.\\n');
-%$     try
-%$         t(2) = dassert(str, expected_str_1);
-%$     catch
-%$         t(2) = false;
-%$     end
-%$ end
-%$
-%$ T = all(t);
-%$ cd(OPATH);
+OPATH = pwd();
+DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
+cd([DSERIES_PATH '/@dates/private']);
+
+try
+    toto = dates();
+    str = dates4display(toto, 'toto', 5);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    expected_str_1 = sprintf('toto is an empty dates object.\\n');
+    try
+        t(2) = dassert(str, expected_str_1);
+    catch
+        t(2) = false;
+    end
+end
+
+T = all(t);
+cd(OPATH);
 %@eof:1
 
 %@test:2
-%$ OPATH = pwd();
-%$ DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
-%$ cd([DSERIES_PATH '/@dates/private']);
-%$
-%$ try
-%$     toto = dates('1950Q1'):dates('1950Q2');
-%$     str = dates4display(toto, 'toto', 5);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     expected_str_2 = sprintf('toto = <dates: 1950Q1, 1950Q2>\\n');
-%$     try
-%$         t(2) = dassert(str, expected_str_2);
-%$     catch
-%$         t(2) = false;
-%$     end
-%$ end
-%$
-%$ T = all(t);
-%$ cd(OPATH);
+OPATH = pwd();
+DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
+cd([DSERIES_PATH '/@dates/private']);
+
+try
+    toto = dates('1950Q1'):dates('1950Q2');
+    str = dates4display(toto, 'toto', 5);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    expected_str_2 = sprintf('toto = <dates: 1950Q1, 1950Q2>\\n');
+    try
+        t(2) = dassert(str, expected_str_2);
+    catch
+        t(2) = false;
+    end
+end
+
+T = all(t);
+cd(OPATH);
 %@eof:2
 
 %@test:3
-%$ OPATH = pwd();
-%$ DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
-%$ cd([DSERIES_PATH '/@dates/private']);
-%$
-%$ try
-%$     toto = dates('1950Q1'):dates('1951Q1');
-%$     str = dates4display(toto, 'toto', 4);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     expected_str_3 = sprintf('toto = <dates: 1950Q1, 1950Q2, ..., 1950Q4, 1951Q1>\\n');
-%$     try
-%$         t(2) = dassert(str, expected_str_3);
-%$     catch
-%$         t(2) = false;
-%$     end
-%$ end
-%$
-%$ T = all(t);
-%$ cd(OPATH);
+OPATH = pwd();
+DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
+cd([DSERIES_PATH '/@dates/private']);
+
+try
+    toto = dates('1950Q1'):dates('1951Q1');
+    str = dates4display(toto, 'toto', 4);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    expected_str_3 = sprintf('toto = <dates: 1950Q1, 1950Q2, ..., 1950Q4, 1951Q1>\\n');
+    try
+        t(2) = dassert(str, expected_str_3);
+    catch
+        t(2) = false;
+    end
+end
+
+T = all(t);
+cd(OPATH);
 %@eof:3
 
 %@test:4
-%$ OPATH = pwd();
-%$ DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
-%$ cd([DSERIES_PATH '/@dates/private']);
-%$
-%$ try
-%$     toto = dates('1950Q1'):dates('1951Q1');
-%$     str = dates4display(toto, 'toto', 6);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     expected_str_4 = sprintf('toto = <dates: 1950Q1, 1950Q2, 1950Q3, 1950Q4, 1951Q1>\\n');
-%$     try
-%$         t(2) = dassert(str, expected_str_4);
-%$     catch
-%$         t(2) = false;
-%$     end
-%$ end
-%$
-%$ T = all(t);
-%$ cd(OPATH);
+OPATH = pwd();
+DSERIES_PATH = strrep(which('initialize_dseries_class'),'/initialize_dseries_class.m','');
+cd([DSERIES_PATH '/@dates/private']);
+
+try
+    toto = dates('1950Q1'):dates('1951Q1');
+    str = dates4display(toto, 'toto', 6);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    expected_str_4 = sprintf('toto = <dates: 1950Q1, 1950Q2, 1950Q3, 1950Q4, 1951Q1>\\n');
+    try
+        t(2) = dassert(str, expected_str_4);
+    catch
+        t(2) = false;
+    end
+end
+
+T = all(t);
+cd(OPATH);
 %@eof:4
diff --git a/src/@dates/remove.m b/src/@dates/remove.m
index cb6aadb..dc67fbe 100644
--- a/src/@dates/remove.m
+++ b/src/@dates/remove.m
@@ -1,4 +1,4 @@
-function o = remove(o, p) % --*-- Unitary tests --*--
+function o = remove(o, p)
 
 % remove method for dates class (removes dates).
 %
@@ -15,7 +15,7 @@ function o = remove(o, p) % --*-- Unitary tests --*--
 %
 % See also pop, setdiff
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -45,7 +45,7 @@ end
 o = copy(o);
 o.remove_(p);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/remove_.m b/src/@dates/remove_.m
index 8e4c530..c938325 100644
--- a/src/@dates/remove_.m
+++ b/src/@dates/remove_.m
@@ -1,4 +1,4 @@
-function o = remove_(o, p) % --*-- Unitary tests --*--
+function o = remove_(o, p)
 
 % remove method for dates class (removes dates).
 %
@@ -15,7 +15,7 @@ function o = remove_(o, p) % --*-- Unitary tests --*--
 %
 % See also pop, setdiff
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -46,7 +46,7 @@ time = setdiff(o.time, p.time, 'legacy');
 
 o.time = time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/setdiff.m b/src/@dates/setdiff.m
index 2856ede..4f0462e 100644
--- a/src/@dates/setdiff.m
+++ b/src/@dates/setdiff.m
@@ -1,4 +1,4 @@
-function [q, io] = setdiff(o, p) % --*-- Unitary tests --*--
+function [q, io] = setdiff(o, p)
 
 % Overloads setdiff function for dates objects.
 %
@@ -12,7 +12,7 @@ function [q, io] = setdiff(o, p) % --*-- Unitary tests --*--
 %
 % See also pop, remove.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -57,7 +57,7 @@ end
 q = dates(o.freq);
 q.time = time;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
diff --git a/src/@dates/sort.m b/src/@dates/sort.m
index 1b8f231..0b15364 100644
--- a/src/@dates/sort.m
+++ b/src/@dates/sort.m
@@ -1,4 +1,4 @@
-function o = sort(o) % --*-- Unitary tests --*--
+function o = sort(o)
 
 % Sort method for dates class (with copy).
 %
@@ -8,7 +8,7 @@ function o = sort(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates] with dates sorted by increasing order.
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2022 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
@@ -26,7 +26,7 @@ function o = sort(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.sort_();
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/sort_.m b/src/@dates/sort_.m
index 25e49cf..c66be81 100644
--- a/src/@dates/sort_.m
+++ b/src/@dates/sort_.m
@@ -1,4 +1,4 @@
-function o = sort_(o) % --*-- Unitary tests --*--
+function o = sort_(o)
 
 % Sort method for dates class (in place modification).
 %
@@ -8,7 +8,7 @@ function o = sort_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates] with dates sorted by increasing order.
 
-% Copyright © 2015-2021 Dynare Team
+% Copyright © 2015-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
@@ -29,7 +29,7 @@ end
 
 o.time = sort(o.time);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -101,4 +101,4 @@ if t(1)
     t(2) = isequal(d,dates('1999-01-02', '2000-01-01', '2001-01-01'));
 end
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dates/strings.m b/src/@dates/strings.m
index 6d3a2d7..b585819 100644
--- a/src/@dates/strings.m
+++ b/src/@dates/strings.m
@@ -1,4 +1,4 @@
-function m = strings(o) % --*-- Unitary tests --*--
+function m = strings(o)
 
 % Returns a cell array of strings containing the dates
 %
@@ -8,7 +8,7 @@ function m = strings(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - m [cell of char] object with n elements.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -29,7 +29,7 @@ for i = 1:o.length()
     m(i) = { date2string(o.time(i), o.freq) };
 end
 
-return
+return  % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates objects
diff --git a/src/@dates/subsasgn.m b/src/@dates/subsasgn.m
index e3b1953..e2f4ac1 100644
--- a/src/@dates/subsasgn.m
+++ b/src/@dates/subsasgn.m
@@ -1,6 +1,6 @@
-function val = subsasgn(val, idx, rhs) % --*-- Unitary tests --*--
+function val = subsasgn(val, idx, rhs)
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -17,7 +17,7 @@ function val = subsasgn(val, idx, rhs) % --*-- Unitary tests --*--
 
 error('dates::subsasgn: Members of dates class are private')
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates objects
@@ -32,4 +32,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dates/subsref.m b/src/@dates/subsref.m
index 9b67144..eb96695 100644
--- a/src/@dates/subsref.m
+++ b/src/@dates/subsref.m
@@ -1,4 +1,4 @@
-function B = subsref(A,S) % --*-- Unitary tests --*--
+function B = subsref(A,S)
 
 % Overloads the subsref method for dates objects.
 %
@@ -13,7 +13,7 @@ function B = subsref(A,S) % --*-- Unitary tests --*--
 % 1. The type of the returned argument depends on the content of S.
 % 2. See the matlab's documentation about the subsref method.
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-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
@@ -256,7 +256,7 @@ if ~isempty(S)
     B = subsref(B, S);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -465,4 +465,4 @@ if t(1)
     t(4) = isequal(d.ndat(),1);
 end
 T = all(t);
-%@eof:9
\ No newline at end of file
+%@eof:9
diff --git a/src/@dates/uminus.m b/src/@dates/uminus.m
index 6d20b3b..814d837 100644
--- a/src/@dates/uminus.m
+++ b/src/@dates/uminus.m
@@ -1,4 +1,4 @@
-function o = uminus(o)  % --*-- Unitary tests --*--
+function o = uminus(o)
 
 % Overloads the unary minus operator for dates objects. Shifts all the elements by one period.
 %
@@ -8,7 +8,7 @@ function o = uminus(o)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -25,7 +25,7 @@ function o = uminus(o)  % --*-- Unitary tests --*--
 
 o.time(:) = o.time(:)-1;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -89,4 +89,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dates/union.m b/src/@dates/union.m
index 523f989..93c162d 100644
--- a/src/@dates/union.m
+++ b/src/@dates/union.m
@@ -1,4 +1,4 @@
-function o = union(varargin) % --*-- Unitary tests --*--
+function o = union(varargin)
 
 % Overloads union function for dates objects (removes repetitions if any).
 %
@@ -11,7 +11,7 @@ function o = union(varargin) % --*-- Unitary tests --*--
 % REMARKS
 % 1. Elements in o are sorted by increasing order.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -32,12 +32,12 @@ end
 
 if isequal(nargin, 1)
     o = sort(unique(varargin{1}));
-    return;
+    return
 end
 
 o = sort(unique(horzcat(varargin{:})));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates objects
@@ -75,4 +75,4 @@ e2 = union(d1,d2);
 t(1) = isequal(e1, d1);
 t(2) = isequal(e2, d3);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dates/unique.m b/src/@dates/unique.m
index 338865c..ecdbc79 100644
--- a/src/@dates/unique.m
+++ b/src/@dates/unique.m
@@ -1,4 +1,4 @@
-function o = unique(o) % --*-- Unitary tests --*--
+function o = unique(o)
 
 % Overloads the unique function for dates objects.
 %
@@ -11,7 +11,7 @@ function o = unique(o) % --*-- Unitary tests --*--
 % REMARKS
 % 1. Only the last occurence of a date is kept.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -33,7 +33,7 @@ end
 o = copy(o);
 o.unique_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -141,4 +141,4 @@ if t(1)
     t(2) = isequal(d2, dates('2000-01-02', '2000-01-01'));
 end
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dates/unique_.m b/src/@dates/unique_.m
index 60d031a..7d59826 100644
--- a/src/@dates/unique_.m
+++ b/src/@dates/unique_.m
@@ -1,4 +1,4 @@
-function o = unique_(o) % --*-- Unitary tests --*--
+function o = unique_(o)
 
 % Overloads the unique function for dates objects (in place modification).
 %
@@ -11,7 +11,7 @@ function o = unique_(o) % --*-- Unitary tests --*--
 % REMARKS
 % 1. Only the last occurence of a date is kept.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -34,7 +34,7 @@ end
 
 o.time = o.time(sort(id));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/uplus.m b/src/@dates/uplus.m
index 15b3a14..d04569f 100644
--- a/src/@dates/uplus.m
+++ b/src/@dates/uplus.m
@@ -1,4 +1,4 @@
-function o = uplus(o)  % --*-- Unitary tests --*--
+function o = uplus(o)
 
 % Overloads the unary plus operator for dates objects. Shifts all the elements by one period.
 %
@@ -8,7 +8,7 @@ function o = uplus(o)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dates]
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -25,7 +25,7 @@ function o = uplus(o)  % --*-- Unitary tests --*--
 
 o.time = o.time+1;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -90,4 +90,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dates/vertcat.m b/src/@dates/vertcat.m
index aa381de..761fd55 100644
--- a/src/@dates/vertcat.m
+++ b/src/@dates/vertcat.m
@@ -1,4 +1,4 @@
-function o = vertcat(varargin) % --*-- Unitary tests --*--
+function o = vertcat(varargin)
 
 % Overloads the vertcat method for dates objects.
 %
@@ -15,7 +15,7 @@ function o = vertcat(varargin) % --*-- Unitary tests --*--
 %
 %  Defines a dates object D containing the dates appearing in A, B and C.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2022 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
@@ -36,7 +36,7 @@ end
 
 o = horzcat(varargin{:});
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
diff --git a/src/@dates/year.m b/src/@dates/year.m
index 657033e..77bce3f 100644
--- a/src/@dates/year.m
+++ b/src/@dates/year.m
@@ -1,8 +1,8 @@
-function y = year(d)  % --*-- Unitary tests --*--
+function y = year(d)
 
 % Returns the year.
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2022 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
@@ -29,7 +29,7 @@ else
     y = d.time;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/abs.m b/src/@dseries/abs.m
index e9cc534..29f8bd9 100644
--- a/src/@dseries/abs.m
+++ b/src/@dseries/abs.m
@@ -1,4 +1,4 @@
-function o = abs(o) % --*-- Unitary tests --*--
+function o = abs(o)
 
 % Apply the absolute value to all the variables in a dseries object (without in place modification).
 %
@@ -8,7 +8,7 @@ function o = abs(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = abs(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.abs_;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
diff --git a/src/@dseries/abs_.m b/src/@dseries/abs_.m
index e0ab59c..592c1fa 100644
--- a/src/@dseries/abs_.m
+++ b/src/@dseries/abs_.m
@@ -1,4 +1,4 @@
-function o = abs_(o) % --*-- Unitary tests --*--
+function o = abs_(o)
 
 % Apply the absolute value to all the variables in a dseries object (in place modification).
 %
@@ -8,7 +8,7 @@ function o = abs_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ end
 
 o.data = abs(o.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -99,4 +99,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/align.m b/src/@dseries/align.m
index ec4ca1c..3b17c2c 100644
--- a/src/@dseries/align.m
+++ b/src/@dseries/align.m
@@ -1,4 +1,4 @@
-function [o, p] = align(o, p) % --*-- Unitary tests --*--
+function [o, p] = align(o, p)
 
 % If necessay completes dseries object o and p so that they are defined on the same time range
 % (in place modification).
@@ -11,7 +11,7 @@ function [o, p] = align(o, p) % --*-- Unitary tests --*--
 % - o [dseries]
 % - p [dseries]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -32,7 +32,7 @@ o = copy(o);
 p = copy(p);
 align_(o, p);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -128,4 +128,4 @@ if t(1)
   t(5) = dassert(ts2.data,[B; NaN(1,2)], 1e-15);
 end
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/align_.m b/src/@dseries/align_.m
index 3fea71d..7e19c43 100644
--- a/src/@dseries/align_.m
+++ b/src/@dseries/align_.m
@@ -1,4 +1,4 @@
-function [o, p] = align_(o, p) % --*-- Unitary tests --*--
+function [o, p] = align_(o, p)
 
 % If necessay completes dseries object o and p so that they are defined on the same time range
 % (in place modification).
@@ -11,7 +11,7 @@ function [o, p] = align_(o, p) % --*-- Unitary tests --*--
 % - o [dseries]
 % - p [dseries]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -62,7 +62,7 @@ end
 o.dates = init:init+(nobs(o)-1);
 p.dates = init:init+(nobs(p)-1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -161,4 +161,4 @@ if t(1)
   t(6) = dassert(ts1.dates, ts2.dates);
 end
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/backcast.m b/src/@dseries/backcast.m
index 62ca613..99d588b 100644
--- a/src/@dseries/backcast.m
+++ b/src/@dseries/backcast.m
@@ -1,4 +1,4 @@
-function q = backcast(o, p, diff)  % --*-- Unitary tests --*--
+function q = backcast(o, p, diff)
 
 % Chains two dseries objects.
 %
@@ -13,7 +13,7 @@ function q = backcast(o, p, diff)  % --*-- Unitary tests --*--
 % The two dseries objects must have common frequency and the same number of variables. Also the
 % two samples must overlap.
 
-% Copyright © 2019-2020 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -38,7 +38,7 @@ else
     q.backcast_(p);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
   y = ones(10,1);
@@ -66,4 +66,4 @@ return
   end
 
   T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/backcast_.m b/src/@dseries/backcast_.m
index 10be86a..352843f 100644
--- a/src/@dseries/backcast_.m
+++ b/src/@dseries/backcast_.m
@@ -1,6 +1,6 @@
-function o = backcast_(o, p, diff)  % --*-- Unitary tests --*--
+function o = backcast_(o, p, diff)
 
-% Copyright © 2019-2020 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -65,7 +65,7 @@ for i=1:o.vobs
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
   y = ones(10,1);
@@ -121,4 +121,4 @@ return
   end
 
   T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/baxter_king_filter.m b/src/@dseries/baxter_king_filter.m
index fb01076..610585f 100644
--- a/src/@dseries/baxter_king_filter.m
+++ b/src/@dseries/baxter_king_filter.m
@@ -1,4 +1,4 @@
-function o = baxter_king_filter(o, high_frequency, low_frequency, K) % --*-- Unitary tests --*--
+function o = baxter_king_filter(o, high_frequency, low_frequency, K)
 
 % Implementation of Baxter and King (1999) band pass filter for dseries objects. The code is adapted from
 % the one provided by Baxter and King. This filter isolates business cycle fluctuations with a period of length
@@ -17,7 +17,7 @@ function o = baxter_king_filter(o, high_frequency, low_frequency, K) % --*-- Uni
 % This filter use a (symmetric) moving average smoother, so that K observations at the beginning and at the end of the
 % sample are lost in the computation of the filter.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -59,7 +59,7 @@ end
 o = copy(o);
 o.baxter_king_filter_(high_frequency, low_frequency, K);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 plot_flag = false;
diff --git a/src/@dseries/baxter_king_filter_.m b/src/@dseries/baxter_king_filter_.m
index b4f52b4..0684224 100644
--- a/src/@dseries/baxter_king_filter_.m
+++ b/src/@dseries/baxter_king_filter_.m
@@ -1,4 +1,4 @@
-function o = baxter_king_filter_(o, high_frequency, low_frequency, K) % --*-- Unitary tests --*--
+function o = baxter_king_filter_(o, high_frequency, low_frequency, K)
 
 % Implementation of Baxter and King (1999) band pass filter for dseries objects. The code is adapted from
 % the one provided by Baxter and King. This filter isolates business cycle fluctuations with a period of length
@@ -17,7 +17,7 @@ function o = baxter_king_filter_(o, high_frequency, low_frequency, K) % --*-- Un
 % This filter use a (symmetric) moving average smoother, so that K observations at the beginning and at the end of the
 % sample are lost in the computation of the filter.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -109,7 +109,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 plot_flag = false;
diff --git a/src/@dseries/center.m b/src/@dseries/center.m
index 317a709..f796e20 100644
--- a/src/@dseries/center.m
+++ b/src/@dseries/center.m
@@ -1,4 +1,4 @@
-function o = center(o, geometric) % --*-- Unitary tests --*--
+function o = center(o, geometric)
 
 % Centers dseries object o around its mean (arithmetic or geometric).
 %
@@ -9,7 +9,7 @@ function o = center(o, geometric) % --*-- Unitary tests --*--
 % OUTPUTS
 %  - o             dseries object.
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -33,7 +33,7 @@ end
 o = copy(o);
 o.center_(geometric);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
diff --git a/src/@dseries/center_.m b/src/@dseries/center_.m
index cc4f130..c8d1602 100644
--- a/src/@dseries/center_.m
+++ b/src/@dseries/center_.m
@@ -1,4 +1,4 @@
-function o = center_(o, geometric) % --*-- Unitary tests --*--
+function o = center_(o, geometric)
 
 % Centers dseries object o around its mean (arithmetic or geometric).
 %
@@ -9,7 +9,7 @@ function o = center_(o, geometric) % --*-- Unitary tests --*--
 % OUTPUTS
 %  - o             dseries object.
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -46,7 +46,7 @@ for i=1:o.vobs
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
diff --git a/src/@dseries/chain.m b/src/@dseries/chain.m
index e07009e..7a42c06 100644
--- a/src/@dseries/chain.m
+++ b/src/@dseries/chain.m
@@ -1,4 +1,4 @@
-function q = chain(o, p)  % --*-- Unitary tests --*--
+function q = chain(o, p)
 
 % Chains two dseries objects.
 %
@@ -13,7 +13,7 @@ function q = chain(o, p)  % --*-- Unitary tests --*--
 % The two dseries objects must have common frequency and the same number of variables. Also the
 % two samples must overlap.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -33,7 +33,7 @@ function q = chain(o, p)  % --*-- Unitary tests --*--
 q = copy(o);
 q.chain_(p);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -57,4 +57,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/chain_.m b/src/@dseries/chain_.m
index d34c18b..b6cec5b 100644
--- a/src/@dseries/chain_.m
+++ b/src/@dseries/chain_.m
@@ -1,6 +1,6 @@
-function o = chain_(o, p)  % --*-- Unitary tests --*--
+function o = chain_(o, p)
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -77,7 +77,7 @@ for i=1:o.vobs
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -124,4 +124,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/copy.m b/src/@dseries/copy.m
index 92e3127..f8488a5 100644
--- a/src/@dseries/copy.m
+++ b/src/@dseries/copy.m
@@ -1,4 +1,4 @@
-function p = copy(o) % --*-- Unitary tests --*--
+function p = copy(o)
 
 % Do a copy of a dseries object.
 %
@@ -8,7 +8,7 @@ function p = copy(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - p [dates]
 
-% Copyright © 2015-2021 Dynare Team
+% Copyright © 2015-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
@@ -31,7 +31,7 @@ p.dates = o.dates;
 p.ops   = o.ops;
 p.tags  = o.tags;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -53,4 +53,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/cumprod.m b/src/@dseries/cumprod.m
index 9365054..a1f9219 100644
--- a/src/@dseries/cumprod.m
+++ b/src/@dseries/cumprod.m
@@ -1,4 +1,4 @@
-function o = cumprod(varargin) % --*-- Unitary tests --*--
+function o = cumprod(varargin)
 
 % Overloads matlab's cumprod function for dseries objects.
 %
@@ -10,7 +10,7 @@ function o = cumprod(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - B     dseries object.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -39,7 +39,7 @@ else
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
diff --git a/src/@dseries/cumprod_.m b/src/@dseries/cumprod_.m
index c2ce917..e2e5ae3 100644
--- a/src/@dseries/cumprod_.m
+++ b/src/@dseries/cumprod_.m
@@ -1,4 +1,4 @@
-function o = cumprod_(varargin) % --*-- Unitary tests --*--
+function o = cumprod_(varargin)
 
 % Overloads matlab's cumprod function for dseries objects.
 %
@@ -10,7 +10,7 @@ function o = cumprod_(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o     dseries object.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -128,7 +128,7 @@ switch nargin
     error('dseries::cumprod: Wrong number of input arguments!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
@@ -291,4 +291,4 @@ if t(1)
   t(2) = dassert(ts0.data, A);
 end
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/@dseries/cumsum.m b/src/@dseries/cumsum.m
index 14f0691..d1575af 100644
--- a/src/@dseries/cumsum.m
+++ b/src/@dseries/cumsum.m
@@ -1,4 +1,4 @@
-function o = cumsum(varargin) % --*-- Unitary tests --*--
+function o = cumsum(varargin)
 
 % Overloads matlab's cumsum function for dseries objects.
 %
@@ -10,7 +10,7 @@ function o = cumsum(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o     dseries object.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -39,7 +39,7 @@ else
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
@@ -98,4 +98,4 @@ if t(1)
   t(3) = dassert(ts1.data, A);
 end
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/cumsum_.m b/src/@dseries/cumsum_.m
index 2e844e7..7782dd0 100644
--- a/src/@dseries/cumsum_.m
+++ b/src/@dseries/cumsum_.m
@@ -1,4 +1,4 @@
-function o = cumsum_(varargin) % --*-- Unitary tests --*--
+function o = cumsum_(varargin)
 
 % Overloads matlab's cumsum function for dseries objects.
 %
@@ -10,7 +10,7 @@ function o = cumsum_(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o     dseries object.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -126,7 +126,7 @@ switch nargin
     error('dseries::cumsum: Wrong number of input arguments!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
@@ -293,4 +293,4 @@ if t(1)
   t(2) = dassert(ts0.data, A);
 end
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/@dseries/detrend.m b/src/@dseries/detrend.m
index 850abcb..ea470c5 100644
--- a/src/@dseries/detrend.m
+++ b/src/@dseries/detrend.m
@@ -1,4 +1,4 @@
-function o = detrend(o, model) % --*-- Unitary tests --*--
+function o = detrend(o, model)
 
 % Detrends a dseries object with a polynomial of order model.
 %
@@ -9,7 +9,7 @@ function o = detrend(o, model) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o       [dseries]   detrended time series.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -34,7 +34,7 @@ end
 o = copy(o);
 o.detrend_(model);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -53,4 +53,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/detrend_.m b/src/@dseries/detrend_.m
index 9b9862b..a71659b 100644
--- a/src/@dseries/detrend_.m
+++ b/src/@dseries/detrend_.m
@@ -1,4 +1,4 @@
-function o = detrend_(o, model) % --*-- Unitary tests --*--
+function o = detrend_(o, model)
 
 % Detrends a dseries object with a polynomial of order model.
 %
@@ -9,7 +9,7 @@ function o = detrend_(o, model) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o       [dseries]   detrended time series.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -60,7 +60,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -146,4 +146,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/dgrowth.m b/src/@dseries/dgrowth.m
index 38f27bd..5d92e1d 100644
--- a/src/@dseries/dgrowth.m
+++ b/src/@dseries/dgrowth.m
@@ -1,4 +1,4 @@
-function o = dgrowth(o) % --*-- Unitary tests --*--
+function o = dgrowth(o)
 
 % Computes daily growth rates.
 %
@@ -8,7 +8,7 @@ function o = dgrowth(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = dgrowth(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.dgrowth_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/dgrowth_.m b/src/@dseries/dgrowth_.m
index 1c22890..2fa6471 100644
--- a/src/@dseries/dgrowth_.m
+++ b/src/@dseries/dgrowth_.m
@@ -1,4 +1,4 @@
-function o = dgrowth_(o) % --*-- Unitary tests --*--
+function o = dgrowth_(o)
 
 % Computes daily growth rates.
 %
@@ -8,7 +8,7 @@ function o = dgrowth_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ switch frequency(o)
   case 12
     error('dseries::dgrowth: I cannot compute daily growth rates from monthly data!')
   case 52
-    error('dseries::dgrowth: I cannot compute daily growth rates from weekly data!')    
+    error('dseries::dgrowth: I cannot compute daily growth rates from weekly data!')
   case 365
     o.data(2:end,:) = o.data(2:end,:)./o.data(1:end-1,:) - 1;
     o.data(1,:) = NaN;
@@ -51,7 +51,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/diff.m b/src/@dseries/diff.m
index 3a44274..781615e 100644
--- a/src/@dseries/diff.m
+++ b/src/@dseries/diff.m
@@ -1,4 +1,4 @@
-function o = diff(o) % --*-- Unitary tests --*--
+function o = diff(o)
 
 % Computes differences.
 %
@@ -8,7 +8,7 @@ function o = diff(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = diff(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.diff_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/diff_.m b/src/@dseries/diff_.m
index 5acec57..8609a99 100644
--- a/src/@dseries/diff_.m
+++ b/src/@dseries/diff_.m
@@ -1,4 +1,4 @@
-function o = diff_(o) % --*-- Unitary tests --*--
+function o = diff_(o)
 
 % Computes differences, independently of the frequency.
 %
@@ -8,7 +8,7 @@ function o = diff_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -36,7 +36,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -55,4 +55,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/dseries.m b/src/@dseries/dseries.m
index bb58198..738566d 100644
--- a/src/@dseries/dseries.m
+++ b/src/@dseries/dseries.m
@@ -1,8 +1,8 @@
-classdef dseries<handle % --*-- Unitary tests --*--
+classdef dseries<handle
 
 % Class for time series.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -178,7 +178,7 @@ methods
         end
         end % dseries
 end % methods
-end % classdef
+end % classdef --*-- Unit tests --*--
 
 %@test:1
 %$ % Test if we can instantiate an empty dseries object.
diff --git a/src/@dseries/eq.m b/src/@dseries/eq.m
index 6ddc1da..d32e485 100644
--- a/src/@dseries/eq.m
+++ b/src/@dseries/eq.m
@@ -1,4 +1,4 @@
-function b = eq(o, p) % --*-- Unitary tests --*--
+function b = eq(o, p)
 
 % Overloads eq (==) operator.
 %
@@ -12,7 +12,7 @@ function b = eq(o, p) % --*-- Unitary tests --*--
 % REMARKS
 %  If the number of variables, the number of observations or the frequencies are different in A and B, the function returns a zero scalar.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -63,7 +63,7 @@ end
 
 b = eq(o.data, p.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -88,4 +88,4 @@ if length(t)>1
    t(2) = dassert(a,logical(ones(10,3)));
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/exist.m b/src/@dseries/exist.m
index c01dbc2..0d764ca 100644
--- a/src/@dseries/exist.m
+++ b/src/@dseries/exist.m
@@ -1,4 +1,4 @@
-function l = exist(o, varname) % --*-- Unitary tests --*--
+function l = exist(o, varname)
 
 % Tests if a variable exists in dseries object o.
 %
@@ -9,7 +9,7 @@ function l = exist(o, varname) % --*-- Unitary tests --*--
 % OUTPUTS
 %  - l       [logical], equal to 1 (true) iff varname is a variable in dseries object o.
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -36,7 +36,7 @@ end
 
 l = ~isempty(strmatch(varname, o.name, 'exact'));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/exp.m b/src/@dseries/exp.m
index ee3d1ad..0fb3487 100644
--- a/src/@dseries/exp.m
+++ b/src/@dseries/exp.m
@@ -1,4 +1,4 @@
-function o = exp(o) % --*-- Unitary tests --*--
+function o = exp(o)
 
 % Apply the exponential to all the variables in a dseries object (without in place modification).
 %
@@ -8,7 +8,7 @@ function o = exp(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = exp(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.exp_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -79,4 +79,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/exp_.m b/src/@dseries/exp_.m
index 83fcd2b..3e97d42 100644
--- a/src/@dseries/exp_.m
+++ b/src/@dseries/exp_.m
@@ -1,4 +1,4 @@
-function o = exp_(o) % --*-- Unitary tests --*--
+function o = exp_(o)
 
 % Apply the exponential to all the variables in a dseries object (in place modification).
 %
@@ -8,7 +8,7 @@ function o = exp_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2015-2021 Dynare Team
+% Copyright © 2015-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -93,4 +93,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/extract.m b/src/@dseries/extract.m
index ca33710..6202c6e 100644
--- a/src/@dseries/extract.m
+++ b/src/@dseries/extract.m
@@ -1,8 +1,8 @@
-function p = extract(o, varargin) % --*-- Unitary tests --*--
+function p = extract(o, varargin)
 
 % Extract some variables from a database.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -99,7 +99,7 @@ p.name = o.name(idVariableName);
 p.tex = o.tex(idVariableName);
 p.ops = o.ops(idVariableName);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
@@ -189,4 +189,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/firstdate.m b/src/@dseries/firstdate.m
index ab83478..3b119f9 100644
--- a/src/@dseries/firstdate.m
+++ b/src/@dseries/firstdate.m
@@ -1,6 +1,6 @@
-function f = firstdate(o) % --*-- Unitary tests --*--
+function f = firstdate(o)
 
-% Copyright © 2014-2021 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -19,7 +19,7 @@ function f = firstdate(o) % --*-- Unitary tests --*--
 
 f = o.dates(1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/firstobservedperiod.m b/src/@dseries/firstobservedperiod.m
index 3b89055..b665075 100644
--- a/src/@dseries/firstobservedperiod.m
+++ b/src/@dseries/firstobservedperiod.m
@@ -1,4 +1,4 @@
-function d = firstobservedperiod(o) % --*-- Unitary tests --*--
+function d = firstobservedperiod(o)
 
 % Returns the first period where all the variables are observed (first period without NaNs).
 %
@@ -8,7 +8,7 @@ function d = firstobservedperiod(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - d [dates]      First period where the N variables are observed (without NaNs).
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -34,7 +34,7 @@ end
 
 d = firstdate(o)+(c(1)-1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/flip.m b/src/@dseries/flip.m
index 104cc0d..22505cc 100644
--- a/src/@dseries/flip.m
+++ b/src/@dseries/flip.m
@@ -1,4 +1,4 @@
-function o = flip(o) % --*-- Unitary tests --*--
+function o = flip(o)
 
 % Flips the rows in the data member (without changing the
 % periods order)
@@ -9,7 +9,7 @@ function o = flip(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -29,7 +29,8 @@ function o = flip(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.flip_;
 
-return
+return % --*-- Unit tests --*--
+
 %@test:1
 % Define a dates object
 data = transpose(1:3);
@@ -50,4 +51,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/flip_.m b/src/@dseries/flip_.m
index b2a1054..5a86666 100644
--- a/src/@dseries/flip_.m
+++ b/src/@dseries/flip_.m
@@ -1,4 +1,4 @@
-function o = flip_(o) % --*-- Unitary tests --*--
+function o = flip_(o)
 
 % Flips the rows in the data member (without changing the
 % periods order)
@@ -9,7 +9,7 @@ function o = flip_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -36,7 +36,8 @@ end
 
 o.data = flipud(o.data);
 
-return
+return % --*-- Unit tests --*--
+
 %@test:1
 % Define a dates object
 data = [transpose(1:3) transpose(4:6)];
@@ -67,35 +68,35 @@ T = all(t);
 %@eof:1
 
 %@test:2
-%$ % Define a dates object
-%$ data = [ones(10,1), -ones(10,1)];
-%$ o = dseries(data);
-%$ q = o;
-%$ r = copy(o);
-%$
-%$ % Call the tested routine.
-%$ try
-%$     o.abs_();
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$      t(2) = dassert(length(o.name), 2);
-%$      t(3) = dassert(o.name{1},'Variable_1');
-%$      t(4) = dassert(o.name{2},'Variable_2');
-%$      t(5) = dassert(q.name{1},'Variable_1');
-%$      t(6) = dassert(q.name{2},'Variable_2');
-%$      t(7) = dassert(r.name{1},'Variable_1');
-%$      t(8) = dassert(r.name{2},'Variable_2');
-%$      t(9) = dassert(o.ops{1},'abs(Variable_1)');
-%$      t(10) = dassert(o.ops{2},'abs(Variable_2)');
-%$      t(11) = dassert(q.ops{1},'abs(Variable_1)');
-%$      t(12) = dassert(q.ops{2},'abs(Variable_2)');
-%$      t(13) = isempty(r.ops{1});
-%$      t(14) = isempty(r.ops{2});
-%$ end
-%$
-%$ T = all(t);
-%@eof:2
\ No newline at end of file
+% Define a dates object
+data = [ones(10,1), -ones(10,1)];
+o = dseries(data);
+q = o;
+r = copy(o);
+
+% Call the tested routine.
+try
+    o.abs_();
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+     t(2) = dassert(length(o.name), 2);
+     t(3) = dassert(o.name{1},'Variable_1');
+     t(4) = dassert(o.name{2},'Variable_2');
+     t(5) = dassert(q.name{1},'Variable_1');
+     t(6) = dassert(q.name{2},'Variable_2');
+     t(7) = dassert(r.name{1},'Variable_1');
+     t(8) = dassert(r.name{2},'Variable_2');
+     t(9) = dassert(o.ops{1},'abs(Variable_1)');
+     t(10) = dassert(o.ops{2},'abs(Variable_2)');
+     t(11) = dassert(q.ops{1},'abs(Variable_1)');
+     t(12) = dassert(q.ops{2},'abs(Variable_2)');
+     t(13) = isempty(r.ops{1});
+     t(14) = isempty(r.ops{2});
+end
+
+T = all(t);
+%@eof:2
diff --git a/src/@dseries/ge.m b/src/@dseries/ge.m
index a82a66b..cf2e68b 100644
--- a/src/@dseries/ge.m
+++ b/src/@dseries/ge.m
@@ -1,8 +1,8 @@
-function l = ge(varargin)  % --*-- Unitary tests --*--
+function l = ge(varargin)
 
 % Overloads the >= operator for dseries objects.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -25,4 +25,4 @@ elseif isdseries(varargin{1}) && isdseries(varargin{2})
     l = varargin{1}.data>=varargin{2}.data;
 else
     error('dseries::ge: wrong types.')
-end
\ No newline at end of file
+end
diff --git a/src/@dseries/gt.m b/src/@dseries/gt.m
index 84d067c..4ee4f56 100644
--- a/src/@dseries/gt.m
+++ b/src/@dseries/gt.m
@@ -1,8 +1,8 @@
-function l = gt(varargin)  % --*-- Unitary tests --*--
+function l = gt(varargin)
 
 % Overloads the > operator for dseries objects.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -25,4 +25,4 @@ elseif isdseries(varargin{1}) && isdseries(varargin{2})
     l = varargin{1}.data>varargin{2}.data;
 else
     error('dseries::gt: wrong types.')
-end
\ No newline at end of file
+end
diff --git a/src/@dseries/hdiff.m b/src/@dseries/hdiff.m
index be154b7..7d6240f 100644
--- a/src/@dseries/hdiff.m
+++ b/src/@dseries/hdiff.m
@@ -1,4 +1,4 @@
-function o = hdiff(o) % --*-- Unitary tests --*--
+function o = hdiff(o)
 
 % Computes bi-annual differences.
 %
@@ -8,7 +8,7 @@ function o = hdiff(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = hdiff(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.hdiff_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/hdiff_.m b/src/@dseries/hdiff_.m
index bc5d557..f31c0be 100644
--- a/src/@dseries/hdiff_.m
+++ b/src/@dseries/hdiff_.m
@@ -1,4 +1,4 @@
-function o = hdiff_(o) % --*-- Unitary tests --*--
+function o = hdiff_(o)
 
 % Computes bi-annual differences.
 %
@@ -8,7 +8,7 @@ function o = hdiff_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -53,7 +53,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -141,4 +141,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/@dseries/hgrowth.m b/src/@dseries/hgrowth.m
index 5604b65..6b64ec6 100644
--- a/src/@dseries/hgrowth.m
+++ b/src/@dseries/hgrowth.m
@@ -1,4 +1,4 @@
-function o = hgrowth(o) % --*-- Unitary tests --*--
+function o = hgrowth(o)
 
 % Computes bi-annual growth rates.
 %
@@ -8,7 +8,7 @@ function o = hgrowth(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = hgrowth(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.hgrowth_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(2,1);
@@ -51,4 +51,4 @@ if length(t)>1
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/hgrowth_.m b/src/@dseries/hgrowth_.m
index b1eefed..8413a96 100644
--- a/src/@dseries/hgrowth_.m
+++ b/src/@dseries/hgrowth_.m
@@ -1,4 +1,4 @@
-function o = hgrowth_(o) % --*-- Unitary tests --*--
+function o = hgrowth_(o)
 
 % Computes bi-annual growth rates.
 %
@@ -8,7 +8,7 @@ function o = hgrowth_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -53,7 +53,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/horzcat.m b/src/@dseries/horzcat.m
index 1f75c3e..aae0830 100644
--- a/src/@dseries/horzcat.m
+++ b/src/@dseries/horzcat.m
@@ -1,4 +1,4 @@
-function o = horzcat(varargin) % --*-- Unitary tests --*--
+function o = horzcat(varargin)
 
 % Overloads horzcat method for dseries objects.
 %
@@ -20,7 +20,7 @@ function o = horzcat(varargin) % --*-- Unitary tests --*--
 % REMARKS
 %  o o1, o2, ... must not have common variables.
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -136,7 +136,7 @@ else
     a.dates = sort(unique(hd));
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
diff --git a/src/@dseries/hpcycle.m b/src/@dseries/hpcycle.m
index 1458c02..5945176 100644
--- a/src/@dseries/hpcycle.m
+++ b/src/@dseries/hpcycle.m
@@ -1,4 +1,4 @@
-function o = hpcycle(o, lambda) % --*-- Unitary tests --*--
+function o = hpcycle(o, lambda)
 
 % Extracts the cycle component from a dseries object using Hodrick Prescott filter.
 %
@@ -9,7 +9,7 @@ function o = hpcycle(o, lambda) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o          [dseries]  Cyclical component of the original time series.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -37,7 +37,7 @@ end
 o = copy(o);
 o.hpcycle_(lambda);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 d = randn(100,1);
@@ -55,4 +55,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/hpcycle_.m b/src/@dseries/hpcycle_.m
index c5c3a2b..7957c25 100644
--- a/src/@dseries/hpcycle_.m
+++ b/src/@dseries/hpcycle_.m
@@ -1,4 +1,4 @@
-function o = hpcycle_(o, lambda) % --*-- Unitary tests --*--
+function o = hpcycle_(o, lambda)
 
 % Extracts the cycle component from a dseries object using Hodrick Prescott filter.
 %
@@ -9,7 +9,7 @@ function o = hpcycle_(o, lambda) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o          [dseries]  Cyclical component of the original time series.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -53,7 +53,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 plot_flag = false;
@@ -104,4 +104,4 @@ if plot_flag
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/hptrend.m b/src/@dseries/hptrend.m
index ba35138..9561cfd 100644
--- a/src/@dseries/hptrend.m
+++ b/src/@dseries/hptrend.m
@@ -1,4 +1,4 @@
-function o = hptrend(o, lambda) % --*-- Unitary tests --*--
+function o = hptrend(o, lambda)
 
 % Extracts the trend component from a dseries object using Hodrick Prescott filter.
 %
@@ -9,7 +9,7 @@ function o = hptrend(o, lambda) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Trend component of the original time series.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -37,7 +37,7 @@ end
 o = copy(o);
 o.hptrend_(lambda);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 d = randn(100,1);
@@ -55,4 +55,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/hptrend_.m b/src/@dseries/hptrend_.m
index 8fcf22d..763c0ef 100644
--- a/src/@dseries/hptrend_.m
+++ b/src/@dseries/hptrend_.m
@@ -1,4 +1,4 @@
-function o = hptrend_(o, lambda) % --*-- Unitary tests --*--
+function o = hptrend_(o, lambda)
 
 % Extracts the trend component from a dseries object using Hodrick Prescott filter.
 %
@@ -9,7 +9,7 @@ function o = hptrend_(o, lambda) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Trend component of the original time series.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -52,7 +52,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 plot_flag = 0;
@@ -102,4 +102,4 @@ if plot_flag
     system('convert -density 300 ../doc/dynare.plots/HPTrend.eps ../doc/dynare.plots/HPTrend.jpg');
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/insert.m b/src/@dseries/insert.m
index 6b616de..aba6cec 100644
--- a/src/@dseries/insert.m
+++ b/src/@dseries/insert.m
@@ -1,4 +1,4 @@
-function o = insert(o, p, id) % --*-- Unitary tests --*--
+function o = insert(o, p, id)
 
 % Adds a variable in a dseries object.
 %
@@ -10,7 +10,7 @@ function o = insert(o, p, id) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o    [dseries]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -109,7 +109,7 @@ for i=1:n
     id = id+1;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -187,4 +187,4 @@ if length(t)>1
    t(5) = dassert(ts2.data, B);
 end
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/isinf.m b/src/@dseries/isinf.m
index 92b047a..1142855 100644
--- a/src/@dseries/isinf.m
+++ b/src/@dseries/isinf.m
@@ -1,4 +1,4 @@
-function b = isinf(o) % --*-- Unitary tests --*--
+function b = isinf(o)
 
 % Returns an array of logicals (true/false).
 % Element (t,i) is true iff the i-th variable at period number t is an Inf.
@@ -10,7 +10,7 @@ function b = isinf(o) % --*-- Unitary tests --*--
 % - b [logical]    T*N array of logicals.
 
 
-% Copyright © 2019-2021 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -29,7 +29,7 @@ function b = isinf(o) % --*-- Unitary tests --*--
 
 b = isinf(o.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/isnan.m b/src/@dseries/isnan.m
index f42c9e6..e037749 100644
--- a/src/@dseries/isnan.m
+++ b/src/@dseries/isnan.m
@@ -1,4 +1,4 @@
-function b = isnan(o) % --*-- Unitary tests --*--
+function b = isnan(o)
 
 % Returns an array of logicals (true/false). Element (t,i) is true iff the i-th variable at
 % period number t is a NaN.
@@ -10,7 +10,7 @@ function b = isnan(o) % --*-- Unitary tests --*--
 % - b [logical]    T*N array of logicals.
 
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -29,7 +29,7 @@ function b = isnan(o) % --*-- Unitary tests --*--
 
 b = isnan(o.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/isreal.m b/src/@dseries/isreal.m
index 121ec46..55d9436 100644
--- a/src/@dseries/isreal.m
+++ b/src/@dseries/isreal.m
@@ -1,4 +1,4 @@
-function b = isreal(o) % --*-- Unitary tests --*--
+function b = isreal(o)
 
 % Returns true if no elements of the dseries have imaginary parts.
 % Returns false otherwise
@@ -10,7 +10,7 @@ function b = isreal(o) % --*-- Unitary tests --*--
 % - b [logical]    true iff all elements are real; false otherwise.
 
 
-% Copyright © 2019-2021 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -29,7 +29,7 @@ function b = isreal(o) % --*-- Unitary tests --*--
 
 b = isreal(o.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/lag.m b/src/@dseries/lag.m
index 9f9f9ad..b5fe32a 100644
--- a/src/@dseries/lag.m
+++ b/src/@dseries/lag.m
@@ -1,4 +1,4 @@
-function o = lag(o, p) % --*-- Unitary tests --*--
+function o = lag(o, p)
 
 % Returns a lagged time series
 %
@@ -23,7 +23,7 @@ function o = lag(o, p) % --*-- Unitary tests --*--
 %    4Y | 3
 %    5Y | 4
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -48,7 +48,7 @@ end
 o = copy(o);
 o.lag_(p);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -71,4 +71,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/lag_.m b/src/@dseries/lag_.m
index c283205..5a1bc09 100644
--- a/src/@dseries/lag_.m
+++ b/src/@dseries/lag_.m
@@ -1,4 +1,4 @@
-function o = lag_(o, p) % --*-- Unitary tests --*--
+function o = lag_(o, p)
 
 % Returns a lagged time series
 %
@@ -23,7 +23,7 @@ function o = lag_(o, p) % --*-- Unitary tests --*--
 %    4Y | 3
 %    5Y | 4
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,7 +75,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -111,4 +111,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/lastobservedperiod.m b/src/@dseries/lastobservedperiod.m
index e2e565e..0f96f0b 100644
--- a/src/@dseries/lastobservedperiod.m
+++ b/src/@dseries/lastobservedperiod.m
@@ -1,4 +1,4 @@
-function d = lastobservedperiod(o) % --*-- Unitary tests --*--
+function d = lastobservedperiod(o)
 
 % Returns the last period where all the variables are observed (last period without NaNs).
 %
@@ -8,7 +8,7 @@ function d = lastobservedperiod(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - b [dates]      Last period where the N variables are observed (without NaNs).
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -32,7 +32,7 @@ if isempty(c)
 end
 d = firstdate(o)+(c(end)-1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -64,4 +64,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/le.m b/src/@dseries/le.m
index 27e2281..f314f07 100644
--- a/src/@dseries/le.m
+++ b/src/@dseries/le.m
@@ -1,8 +1,8 @@
-function l = le(varargin)  % --*-- Unitary tests --*--
+function l = le(varargin)
 
 % Overloads the <= operator for dseries objects.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -25,4 +25,4 @@ elseif isdseries(varargin{1}) && isdseries(varargin{2})
     l = varargin{1}.data<=varargin{2}.data;
 else
     error('dseries::le: wrong types.')
-end
\ No newline at end of file
+end
diff --git a/src/@dseries/lead.m b/src/@dseries/lead.m
index 33bf135..ab64af2 100644
--- a/src/@dseries/lead.m
+++ b/src/@dseries/lead.m
@@ -1,4 +1,4 @@
-function o = lead(o, p) % --*-- Unitary tests --*--
+function o = lead(o, p)
 
 % Returns a lagged time series
 %
@@ -23,7 +23,7 @@ function o = lead(o, p) % --*-- Unitary tests --*--
 %    4Y | 5
 %    5Y | NaN
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -48,7 +48,7 @@ end
 o = copy(o);
 o.lead_(p);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -72,4 +72,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/lead_.m b/src/@dseries/lead_.m
index bc3abfe..d463441 100644
--- a/src/@dseries/lead_.m
+++ b/src/@dseries/lead_.m
@@ -1,4 +1,4 @@
-function o = lead_(o, p) % --*-- Unitary tests --*--
+function o = lead_(o, p)
 
 % Returns a leaded time series
 %
@@ -23,7 +23,7 @@ function o = lead_(o, p) % --*-- Unitary tests --*--
 %    4Y | 5
 %    5Y | NaN
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,7 +75,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -128,4 +128,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/length.m b/src/@dseries/length.m
index 21f386a..86d89f3 100644
--- a/src/@dseries/length.m
+++ b/src/@dseries/length.m
@@ -1,8 +1,8 @@
-function length(o) % --*-- Unitary tests --*--
+function length(o)
 
 % Overloads size function.
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -22,7 +22,7 @@ function length(o) % --*-- Unitary tests --*--
 error(['dseries::length: we do not support the length operator on ' ...
        'dseries. Please use ''nobs'' or ''vobs''']);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -35,4 +35,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/lineartrend.m b/src/@dseries/lineartrend.m
index 7db400a..336161b 100644
--- a/src/@dseries/lineartrend.m
+++ b/src/@dseries/lineartrend.m
@@ -1,4 +1,4 @@
-function trend = lineartrend(o) % --*-- Unitary tests --*--
+function trend = lineartrend(o)
 
 % Returns a trend centered on zero.
 %
@@ -6,12 +6,12 @@ function trend = lineartrend(o) % --*-- Unitary tests --*--
 % - o       [dseries]   time series with T observations.
 %
 % OUTPUTS
-% - trend   [double]    T*1 vector, linear trend with mean zero. 
+% - trend   [double]    T*1 vector, linear trend with mean zero.
 %
 % REMARKS
-% the generic element of the returned argument is trend(i) = i - (n+1)/2.    
+% the generic element of the returned argument is trend(i) = i - (n+1)/2.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -30,7 +30,7 @@ function trend = lineartrend(o) % --*-- Unitary tests --*--
 
 trend = transpose(1:o.nobs)-.5*(o.nobs+1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 a = dseries(randn(100,3));
@@ -45,4 +45,4 @@ if t(1)
     t(2) = abs(mean(trend))<1e-12;
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/log.m b/src/@dseries/log.m
index d6169c8..c9ef2f9 100644
--- a/src/@dseries/log.m
+++ b/src/@dseries/log.m
@@ -1,4 +1,4 @@
-function o = log(o) % --*-- Unitary tests --*--
+function o = log(o)
 
 % Apply the logarithm to all the variables in a dseries object (without in place modification).
 %
@@ -8,7 +8,7 @@ function o = log(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = log(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.log_;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -79,4 +79,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/log_.m b/src/@dseries/log_.m
index ad35952..6379c69 100644
--- a/src/@dseries/log_.m
+++ b/src/@dseries/log_.m
@@ -1,4 +1,4 @@
-function o = log_(o) % --*-- Unitary tests --*--
+function o = log_(o)
 
 % Apply the logarithm to all the variables in a dseries object (in place modification).
 %
@@ -8,7 +8,7 @@ function o = log_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -38,7 +38,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -96,4 +96,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/lt.m b/src/@dseries/lt.m
index e7fe511..c42abab 100644
--- a/src/@dseries/lt.m
+++ b/src/@dseries/lt.m
@@ -1,8 +1,8 @@
-function l = lt(varargin)  % --*-- Unitary tests --*--
+function l = lt(varargin)
 
 % Overloads the < operator for dseries objects.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -25,4 +25,4 @@ elseif isdseries(varargin{1}) && isdseries(varargin{2})
     l = varargin{1}.data<varargin{2}.data;
 else
     error('dseries::lt: wrong types.')
-end
\ No newline at end of file
+end
diff --git a/src/@dseries/mdiff.m b/src/@dseries/mdiff.m
index ccbb125..ec26df4 100644
--- a/src/@dseries/mdiff.m
+++ b/src/@dseries/mdiff.m
@@ -1,4 +1,4 @@
-function o = mdiff(o) % --*-- Unitary tests --*--
+function o = mdiff(o)
 
 % Computes monthly differences.
 %
@@ -8,7 +8,7 @@ function o = mdiff(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = mdiff(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.mdiff_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/mdiff_.m b/src/@dseries/mdiff_.m
index 2c62147..381e21c 100644
--- a/src/@dseries/mdiff_.m
+++ b/src/@dseries/mdiff_.m
@@ -1,4 +1,4 @@
-function o = mdiff_(o) % --*-- Unitary tests --*--
+function o = mdiff_(o)
 
 % Computes monthly differences.
 %
@@ -8,7 +8,7 @@ function o = mdiff_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -51,7 +51,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/mean.m b/src/@dseries/mean.m
index e4ecb87..520f54f 100644
--- a/src/@dseries/mean.m
+++ b/src/@dseries/mean.m
@@ -1,4 +1,4 @@
-function m = mean(o, geometric) % --*-- Unitary tests --*--
+function m = mean(o, geometric)
 
 % Returns the mean of the variables in a @dseries object o.
 %
@@ -9,7 +9,7 @@ function m = mean(o, geometric) % --*-- Unitary tests --*--
 % OUTPUTS
 %  o m             1*vobs(o) vector of doubles.
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -36,7 +36,7 @@ else
     m = mean(o.data);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -126,4 +126,4 @@ if t(1)
    t(7) = isequal(m1, m3);
 end
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dseries/merge.m b/src/@dseries/merge.m
index dab0d99..566946d 100644
--- a/src/@dseries/merge.m
+++ b/src/@dseries/merge.m
@@ -1,11 +1,11 @@
-function q = merge(o, p, legacy) % --*-- Unitary tests --*--
+function q = merge(o, p, legacy)
 
 % Merge method for dseries objects.
 %
 % INPUTS
 % - o                [dseries]
 % - p                [dseries]
-% - legacy           [logical]      revert to legacy behaviour if `true` (default is `false`), 
+% - legacy           [logical]      revert to legacy behaviour if `true` (default is `false`),
 %
 % OUTPUTS
 % - q                [dseries]
@@ -15,7 +15,7 @@ function q = merge(o, p, legacy) % --*-- Unitary tests --*--
 % in p take precedence except if p has NaNs (the exception can be
 % removed by setting the third argument, legacy, equal to true).
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -153,7 +153,7 @@ end
 
 q.dates = q_init:q_init+(nobs(q)-1);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -283,4 +283,4 @@ if t(1)
    t(6) = dassert(y.dates(1), z.dates(1));
 end
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/@dseries/mgrowth.m b/src/@dseries/mgrowth.m
index e0ecf8c..82a87b3 100644
--- a/src/@dseries/mgrowth.m
+++ b/src/@dseries/mgrowth.m
@@ -1,4 +1,4 @@
-function o = mgrowth(o) % --*-- Unitary tests --*--
+function o = mgrowth(o)
 
 % Computes monthly growth rates.
 %
@@ -8,7 +8,7 @@ function o = mgrowth(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = mgrowth(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.mgrowth_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/mgrowth_.m b/src/@dseries/mgrowth_.m
index e11c6d3..9d5daa7 100644
--- a/src/@dseries/mgrowth_.m
+++ b/src/@dseries/mgrowth_.m
@@ -1,4 +1,4 @@
-function o = mgrowth_(o) % --*-- Unitary tests --*--
+function o = mgrowth_(o)
 
 % Computes monthly growth rates.
 %
@@ -8,7 +8,7 @@ function o = mgrowth_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -51,7 +51,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/minus.m b/src/@dseries/minus.m
index 1c88fc6..b7412fb 100644
--- a/src/@dseries/minus.m
+++ b/src/@dseries/minus.m
@@ -1,4 +1,4 @@
-function q = minus(o, p) % --*-- Unitary tests --*--
+function q = minus(o, p)
 
 % Overloads the minus operator for dseries objects.
 %
@@ -27,7 +27,7 @@ function q = minus(o, p) % --*-- Unitary tests --*--
 %  4Y | 0
 %  5Y | 0
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -169,7 +169,7 @@ end
 
 q.data = bsxfun(@minus, o.data, p.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/mpower.m b/src/@dseries/mpower.m
index 1b05610..43b1f8a 100644
--- a/src/@dseries/mpower.m
+++ b/src/@dseries/mpower.m
@@ -1,4 +1,4 @@
-function q = mpower(o, p) % --*-- Unitary tests --*--
+function q = mpower(o, p)
 
 % Overloads the power (^) operator for dseries objects.
 %
@@ -9,7 +9,7 @@ function q = mpower(o, p) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q [dseries]           T observations and N variables.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -88,7 +88,7 @@ end
 
 error('dseries:WrongInputArguments', 'Wrong calling sequence! Please check the manual.')
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -168,4 +168,4 @@ if t(1)
    t(6) = dassert(ts2.tex,{'MyVar1';'MyVar2'});
 end
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/mrdivide.m b/src/@dseries/mrdivide.m
index 972c906..ca29a1b 100644
--- a/src/@dseries/mrdivide.m
+++ b/src/@dseries/mrdivide.m
@@ -1,4 +1,4 @@
-function q = mrdivide(o, p) % --*-- Unitary tests --*--
+function q = mrdivide(o, p)
 
 % Overloads the mrdivde (/) operator for dseries objects.
 %
@@ -9,7 +9,7 @@ function q = mrdivide(o, p) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q [dseries]           T observations and N variables.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -121,7 +121,7 @@ else
     error()
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/mtimes.m b/src/@dseries/mtimes.m
index 43307db..d1edade 100644
--- a/src/@dseries/mtimes.m
+++ b/src/@dseries/mtimes.m
@@ -1,4 +1,4 @@
-function q = mtimes(o, p) % --*-- Unitary tests --*--
+function q = mtimes(o, p)
 
 % Overloads the mtimes (*) operator for dseries objects.
 %
@@ -9,7 +9,7 @@ function q = mtimes(o, p) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q [dseries]           T observations and N variables.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -133,7 +133,7 @@ else
     error()
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/nanmean.m b/src/@dseries/nanmean.m
index dd21abf..368f6bc 100644
--- a/src/@dseries/nanmean.m
+++ b/src/@dseries/nanmean.m
@@ -1,4 +1,4 @@
-function m = nanmean(o, geometric) % --*-- Unitary tests --*--
+function m = nanmean(o, geometric)
 
 % Returns the mean of the variables in a @dseries object o (robust
 % to the presence of NaNs).
@@ -10,7 +10,7 @@ function m = nanmean(o, geometric) % --*-- Unitary tests --*--
 % OUTPUTS
 % - m             1*vobs(o) vector of doubles.
 
-% Copyright © 2019-2021 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -41,7 +41,7 @@ else
     m = nanmean(o.data);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -83,4 +83,4 @@ if t(1)
    t(3) = dassert(m, [sqrt(3), 0]);
 end
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/nanstd.m b/src/@dseries/nanstd.m
index 757417b..31042f2 100644
--- a/src/@dseries/nanstd.m
+++ b/src/@dseries/nanstd.m
@@ -1,4 +1,4 @@
-function s = nanstd(o, geometric) % --*-- Unitary tests --*--
+function s = nanstd(o, geometric)
 
 % Returns the standard deviation of the variables in a @dseries object o (robust
 % to the presence of NaNs).
@@ -13,7 +13,7 @@ function s = nanstd(o, geometric) % --*-- Unitary tests --*--
 % REMARKS
 % See https://en.wikipedia.org/wiki/Geometric_standard_deviation for the definition of the geometric standard deviation.
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -49,7 +49,7 @@ else
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -71,4 +71,4 @@ if t(1)
    t(3) = all(abs(s-ones(1, 2))<1e-3);
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/ne.m b/src/@dseries/ne.m
index 9d14706..1bad95a 100644
--- a/src/@dseries/ne.m
+++ b/src/@dseries/ne.m
@@ -1,4 +1,4 @@
-function b = ne(o, p) % --*-- Unitary tests --*--
+function b = ne(o, p)
 
 % Overloads ne (~=) operator.
 %
@@ -12,7 +12,7 @@ function b = ne(o, p) % --*-- Unitary tests --*--
 % REMARKS
 %  If the number of variables, the number of observations or the frequencies are different in A and B, the function returns one.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,7 +75,7 @@ end
 
 b = ne(o.data, p.data);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -103,4 +103,4 @@ if t(1)
    t(2) = dassert(a,logical([ones(10,2), ones(10,1)]));
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/onesidedhpcycle.m b/src/@dseries/onesidedhpcycle.m
index 43daf89..51d80f0 100644
--- a/src/@dseries/onesidedhpcycle.m
+++ b/src/@dseries/onesidedhpcycle.m
@@ -1,4 +1,4 @@
-function o = onesidedhpcycle(o, lambda, init) % --*-- Unitary tests --*--
+function o = onesidedhpcycle(o, lambda, init)
 
 % Extracts the cycle component from a dseries object using the one sided HP filter.
 %
@@ -9,7 +9,7 @@ function o = onesidedhpcycle(o, lambda, init) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Cycle component of the original time series.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -46,7 +46,7 @@ else
     o.onesidedhpcycle_(lambda, init);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 d = randn(100,1);
@@ -64,4 +64,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/onesidedhpcycle_.m b/src/@dseries/onesidedhpcycle_.m
index 0cc1dab..3194ac2 100644
--- a/src/@dseries/onesidedhpcycle_.m
+++ b/src/@dseries/onesidedhpcycle_.m
@@ -1,4 +1,4 @@
-function o = onesidedhpcycle_(o, lambda, init) % --*-- Unitary tests --*--
+function o = onesidedhpcycle_(o, lambda, init)
 
 % Extracts the cycle component from a dseries object using a one sided HP filter.
 %
@@ -9,7 +9,7 @@ function o = onesidedhpcycle_(o, lambda, init) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Cycle component of the original time series.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,11 +75,11 @@ if nargin>2
     trend = o.hptrend(lambda);
     x0 = trend.data(1:2,:);
     [~, o.data] = one_sided_hp_filter(o.data, lambda, x0);
-else    
+else
     [~, o.data] = one_sided_hp_filter(o.data, lambda);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 e = .2*randn(200,1);
@@ -124,4 +124,4 @@ catch
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/onesidedhptrend.m b/src/@dseries/onesidedhptrend.m
index 75a82ca..3c38a77 100644
--- a/src/@dseries/onesidedhptrend.m
+++ b/src/@dseries/onesidedhptrend.m
@@ -1,4 +1,4 @@
-function o = onesidedhptrend(o, lambda, init) % --*-- Unitary tests --*--
+function o = onesidedhptrend(o, lambda, init)
 
 % Extracts the trend component from a dseries object using the one sided HP filter.
 %
@@ -9,7 +9,7 @@ function o = onesidedhptrend(o, lambda, init) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Trend component of the original time series.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -46,7 +46,7 @@ else
     o.onesidedhptrend_(lambda, init);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 d = randn(100,1);
@@ -64,4 +64,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/onesidedhptrend_.m b/src/@dseries/onesidedhptrend_.m
index 38458c6..71890c0 100644
--- a/src/@dseries/onesidedhptrend_.m
+++ b/src/@dseries/onesidedhptrend_.m
@@ -1,4 +1,4 @@
-function o = onesidedhptrend_(o, lambda) % --*-- Unitary tests --*--
+function o = onesidedhptrend_(o, lambda)
 
 % Extracts the trend component from a dseries object using a one sided HP filter.
 %
@@ -9,7 +9,7 @@ function o = onesidedhptrend_(o, lambda) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o         [dseries]   Trend component of the original time series.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,11 +75,11 @@ if nargin>2
     trend = o.hptrend(lambda);
     x0 = trend.data(1:2,:);
     o.data = one_sided_hp_filter(o.data, lambda, x0);
-else    
+else
     o.data = one_sided_hp_filter(o.data, lambda);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 e = .2*randn(200,1);
@@ -124,4 +124,4 @@ catch
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/plus.m b/src/@dseries/plus.m
index 5f65892..0b995df 100644
--- a/src/@dseries/plus.m
+++ b/src/@dseries/plus.m
@@ -1,4 +1,4 @@
-function q = plus(o, p) % --*-- Unitary tests --*--
+function q = plus(o, p)
 
 % Overloads the plus (+) operator for dseries objects.
 %
@@ -9,7 +9,7 @@ function q = plus(o, p) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q [dseries]
 
-% Copyright © 2011-2021 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -109,11 +109,11 @@ if isdseries(o) && isdseries(p)
             idp = ido;
         end
     end
-    
+
     if ~isequal(frequency(o),frequency(p))
         error(['dseries::plus: Cannot add ' inputname(1) ' and ' inputname(2) ' (frequencies are different)!'])
     end
-    
+
     if ~isequal(nobs(o), nobs(p)) || ~isequal(firstdate(o),firstdate(p))
         [o, p] = align(o, p);
     end
@@ -142,7 +142,7 @@ else
     error()
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/pop.m b/src/@dseries/pop.m
index e6ab629..6b616b4 100644
--- a/src/@dseries/pop.m
+++ b/src/@dseries/pop.m
@@ -1,4 +1,4 @@
-function [o, id] = pop(o, a) % --*-- Unitary tests --*--
+function [o, id] = pop(o, a)
 
 % Removes a variable from a dseries object.
 %
@@ -9,7 +9,7 @@ function [o, id] = pop(o, a) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]  T observations and N-1 variables.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -34,7 +34,7 @@ else
     [o, id] = pop_(o, a);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/pop_.m b/src/@dseries/pop_.m
index 68f4151..817f250 100644
--- a/src/@dseries/pop_.m
+++ b/src/@dseries/pop_.m
@@ -1,4 +1,4 @@
-function [o, id] = pop_(o, a) % --*-- Unitary tests --*--
+function [o, id] = pop_(o, a)
 
 % Removes a variable from a dseries object.
 %
@@ -9,7 +9,7 @@ function [o, id] = pop_(o, a) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]  T observations and N-1 variables.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -44,10 +44,10 @@ o.tex(id) = [];
 o.ops(id) = [];
 otagnames = fieldnames(o.tags);
 for i=1:length(otagnames)
-    o.tags.(otagnames{i})(id) = []; 
+    o.tags.(otagnames{i})(id) = [];
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -76,4 +76,4 @@ if t(1)
    t(5) = dassert(ts1.tags.type, {1;3});
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/projection.m b/src/@dseries/projection.m
index ef038b7..520361e 100644
--- a/src/@dseries/projection.m
+++ b/src/@dseries/projection.m
@@ -1,4 +1,4 @@
-function o = projection(o, info, periods)  % --*-- Unitary tests --*--
+function o = projection(o, info, periods)
 
 % Projects variables in dseries object o.
 %
@@ -31,6 +31,23 @@ function o = projection(o, info, periods)  % --*-- Unitary tests --*--
 %     or 'AR' with an autoregressive parameter equal to one (random walk).
 % [4] This projection routine only deals with exponential trends.
 
+% Copyright © 2022-2023 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/>.
+
 % Fetch useful lines in info
 [~, ~, i2] = intersect(o.name, info(:,1));
 INFO = info(i2,:);
@@ -70,7 +87,7 @@ for i=1:rows(INFO)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -91,4 +108,4 @@ if t(1)
     t(6) = dassert(ts.A4.data(11), 0.2, 1e-15) & dassert(ts.A4.data(12), 0.2*0.4-0.2, 1e-15);
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/qdiff.m b/src/@dseries/qdiff.m
index d6d5074..bf060c8 100644
--- a/src/@dseries/qdiff.m
+++ b/src/@dseries/qdiff.m
@@ -1,4 +1,4 @@
-function o = qdiff(o) % --*-- Unitary tests --*--
+function o = qdiff(o)
 
 % Computes quaterly differences.
 %
@@ -8,7 +8,7 @@ function o = qdiff(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = qdiff(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.qdiff_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/qdiff_.m b/src/@dseries/qdiff_.m
index a65fafc..10d1543 100644
--- a/src/@dseries/qdiff_.m
+++ b/src/@dseries/qdiff_.m
@@ -1,4 +1,4 @@
-function o = qdiff_(o) % --*-- Unitary tests --*--
+function o = qdiff_(o)
 
 % Computes quaterly differences.
 %
@@ -8,7 +8,7 @@ function o = qdiff_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -52,7 +52,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/qgrowth.m b/src/@dseries/qgrowth.m
index c4c380b..9f8e79a 100644
--- a/src/@dseries/qgrowth.m
+++ b/src/@dseries/qgrowth.m
@@ -1,4 +1,4 @@
-function o = qgrowth(o) % --*-- Unitary tests --*--
+function o = qgrowth(o)
 
 % Computes quaterly growth rates.
 %
@@ -8,7 +8,7 @@ function o = qgrowth(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = qgrowth(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.qgrowth_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -48,4 +48,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/qgrowth_.m b/src/@dseries/qgrowth_.m
index 6aa03ba..78aff69 100644
--- a/src/@dseries/qgrowth_.m
+++ b/src/@dseries/qgrowth_.m
@@ -1,4 +1,4 @@
-function o = qgrowth_(o) % --*-- Unitary tests --*--
+function o = qgrowth_(o)
 
 % Computes quaterly growth rates.
 %
@@ -8,7 +8,7 @@ function o = qgrowth_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -52,7 +52,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/remove.m b/src/@dseries/remove.m
index 0739e06..738c9f0 100644
--- a/src/@dseries/remove.m
+++ b/src/@dseries/remove.m
@@ -1,8 +1,8 @@
-function o = remove(o, a) % --*-- Unitary tests --*--
+function o = remove(o, a)
 
 % Removes a variable from a dseries object (alias for the pop method).
 
-% Copyright © 2014-2017 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -19,4 +19,4 @@ function o = remove(o, a) % --*-- Unitary tests --*--
 % You should have received a copy of the GNU General Public License
 % along with Dynare.  If not, see <https://www.gnu.org/licenses/>.
 
-o = pop(o, a);
\ No newline at end of file
+o = pop(o, a);
diff --git a/src/@dseries/remove_.m b/src/@dseries/remove_.m
index 474ccec..50e9c8a 100644
--- a/src/@dseries/remove_.m
+++ b/src/@dseries/remove_.m
@@ -1,8 +1,8 @@
-function o = remove_(o, a) % --*-- Unitary tests --*--
+function o = remove_(o, a)
 
 % Removes a variable from a dseries object (alias for the pop_ method).
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -19,4 +19,4 @@ function o = remove_(o, a) % --*-- Unitary tests --*--
 % You should have received a copy of the GNU General Public License
 % along with Dynare.  If not, see <https://www.gnu.org/licenses/>.
 
-o = pop_(o, a);
\ No newline at end of file
+o = pop_(o, a);
diff --git a/src/@dseries/rename.m b/src/@dseries/rename.m
index d0f71a9..b541d94 100644
--- a/src/@dseries/rename.m
+++ b/src/@dseries/rename.m
@@ -1,4 +1,4 @@
-function o = rename(o, old, new) % --*-- Unitary tests --*--
+function o = rename(o, old, new)
 
 % Renames variables in a dseries object.
 %
@@ -10,7 +10,7 @@ function o = rename(o, old, new) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o     [dseries]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -34,7 +34,7 @@ else
     o.rename_(old, new);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 ts = dseries([transpose(1:5), transpose(6:10)],'1950q1',{'Output'; 'Consumption'}, {'Y_t'; 'C_t'});
@@ -72,4 +72,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/@dseries/rename_.m b/src/@dseries/rename_.m
index 794d4a6..96f6301 100644
--- a/src/@dseries/rename_.m
+++ b/src/@dseries/rename_.m
@@ -1,4 +1,4 @@
-function o = rename_(o, old, new) % --*-- Unitary tests --*--
+function o = rename_(o, old, new)
 
 % Renames variables in a dseries object.
 %
@@ -10,7 +10,7 @@ function o = rename_(o, old, new) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o     [dseries]
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -83,7 +83,7 @@ else
     o.name(idname) = {new};
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 ts = dseries([transpose(1:5), transpose(6:10)],'1950q1',{'Output'; 'Consumption'}, {'Y_t'; 'C_t'});
diff --git a/src/@dseries/round.m b/src/@dseries/round.m
index 509bd1f..d1b2411 100644
--- a/src/@dseries/round.m
+++ b/src/@dseries/round.m
@@ -1,4 +1,4 @@
-function o = round(o, n) % --*-- Unitary tests --*--
+function o = round(o, n)
 
 % Rounds to nearest decimal or integer.
 %
@@ -9,7 +9,7 @@ function o = round(o, n) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -33,7 +33,7 @@ end
 o = copy(o);
 o.round_(n);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -54,4 +54,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/round_.m b/src/@dseries/round_.m
index d6389e2..1a6be10 100644
--- a/src/@dseries/round_.m
+++ b/src/@dseries/round_.m
@@ -1,4 +1,4 @@
-function o = round_(o, n) % --*-- Unitary tests --*--
+function o = round_(o, n)
 
 % Rounds to nearest decimal or integer.
 %
@@ -9,7 +9,7 @@ function o = round_(o, n) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o [dseries]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -44,7 +44,7 @@ else
     o.data = round(o.data, n);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dates object
@@ -64,4 +64,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/save.m b/src/@dseries/save.m
index 253e372..5977476 100644
--- a/src/@dseries/save.m
+++ b/src/@dseries/save.m
@@ -1,8 +1,8 @@
-function save(o, basename, format) % --*-- Unitary tests --*--
+function save(o, basename, format)
 
 % Saves a dseries object on disk.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -128,7 +128,7 @@ switch format
     error('Unknown format!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
@@ -212,4 +212,4 @@ end
 delete('dynare_series.mat');
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dseries/set_names.m b/src/@dseries/set_names.m
index 8ac5b80..db2f702 100644
--- a/src/@dseries/set_names.m
+++ b/src/@dseries/set_names.m
@@ -1,4 +1,4 @@
-function o = set_names(o, varargin) % --*-- Unitary tests --*--
+function o = set_names(o, varargin)
 
 % Specifies names of the variables in a dseries object (in place modification).
 %
@@ -9,7 +9,7 @@ function o = set_names(o, varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o                 [dseries]
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -50,7 +50,7 @@ for i=1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -76,4 +76,4 @@ if length(t)>1
    t(4) = dassert(ts1.name,A_name);
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/size.m b/src/@dseries/size.m
index caa5a85..197e40a 100644
--- a/src/@dseries/size.m
+++ b/src/@dseries/size.m
@@ -1,8 +1,8 @@
-function varargout = size(o, varargin) % --*-- Unitary tests --*--
+function varargout = size(o, varargin)
 
 % Overloads size function.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,9 +35,9 @@ switch nargout
     error('dseries::size: Wrong calling sequence! Cannot return more than two arguments.')
 end
 
-return
+return % --*-- Unit tests --*--
 
-%@test:2
+%@test:1
 ts = dseries(randn(10,1));
 try
     dimension = ts.size();
@@ -52,4 +52,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/std.m b/src/@dseries/std.m
index 85692ef..695013f 100644
--- a/src/@dseries/std.m
+++ b/src/@dseries/std.m
@@ -1,4 +1,4 @@
-function s = std(o, geometric) % --*-- Unitary tests --*--
+function s = std(o, geometric)
 
 % Returns the standard deviation of the variables in a @dseries object o.
 % See https://en.wikipedia.org/wiki/Geometric_standard_deviation
@@ -10,7 +10,7 @@ function s = std(o, geometric) % --*-- Unitary tests --*--
 % OUTPUTS
 % - s             [double]    1*N vector.
 
-% Copyright © 2016-2021 Dynare Team
+% Copyright © 2016-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -38,7 +38,7 @@ else
     s = std(o.data);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a dataset.
@@ -124,4 +124,4 @@ if t(1)
    t(3) = dassert(max(abs(s-[.1, .1]))<.0001, true);
 end
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dseries/subsample.m b/src/@dseries/subsample.m
index 4f0d5d0..a373ef2 100644
--- a/src/@dseries/subsample.m
+++ b/src/@dseries/subsample.m
@@ -1,4 +1,4 @@
-function p = subsample(o, d1, d2)  % --*-- Unitary tests --*--
+function p = subsample(o, d1, d2)
 
 % function p = subsample(o, d1, d2)
 
@@ -24,7 +24,7 @@ function p = subsample(o, d1, d2)  % --*-- Unitary tests --*--
 % 3Y | 3
 % 4Y | 4
 
-% Copyright © 2019-2021 Dynare Team
+% Copyright © 2019-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -93,7 +93,7 @@ p.dates = o.dates(idx);
 p.ops   = o.ops;
 p.tags  = o.tags;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/@dseries/subsasgn.m b/src/@dseries/subsasgn.m
index 886da17..a4eb972 100644
--- a/src/@dseries/subsasgn.m
+++ b/src/@dseries/subsasgn.m
@@ -1,4 +1,4 @@
-function A = subsasgn(A,S,B) % --*-- Unitary tests --*--
+function A = subsasgn(A,S,B)
 
 %@info:
 %! @deftypefn {Function File} {@var{A} =} subsasgn (@var{A}, @var{S}, @var{B})
@@ -8,7 +8,7 @@ function A = subsasgn(A,S,B) % --*-- Unitary tests --*--
 %! @end deftypefn
 %@eod:
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -234,7 +234,7 @@ if merge_dseries_objects
     A = merge(A, B, true);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -859,4 +859,4 @@ if t(1)
    end
 end
 T = all(t);
-%@eof:22
\ No newline at end of file
+%@eof:22
diff --git a/src/@dseries/subsref.m b/src/@dseries/subsref.m
index b5c101b..e2c4c8f 100644
--- a/src/@dseries/subsref.m
+++ b/src/@dseries/subsref.m
@@ -1,4 +1,4 @@
-function r = subsref(o, S) % --*-- Unitary tests --*--
+function r = subsref(o, S)
 
 % Overloads the subsref method for dseries class.
 %
@@ -15,7 +15,7 @@ function r = subsref(o, S) % --*-- Unitary tests --*--
 %                     by applying a public method on `o`, or a dseries object built by extracting
 %                     a variable from `o`, or a dseries object containing a subsample.
 
-% Copyright © 2011-2022 Dynare Team
+% Copyright © 2011-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -319,7 +319,7 @@ if ~isempty(S)
     r = subsref(r, S);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
  % Define a data set.
diff --git a/src/@dseries/tag.m b/src/@dseries/tag.m
index 09bdd75..e20f382 100644
--- a/src/@dseries/tag.m
+++ b/src/@dseries/tag.m
@@ -1,4 +1,4 @@
-function o = tag(o, a, b, c) % --*-- Unitary tests --*--
+function o = tag(o, a, b, c)
 
 % Add tag to a dseries oject (in place modification).
 
@@ -10,8 +10,8 @@ function o = tag(o, a, b, c) % --*-- Unitary tests --*--
 %
 % OUTPUT
 % - o        [dseries]    Updated with tag
-    
-% Copyright © 2017-2021 Dynare Team
+
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -47,7 +47,7 @@ else
     o.tags.(a)(id) = {c};
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 ts = dseries(randn(10, 3));
@@ -121,4 +121,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/@dseries/tex_rename.m b/src/@dseries/tex_rename.m
index 0c20470..bf0e670 100644
--- a/src/@dseries/tex_rename.m
+++ b/src/@dseries/tex_rename.m
@@ -1,6 +1,6 @@
-function o = tex_rename(o, varargin) % --*-- Unitary tests --*--
+function o = tex_rename(o, varargin)
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -20,7 +20,7 @@ function o = tex_rename(o, varargin) % --*-- Unitary tests --*--
 o = copy(o);
 o.tex_rename_(varargin{:});
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 ts = dseries([transpose(1:5), transpose(6:10)],'1950q1',{'Output'; 'Consumption'}, {'Y_t'; 'C_t'});
@@ -43,4 +43,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/tex_rename_.m b/src/@dseries/tex_rename_.m
index aa147c1..7157a7d 100644
--- a/src/@dseries/tex_rename_.m
+++ b/src/@dseries/tex_rename_.m
@@ -1,4 +1,4 @@
-function o = tex_rename_(o, varargin) % --*-- Unitary tests --*--
+function o = tex_rename_(o, varargin)
 
 % Changes TeX name of variables in a dseries object.
 %
@@ -16,7 +16,7 @@ function o = tex_rename_(o, varargin) % --*-- Unitary tests --*--
 % first element is the name of the variable and the second argument is the new texname for this
 % variable.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -64,7 +64,7 @@ else
     o.tex(idname) = {newtexname};
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 ts = dseries([transpose(1:5), transpose(6:10)],'1950q1',{'Output'; 'Consumption'}, {'Y_t'; 'C_t'});
@@ -102,4 +102,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/@dseries/uminus.m b/src/@dseries/uminus.m
index 3b3d10e..11c8db0 100644
--- a/src/@dseries/uminus.m
+++ b/src/@dseries/uminus.m
@@ -1,4 +1,4 @@
-function o = uminus(o) % --*-- Unitary tests --*--
+function o = uminus(o)
 
 % Overloads the uminus method for dseries objects.
 %
@@ -8,7 +8,7 @@ function o = uminus(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -37,7 +37,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
diff --git a/src/@dseries/vertcat.m b/src/@dseries/vertcat.m
index 9bcebbc..6551366 100644
--- a/src/@dseries/vertcat.m
+++ b/src/@dseries/vertcat.m
@@ -1,4 +1,4 @@
-function o = vertcat(varargin) % --*-- Unitary tests --*--
+function o = vertcat(varargin)
 
 % Overloads vertcat method for dseries objects.
 %
@@ -17,7 +17,7 @@ function o = vertcat(varargin) % --*-- Unitary tests --*--
 %
 %  merges the samples.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -70,7 +70,7 @@ else
 end
 d.dates = [b.dates; c.dates];
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a data set.
diff --git a/src/@dseries/ydiff.m b/src/@dseries/ydiff.m
index 3c4fdcb..4c19282 100644
--- a/src/@dseries/ydiff.m
+++ b/src/@dseries/ydiff.m
@@ -1,4 +1,4 @@
-function o = ydiff(o) % --*-- Unitary tests --*--
+function o = ydiff(o)
 
 % Computes yearly differences.
 %
@@ -8,7 +8,7 @@ function o = ydiff(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2020 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = ydiff(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.ydiff_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -50,4 +50,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/ydiff_.m b/src/@dseries/ydiff_.m
index 582370a..710d3cd 100644
--- a/src/@dseries/ydiff_.m
+++ b/src/@dseries/ydiff_.m
@@ -1,4 +1,4 @@
-function o = ydiff_(o) % --*-- Unitary tests --*--
+function o = ydiff_(o)
 
 % Computes yearly differences.
 %
@@ -8,7 +8,7 @@ function o = ydiff_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2020 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -56,7 +56,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -159,4 +159,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/@dseries/ygrowth.m b/src/@dseries/ygrowth.m
index 718b7e2..cefb86b 100644
--- a/src/@dseries/ygrowth.m
+++ b/src/@dseries/ygrowth.m
@@ -1,4 +1,4 @@
-function o = ygrowth(o) % --*-- Unitary tests --*--
+function o = ygrowth(o)
 
 % Computes yearly growth rates.
 %
@@ -8,7 +8,7 @@ function o = ygrowth(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function o = ygrowth(o) % --*-- Unitary tests --*--
 o = copy(o);
 o.ygrowth_();
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(2,1);
@@ -51,4 +51,4 @@ if length(t)>1
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@dseries/ygrowth_.m b/src/@dseries/ygrowth_.m
index d356455..1e2e9e9 100644
--- a/src/@dseries/ygrowth_.m
+++ b/src/@dseries/ygrowth_.m
@@ -1,4 +1,4 @@
-function o = ygrowth_(o) % --*-- Unitary tests --*--
+function o = ygrowth_(o)
 
 % Computes yearly growth rates.
 %
@@ -8,7 +8,7 @@ function o = ygrowth_(o) % --*-- Unitary tests --*--
 % OUTPUTS
 % - o   [dseries]
 
-% Copyright © 2012-2020 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -31,7 +31,7 @@ switch frequency(o)
     o.data(1,:) = NaN;
   case 2
     o.data(3:end,:) = o.data(3:end,:)./o.data(1:end-2,:) - 1;
-    o.data(1:2,:) = NaN;  
+    o.data(1:2,:) = NaN;
   case 4
     o.data(5:end,:) = o.data(5:end,:)./o.data(1:end-4,:) - 1;
     o.data(1:4,:) = NaN;
@@ -56,7 +56,7 @@ for i = 1:vobs(o)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -135,4 +135,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/@x13/clean.m b/src/@x13/clean.m
index fe58c04..b6deb0b 100644
--- a/src/@x13/clean.m
+++ b/src/@x13/clean.m
@@ -1,8 +1,8 @@
-function clean(o) % --*-- Unitary tests --*--
+function clean(o)
 
 % Erase generated files if any.
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -22,7 +22,7 @@ if ~isempty(o.results)
     delete(sprintf('%s.*', basename))
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -37,4 +37,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@x13/print.m b/src/@x13/print.m
index a733752..05216a6 100644
--- a/src/@x13/print.m
+++ b/src/@x13/print.m
@@ -1,8 +1,8 @@
-function basename = print(o, basename) % --*-- Unitary tests --*--
+function basename = print(o, basename)
 
 % Prints spc file.
 
-% Copyright © 2017 Dynare Team
+% Copyright © 2017-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -392,7 +392,7 @@ end
 
 fclose(fid);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -423,4 +423,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@x13/subsasgn.m b/src/@x13/subsasgn.m
index 48a4aea..aa5980b 100644
--- a/src/@x13/subsasgn.m
+++ b/src/@x13/subsasgn.m
@@ -1,6 +1,6 @@
-function val = subsasgn(val, idx, rhs) % --*-- Unitary tests --*--
+function val = subsasgn(val, idx, rhs)
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -17,7 +17,7 @@ function val = subsasgn(val, idx, rhs) % --*-- Unitary tests --*--
 
 error('Members of x13 class are private')
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(3,1);
@@ -47,4 +47,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@x13/subsref.m b/src/@x13/subsref.m
index c5373f6..80ae417 100644
--- a/src/@x13/subsref.m
+++ b/src/@x13/subsref.m
@@ -1,8 +1,8 @@
-function o = subsref(o, S) % --*-- Unitary tests --*--
+function o = subsref(o, S)
 
 % Overloads the subsref method.
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -95,7 +95,7 @@ switch S(1).type
     error('x13:: I do not understand what you are asking for!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(4,1);
@@ -133,4 +133,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/@x13/x13.m b/src/@x13/x13.m
index 2d1773f..755c53d 100644
--- a/src/@x13/x13.m
+++ b/src/@x13/x13.m
@@ -1,4 +1,4 @@
-classdef x13<handle % --*-- Unitary tests --*--
+classdef x13<handle
 
 % Class for X13 toolbox.
 
@@ -116,9 +116,8 @@ methods
         o.results = struct();
         o.commands = {};
     end
-end
-
-end
+end % methods
+end % classdef --*-- Unit tests --*--
 
 %@test:1
 %$
@@ -156,4 +155,4 @@ end
 %$ end
 %$
 %$ T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/mdbnomics2dseries/mdbnomics2dseries.m b/src/mdbnomics2dseries/mdbnomics2dseries.m
index 12712b3..ba10e49 100644
--- a/src/mdbnomics2dseries/mdbnomics2dseries.m
+++ b/src/mdbnomics2dseries/mdbnomics2dseries.m
@@ -1,4 +1,4 @@
-function ds = mdbnomics2dseries(varargin)   % --*-- Unitary tests --*--
+function ds = mdbnomics2dseries(varargin)
 
 % Given cell array from the mdbnomics library, it returns a dseries object.
 %
@@ -10,7 +10,7 @@ function ds = mdbnomics2dseries(varargin)   % --*-- Unitary tests --*--
 % OUTPUTS
 % - ds [dseries]
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-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
@@ -26,107 +26,108 @@ function ds = mdbnomics2dseries(varargin)   % --*-- Unitary tests --*--
 % along with Dynare.  If not, see <https://www.gnu.org/licenses/>.
 
 switch nargin
-    case 0
-        % Return empty object.
-        error('mdbnomics2dseries:WrongInputArguments', 'Input must be non empty!');
-    case 1
-        if iscell(varargin{1})
-            switch length(varargin{1})
-                case 0
-                    error('mdbnomics2dseries:WrongInputArguments', 'Input must be non empty!');
-                otherwise
-                    o.data = varargin{1}(2:end,:);
-                    o.cols = varargin{1}(1,:);
-                    col_idx = {'x_frequency', 'dataset_code', 'series_code', 'original_period', 'period', 'value'};
-                    for ii = 1:size(col_idx,2)
-                        o.col_idx.(col_idx{ii}) = find(strcmp(col_idx{ii}, o.cols));
-                    end
-                    % Check if database has multiple frequencies
-                    if size(unique(o.data(:, o.col_idx.x_frequency)),1) > 1
-                        error('mdbnomics2dseries:DatabaseCheck: The database, that you are trying to convert, contains multiple frequencies. Currently, this type of dseries conversion is not supported. Please select a section of your database with uniform frequency.');
-                    end
-                    ds = convert_mdbnomics(o);
+  case 0
+    % Return empty object.
+    error('mdbnomics2dseries:WrongInputArguments', 'Input must be non empty!');
+  case 1
+    if iscell(varargin{1})
+        switch length(varargin{1})
+          case 0
+            error('mdbnomics2dseries:WrongInputArguments', 'Input must be non empty!');
+          otherwise
+            o.data = varargin{1}(2:end,:);
+            o.cols = varargin{1}(1,:);
+            col_idx = {'x_frequency', 'dataset_code', 'series_code', 'original_period', 'period', 'value'};
+            for ii = 1:size(col_idx,2)
+                o.col_idx.(col_idx{ii}) = find(strcmp(col_idx{ii}, o.cols));
             end
+            % Check if database has multiple frequencies
+            if size(unique(o.data(:, o.col_idx.x_frequency)),1) > 1
+                error('mdbnomics2dseries:DatabaseCheck: The database, that you are trying to convert, contains multiple frequencies. Currently, this type of dseries conversion is not supported. Please select a section of your database with uniform frequency.');
+            end
+            ds = convert_mdbnomics(o);
         end
-    otherwise
-        error('mdbnomics2dseries:WrongInputArguments', 'Too many input arguments! Please check the manual.')
-end
+    end
+  otherwise
+    error('mdbnomics2dseries:WrongInputArguments', 'Too many input arguments! Please check the manual.')
 end
 
+return % --*-- Unit tests --*--
+
 %@test:1
-%$ try
-%$     dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
-%$     df_id = load([dseries_src_root '../tests/data/mdbnomics/df_id.mat']);
-%$     df = df_id.df;
-%$     ds = mdbnomics2dseries(df);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     t(2) = dassert(ds.name, {'EA19_1_0_0_0_ZUTN'});
-%$     t(3) = dassert(ds.vobs, 1);
-%$     t(4) = dassert(ds.tags.provider_code, {'AMECO'});
-%$     t(5:10) = isfield(ds.tags, {'freq', 'unit', 'geo', 'Frequency', 'Unit', 'Country'});
-%$ end
-%$
-%$ T = all(t);
+try
+    dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
+    df_id = load([dseries_src_root '../tests/data/mdbnomics/df_id.mat']);
+    df = df_id.df;
+    ds = mdbnomics2dseries(df);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    t(2) = dassert(ds.name, {'EA19_1_0_0_0_ZUTN'});
+    t(3) = dassert(ds.vobs, 1);
+    t(4) = dassert(ds.tags.provider_code, {'AMECO'});
+    t(5:10) = isfield(ds.tags, {'freq', 'unit', 'geo', 'Frequency', 'Unit', 'Country'});
+end
+
+T = all(t);
 %@eof:1
 
 %@test:2
-%$ try
-%$     dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
-%$     df_bi = load([dseries_src_root '../tests/data/mdbnomics/df_bi-annual.mat']);
-%$     df = df_bi.df;
-%$     ds = mdbnomics2dseries(df);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     t(2) = dassert(ds.name, {'S_EUR_AL'});
-%$     t(3) = dassert(ds.vobs, 1);
-%$     t(4) = dassert(ds.dates.freq, 2);
-%$     t(4) = dassert(ds.tags.provider_code, {'Eurostat'});
-%$     t(5:10) = isfield(ds.tags, {'FREQ', 'currency', 'geo', 'Frequency', 'Currency', 'Geopolitical_entity__reporting_'});
-%$ end
-%$
-%$ T = all(t);
+try
+    dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
+    df_bi = load([dseries_src_root '../tests/data/mdbnomics/df_bi-annual.mat']);
+    df = df_bi.df;
+    ds = mdbnomics2dseries(df);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    t(2) = dassert(ds.name, {'S_EUR_AL'});
+    t(3) = dassert(ds.vobs, 1);
+    t(4) = dassert(ds.dates.freq, 2);
+    t(4) = dassert(ds.tags.provider_code, {'Eurostat'});
+    t(5:10) = isfield(ds.tags, {'FREQ', 'currency', 'geo', 'Frequency', 'Currency', 'Geopolitical_entity__reporting_'});
+end
+
+T = all(t);
 %@eof:2
 
 %@test:3
-%$ try
-%$     dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
-%$     df_dataset = load([dseries_src_root '../tests/data/mdbnomics/df_dataset.mat']);
-%$     df = df_dataset.df;
-%$     ds = mdbnomics2dseries(df);
-%$     t(1) = true;
-%$ catch
-%$     t(1) = false;
-%$ end
-%$
-%$ if t(1)
-%$     t(2) = dassert(ds.vobs, 49);
-%$     t(3) = dassert(ds.dates.freq, 1);
-%$     t(4) = dassert(length(unique(ds.tags.provider_code)), 1);
-%$     t(5) = dassert(size(ds.tags.provider_code,1), 49);
-%$ end
-%$
-%$ T = all(t);
+try
+    dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
+    df_dataset = load([dseries_src_root '../tests/data/mdbnomics/df_dataset.mat']);
+    df = df_dataset.df;
+    ds = mdbnomics2dseries(df);
+    t(1) = true;
+catch
+    t(1) = false;
+end
+
+if t(1)
+    t(2) = dassert(ds.vobs, 49);
+    t(3) = dassert(ds.dates.freq, 1);
+    t(4) = dassert(length(unique(ds.tags.provider_code)), 1);
+    t(5) = dassert(size(ds.tags.provider_code,1), 49);
+end
+
+T = all(t);
 %@eof:3
 
 %@test:4
-%$ try
-%$     dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
-%$     df_multi = load([dseries_src_root '../tests/data/mdbnomics/df_multi-freq.mat']);
-%$     df = df_multi.df;
-%$     ds = mdbnomics2dseries(df);
-%$     t(1) = false;
-%$ catch
-%$     t(1) = true;
-%$ end
-%$
-%$ T = all(t);
+try
+    dseries_src_root = strrep(which('initialize_dseries_class'),'initialize_dseries_class.m','');
+    df_multi = load([dseries_src_root '../tests/data/mdbnomics/df_multi-freq.mat']);
+    df = df_multi.df;
+    ds = mdbnomics2dseries(df);
+    t(1) = false;
+catch
+    t(1) = true;
+end
+
+T = all(t);
 %@eof:4
diff --git a/src/read/load_csv_file_data.m b/src/read/load_csv_file_data.m
index 0429105..bfddf0c 100644
--- a/src/read/load_csv_file_data.m
+++ b/src/read/load_csv_file_data.m
@@ -1,4 +1,4 @@
-function [freq, init, data, varlist] = load_csv_file_data(file) % --*-- Unitary tests --*--
+function [freq, init, data, varlist] = load_csv_file_data(file)
 %function [freq, init, data, varlist] = load_csv_file_data(file)
 % Loads data in a csv file.
 %
@@ -16,7 +16,7 @@ function [freq, init, data, varlist] = load_csv_file_data(file) % --*-- Unitary
 %  names. Similarly, if the first column does not contain dates, then
 %  freq will be 1 and init will be year 1.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -193,7 +193,7 @@ varlist = transpose(varlist);
 % Remove double quotes if any
 varlist = strrep(varlist,'"','');
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Download csv file with data.
diff --git a/src/read/load_m_file_data.m b/src/read/load_m_file_data.m
index ae6efc7..4c6017a 100644
--- a/src/read/load_m_file_data.m
+++ b/src/read/load_m_file_data.m
@@ -1,4 +1,4 @@
-function [freq, init, data, varlist, tex, ops, tags] = load_m_file_data(file) % --*-- Unitary tests --*--
+function [freq, init, data, varlist, tex, ops, tags] = load_m_file_data(file)
 
 % Loads data in a matlab/octave script.
 %
@@ -16,7 +16,7 @@ function [freq, init, data, varlist, tex, ops, tags] = load_m_file_data(file) %
 % be a scalar integer and INIT__ a string like '1938M11', '1945Q3', '1973W3' or '2009A'. If these variables are not specified
 % default values for freq and init are 1 and dates(1,1).
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -55,7 +55,7 @@ end
 
 [freq, init, data, varlist, tex, ops, tags] = load_mat_file_data(tmp__dataSet__struct);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
  % Create a data m-file
diff --git a/src/read/load_mat_file_data.m b/src/read/load_mat_file_data.m
index 74753d5..d4bb714 100644
--- a/src/read/load_mat_file_data.m
+++ b/src/read/load_mat_file_data.m
@@ -1,4 +1,4 @@
-function [freq, init, data, varlist, tex, ops, tags] = load_mat_file_data(file)  % --*-- Unitary tests --*--
+function [freq, init, data, varlist, tex, ops, tags] = load_mat_file_data(file)
 
 % Loads data in a matlab/octave mat-file.
 %
@@ -17,7 +17,7 @@ function [freq, init, data, varlist, tex, ops, tags] = load_mat_file_data(file)
 % be a scalar integer and INIT__ a string like '1938M11', '1945Q3', '1973W3' or '2009A'. If these variables are not specified
 % default values for freq and init are 1 and dates(1,1).
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -120,7 +120,7 @@ if isempty(ops)
     ops = cell(length(varlist), 1);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Create a data mat-file
@@ -199,4 +199,4 @@ t(8) = dassert(ops,{'grandad(hagop)';'dad(bedros)'});
 t(9) = dassert(data(:,1),[1;2;3;4;5]);
 t(10) = dassert(data(:,2),[2;3;4;5;6]);
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/convert/date2string.m b/src/utilities/convert/date2string.m
index 47a40e2..851c966 100644
--- a/src/utilities/convert/date2string.m
+++ b/src/utilities/convert/date2string.m
@@ -1,4 +1,4 @@
-function s = date2string(varargin) % --*-- Unitary tests --*--
+function s = date2string(varargin)
 
 % Returns date as a string.
 %
@@ -10,7 +10,7 @@ function s = date2string(varargin) % --*-- Unitary tests --*--
 % OUTPUTS
 % - s               [char]        1×n array.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -75,7 +75,7 @@ else
     s = sprintf('%i%s%i', year, freq2string(freq), subperiod);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -161,4 +161,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/utilities/convert/dates2M.m b/src/utilities/convert/dates2M.m
index 4591f09..c8d2602 100644
--- a/src/utilities/convert/dates2M.m
+++ b/src/utilities/convert/dates2M.m
@@ -1,4 +1,4 @@
-function m = dates2M(d) % --*-- Unitary tests --*--
+function m = dates2M(d)
 
 % INPUTS
 % - d   [dates]    daily frequency object with n elements.
@@ -6,7 +6,7 @@ function m = dates2M(d) % --*-- Unitary tests --*--
 % OUTPUTS
 % - m   [dates]    monthly frequency object with n elements (with repetitions).
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -32,7 +32,7 @@ end
 
 m = dates(12, time(:,1:2));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -73,4 +73,4 @@ catch
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/utilities/convert/dates2Q.m b/src/utilities/convert/dates2Q.m
index c29daee..1ba673c 100644
--- a/src/utilities/convert/dates2Q.m
+++ b/src/utilities/convert/dates2Q.m
@@ -1,4 +1,4 @@
-function m = dates2Q(d)  % --*-- Unitary tests --*--
+function m = dates2Q(d)
 
 % INPUTS
 % - d   [dates]    daily, monthly or quarterly frequency object with n elements.
@@ -6,7 +6,7 @@ function m = dates2Q(d)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - q   [dates]    quaterly frequency object with n elements (with repetitions).
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-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
@@ -46,7 +46,7 @@ time(isq4,2) = 4;
 
 m = dates(4, time(:,1:2));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -102,4 +102,4 @@ catch
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/utilities/convert/dates2S.m b/src/utilities/convert/dates2S.m
index 2ff307d..1453759 100644
--- a/src/utilities/convert/dates2S.m
+++ b/src/utilities/convert/dates2S.m
@@ -1,4 +1,4 @@
-function m = dates2S(d) % --*-- Unitary tests --*--
+function m = dates2S(d)
 
 % INPUTS
 % - d   [dates]    daily, monthly, quarterly or bi-annual frequency object with n elements.
@@ -6,7 +6,7 @@ function m = dates2S(d) % --*-- Unitary tests --*--
 % OUTPUTS
 % - q   [dates]    bi-annual frequency object with n elements (with repetitions).
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -47,7 +47,7 @@ time(ish2,2) = 2;
 
 m = dates(2, time(:,1:2));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -107,4 +107,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/utilities/convert/dates2Y.m b/src/utilities/convert/dates2Y.m
index 03314cc..6057456 100644
--- a/src/utilities/convert/dates2Y.m
+++ b/src/utilities/convert/dates2Y.m
@@ -1,4 +1,4 @@
-function y = dates2Y(d) % --*-- Unitary tests --*--
+function y = dates2Y(d)
 
 % INPUTS
 % - d   [dates]    daily, monthly, qarterly or bi-annual frequency object with n elements.
@@ -6,7 +6,7 @@ function y = dates2Y(d) % --*-- Unitary tests --*--
 % OUTPUTS
 % - y   [dates]    annual frequency object with n elements (with repetitions).
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-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
@@ -34,7 +34,7 @@ end
 
 y = dates(1, time(:,1));
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -94,4 +94,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/utilities/convert/dseries2M.m b/src/utilities/convert/dseries2M.m
index dfae0a4..27bb1db 100644
--- a/src/utilities/convert/dseries2M.m
+++ b/src/utilities/convert/dseries2M.m
@@ -1,4 +1,4 @@
-function ts = dseries2M(ds, method) % --*-- Unitary tests --*--
+function ts = dseries2M(ds, method)
 
 % INPUTS
 % - ds       [dseries]    daily frequency object with n elements.
@@ -11,7 +11,7 @@ function ts = dseries2M(ds, method) % --*-- Unitary tests --*--
 % OUTPUTS
 % - ts       [dseries]    monthly frequency object with m<n elements.
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-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
@@ -80,7 +80,7 @@ end
 
 ts = dseries(tdata, dM, ds.name, ds.tex);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/utilities/convert/dseries2Q.m b/src/utilities/convert/dseries2Q.m
index 5754ca3..d34aec3 100644
--- a/src/utilities/convert/dseries2Q.m
+++ b/src/utilities/convert/dseries2Q.m
@@ -1,4 +1,4 @@
-function ts = dseries2Q(ds, method)  % --*-- Unitary tests --*--
+function ts = dseries2Q(ds, method)
 
 % INPUTS
 % - ds       [dseries]    daily or monthly frequency object with n elements.
@@ -11,7 +11,7 @@ function ts = dseries2Q(ds, method)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - ts       [dseries]    quaterly frequency object with m<n elements.
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-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
@@ -97,7 +97,7 @@ end
 
 ts = dseries(tdata, dQ, ds.name, ds.tex);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -233,4 +233,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:8
\ No newline at end of file
+%@eof:8
diff --git a/src/utilities/convert/dseries2S.m b/src/utilities/convert/dseries2S.m
index 9513f75..7b54166 100644
--- a/src/utilities/convert/dseries2S.m
+++ b/src/utilities/convert/dseries2S.m
@@ -1,4 +1,4 @@
-function ts = dseries2S(ds, method)  % --*-- Unitary tests --*--
+function ts = dseries2S(ds, method)
 
 % INPUTS
 % - ds       [dseries]    daily, monthly or quarterly frequency object with n elements.
@@ -11,7 +11,7 @@ function ts = dseries2S(ds, method)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - ts       [dseries]    bi-annual frequency object with m<n elements.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -112,7 +112,7 @@ end
 
 ts = dseries(tdata, dH, ds.name, ds.tex);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -316,4 +316,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:12
\ No newline at end of file
+%@eof:12
diff --git a/src/utilities/convert/dseries2Y.m b/src/utilities/convert/dseries2Y.m
index 059f8ff..7ee6b4f 100644
--- a/src/utilities/convert/dseries2Y.m
+++ b/src/utilities/convert/dseries2Y.m
@@ -1,4 +1,4 @@
-function ts = dseries2Y(ds, method)  % --*-- Unitary tests --*--
+function ts = dseries2Y(ds, method)
 
 % INPUTS
 % - ds       [dseries]    daily, monthly, quarterly, bi-annual or yearly frequency object with n elements.
@@ -11,7 +11,7 @@ function ts = dseries2Y(ds, method)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - ts       [dseries]    annual frequency object with m<n elements.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-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
@@ -123,7 +123,7 @@ end
 
 ts = dseries(tdata, dY, ds.name, ds.tex);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -327,4 +327,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:12
\ No newline at end of file
+%@eof:12
diff --git a/src/utilities/convert/freq2string.m b/src/utilities/convert/freq2string.m
index 249abde..505fb47 100644
--- a/src/utilities/convert/freq2string.m
+++ b/src/utilities/convert/freq2string.m
@@ -1,4 +1,4 @@
-function s = freq2string(freq) % --*-- Unitary tests --*--
+function s = freq2string(freq)
 
 % INPUTS
 % - freq  [integer]   scalar equal to 1, 2, 4, 12, 52, or 365 (resp. annual, bi-annual, quaterly, monthly, weekly, or daily)
@@ -6,7 +6,7 @@ function s = freq2string(freq) % --*-- Unitary tests --*--
 % OUTPUTS
 % - s     [char]      scalar equal to Y, S, Q, M, W, or D (resp. annual, bi-annual, quaterly, monthly, weekly, or daily)
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -40,7 +40,7 @@ switch freq
     error('dates::freq2string: Unknown frequency!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/utilities/convert/string2date.m b/src/utilities/convert/string2date.m
index c3ce5de..a39aff9 100644
--- a/src/utilities/convert/string2date.m
+++ b/src/utilities/convert/string2date.m
@@ -1,6 +1,6 @@
-function date = string2date(a) % --*-- Unitary tests --*--
+function date = string2date(a)
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -64,7 +64,7 @@ if isdaily(a)
     return
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define some dates
@@ -123,4 +123,4 @@ t(14)= isequal(d8.freq,e_freq_8);
 t(15)= isequal(d9.freq,e_freq_9);
 t(16)= isequal(d9.time(1),e_date_9(1));
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/convert/string2freq.m b/src/utilities/convert/string2freq.m
index 1df5f96..66d7780 100644
--- a/src/utilities/convert/string2freq.m
+++ b/src/utilities/convert/string2freq.m
@@ -1,4 +1,4 @@
-function freq = string2freq(s) % --*-- Unitary tests --*--
+function freq = string2freq(s)
 
 % INPUTS
 % - s        [char]     scalar equal to Y, H, Q, M, or D (resp. annual, bi-annual, quaterly, monthly, or daily)
@@ -6,7 +6,7 @@ function freq = string2freq(s) % --*-- Unitary tests --*--
 % OUTPUTS
 % - freq     [integer]  scalar equal to 1, 2, 4, 12, or 365 (resp. annual, bi-annual, quaterly, monthly, or daily)
 
-% Copyright © 2013-2020 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -38,7 +38,7 @@ switch upper(s)
     error('dates::freq2string: Unknown frequency!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/utilities/from/from.m b/src/utilities/from/from.m
index 6a6a761..03e00dc 100644
--- a/src/utilities/from/from.m
+++ b/src/utilities/from/from.m
@@ -1,6 +1,6 @@
-function from(varargin)   % --*-- Unitary tests --*--
+function from(varargin)
 
-% Copyright © 2014-2018 Dynare Team
+% Copyright © 2014-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -423,7 +423,7 @@ if ~isempty(idv)
 end
 i = 0;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
diff --git a/src/utilities/get/get_series_start_date.m b/src/utilities/get/get_series_start_date.m
index 31bfd67..d05d832 100644
--- a/src/utilities/get/get_series_start_date.m
+++ b/src/utilities/get/get_series_start_date.m
@@ -1,4 +1,4 @@
-function series_start_date = get_series_start_date(frequency, original_period)   % --*-- Unitary tests --*--
+function series_start_date = get_series_start_date(frequency, original_period)
 % Given cell array obtained using from the mdbnomics library,
 % it returns a cell array of metadata ot be appended to a dseries object.
 %
@@ -9,7 +9,7 @@ function series_start_date = get_series_start_date(frequency, original_period)
 % OUTPUTS
 % - series_start_date [string]
 
-% Copyright © 2020-2021 Dynare Team
+% Copyright © 2020-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
@@ -40,7 +40,7 @@ switch frequency
         error('mdbnomics2dseries::get_series_start_date: The frequency of the dataset is currently unsupported!');
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -115,4 +115,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/utilities/insert/insert_column_vector_in_a_matrix.m b/src/utilities/insert/insert_column_vector_in_a_matrix.m
index 068241c..bbb4903 100644
--- a/src/utilities/insert/insert_column_vector_in_a_matrix.m
+++ b/src/utilities/insert/insert_column_vector_in_a_matrix.m
@@ -1,4 +1,4 @@
-function matrix = insert_column_vector_in_a_matrix(matrix, vector, i)  % --*-- Unitary tests --*--
+function matrix = insert_column_vector_in_a_matrix(matrix, vector, i)
 
 % Inserts a vector in a matrix in a column specified by i.
 %
@@ -10,7 +10,7 @@ function matrix = insert_column_vector_in_a_matrix(matrix, vector, i)  % --*-- U
 % OUTPUTS
 % - matrix [float]   n*(m+1) matrix.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -52,7 +52,7 @@ switch i
     matrix = [matrix(:,1:i-1), vector, matrix(:,i:m)];
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 % Define a datasets.
@@ -73,4 +73,4 @@ if t(1)
   t(4) = dassert(A2,[A(:,1), b, A(:,2:end)],1e-15);
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/insert/insert_object_in_a_one_dimensional_cell_array.m b/src/utilities/insert/insert_object_in_a_one_dimensional_cell_array.m
index 5ec08f4..f82dcf4 100644
--- a/src/utilities/insert/insert_object_in_a_one_dimensional_cell_array.m
+++ b/src/utilities/insert/insert_object_in_a_one_dimensional_cell_array.m
@@ -1,4 +1,4 @@
-function one_dimensional_cell_array = insert_object_in_a_one_dimensional_cell_array(one_dimensional_cell_array, object, i) % --*-- Unitary tests --*--
+function one_dimensional_cell_array = insert_object_in_a_one_dimensional_cell_array(one_dimensional_cell_array, object, i)
 
 % Inserts an object in a one dimensional cell array in a position specified by i.
 %
@@ -10,7 +10,7 @@ function one_dimensional_cell_array = insert_object_in_a_one_dimensional_cell_ar
 % OUTPUTS
 % - one_dimensional_cell_array [any]     cell array with n+1 elements.
 
-% Copyright © 2013-2021 Dynare Team
+% 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
@@ -60,7 +60,7 @@ if nc>nr
     one_dimensional_cell_array = transpose(one_dimensional_cell_array);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 A = {'A1'; 'A2'; 'A3'}; b = [1, pi];
@@ -80,4 +80,4 @@ if t(1)
   t(4) = dassert(C2,[A(1); b; A(2:3)]);
 end
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isbiannual.m b/src/utilities/is/isbiannual.m
index ec206f9..3fa2b5d 100644
--- a/src/utilities/is/isbiannual.m
+++ b/src/utilities/is/isbiannual.m
@@ -1,4 +1,4 @@
-function b = isbiannual(str)  % --*-- Unitary tests --*--
+function b = isbiannual(str)
 
 % Tests if the input can be interpreted as a bi-annual date (semester).
 %
@@ -8,7 +8,7 @@ function b = isbiannual(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       integer scalar, equal to 1 if str can be interpreted as a bi-annual date or 0 otherwise.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950H2';
@@ -58,4 +58,4 @@ t(7) = dassert(isbiannual(date_7),false);
 t(8) = dassert(isbiannual(date_8),true);
 t(9) = dassert(isbiannual(date_9),true);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isdaily.m b/src/utilities/is/isdaily.m
index 9670dcd..894be69 100644
--- a/src/utilities/is/isdaily.m
+++ b/src/utilities/is/isdaily.m
@@ -1,4 +1,4 @@
-function b = isdaily(str)  % --*-- Unitary tests --*--
+function b = isdaily(str)
 
 % Tests if the input can be interpreted as a daily date.
 %
@@ -8,7 +8,7 @@ function b = isdaily(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - b       [logical]  scalar, equal to true iff str can be interpreted as a daily date.
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -41,7 +41,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t(1) = isdaily('1950-03-29');
@@ -56,4 +56,4 @@ t(9) = ~isdaily('2000-12-00');
 t(10) = ~isdaily('1999-02-29');
 t(11) = isdaily('2000-02-29');
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isdate.m b/src/utilities/is/isdate.m
index 41ba5e9..78e967e 100644
--- a/src/utilities/is/isdate.m
+++ b/src/utilities/is/isdate.m
@@ -1,4 +1,4 @@
-function b = isdate(str)  % --*-- Unitary tests --*--
+function b = isdate(str)
 
 % Tests if the input string can be interpreted as a date.
 %
@@ -8,7 +8,7 @@ function b = isdate(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - b       [logical]   scalar equal to true iff str can be interpreted as a date.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -32,7 +32,7 @@ end
 
 b = isstringdate(str);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = 1950;
@@ -59,4 +59,4 @@ t(9) = isdate(date_9);
 t(10) = isdate(date_10);
 t(11) = ~isdate(date_11);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isdates.m b/src/utilities/is/isdates.m
index e781ddd..27887ca 100644
--- a/src/utilities/is/isdates.m
+++ b/src/utilities/is/isdates.m
@@ -1,8 +1,8 @@
-function B = isdates(A) % --*-- Unitary tests --*--
+function B = isdates(A)
 
 % Tests if the input A is a dates object.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -21,7 +21,7 @@ function B = isdates(A) % --*-- Unitary tests --*--
 
 B = isa(A,'dates');
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -111,4 +111,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/utilities/is/isfreq.m b/src/utilities/is/isfreq.m
index af117ad..9cfe50c 100644
--- a/src/utilities/is/isfreq.m
+++ b/src/utilities/is/isfreq.m
@@ -1,4 +1,4 @@
-function B = isfreq(A) % --*-- Unitary tests --*--
+function B = isfreq(A)
 
 % Tests if A can be interpreted as a frequency.
 %
@@ -8,7 +8,7 @@ function B = isfreq(A) % --*-- Unitary tests --*--
 % OUTPUTS
 % - B     [logical]        scalar equal to true iff A can be interpreted as a frequency.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -38,7 +38,7 @@ if isnumeric(A) && isequal(length(A),1) && ismember(A,[1 2 4 12 52 365])
     B = true;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -128,4 +128,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/utilities/is/isleapyear.m b/src/utilities/is/isleapyear.m
index ae3a984..98ea804 100644
--- a/src/utilities/is/isleapyear.m
+++ b/src/utilities/is/isleapyear.m
@@ -1,4 +1,4 @@
-function b = isleapyear(y) % --*-- Unitary tests --*--
+function b = isleapyear(y)
 
 % Returns true iff y is a leap year.
 %
@@ -9,7 +9,7 @@ function b = isleapyear(y) % --*-- Unitary tests --*--
 % - b     [logical]     scalar or vector, equal to true iff y is a leap year.
 
 
-% Copyright © 2015-2021 Dynare Team
+% Copyright © 2015-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -37,7 +37,7 @@ else
     error('Input must be an integer scalar or vector of integers.')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 LeapYears = [1904; 	1908; 	1912; 	1916; 	1920; 	1924; 	1928; 	1932; 	1936;
@@ -130,4 +130,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/is/islongyear.m b/src/utilities/is/islongyear.m
index 22de544..2d2c2c7 100644
--- a/src/utilities/is/islongyear.m
+++ b/src/utilities/is/islongyear.m
@@ -1,4 +1,4 @@
-function b = islongyear(y) % --*-- Unitary tests --*--
+function b = islongyear(y)
 
 % Returns true iff y is a year with 53 weeks.
 %
@@ -9,7 +9,7 @@ function b = islongyear(y) % --*-- Unitary tests --*--
 % - b     [logical]     scalar or vector, equal to true iff y is a leap year.
 
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -28,7 +28,7 @@ function b = islongyear(y) % --*-- Unitary tests --*--
 
 b = mod(y+fix(y/4)-fix(y/100)+fix(y/400),7)==4 | mod(y-1+fix((y-1)/4)-fix((y-1)/100)+fix((y-1)/400),7)==3;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 Years = [4; 9 ; 15 ; 20 ; 26 ; 32 ; 37 ; 43 ; 48 ; 54 ; 60 ; 65 ; 71 ; 76 ; 82 ; 88 ; 93 ; 99 ;
@@ -69,4 +69,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/is/ismonthly.m b/src/utilities/is/ismonthly.m
index 9e3c1f5..fd1ade7 100644
--- a/src/utilities/is/ismonthly.m
+++ b/src/utilities/is/ismonthly.m
@@ -1,4 +1,4 @@
-function b = ismonthly(str)  % --*-- Unitary tests --*--
+function b = ismonthly(str)
 
 % Tests if the input can be interpreted as a monthly date.
 %
@@ -8,7 +8,7 @@ function b = ismonthly(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       integer scalar, equal to 1 if str can be interpreted as a monthly date or 0 otherwise.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950M2';
@@ -56,4 +56,4 @@ t(6) = dassert(ismonthly(date_6),false);
 t(7) = dassert(ismonthly(date_7),false);
 t(8) = dassert(ismonthly(date_8),false);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isquarterly.m b/src/utilities/is/isquarterly.m
index fc46cb5..4b4e3cb 100644
--- a/src/utilities/is/isquarterly.m
+++ b/src/utilities/is/isquarterly.m
@@ -1,4 +1,4 @@
-function b = isquarterly(str)  % --*-- Unitary tests --*--
+function b = isquarterly(str)
 
 % Tests if the input can be interpreted as a quarterly date.
 %
@@ -8,7 +8,7 @@ function b = isquarterly(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       integer scalar, equal to 1 if str can be interpreted as a quarterly date or 0 otherwise.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950Q2';
@@ -54,4 +54,4 @@ t(5) = dassert(isquarterly(date_5),false);
 t(6) = dassert(isquarterly(date_6),false);
 t(7) = dassert(isquarterly(date_7),false);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isstringdate.m b/src/utilities/is/isstringdate.m
index 53b46d9..4f21c91 100644
--- a/src/utilities/is/isstringdate.m
+++ b/src/utilities/is/isstringdate.m
@@ -1,4 +1,4 @@
-function b = isstringdate(str)  % --*-- Unitary tests --*--
+function b = isstringdate(str)
 
 % Tests if the input string can be interpreted as a date.
 %
@@ -8,7 +8,7 @@ function b = isstringdate(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       integer scalar, equal to 1 if str can be interpreted as a date or 0 otherwise.
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -31,7 +31,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950M2';
@@ -77,4 +77,4 @@ t(19) = ~isstringdate(date_19);
 t(20) = isstringdate(date_20);
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/issubperiod.m b/src/utilities/is/issubperiod.m
index ec38a97..9a85e65 100644
--- a/src/utilities/is/issubperiod.m
+++ b/src/utilities/is/issubperiod.m
@@ -1,6 +1,6 @@
-function C = issubperiod(A,B) % --*-- Unitary tests --*--
+function C = issubperiod(A,B)
 
-% Copyright © 2013-2021 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -27,7 +27,7 @@ else
     error('issubperiod:: Second input argument must be equal to 1, 2, 4, 12 or 52 (frequency)!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -87,4 +87,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:4
\ No newline at end of file
+%@eof:4
diff --git a/src/utilities/is/isvalidday.m b/src/utilities/is/isvalidday.m
index 717642b..8c702dc 100644
--- a/src/utilities/is/isvalidday.m
+++ b/src/utilities/is/isvalidday.m
@@ -1,4 +1,4 @@
-function b = isvalidday(day)  % --*-- Unitary tests --*--
+function b = isvalidday(day)
 
 % Tests if the input can be interpreted as a daily date.
 %
@@ -8,7 +8,7 @@ function b = isvalidday(day)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - b       [logical]  scalar, equal to true iff day can be interpreted as a daily date.
 
-% Copyright © 2020 Dynare Team
+% Copyright © 2020-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -61,7 +61,7 @@ end
 
 b = true;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t(1) = isvalidday([2000,1,1]) && isvalidday([2000,1,31]) && ~isvalidday([2000, 1, 32]);
@@ -81,4 +81,4 @@ t(14) = ~isvalidday([2000, 1, 0]);
 t(15) = ~isvalidday([2000, 0, 1]);
 t(16) = isvalidday([0, 1, 1]);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isweekly.m b/src/utilities/is/isweekly.m
index 75f31de..c9da49e 100644
--- a/src/utilities/is/isweekly.m
+++ b/src/utilities/is/isweekly.m
@@ -1,4 +1,4 @@
-function b = isweekly(str)  % --*-- Unitary tests --*--
+function b = isweekly(str)
 
 % Tests if the input can be interpreted as a monthly date.
 %
@@ -8,7 +8,7 @@ function b = isweekly(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       logical scalar, equal to true if str can be interpreted as a weekly date or false otherwise.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -46,7 +46,7 @@ if b
 end
 
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950M2';
@@ -71,4 +71,4 @@ t(8) = ~isweekly(date_8);
 t(9) = ~isweekly(date_9);
 t(10) = isweekly(date_10);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/is/isyearly.m b/src/utilities/is/isyearly.m
index 812a1ef..d253c1e 100644
--- a/src/utilities/is/isyearly.m
+++ b/src/utilities/is/isyearly.m
@@ -1,4 +1,4 @@
-function b = isyearly(str)  % --*-- Unitary tests --*--
+function b = isyearly(str)
 
 % Tests if the input can be interpreted as a yearly date.
 %
@@ -8,7 +8,7 @@ function b = isyearly(str)  % --*-- Unitary tests --*--
 % OUTPUTS
 %  o b       integer scalar, equal to 1 if str can be interpreted as a yearly date or 0 otherwise.
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -35,7 +35,7 @@ else
     b = false;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 date_1 = '1950M2';
@@ -56,4 +56,4 @@ t(6) = dassert(isyearly(date_6),true);
 t(7) = dassert(isyearly(date_7),true);
 t(8) = dassert(isyearly(date_8),false);
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/misc/dummy.m b/src/utilities/misc/dummy.m
index ea642a5..0eb0399 100644
--- a/src/utilities/misc/dummy.m
+++ b/src/utilities/misc/dummy.m
@@ -1,4 +1,4 @@
-function o = dummy(iperiod, len, subperiod, name)  % --*-- Unitary tests --*--
+function o = dummy(iperiod, len, subperiod, name)
 
 % Returns a dummy in a dseries object.
 %
@@ -11,7 +11,7 @@ function o = dummy(iperiod, len, subperiod, name)  % --*-- Unitary tests --*--
 % OUTPUTS
 % - o           [dseries]
 
-% Copyright © 2022 Dynare Team
+% Copyright © 2022-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -53,7 +53,7 @@ end
 
 o = dseries(d, iperiod, name);
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -98,4 +98,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/utilities/misc/firstdayofyear.m b/src/utilities/misc/firstdayofyear.m
index 6eff5b2..3bda947 100644
--- a/src/utilities/misc/firstdayofyear.m
+++ b/src/utilities/misc/firstdayofyear.m
@@ -1,4 +1,4 @@
-function daycode = firstdayofyear(y) % --*-- Unitary tests --*--
+function daycode = firstdayofyear(y)
 
 % Returns the first day of year y.
 %
@@ -19,7 +19,7 @@ function daycode = firstdayofyear(y) % --*-- Unitary tests --*--
 %   6 -> Saturday
 %   7 -> Sunday
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -42,7 +42,7 @@ daycode = mod(1+5*mod(y-1,4)+4*mod(y-1,100)+6*mod(y-1,400), 7);
 id = (daycode==0);
 daycode(id) = 7;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -72,4 +72,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/misc/firstweekofyear.m b/src/utilities/misc/firstweekofyear.m
index 3455144..fe7a1b8 100644
--- a/src/utilities/misc/firstweekofyear.m
+++ b/src/utilities/misc/firstweekofyear.m
@@ -1,16 +1,16 @@
-function [d, m, y] = firstweekofyear(y) % --*-- Unitary tests --*--
+function [d, m, y] = firstweekofyear(y)
 
 % Returns the month and day the first week of the year.
 %
-% INPUTS 
+% INPUTS
 % - y   [integer] the year.
 %
-% OUTPUS 
+% OUTPUS
 % - y   [integer] the year of the first iso-8601 week of year y (may be equal to y-1)
 % - m   [integer] the month number (1 or 12)
 % - d   [integer] the day of the month.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -45,7 +45,7 @@ if any(id)
     d(id) = 9-d(id);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = false(2,1);
@@ -84,8 +84,8 @@ if t(1)
     t(7) = m(3)==12;
     t(8) = y(1)==1972;
     t(9) = y(2)==1973;
-    t(10) = y(3)==1973;    
+    t(10) = y(3)==1973;
 end
 
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/misc/weeksince.m b/src/utilities/misc/weeksince.m
index 830d347..1723bde 100644
--- a/src/utilities/misc/weeksince.m
+++ b/src/utilities/misc/weeksince.m
@@ -1,16 +1,16 @@
-function o = weeksince(y, w, z) % --*-- Unitary tests --*--
+function o = weeksince(y, w, z)
 
 % Returns the number of ISO weeks since base year z.
 %
-% INPUTS 
+% INPUTS
 % - y   [integer] scalar or vector, the year.
 % - w   [integer] scalar or vector, the week.
 % - z   [integer] scalar, base year.
 %
-% OUTPUTS 
+% OUTPUTS
 % - o   [integer] scalar or vector, number of weeks.
 
-% Copyright © 2021 Dynare Team
+% Copyright © 2021-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -41,7 +41,7 @@ end
 
 o = 52*(y-z)+n+w;
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -116,4 +116,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/utilities/missing/isint/isint.m b/src/utilities/missing/isint/isint.m
index 7b550cc..40a0d2f 100644
--- a/src/utilities/missing/isint/isint.m
+++ b/src/utilities/missing/isint/isint.m
@@ -1,4 +1,4 @@
-function [l,c,d] = isint(a) % --*-- Unitary tests --*--
+function [l,c,d] = isint(a)
 
 %  This function tests if the input argument is an integer.
 %
@@ -13,7 +13,7 @@ function [l,c,d] = isint(a) % --*-- Unitary tests --*--
 %  REMARKS
 %  p+q is equal to the product of m by n.
 
-% Copyright © 2009-2021 Dynare Team
+% Copyright © 2009-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
@@ -44,7 +44,7 @@ if nargout>1
     d = find(l==false);
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 a = 1938;
@@ -128,4 +128,4 @@ if t(1)
 end
 
 T = all(t);
-%@eof:5
\ No newline at end of file
+%@eof:5
diff --git a/src/utilities/missing/shiftS/shiftS.m b/src/utilities/missing/shiftS/shiftS.m
index 2f8fd1d..349355f 100644
--- a/src/utilities/missing/shiftS/shiftS.m
+++ b/src/utilities/missing/shiftS/shiftS.m
@@ -1,8 +1,8 @@
-function S = shiftS(S,n) % --*-- Unitary tests --*--
+function S = shiftS(S,n)
 
 % Removes the first n elements of a one dimensional cell array.
 
-% Copyright © 2009-2019 Dynare Team
+% Copyright © 2009-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
@@ -23,7 +23,7 @@ else
     S = {};
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 Cell = {'1', '2', '3'};
diff --git a/src/utilities/op/add_periods_to_array_of_dates.m b/src/utilities/op/add_periods_to_array_of_dates.m
index ae43ad2..eba7cd9 100644
--- a/src/utilities/op/add_periods_to_array_of_dates.m
+++ b/src/utilities/op/add_periods_to_array_of_dates.m
@@ -1,8 +1,8 @@
-function time = add_periods_to_array_of_dates(time, freq, p)  % --*-- Unitary tests --*--
+function time = add_periods_to_array_of_dates(time, freq, p)
 
 % Adds a p periods (p can be negative) to a date (or a set of dates) characterized by array time and frequency freq.
 
-% Copyright © 2013-2020 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -43,7 +43,7 @@ if ~isempty(id2)
     time(id2,2) = time(id2,2) + freq;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t(1) = isequal(add_periods_to_array_of_dates([1950 1], 4, 1),[1950 2]);
@@ -87,4 +87,4 @@ t(6) = isequal(add_periods_to_array_of_dates([1950 1], 4, -6),[1948 3]);
 t(7) = isequal(add_periods_to_array_of_dates([1950 1], 4, -7),[1948 2]);
 t(8) = isequal(add_periods_to_array_of_dates([1950 1], 4, -8),[1948 1]);
 T = all(t);
-%@eof:6
\ No newline at end of file
+%@eof:6
diff --git a/src/utilities/op/add_periods_to_date.m b/src/utilities/op/add_periods_to_date.m
index 57887d5..f549ec8 100644
--- a/src/utilities/op/add_periods_to_date.m
+++ b/src/utilities/op/add_periods_to_date.m
@@ -1,8 +1,8 @@
-function time = add_periods_to_date(time, freq, p)  % --*-- Unitary tests --*--
+function time = add_periods_to_date(time, freq, p)
 
 % Adds a p periods (p can be negative) to a date (or a set of dates) characterized by array time and frequency freq.
 
-% Copyright © 2013-2017 Dynare Team
+% Copyright © 2013-2023 Dynare Team
 %
 % This file is part of Dynare.
 %
@@ -36,7 +36,7 @@ if time(2)<1
     time(2) = time(2) + freq;
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t(1) = dassert(add_periods_to_date([1950 1], 4, 1),[1950 2]);
@@ -60,4 +60,4 @@ t(6) = dassert(add_periods_to_date([1950 1], 4, -6),[1948 3]);
 t(7) = dassert(add_periods_to_date([1950 1], 4, -7),[1948 2]);
 t(8) = dassert(add_periods_to_date([1950 1], 4, -8),[1948 1]);
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/str/name2tex.m b/src/utilities/str/name2tex.m
index 2fd60c0..be4316a 100644
--- a/src/utilities/str/name2tex.m
+++ b/src/utilities/str/name2tex.m
@@ -1,4 +1,4 @@
-function tex = name2tex(name, info) % --*-- Unitary tests --*--
+function tex = name2tex(name, info)
 
 % Converts plain text name into tex name.
 %
@@ -11,7 +11,7 @@ function tex = name2tex(name, info) % --*-- Unitary tests --*--
 % OUTPUTS
 % - tex  [string or cell of strings]
 
-% Copyright © 2012-2021 Dynare Team
+% Copyright © 2012-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
@@ -80,7 +80,7 @@ if info
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(16,1);
diff --git a/src/utilities/x13/checkcommandcompatibility.m b/src/utilities/x13/checkcommandcompatibility.m
index 9c5aa1a..fe395d3 100644
--- a/src/utilities/x13/checkcommandcompatibility.m
+++ b/src/utilities/x13/checkcommandcompatibility.m
@@ -1,8 +1,8 @@
-function checkcommandcompatibility(o, comm) % --*-- Unitary tests --*--
+function checkcommandcompatibility(o, comm)
 
 % Checks for compatibility of X13 commands.
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -39,7 +39,7 @@ switch comm
   otherwise
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(2,1);
@@ -63,4 +63,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/x13/checkoptioncompatibility.m b/src/utilities/x13/checkoptioncompatibility.m
index 155e510..3f0f9c2 100644
--- a/src/utilities/x13/checkoptioncompatibility.m
+++ b/src/utilities/x13/checkoptioncompatibility.m
@@ -1,8 +1,8 @@
-function checkoptioncompatibility(o) % --*-- Unitary tests --*--
+function checkoptioncompatibility(o)
 
 % Checks for compatibility of options in different X13 commands.
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -37,7 +37,7 @@ if ~isempty(o.estimate.file)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 t = zeros(3,1);
@@ -68,4 +68,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/x13/isoption.m b/src/utilities/x13/isoption.m
index 4568db8..79637c7 100644
--- a/src/utilities/x13/isoption.m
+++ b/src/utilities/x13/isoption.m
@@ -1,6 +1,6 @@
-function b = isoption(command, option) % --*-- Unitary tests --*--
+function b = isoption(command, option)
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -79,7 +79,7 @@ switch command
     error('x13:isoption: Unknown block!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -97,4 +97,4 @@ b = isoption('arima','ar');
 c = ~isoption('arima','unicorn');
 t = [b;c];
 T = all(t);
-%@eof:2
\ No newline at end of file
+%@eof:2
diff --git a/src/utilities/x13/printoption.m b/src/utilities/x13/printoption.m
index e8f4055..af5760b 100644
--- a/src/utilities/x13/printoption.m
+++ b/src/utilities/x13/printoption.m
@@ -1,6 +1,6 @@
-function printoption(fid, optname, optvalue) % --*-- Unitary tests --*--
+function printoption(fid, optname, optvalue)
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -43,7 +43,7 @@ elseif isreal(optvalue)
     end
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 fid = fopen('test.spc', 'w');
@@ -60,4 +60,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/x13/printspan.m b/src/utilities/x13/printspan.m
index 8c5a794..d5ffabf 100644
--- a/src/utilities/x13/printspan.m
+++ b/src/utilities/x13/printspan.m
@@ -1,6 +1,6 @@
-function printspan(fid, period1, period2) % --*-- Unitary tests --*--
+function printspan(fid, period1, period2)
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -39,7 +39,7 @@ switch period1.freq
     error('x13:regression: This is a bug! Please contact the authors.')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -81,4 +81,4 @@ catch
 end
 
 T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
diff --git a/src/utilities/x13/printstart.m b/src/utilities/x13/printstart.m
index 017c954..e94261a 100644
--- a/src/utilities/x13/printstart.m
+++ b/src/utilities/x13/printstart.m
@@ -1,4 +1,4 @@
-function printstart(fid, period) % --*-- Unitary tests --*--
+function printstart(fid, period)
 
 % Copyright © 2017-2021 Dynare Team
 %
@@ -31,7 +31,7 @@ switch period.freq
     error('x13:regression: This frequency is not available in X-13 class.')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -44,4 +44,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/src/utilities/x13/setdefaultmember.m b/src/utilities/x13/setdefaultmember.m
index 07915d4..1258913 100644
--- a/src/utilities/x13/setdefaultmember.m
+++ b/src/utilities/x13/setdefaultmember.m
@@ -1,8 +1,8 @@
-function s = setdefaultmember(name) % --*-- Unitary tests --*--
+function s = setdefaultmember(name)
 
 % Sets members of X13 object to default values (empty).
 
-% Copyright © 2017-2021 Dynare Team
+% Copyright © 2017-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
@@ -83,7 +83,7 @@ switch name
     error('x13:setdefaultmember: Unknown member!')
 end
 
-return
+return % --*-- Unit tests --*--
 
 %@test:1
 try
@@ -94,4 +94,4 @@ catch
 end
 
 T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1
diff --git a/tests/runalltests.m b/tests/runalltests.m
index d29ab02..5231a87 100644
--- a/tests/runalltests.m
+++ b/tests/runalltests.m
@@ -1,6 +1,6 @@
 function runalltests()
 
-% Copyright © 2015-2022 Dynare Team
+% Copyright © 2015-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
@@ -50,7 +50,7 @@ if isoctave()
     addpath([unit_tests_root 'fake']);
 end
 
-r = run_unitary_tests_in_directory(dseries_src_root(1:end-1));
+r = run_unit_tests_in_directory(dseries_src_root(1:end-1));
 
 delete('*.log');
 
-- 
GitLab