Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Replacing `find` with logical indexing
Replies: 6   Last Post: Jul 10, 2013 1:47 PM

 Messages: [ Previous | Next ]
 dpb Posts: 9,850 Registered: 6/7/07
Re: Replacing `find` with logical indexing
Posted: Jul 9, 2013 9:20 PM

On 7/9/2013 7:44 PM, Jeff wrote:
> I am trying to find the smallest nonzero magnitude element of a vector,
> D. Note that D is not integers. I am eliminating nonzero entries using
> the `find` command (first line below). The editor says to replace `find`
> with logical indexing, but replacing `D(find(D))` with `D(D)` returns an
> error. This is not slowing my code very much, but the editor's warning
> is mildly bothersome. How do I use logical indexing to fix this?
>
> nonZeroD = abs(D(find(D)));

nonZeroD = abs(D(D~-0));

> lambda = max(nonZeroD(nonZeroD == min(nonZeroD)));
> ICmode = find(D==lambda, 1, 'first');
>
> I also note that the editor doesn't complain about the second `find`.

Depends on what's done w/ the result. In the second you're using the
truncated version of returning only the first position so logical

I'm not absolutely positive I follow what you're trying to do but why
isn't it just

DNZmin=min(abs(D(D~-0))); % ?

What's the point of the max(DNZ(DNZ==min(DNZ)))? If you've only
collected terms in the vector ==min(vector), then the max is the min by
that selection process.

--

Date Subject Author
7/9/13 Jeff
7/9/13 dpb
7/9/13 Jeff
7/9/13 Derek Goring
7/9/13 dpb
7/10/13 Jeff
7/10/13 dpb