Does anybody know a fast algorithm for this?

L

Luna Moon

In a 2D array of double(floating point) numbers,

ideally it should have an ordering on two dimensions:

for left to right, the numbers should increase,

from bottom to up, the numbers should increase,

if in the middle of some row and some column, there are valleys and
peaks, which means that the ordering is violated,

I want to detect these violations and set the peaks or valleys to some
predefined constants, say -1.

my question is:

how to do that efficiently in C/C++, also also Matlab?

Basically, the question is how to detect peaks and valleys which
violate the ordering and set them to an alert number, say -1. This is
like "highlighting" ...

Thanks!
 
G

Gianni Mariani

Luna said:
In a 2D array of double(floating point) numbers,

ideally it should have an ordering on two dimensions:

for left to right, the numbers should increase,

from bottom to up, the numbers should increase,

if in the middle of some row and some column, there are valleys and
peaks, which means that the ordering is violated,

I want to detect these violations and set the peaks or valleys to some
predefined constants, say -1.

my question is:

how to do that efficiently in C/C++, also also Matlab?

Basically, the question is how to detect peaks and valleys which
violate the ordering and set them to an alert number, say -1. This is
like "highlighting" ...


Does the 2D array fit into the primary cache ? If it does, then a
brainless for loop for each array will work fine. If however it does
not fit into cache you would benefit greatly if you tiled your algorithm
(i.e. did the analysis in tiles). Doing it in tiles allows you to
multi-thread the code so you can run analysis of multiple tiles
simultaneously.
 
K

Kai-Uwe Bux

Luna said:
In a 2D array of double(floating point) numbers,

ideally it should have an ordering on two dimensions:

for left to right, the numbers should increase,

from bottom to up, the numbers should increase,

if in the middle of some row and some column, there are valleys and
peaks, which means that the ordering is violated,

I want to detect these violations and set the peaks or valleys to some
predefined constants, say -1.

my question is:

how to do that efficiently in C/C++, also also Matlab?
[snip]

I don't think that this is really language dependent.

Note that not a single entry is in violation of the order but a pair of
adjacent (horizontally or vertically) entries (this raises the question,
where you want to put the -1).

The entries form a pattern like this:

x--x--x--x
| | | |
x--x--x--x
| | | |
x--x--x--x

You can assign values to the entries so that there is exactly one violating
edge and you can even make it so that this edge can be anywhere. For this
reasin, just verifying that no violations happen requires inspection of all
horizontal and vertical edges in the pattern. The number of edges is
approximately twice the number of entries in the matrix. Therefore, no
algorithm can have better than linear (in the number of matrix entries)
runtime in the worst case. Thus, I would go with the simple, straight
forward implementation.


Best

Kai-Uwe Bux
 

Ask a Question

Want to reply to this thread or ask your own question?

You'll need to choose a username for the site, which only take a couple of moments. After that, you can post your question and our members will help you out.

Ask a Question

Members online

No members online now.

Forum statistics

Threads
473,744
Messages
2,569,484
Members
44,903
Latest member
orderPeak8CBDGummies

Latest Threads

Top