[Obsolete] st_deviation 37/15337/3
Bruno JOFRET [Tue, 7 Oct 2014 12:33:12 +0000 (14:33 +0200)]
Change-Id: I1812ffb1994a0f62b1275d8931c2122e09477bc9

14 files changed:
scilab/CHANGES_5.5.X
scilab/modules/helptools/data/configuration/scilab_macros.txt
scilab/modules/statistics/help/en_US/central_tendency/mean.xml
scilab/modules/statistics/help/en_US/descriptive_statistics/st_deviation.xml [deleted file]
scilab/modules/statistics/help/en_US/measures_shape/cmoment.xml
scilab/modules/statistics/help/fr_FR/descriptive_statistics/median.xml
scilab/modules/statistics/help/fr_FR/descriptive_statistics/st_deviation.xml [deleted file]
scilab/modules/statistics/help/ja_JP/central_tendency/mean.xml
scilab/modules/statistics/help/ja_JP/descriptive_statistics/st_deviation.xml [deleted file]
scilab/modules/statistics/help/ja_JP/measures_shape/cmoment.xml
scilab/modules/statistics/locales/statistics.pot
scilab/modules/statistics/macros/st_deviation.sci [deleted file]
scilab/modules/umfpack/tests/unit_tests/utm300.dia.ref
scilab/modules/umfpack/tests/unit_tests/utm300.tst

index 812bab8..f0c1434 100644 (file)
@@ -28,6 +28,8 @@ Removed Functions
 
 * regress removed. Please use reglin instead.
 
+* st_deviation removed. Please use stdev instead.
+
 
 Compilation
 ===========
index 8485544..8455e0a 100644 (file)
@@ -88,7 +88,7 @@ y=mean(A,'m')
                 <link linkend="median">median</link>
             </member>
             <member>
-                <link linkend="st_deviation">st_deviation</link>
+                <link linkend="stdev">stdev</link>
             </member>
         </simplelist>
     </refsection>
