From octave-sources-request at bevo dot che dot wisc dot edu Thu Mar 28 02:04:51 2002
Subject: mat2str
From: Andreas Helms <helms at puppis dot astro dot physik dot uni-potsdam dot de>
To: octave-sources at bevo dot che dot wisc dot edu
Date: Thu, 28 Mar 2002 01:28:16 -0600



I enhanced the function mat2str from octave-forge a little bit so it can
format a matrix with arbitrary precision.

Andreas

 
-------------------------------------------------


## Copyright (C) 1996 John W. Eaton
##
## This program 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.
##
## This program 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 this program; see the file COPYING.  If not, write to the Free
## Software Foundation, 59 Temple Place - Suite 330, Boston, MA
## 02111-1307, USA.

## usage: mat2str (x)
##        mat2str (x,n)
## mat2str (x) format x as a string suitable for use in eval.
## mat2str(A,n) converts matrix x using n digits of precision.
## 
##
## See also: sprintf, int2str

## Author: jwe
## Modified by: Ariel Tankus, 15.6.98 .
## Modified by: Paul Kienzle, 15.7.00, to handle matrices.
## Modified by: Andreas Helms, 21.3.02, to handle arbitrary precision. 

function retval = mat2str (x,n)

  if (nargin == 1)
    [nr, nc] = size(x);
    if (nr*nc == 1)
      retval = sprintf ("%.100g", x);
    else
      retval = sprintf (" %.100g,", x.');
      retval(1) = "[";
      retval(length(retval)) = "]";
      idx = find (retval == ",");
      retval(idx(nc:nc:length(idx))) = ";";
    endif
  elseif (nargin == 2)
    [nr, nc] = size(x);
    format = ['%.' int2str(n) 'g'];
    if (nr*nc == 1)
      retval = sprintf (format, x);
    else
      retval = sprintf ([" " ,format , ","] , x.');
      retval(1) = "[";
      retval(length(retval)) = "]";
      idx = find (retval == ",");
      retval(idx(nc:nc:length(idx))) = ";";
    endif	
			
    else
    usage ("mat2str (x,n)");
  endif

endfunction

-------------------------------------------------


Lehrstuhl Astrophysik   mailto:helms at astro dot physik dot uni-potsdam dot de
Universitaet Potsdam    http://www.astro.physik.uni-potsdam.de/~helms
D-14469 Potsdam         ftp://ftp.astro.physik.uni-potsdam.de/pub/helms
Phone: +49 331 977 1556  FAX   : +49 331 977 1107