# HG changeset patch
# User jwe
# Date 1020134404 0
# Node ID 9eb96199e0f186ef64bd8b30461cf942ef1175d3
# Parent  f11cfbfc84b5f959e7f452e1b7e5a88e8c36ab3a
[project @ 2002-04-30 02:40:04 by jwe]

diff --git a/scripts/ChangeLog b/scripts/ChangeLog
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,3 +1,7 @@
+2002-04-29  Paul Kienzle <pkienzle@users.sf.net>
+
+	* general/repmat.m: New function.
+
 2002-04-25  Paul Kienzle <pkienzle@users.sf.net>
 
 	* audio/lin2mu.m: Accept matrices and [-1,1] normalized audio.
diff --git a/scripts/general/repmat.m b/scripts/general/repmat.m
new file mode 100644
--- /dev/null
+++ b/scripts/general/repmat.m
@@ -0,0 +1,53 @@
+## Copyright (C) 2000 Paul Kienzle
+##
+## This file is part of Octave.
+##
+## Octave 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 2, or (at your option)
+## any later version.
+##
+## Octave 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 Octave; see the file COPYING.  If not, write to the Free
+## Software Foundation, 59 Temple Place - Suite 330, Boston, MA
+## 02111-1307, USA.
+
+## -*- texinfo -*-
+## @deftypefn {Function File} {} repmat (@var{A}, @var{m}, @var{n})
+## @deftypefnx {Function File} {} repmat (@var{A}, [@var{m} @var{n}])
+## Form a block matrix of size @var{m} by @var{n}, with a copy of matrix
+## @var{A} as each element.  If @var{n} is not specified, form an 
+## @var{m} by @var{m} block matrix.
+## @end deftypefn
+
+## Author: Paul Kienzle <pkienzle@kienzle.powernet.co.uk>
+## Created: July 2000
+
+function x = repmat (b, m, n)
+  if (nargin < 2 || nargin > 3)
+    usage ("repmat (a, m, n)");
+  endif
+
+  if (nargin == 2)
+    if (is_scalar (m))
+      n = m;
+    elseif (is_vector (m) && length (m) == 2)
+      n = m(2);
+      m = m(1);
+    else
+      error ("repmat: only builds 2D matrices")
+    endif
+  endif
+
+  if (isstr (a))
+    x = setstr (kron (ones (m, n), toascii (a)));
+  else
+    x = kron (ones (m, n), a);
+  endif
+
+endfunction