diff --git a/scilab/modules/statistics/help/en_US/descriptive_statistics/st_deviation.xml b/scilab/modules/statistics/help/en_US/descriptive_statistics/st_deviation.xml
deleted file mode 100644 (file)
index e573d29..0000000
+++ /dev/null
@@ -1,117 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<!--
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2013 - Scilab Enterprises - Paul Bignier: made obsolete, merged in stdev
- * Copyright (C) 2000 - INRIA -
- *
- * This file must be used under the terms of the CeCILL.
- * This source file is licensed as described in the file COPYING, which
- * you should have received as part of this distribution.  The terms
- * are also available at
- * http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
- *
- -->
-<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" xmlns:scilab="http://www.scilab.org" xml:lang="en" xml:id="st_deviation">
-    <refnamediv>
-        <refname>st_deviation</refname>
-        <refpurpose>standard deviation (row or
-            column-wise) of vector/matrix entries
-            <emphasis role="bold">This function is obsolete.</emphasis>
-        </refpurpose>
-    </refnamediv>
-    <refnamediv xml:id="stdev">
-        <refname>stdev</refname>
-        <refpurpose>standard deviation (row or column-wise) of vector/matrix entries</refpurpose>
-    </refnamediv>
-    <refsynopsisdiv>
-        <title>Calling Sequence</title>
-        <synopsis>y=st_deviation(x)
-            y=st_deviation(x,'r')
-            y=st_deviation(x,'c')
-            y=stdev(x)
-            y=stdev(x,'r')
-            y=stdev(x,'c')
-        </synopsis>
-    </refsynopsisdiv>
-    <refsection>
-        <title>Arguments</title>
-        <variablelist>
-            <varlistentry>
-                <term>x</term>
-                <listitem>
-                    <para>real vector or matrix</para>
-                </listitem>
-            </varlistentry>
-            <varlistentry>
-                <term>y</term>
-                <listitem>
-                    <para>scalar or vector</para>
-                </listitem>
-            </varlistentry>
-        </variablelist>
-    </refsection>
-    <refsection>
-        <title>Description</title>
-        <warning>
-            This function is obsolete.
-            It is better to use <link linkend="stdev">stdev</link> instead.
-        </warning>
-        <para>
-            st_deviation computes the "sample" standard deviation, that
-            is, it is normalized by N-1, where N is the sequence length.
-        </para>
-        <para>
-            For a vector or a matrix <literal>x</literal>, <literal>y=st_deviation(x)</literal>  returns in the
-            scalar <literal>y</literal> the standard deviation of all the entries of <literal>x</literal>.
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'r')</literal> (or, equivalently,
-            <literal>y=st_deviation(x,1)</literal>) is the rowwise standard deviation. It returns in each
-            entry of the column vector <literal>y</literal> the standard deviation of each row of <literal>x</literal>.
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'c')</literal> (or, equivalently, <literal>y=st_deviation(x,2)</literal>) is the columnwise st_deviation. It returns in each entry of the row vector
-            <literal>y</literal> the standard deviation of each column  of <literal>x</literal>.
-        </para>
-    </refsection>
-    <refsection>
-        <title>Examples</title>
-        <programlisting role="example"><![CDATA[
-A=[1,2,10;7,7.1,7.01];
-st_deviation(A)
-st_deviation(A,'r')
-st_deviation(A,'c')
- ]]></programlisting>
-    </refsection>
-    <refsection role="see also">
-        <title>See Also</title>
-        <simplelist type="inline">
-            <member>
-                <link linkend="sum">sum</link>
-            </member>
-            <member>
-                <link linkend="median">median</link>
-            </member>
-            <member>
-                <link linkend="mean">mean</link>
-            </member>
-            <member>
-                <link linkend="nanstdev">nanstdev</link>
-            </member>
-            <member>
-                <link linkend="stdevf">stdevf</link>
-            </member>
-        </simplelist>
-    </refsection>
-    <refsection>
-        <title>History</title>
-        <revhistory>
-            <revision>
-                <revnumber>5.5.0</revnumber>
-                <revdescription>
-                    Function tagged as obsolete. Will be removed in 5.5.1. Please use <link linkend="stdev">stdev</link> instead.
-                </revdescription>
-            </revision>
-        </revhistory>
-    </refsection>
-</refentry>
index 55f4996..e73e087 100644 (file)
@@ -79,7 +79,7 @@ mom=cmoment(x,3,'c')
                 <link linkend="median">median</link>
             </member>
             <member>
-                <link linkend="st_deviation">st_deviation</link>
+                <link linkend="stdev">stdev</link>
             </member>
             <member>
                 <link linkend="mean">mean</link>
index e93634f..a70b317 100644 (file)
@@ -106,7 +106,7 @@ median(A,'m')
                 <link linkend="mean">mean</link>
             </member>
             <member>
-                <link linkend="st_deviation">st_deviation</link>
+                <link linkend="stdev">stdev</link>
             </member>
         </simplelist>
     </refsection>
