[Octave-patch-tracker] [patch #10067] [octave forge] (linear-algebra) Improvement of funm

[Simon]

Follow-up Comment #7, patch #10067 (project octave):

I included file funm.m from the RickOne-Repo
<https://github.com/RickOne16/matrix/blob/master/funm/funm.m> into the
linear-algebra package. Some minor changes to funm.m were necessary. List of
changes: 

1) funm is renamed to schur_parlett_funm since the RickOne-Repo and the
linear-algebra package use different FUN parameters. In the RickOne-Repo
parameter FUN has to be a function handle with two parameters x and k (for
exp, sin, cos, log, sinh and cosh exceptions are possible). In the
linear-algebra package the parameter FUN is a string or a function handle with
only one parameter x.  
2) Help text of funm is converted into an octave docstring.
3) ordschur is replaced by LAPACK-Routine xTREXC (see trexc.cc). 
4) In function swapping the line  

ind = cat (2,ind,{mdone:mdone+h(i)-1}) 
 
is replaced by 

ind{i} = mdone:mdone+h (i)-1
 
since the original line causes a crash (zero indices are stored in cell array
ind). The crash can be reproduced with

A = [ 4, 2, 0; 1, 4, 1; 1, 1, 4];
C = funm(A,@exp);

5) Additional BISTs.
6) Additional references in docstring.  

(file #52564)
    _______________________________________________________

Additional Item Attachment:

File name: linear-algebra-2.2.4.tar.gz    Size:397 KB
   
<https://file.savannah.gnu.org/file/linear-algebra-2.2.4.tar.gz?file_id=52564>



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/patch/?10067>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/