[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Indirect Addressing Question
From: |
adler |
Subject: |
Re: Indirect Addressing Question |
Date: |
Thu, 15 Nov 2001 22:16:10 -0500 (EST) |
On Wed, 14 Nov 2001, Paul Kienzle wrote:
> Here is yet another solution, but it requires Andy Adler's sparse functions
> from
> octave-forge:
>
> q = [1:10];
> idx = [5,1,1,1,2,2];
> projection = sum(sparse(1:length(idx),idx,q(idx)+10));
> used = find(projection);
> q(used) = projection(used);
> q
> q =
>
> 33 24 3 4 15 6 7 8 9 10
There's a more direct way of doing this with sparse:
q = [1:10];
idx = [5,1,1,1,2,2];
[i,j,v]=spfind( sparse(1,idx,q(idx)+10,1,5,'sum') );
q(j)=v
q =
33 24 3 4 15 6 7 8 9 10
I introduced the "sum" mode on the suggestion of
Andreas Stahel <address@hidden>, as
a nice way to simplify some FEM problems.
_______________________________________
Andy Adler, address@hidden
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------
- Indirect Addressing Question, Etienne Grossmann, 2001/11/14
- Re: Indirect Addressing Question, Mirek Kwasniak, 2001/11/14
- Re: Indirect Addressing Question, Paul Kienzle, 2001/11/14
- overloading functions (Was: Indirect Addressing Question), adler, 2001/11/15
- Re: overloading functions (Was: Indirect Addressing Question), Paul Kienzle, 2001/11/16
- Re: overloading functions (Was: Indirect Addressing Question), John W. Eaton, 2001/11/16
- Re: overloading functions (Was: Indirect Addressing Question), adler, 2001/11/16
- Re: overloading functions, Paul Kienzle, 2001/11/20
- semantics of sum, prod, cumsum etc., adler, 2001/11/22
- Re: Indirect Addressing Question,
adler <=