diff --git a/scilab/modules/statistics/help/fr_FR/descriptive_statistics/st_deviation.xml b/scilab/modules/statistics/help/fr_FR/descriptive_statistics/st_deviation.xml
deleted file mode 100644 (file)
index faada8f..0000000
+++ /dev/null
@@ -1,115 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<!--
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2000 - INRIA
- *
- * This file must be used under the terms of the CeCILL.
- * This source file is licensed as described in the file COPYING, which
- * you should have received as part of this distribution.  The terms
- * are also available at
- * http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
- *
- -->
-<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" xmlns:scilab="http://www.scilab.org" xml:lang="fr" xml:id="st_deviation">
-    <refnamediv>
-        <refname>st_deviation</refname>
-        <refpurpose>
-            écart-type des termes d'une matrice
-            <emphasis role="bold">Cette fonction est obsolète.</emphasis>
-        </refpurpose>
-    </refnamediv>
-    <refnamediv xml:id="stdev">
-        <refname>stdev</refname>
-        <refpurpose>écart-type des termes d'une matrice  </refpurpose>
-    </refnamediv>
-    <refsynopsisdiv>
-        <title>Séquence d'appel</title>
-        <synopsis>y=st_deviation(x)
-            y=st_deviation(x,'r')
-            y=st_deviation(x,'c')
-            y=stdev(x)
-            y=stdev(x,'r')
-            y=stdev(x,'c')
-        </synopsis>
-    </refsynopsisdiv>
-    <refsection>
-        <title>Paramètres</title>
-        <variablelist>
-            <varlistentry>
-                <term>x  </term>
-                <listitem>
-                    <para>vecteur ou matrice réelle
-                    </para>
-                </listitem>
-            </varlistentry>
-            <varlistentry>
-                <term>y  </term>
-                <listitem>
-                    <para>scalaire ou vecteur
-                    </para>
-                </listitem>
-            </varlistentry>
-        </variablelist>
-    </refsection>
-    <refsection>
-        <title>Description</title>
-        <warning>
-            Cette fonction est <emphasis role="bold">obsolète</emphasis>,
-            utilisez <link linkend="stdev">stdev</link> à la place.
-        </warning>
-        <para>
-            st_deviation calcule l'estimateur sans biais de l'écart-type (normalisé par N-1, ou N est le nombre d'échantillons).
-        </para>
-        <para>
-            Pour un vecteur ou une matrice <literal>x</literal>, <literal>y=st_deviation(x)</literal> renvoie dans <literal>y</literal> l'écart-type des termes de <literal>x</literal>.
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'r')</literal> (ou <literal>y=st_deviation(x,1)</literal>) est l'écart-type suivant l'indice de ligne. Chaque composante du vecteur ligne <literal>y</literal> contient l'écart-type de chaque colonne de <literal>x</literal>.
-        </para>
-        <para>
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'c')</literal> (ou <literal>y=st_deviation(x,2)</literal>) est l'écart-type selon l'indice de colonne. Chaque composante du vecteur colonne <literal>y</literal> contient l'écart-type de chaque ligne de <literal>x</literal>.
-        </para>
-    </refsection>
-    <refsection>
-        <title>Exemples</title>
-        <programlisting role="example"><![CDATA[
-A=[1,2,10;7,7.1,7.01];
-st_deviation(A)
-st_deviation(A,'r')
-st_deviation(A,'c')
- ]]></programlisting>
-    </refsection>
-    <refsection role="see also">
-        <title>Voir aussi</title>
-        <simplelist type="inline">
-            <member>
-                <link linkend="sum">sum</link>
-            </member>
-            <member>
-                <link linkend="median">median</link>
-            </member>
-            <member>
-                <link linkend="mean">mean</link>
-            </member>
-            <member>
-                <link linkend="nanstdev">nanstdev</link>
-            </member>
-            <member>
-                <link linkend="stdevf">stdevf</link>
-            </member>
-        </simplelist>
-    </refsection>
-    <refsection>
-        <title>History</title>
-        <revhistory>
-            <revision>
-                <revnumber>5.5.0</revnumber>
-                <revdescription>
-                    Fonction obsolète. Sera supprimée dans la version 5.5.1. Veuillez utiliser <link linkend="stdev">stdev</link> à la place.
-                </revdescription>
-            </revision>
-        </revhistory>
-    </refsection>
-</refentry>
index c80d039..2230a18 100644 (file)
@@ -86,7 +86,7 @@ y=mean(A,'m')
                 <link linkend="median">median</link>
             </member>
             <member>
-                <link linkend="st_deviation">st_deviation</link>
+                <link linkend="stdev">stdev</link>
             </member>
         </simplelist>
     </refsection>
