Mercurial > octave-antonio
annotate scripts/sparse/spfun.m @ 9245:16f53d29049f
update copyright notices
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Fri, 22 May 2009 10:46:00 -0400 |
| parents | 1bf0ce0930be |
| children | a8ce6bdecce5 |
| rev | line source |
|---|---|
| 9245 | 1 ## Copyright (C) 2004, 2005, 2006, 2007, 2008, 2009 |
| 2 ## David Bateman and Andy Adler | |
| 5164 | 3 ## |
| 7016 | 4 ## This file is part of Octave. |
| 5164 | 5 ## |
| 7016 | 6 ## Octave is free software; you can redistribute it and/or modify it |
| 7 ## under the terms of the GNU General Public License as published by | |
| 8 ## the Free Software Foundation; either version 3 of the License, or (at | |
| 9 ## your option) any later version. | |
| 10 ## | |
| 11 ## Octave is distributed in the hope that it will be useful, but | |
| 12 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 14 ## General Public License for more details. | |
| 5164 | 15 ## |
| 16 ## You should have received a copy of the GNU General Public License | |
| 7016 | 17 ## along with Octave; see the file COPYING. If not, see |
| 18 ## <http://www.gnu.org/licenses/>. | |
| 5164 | 19 |
| 20 ## -*- texinfo -*- | |
| 21 ## @deftypefn {Function File} {@var{y} =} spfun (@var{f},@var{x}) | |
| 22 ## Compute @code{f(@var{x})} for the non-zero values of @var{x}. | |
| 23 ## This results in a sparse matrix with the same structure as | |
|
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
24 ## @var{x}. The function @var{f} can be passed as a string, a |
| 5164 | 25 ## function handle or an inline function. |
| 26 ## @end deftypefn | |
| 27 | |
| 6498 | 28 function t = spfun (f, s) |
| 29 | |
| 5164 | 30 if (nargin != 2) |
| 6046 | 31 print_usage (); |
| 5164 | 32 endif |
| 33 | |
|
7505
f5005d9510f4
Remove dispatched sparse functions and treat in the generic versions of the functions
David Bateman <dbateman@free.fr>
parents:
7151
diff
changeset
|
34 [i, j, v] = find (s); |
|
f5005d9510f4
Remove dispatched sparse functions and treat in the generic versions of the functions
David Bateman <dbateman@free.fr>
parents:
7151
diff
changeset
|
35 [m, n] = size (s); |
| 5164 | 36 |
| 6220 | 37 if (isa (f, "function_handle") || isa (f, "inline function")) |
| 6498 | 38 t = sparse (i, j, f(v), m, n); |
| 5164 | 39 else |
| 6498 | 40 t = sparse(i, j, feval (f, v), m, n); |
| 5164 | 41 endif |
| 42 | |
| 43 endfunction | |
| 44 | |
| 45 %!assert(spfun('exp',[1,2;3,0]),sparse([exp(1),exp(2);exp(3),0])) | |
| 46 %!assert(spfun('exp',sparse([1,2;3,0])),sparse([exp(1),exp(2);exp(3),0])) | |
| 47 %!assert(spfun(@exp,[1,2;3,0]),sparse([exp(1),exp(2);exp(3),0])) | |
| 48 %!assert(spfun(@exp,sparse([1,2;3,0])),sparse([exp(1),exp(2);exp(3),0])) | |
| 49 |