diff --git a/scilab/modules/statistics/help/ja_JP/descriptive_statistics/st_deviation.xml b/scilab/modules/statistics/help/ja_JP/descriptive_statistics/st_deviation.xml
deleted file mode 100644 (file)
index e8f81f8..0000000
+++ /dev/null
@@ -1,100 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<!--
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2000 - INRIA - 
- * 
- * This file must be used under the terms of the CeCILL.
- * This source file is licensed as described in the file COPYING, which
- * you should have received as part of this distribution.  The terms
- * are also available at    
- * http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
- *
- -->
-<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" xmlns:scilab="http://www.scilab.org" xml:lang="ja" xml:id="st_deviation">
-    <refnamediv>
-        <refname>st_deviation</refname>
-        <refpurpose>ベクトル/行列の要素の(行方向または列方向の)標準偏差</refpurpose>
-    </refnamediv>
-    <refnamediv xml:id="stdev">
-        <refname>stdev</refname>
-        <refpurpose>ベクトル/行列の要素の(行方向または列方向の)標準偏差</refpurpose>
-    </refnamediv>
-    <refsynopsisdiv>
-        <title>呼び出し手順</title>
-        <synopsis>y=st_deviation(x)
-            y=st_deviation(x,'r')
-            y=st_deviation(x,'c')
-            y=stdev(x)
-            y=stdev(x,'r')
-            y=stdev(x,'c')
-        </synopsis>
-    </refsynopsisdiv>
-    <refsection>
-        <title>パラメータ</title>
-        <variablelist>
-            <varlistentry>
-                <term>x</term>
-                <listitem>
-                    <para>実数ベクトルまたは行列</para>
-                </listitem>
-            </varlistentry>
-            <varlistentry>
-                <term>y</term>
-                <listitem>
-                    <para>スカラーまたはベクトル</para>
-                </listitem>
-            </varlistentry>
-        </variablelist>
-    </refsection>
-    <refsection>
-        <title>説明</title>
-        <para>
-            st_deviation は, N-1で正規化された"標本"標準偏差を計算します.
-            ただし, Nはデータ長です.
-        </para>
-        <para>
-            <literal>x</literal>がベクトルまたは行列の場合, <literal>y=st_deviation(x)</literal>は
-            <literal>x</literal>の全ての要素の標準偏差をスカラー<literal>y</literal>として返します.
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'r')</literal> (または等価的に,
-            <literal>y=st_deviation(x,1)</literal>) は行毎の標準偏差です. 
-            列ベクトル<literal>y</literal>の各要素について<literal>x</literal>の各行の標準偏差を
-            返します.
-        </para>
-        <para>
-            <literal>y=st_deviation(x,'c')</literal> (または等価的に, <literal>y=st_deviation(x,2)</literal>)
-            は列毎の標準偏差です.行ベクトル<literal>y</literal>の各要素について<literal>x</literal>の各列の
-            標準偏差を返します.
-        </para>
-    </refsection>
-    <refsection>
-        <title>例</title>
-        <programlisting role="example"><![CDATA[ 
-A=[1,2,10;7,7.1,7.01];
-st_deviation(A)
-st_deviation(A,'r')
-st_deviation(A,'c')
- ]]></programlisting>
-    </refsection>
-    <refsection role="see also">
-        <title>参照</title>
-        <simplelist type="inline">
-            <member>
-                <link linkend="sum">sum</link>
-            </member>
-            <member>
-                <link linkend="median">median</link>
-            </member>
-            <member>
-                <link linkend="mean">mean</link>
-            </member>
-            <member>
-                <link linkend="nanstdev">nanstdev</link>
-            </member>
-            <member>
-                <link linkend="stdevf">stdevf</link>
-            </member>
-        </simplelist>
-    </refsection>
-</refentry>
index 67c4c35..87fbb19 100644 (file)
@@ -80,7 +80,7 @@ mom=cmoment(x,3,'c')
                 <link linkend="median">median</link>
             </member>
             <member>
-                <link linkend="st_deviation">st_deviation</link>
+                <link linkend="stdev">stdev</link>
             </member>
             <member>
                 <link linkend="mean">mean</link>
index bab4579..b4db32d 100644 (file)
@@ -698,7 +698,6 @@ msgstr ""
 # File: macros/nanmeanf.sci, line: 48
 # File: macros/sample.sci, line: 51
 # File: macros/samplef.sci, line: 100
-# File: macros/st_deviation.sci, line: 40
 # File: macros/variancef.sci, line: 132
 #, c-format
 msgid ""
diff --git a/scilab/modules/statistics/macros/st_deviation.sci b/scilab/modules/statistics/macros/st_deviation.sci
deleted file mode 100644 (file)
index d79c16d..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2013 - Scilab Enterprises - Paul Bignier: made obsolete, merged in stdev
-// Copyright (C) XXXX - INRIA
-// Copyright (C) XXXX - ENPC
-// Copyright (C) XXXX - Ecole des Mines de Nantes - Philippe Castagliola
-// This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at
-// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
-
-function sd=st_deviation(x,cr)
-    //
-    warnobsolete("stdev", "5.5.1");
-    if argn(2)<2 then cr="*",end
-    if x == [] then sd=%nan;return ;end
-    if typeof(x)=="hypermat" then sd=%hm_stdev(x,cr),return,end
-    [m,n]=size(x);
-    if cr=="*" then
-        n=m*n
-        select n
-        case 0 then sd=%nan
-        case 1 then sd=0
-        else
-            sd=sqrt(sum((x-mean(x)).^2)/(n-1));
-        end
-    elseif cr=="c"|cr==2
-        if n==1 then
-            sd=zeros(m,1)
-        else
-            sd=sqrt(sum((x-mean(x,"c")*ones(x(1,:))).^2,"c")/(n-1));
-        end
-    elseif cr=="r"|cr==1
-        if m==1 then
-            sd=zeros(1,n)
-        else
-            sd=sqrt(sum((x-ones(x(:,1))*mean(x,"r")).^2,"r")/(m-1));
-        end
-    else
-        error(msprintf(gettext("%s: Wrong value for input argument #%d: ''%s'', ''%s'', ''%s'', %d or %d expected.\n"),"st_deviation",2,"*","r","c",1,2));
-    end
-endfunction
index 4aeea50..c635120 100644 (file)
@@ -4,6 +4,7 @@
 //
 //  This file is distributed under the same license as the Scilab package.
 // ============================================================================
+// <-- CLI SHELL MODE -->
 // Check value provided here:
 // http://math.nist.gov/MatrixMarket/data/SPARSKIT/tokamak/utm300.html
 warning("off");
@@ -27,7 +28,7 @@ if nnz(A-A') <> 4382 then bugmes();quit;end
 offsets = [ 0 -50 50 1 -1 -5 5 2 -2 -4 ];
 nonzerosExpected = [ 300 190 186 180 166 159 158 150 134 113 ];
 for j=1:length(offsets)
-       if nnz(diag(A, offsets(j))) <> nonzerosExpected(j) then bugmes();quit;end
+    if nnz(diag(A, offsets(j))) <> nonzerosExpected(j) then bugmes();quit;end
 end
 // TODO add the rest of information
 [ij]=spget(A);
@@ -36,7 +37,7 @@ e=ij(:,1)-ij(:,2);
 m=ceil(mean(abs(e)));
 if m <> 21 then bugmes();quit;end
 // Bandwidths std.dev.
-d=ceil(st_deviation(abs(e)));
+d=ceil(stdev(abs(e)));
 if d <> 23 then bugmes();quit;end
 lastinrow=[find(diff(ij(:,1))>0) size(ij,1)];
 rowlength=diff([0 lastinrow]);
@@ -50,10 +51,10 @@ if indshortest <> 3 then bugmes();quit;end
 // Average nonzeros per row
 m=ceil(mean(rowlength));
 if m <> 11 then bugmes();quit;end
-d=round(st_deviation(rowlength)*10)/10;
+d=round(stdev(rowlength)*10)/10;
 if d <> 7.7 then bugmes();quit;end
 // Column Data
-ji=gsort(ij(:,[2 1]),'lr','i');
+ji=gsort(ij(:,[2 1]),"lr","i");
 lastincol=[find(diff(ji(:,1))>0) size(ji,1)];
 collength=diff([0 lastincol]);
 [longest,indlongest]=max(collength);
@@ -64,7 +65,7 @@ if shortest <> 1 then bugmes();quit;end
 if indshortest <> 16 then bugmes();quit;end
 m=ceil(mean(collength));
 if m <> 11 then bugmes();quit;end
-d=round(st_deviation(collength)*10)/10;
+d=round(stdev(collength)*10)/10;
 if d <> 7.1 then bugmes();quit;end
 // Profile Storage
 lastinrow=[find(diff(ij(:,1))>0) size(ij,1)];
@@ -72,9 +73,9 @@ rowbounds= [ij([1 lastinrow(1:$-1)+1],2) ij(lastinrow,2)];
 rowindex=ij(lastinrow,1);
 rowbandwidth=rowbounds-rowindex*ones(1,2);
 lowband=[-rowbandwidth(find(rowbandwidth(:,1)<=0),1)
-         -rowbandwidth(find(rowbandwidth(:,2)<=0),2)];
+-rowbandwidth(find(rowbandwidth(:,2)<=0),2)];
 uppband=[rowbandwidth(find(rowbandwidth(:,1)>=0),1)
-         rowbandwidth(find(rowbandwidth(:,2)>=0),2)];
+rowbandwidth(find(rowbandwidth(:,2)>=0),2)];
 // lower bandwidth => max
 minlow=max(lowband);
 if minlow <> 74 then bugmes();quit;end
index c570502..7b46105 100644 (file)
@@ -30,7 +30,7 @@ offsets = [ 0 -50 50 1 -1 -5 5 2 -2 -4 ];
 nonzerosExpected = [ 300 190 186 180 166 159 158 150 134 113 ];
 
 for j=1:length(offsets)
-       if nnz(diag(A, offsets(j))) <> nonzerosExpected(j) then pause, end
+    if nnz(diag(A, offsets(j))) <> nonzerosExpected(j) then pause, end
 end
 
 // TODO add the rest of information
@@ -43,7 +43,7 @@ m=ceil(mean(abs(e)));
 if m <> 21 then pause, end
 
 // Bandwidths std.dev.
-d=ceil(st_deviation(abs(e)));
+d=ceil(stdev(abs(e)));
 if d <> 23 then pause, end
 
 lastinrow=[find(diff(ij(:,1))>0) size(ij,1)];
@@ -58,15 +58,15 @@ if indlongest <> 116 then pause, end
 if shortest <> 1 then pause, end
 if indshortest <> 3 then pause, end
 
-// Average nonzeros per row 
+// Average nonzeros per row
 m=ceil(mean(rowlength));
 if m <> 11 then pause, end
 
-d=round(st_deviation(rowlength)*10)/10;
+d=round(stdev(rowlength)*10)/10;
 if d <> 7.7 then pause, end
 
 // Column Data
-ji=gsort(ij(:,[2 1]),'lr','i');
+ji=gsort(ij(:,[2 1]),"lr","i");
 lastincol=[find(diff(ji(:,1))>0) size(ji,1)];
 collength=diff([0 lastincol]);
 
@@ -81,7 +81,7 @@ if indshortest <> 16 then pause, end
 m=ceil(mean(collength));
 if m <> 11 then pause, end
 
-d=round(st_deviation(collength)*10)/10;
+d=round(stdev(collength)*10)/10;
 if d <> 7.1 then pause, end
 
 // Profile Storage
@@ -90,9 +90,9 @@ rowbounds= [ij([1 lastinrow(1:$-1)+1],2) ij(lastinrow,2)];
 rowindex=ij(lastinrow,1);
 rowbandwidth=rowbounds-rowindex*ones(1,2);
 lowband=[-rowbandwidth(find(rowbandwidth(:,1)<=0),1)
-         -rowbandwidth(find(rowbandwidth(:,2)<=0),2)];
+-rowbandwidth(find(rowbandwidth(:,2)<=0),2)];
 uppband=[rowbandwidth(find(rowbandwidth(:,1)>=0),1)
-         rowbandwidth(find(rowbandwidth(:,2)>=0),2)];   
+rowbandwidth(find(rowbandwidth(:,2)>=0),2)];
 
 // lower bandwidth => max
 minlow=max(lowband);