Blitz++ Library

BZ_DECLARE_FUNCTION is elegant for allowing functions to play nicely 
with array expressions (and first/secondIndex variables, etc), but is 
there a straightforward way to do the same with methods?

Presuming the uber-magic solution is impossible, is there a short-ish 
way of writing a templated function to accomplish the same thing?

As a prototypical example:

#include <blitz/array.h>

class foo {
public:
int twotimes(int i) {return 2*i;}
};

BZ_DECLARE_MAGIC(foo::twotimes) // I wish

int threetimes(int i) {return 3*i;}
BZ_DECLARE_FUNCTION(threetimes);

int main() {
Array<int,1> myarray(10);
foo bar; // declare a foo
firstIndex ii;
myarray = threetimes(ii); // works
myarray = bar.twotimes(ii); // doesn't work at all
return 0;
}

Hi Everyone

I was wondering if uneven distributed indexing is possible ?
For example can I have something similar to this:

    TinyVector<int,1> index = 0, 3, 4, 10;
    Array<double,2> data;
    Array<double,1> slice = data(index,0);

I have unsuccessfully tried the 'find' function and also indirections!

Any input is appreciated.
Thanks in advance,
Hamid

[I am forwarding this message from someone who is having trouble posting to
the list.  -- JCC]


From: Paul Sheridan <sherida6@...>
Date: April 22, 2008 9:14:30 AM JST
To: blitz-support@...
Subject: Random Number Generators




Dear all,


I am running blitz 0.9 and g++ (GCC) 3.3.3 (SuSE Linux).


I wrote a test program 'beta.cpp' to generate a random number from the beta
distribution:


#include <iostream>#include <random/beta.h>
using namespace std;using namespace ranlib;
int main(){ Beta<double> betaGen; betaGen.setParameters( 1.0, 5.0 ); double
x = betaGen.random();
cout << x << endl; 
return 0;}


When I compile it, I get the following errors:


Beta_test> g++ beta.cpp -I/home1/usr2/psherida/blitz/include -o beta
beta.cpp: In function `int main()':beta.cpp:9: error: no matching function
for call to `ranlib::Beta<double, ranlib::defaultIRNG,
ranlib::defaultState>::Beta()'/home1/usr2/psherida/blitz/include/random/beta
.h:52: error: candidates are: ranlib::Beta<double, ranlib::defaultIRNG,
ranlib::defaultState>::Beta(const ranlib::Beta<double, ranlib::defaultIRNG,
ranlib::defaultState>&)/home1/usr2/psherida/blitz/include/random/beta.h:57:
error: ranlib::Beta<T, IRNG, stateTag>::Beta(T, T) [with T = double, IRNG =
ranlib::defaultIRNG, stateTag = ranlib::defaultState


I believe I have followed the instructions in Chapter 9 of the blitz manual
precisely. Any suggestions would be much appreciated.


Regards,


Paul Sheridan

Forwarding this reply on behalf of Patrik Jonsson (the list did not like the
original sender's address) ...
-- Julian C.
 

-----Original Message-----
From: Patrik Jonsson [mailto:patrik@...] 
Sent: Tuesday, April 01, 2008 5:28 PM
To: blitz-support@...
Cc: blitz-support@...
Subject: Re: [Blitz-support] (no subject)

Hi Ying,

Short answer is no, there is no sort function in blitz. It's possible the
standard STL sort algorithm would work, but the array iterators were not
really STL-compliant last I tried doing something with them. Worth a try
though.

I don't understand your example. Why should the output not be what you put
in? You are just assigning to an array and then printing it. (And unless you
care about the storage order, there's no need to supply a specific storage
object. You'll get C-storage order by default.)

Regards,

/Patrik

Ying Wu wrote:
> Hi,
>
> My area is Data Mining and recently, I am trying to move from matlab 
> to
> c++ platform, so I find blitz++. Of course, I am a new one using this
> excellent c++ library. 
>
> I have a problem. I try to sort a 1*5 vector like
>
>     GeneralArrayStorage<2> storage;
>     storage.ordering() = firstRank;
>     storage.base() = 0;
>     storage.ascendingFlag() = true;
>     Array<int,2> D(1,5,storage);
>     D=4,6,8,2,3;
>     cout<< "D = " <<D<<endl;
>
> The output should be 2,3,4,6,8
> But why doesn't it work?
>
> Could I ask whether there is some function like sort() in matlab?
>
> Thanks very much,
>
> Best wishes,
>
> Ying Wu
>
> **********************************************************************
> *****************************************
>
> Following the successful merger between the University of Paisley and Bell
College, the University has been renamed as University of the West of
Scotland.
>
> Please note that as part of the next stage in the University's
development, we have changed name to University of the West of Scotland. 
>
> As a result, our web and email addresses have changed and I would be
obliged if you would update your records accordingly. 
>
> New University Web address: 
>
> http://www.uws.ac.uk
>
> New University e-mail address format: 
>
> firstname.surname@...
>
> The University of the West of Scotland is a registered Scottish charity.
Charity number SC002520.
>
> **********************************************************************
> *****************************************
>
> Legal disclaimer
> --------------------------
>
> The information transmitted is the property of the University of the 
> West of Scotland and is intended only for the person or entity to which it
is addressed and may contain confidential and/or privileged material.
Statements and opinions expressed in this e-mail may not represent those of
the company.  Any review, retransmission, dissemination and other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipient is prohibited.
> If you received this in error, please contact the sender immediately and
delete the material from any computer. 
>
> --------------------------
>
> ----------------------------------------------------------------------
> --- Check out the new SourceForge.net Marketplace.
> It's the best place to buy or sell services for just about anything 
> Open Source.
> http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marke
> tplace _______________________________________________
> Blitz-support mailing list
> Blitz-support@...
> https://lists.sourceforge.net/lists/listinfo/blitz-support
>

I want to convert an NxN blitz::Array to an NxNx1 blitz::Array without
making a copy. Is there an easy way to do that? I'm a bit miffed looking at
the docs.

Thanks,

Joe

Hello,

 

I am currently new to Blitz++ (thank you for this really nice library!).

 

I am trying to initialize a large 2D Array as follow, but MinGW (g++ for windows) fails to compile.  Curiously, the initialization works if I include only the first 943 elements (from beginning to #943, excluding value 16.864300 highlighter below).

 

1)       Is it a compiler issue?

2)       Is it a programming issues?

3)       Is it a bug?

4)       How can I get around this?

 

Yan

 

See following code:

 

int main()

{

    blitz::Array<double,2> M(41,45);

    M=

    0.000000, 0.317440,     0.634880,     0.952320,     1.269760,     1.295360,     1.320960,     1.346560,       1.372160,     1.136430,     0.900696,     0.664963,     0.429230,     1.228070,     2.026910,       2.825750,     3.624590,     3.129840,     2.635100,     2.140360,     1.645620,     2.535100,       3.424590,     4.314070,     5.203560,     4.630760,     4.057960,     3.485160,     2.912360,       5.906600,     8.900840,     11.895100,    14.889300,    11.376800,    7.864350,     4.351870,       0.839392,     8.067820,     15.296300,    22.524700,    29.753100,    22.314800,    14.876600,       7.438280,     0.000000,

    0.000000, 0.306013,     0.612026,     0.918039,     1.224050,     1.297180,     1.370310,     1.443440,       1.516570,     1.255360,     0.994140,     0.732923,     0.471707,     1.198290,     1.924860,       2.651440,     3.378020,     3.022150,     2.666290,     2.310420,     1.954550,     2.790080,       3.625600,     4.461130,     5.296650,     4.704710,     4.112770,     3.520830,     2.928890,       5.832990,     8.737080,     11.641200,    14.545300,    11.185000,    7.824790,     4.464550,       1.104310,     8.363060,     15.621800,    22.880600,    30.139300,    22.604500,    15.069700,       7.534830,     0.000000,

    0.000000, 0.293976,     0.587952,     0.881927,     1.175900,     1.299620,     1.423330,     1.547040,       1.670760,     1.378390,     1.086020,     0.793648,     0.501279,     1.159300,     1.817320,       2.475340,     3.133360,     2.925040,     2.716710,     2.508390,     2.300070,     3.059300,       3.818540,     4.577770,     5.337010,     4.751170,     4.165320,     3.579480,     2.993640,       5.783320,     8.572990,     11.362700,    14.152300,    10.961100,    7.769910,     4.578680,       1.387460,     8.698280,     16.009100,    23.319900,    30.630700,    22.973000,    15.315400,       7.657680,     0.000000,

    0.000000, 0.280718,     0.561436,     0.842155,     1.122870,     1.303280,     1.483680,     1.664080,       1.844480,     1.509620,     1.174760,     0.839902,     0.505042,     1.101910,     1.698770,       2.295640,     2.892510,     2.849070,     2.805630,     2.762200,     2.718760,     3.357040,       3.995320,     4.633600,     5.271880,     4.742610,     4.213350,     3.684080,     3.154810,       5.781540,     8.408270,     11.035000,    13.661700,    10.673100,    7.684400,     4.695730,       1.707070,     9.113440,     16.519800,    23.926200,    31.332600,    23.499400,    15.666300,       7.833140,     0.000000,

    0.000000, 0.265630,     0.531260,     0.796890,     1.062520,     1.308770,     1.555020,     1.801270,       2.047520,     1.653160,     1.258810,     0.864449,     0.470092,     1.016910,     1.563730,       2.110550,     2.657360,     2.804820,     2.952290,     3.099750,     3.247210,     3.697540,       4.147870,     4.598200,     5.048520,     4.651550,     4.254580,     3.857600,     3.460630,       5.851610,     8.242600,     10.633600,    13.024600,    10.288800,    7.552960,     4.817160,       2.081350,     9.648520,     17.215700,    24.782900,    32.350000,    24.262500,    16.175000,       8.087510,     0.000000,

    0.000000, 0.248676,     0.497352,     0.746029,     0.994705,     1.316740,     1.638770,     1.960800,       2.282830,     1.809830,     1.336830,     0.863835,     0.390836,     0.901859,     1.412880,       1.923900,     2.434930,     2.801040,     3.167160,     3.533280,     3.899390,     4.084720,       4.270040,     4.455360,     4.640680,     4.463960,     4.287230,     4.110500,     3.933780,       6.006250,     8.078720,     10.151200,    12.223700,    9.790650,     7.357640,     4.924630,       2.491620,     10.315000,    18.138500,    25.961900,    33.785300,    25.339000,    16.892700,       8.446330,     0.000000,

    0.000000, 0.232123,     0.464247,     0.696370,     0.928494,     1.327900,     1.727300,     2.126710,       2.526120,     1.967320,     1.408520,     0.849726,     0.290929,     0.781337,     1.271740,       1.762150,     2.252560,     2.839140,     3.425710,     4.012290,     4.598870,     4.481160,       4.363450,     4.245740,     4.128030,     4.219720,     4.311410,     4.403110,     4.494800,       6.213230,     7.931670,     9.650100,     11.368500,    9.219270,     7.070000,     4.920740,       2.771470,     11.010700,    19.250000,    27.489300,    35.728500,    26.796400,    17.864300,       8.932140,     0.000000,

    0.000000, 0.218814,     0.437628,     0.656442,     0.875255,     1.343010,     1.810760,     2.278510,       2.746270,     2.110030,     1.473800,     0.837571,     0.201339,     0.686684,     1.172030,       1.657370,     2.142720,     2.918680,     3.694650,     4.470610,     5.246580,     4.839110,       4.431640,     4.024170,     3.616690,     3.972200,     4.327710,     4.683220,     5.038730,       6.429140,     7.819550,     9.209960,     10.600400,    8.629670,     6.658970,     4.688270,       2.717580,     11.604900,    20.492200,    29.379500,    38.266800,    28.700100,    19.133400,       9.566700,     0.000000,

    0.000000, 0.211589,     0.423178,     0.634768,     0.846357,     1.362810,     1.879260,     2.395720,       2.912170,     2.222390,     1.532600,     0.842818,     0.153033,     0.649241,     1.145450,       1.641660,     2.137860,     3.039270,     3.940670,     4.842070,     5.743470,     5.110810,       4.478160,     3.845500,     3.212840,     3.774780,     4.336710,     4.898650,     5.460580,       6.610530,     7.760480,     8.910420,     10.060400,    8.076930,     6.093490,     4.110040,       2.126600,     11.966800,    21.806900,    31.647100,    41.487200,    31.115400,    20.743600,       10.371800,    0.000000,

    0.000000, 0.211180,     0.422361,     0.633541,     0.844722,     1.387770,     1.930830,     2.473880,       3.016930,     2.301360,     1.585780,     0.870199,     0.154620,     0.682164,     1.209710,       1.737250,     2.264800,     3.202800,     4.140800,     5.078800,     6.016800,     5.260430,       4.504060,     3.747690,     2.991320,     3.665850,     4.340390,     5.014920,     5.689450,       6.728360,     7.767280,     8.806190,     9.845100,     7.601610,     5.358110,     3.114610,       0.871111,     12.001000,    23.130900,    34.260800,    45.390800,    34.043100,    22.695400,       11.347700,    0.000000,

    0.000000, 0.209874,     0.419747,     0.629621,     0.839494,     1.417240,     1.994990,     2.572740,       3.150490,     2.394180,     1.637880,     0.881579,     0.125276,     0.725883,     1.326490,       1.927100,     2.527710,     3.420540,     4.313380,     5.206220,     6.099050,     5.299740,       4.500430,     3.701110,     2.901800,     3.624030,     4.346250,     5.068480,     5.790700,       6.811190,     7.831690,     8.852180,     9.872670,     7.186310,     4.499940,     1.813570,     -0.872794,     11.753900,    24.380600,    37.007300,    49.634000,    37.225500,    24.817000,    12.408500,       0.000000,

    0.000000, 0.197844,     0.395687,     0.593531,     0.791375,     1.450280,     2.109180,     2.768080,       3.426980,     2.560690,     1.694400,     0.828108,     -0.038182,    0.702645,     1.443470,       2.184300,     2.925120,     3.706100,     4.487070,     5.268040,     6.049010,     5.252420,       4.455830,     3.659250,     2.862660,     3.612940,     4.363210,     5.113490,     5.863770,       6.901970,     7.940170,     8.978370,     10.016600,    6.799150,     3.581730,     0.364302,     -2.853120,     11.307000,    25.467100,    39.627200,    53.787300,    40.340500,    26.893700,    13.446800,       0.000000,

    0.000000, 0.165266,     0.330531,     0.495797,     0.661062,     1.485930,     2.310810,     3.135680,       3.960550,     2.860680,     1.760810,     0.660935,     -0.438937,    0.534695,     1.508330,       2.481960,     3.455590,     4.073060,     4.690520,     5.307990,     5.925450,     5.142160,       4.358870,     3.575580,     2.792290,     3.596230,     4.400180,     5.204130,     6.008080,       7.043630,     8.079180,     9.114730,     10.150300,    6.408250,     2.666200,     -1.075840,    -4.817880,     10.741900,    26.301700,    41.861400,    57.421200,    43.065900,    28.710600,    14.355300,       0.000000,

    0.000000, 0.108004,     0.216008,     0.324012,     0.432016,     1.524760,     2.617510,     3.710250,       4.803000,     3.320240,     1.837480,     0.354714,     -1.128050,    0.180402,     1.488850,       2.797300,     4.105750,     4.523580,     4.941400,     5.359220,     5.777050,     4.992980,       4.208910,     3.424840,     2.640770,     3.553410,     4.466050,     5.378680,     6.291320,       7.261980,     8.232640,     9.203310,     10.174000,    5.990450,     1.806920,     -2.376600,    -6.560130,     10.130500,    26.821200,    43.511900,    60.202500,    45.151900,    30.101300,    15.050600,       0.000000,

    0.000000, 0.044682,     0.089364,     0.134046,     0.178728,     1.573230,     2.967730,     4.362230,       5.756740,     3.830530,     1.904320,     -0.021884,    -1.948090,    -0.257385,    1.433320,       3.124030,     4.814730,     5.014040,     5.213340,     5.412650,     5.611950,     4.830250,       4.048550,     3.266840,     2.485140,     3.527400,     4.569660,     5.611920,     6.654180,       7.514240,     8.374290,     9.234350,     10.094400,    5.557520,     1.020640,     -3.516240,    -8.053120,     9.506040,     27.065200,    44.624400,    62.183500,    46.637600,    31.091800,    15.545900,       0.000000,

    0.000000, -0.000388,    -0.000776,    -0.001164,    -0.001552,    1.639290,     3.280130,     4.920970,       6.561810,     4.248980,     1.936150,     -0.376684,    -2.689520,    -0.639695,    1.410130,       3.459950,     5.509770,     5.489380,     5.468990,     5.448590,     5.428200,     4.679680,       3.931160,     3.182640,     2.434110,     3.576990,     4.719870,     5.862740,     7.005620,       7.740470,     8.475310,     9.210160,     9.945010,     5.129960,     0.314917,     -4.500130,    -9.315180,     8.891830,     27.098800,    45.305800,    63.512800,    47.634600,    31.756400,    15.878200,       0.000000,

    0.000000, -0.002893,    -0.005786,    -0.008680,    -0.011573,    1.730890,     3.473360,     5.215830,       6.958290,     4.433020,     1.907760,     -0.617509,    -3.142780,    -0.827560,    1.487660,       3.802870,     6.118090,     5.894520,     5.670960,     5.447390,     5.223830,     4.566980,       3.910130,     3.253280,     2.596440,     3.760970,     4.925500,     6.090040,     7.254570,       7.880740,     8.506910,     9.133080,     9.759240,     4.728280,     -0.302680,    -5.333640,    -10.364600,    8.311340,     26.987300,    45.663200,    64.339200,    48.254400,    32.169600,    16.084800,       0.000000,

    0.000000, 0.053393,     0.106787,     0.160180,     0.213573,     1.852190,     3.490800,     5.129410,       6.768020,     4.286710,     1.805400,     -0.675906,    -3.157210,    -0.722269,    1.712680,       4.147620,     6.582570,     6.190020,     5.797470,     5.404920,     5.012370,     4.515060,       4.017750,     3.520440,     3.023130,     4.102520,     5.181900,     6.261280,     7.340660,       7.896280,     8.451910,     9.007530,     9.563150,     4.366950,     -0.829251,    -6.025450,    -11.221700,    7.793860,     26.809400,    45.824900,    64.840400,    48.630300,    32.420200,    16.210100,       0.000000,

    0.000000, 0.152356,     0.304713,     0.457069,     0.609426,     1.992080,     3.374730,     4.757390,       6.140040,     3.900590,     1.661140,     -0.578310,    -2.817760,    -0.386161,    2.045440,       4.477040,     6.908630,     6.398840,     5.889040,     5.379240,     4.869450,     4.535680,       4.201920,     3.868150,     3.534380,     4.482350,     5.430320,     6.378280,     7.326250,       7.832990,     8.339740,     8.846480,     9.353230,     4.036340,     -1.280550,    -6.597440,    -11.914300,    7.392050,     26.698400,    46.004800,    65.311200,    48.983400,    32.655600,    16.327800,       0.000000,

    0.000000, 0.269795,     0.539591,     0.809386,     1.079180,     2.135680,     3.192180,     4.248670,       5.305170,     3.411820,     1.518470,     -0.374876,    -2.268230,    0.078168,     2.424560,       4.770960,     7.117350,     6.559560,     6.001770,     5.443980,     4.886190,     4.637810,       4.389420,     4.141040,     3.892660,     4.745580,     5.598500,     6.451430,     7.304350,       7.757920,     8.211480,     8.665040,     9.118610,     3.720800,     -1.677000,    -7.074810,    -12.472600,    7.164420,     26.801500,    46.438500,    66.075500,    49.556600,    33.037800,    16.518900,       0.000000,

    0.000000, 0.381509,     0.763018,     1.144530,     1.526040,     2.268090,     3.010150,     3.752200,       4.494250,     2.957580,     1.420910,     -0.115762,    -1.652430,    0.568119,     2.788670,       5.009220,     7.229780,     6.710760,     6.191750,     5.672740,     5.153720,     4.830400,       4.507070,     4.183740,     3.860420,     4.737310,     5.614210,     6.491110,     7.368010,       7.738110,     8.108210,     8.478310,     8.848400,     3.404690,     -2.039020,    -7.482730,    -12.926400,    7.169510,     27.265500,    47.361400,    67.457400,    50.593000,    33.728700,       16.864300,     0.000000,

    0.000000, 0.472276,     0.944552,     1.416830,     1.889100,     2.386040,     2.882980,     3.379920,       3.876860,     2.641660,     1.406460,     0.171264,     -1.063930,    1.012920,     3.089780,       5.166640,     7.243500,     6.866470,     6.489430,     6.112400,     5.735370,     5.119250,       4.503120,     3.887000,     3.270880,     4.351350,     5.431820,     6.512290,     7.592770,       7.826840,     8.060910,     8.294980,     8.529050,     3.077070,     -2.374910,    -7.826900,    -13.278900,    7.437710,     28.154300,    48.870900,    69.587500,    52.190600,    34.793700,    17.396900,       0.000000,

    0.000000, 0.562790,     1.125580,     1.688370,     2.251160,     2.532760,     2.814350,     3.095950,       3.377550,     2.434360,     1.491170,     0.547978,     -0.395212,    1.469140,     3.333480,       5.197830,     7.062170,     6.942420,     6.822670,     6.702920,     6.583170,     5.497450,       4.411720,     3.325990,     2.240270,     3.676270,     5.112270,     6.548280,     7.984280,       8.022270,     8.060270,     8.098270,     8.136260,     2.745840,     -2.644580,    -8.035010,    -13.425400,    7.886980,     29.199400,    50.511800,    71.824200,    53.868200,    35.912100,    17.956100,       0.000000,

    0.000000, 0.682724,     1.365450,     2.048170,     2.730900,     2.763080,     2.795260,     2.827440,       2.859620,     2.272600,     1.685570,     1.098540,     0.511518,     2.025130,     3.538750,       5.052360,     6.565980,     6.829830,     7.093680,     7.357540,     7.621390,     5.954930,       4.288460,     2.621990,     0.955527,     2.849330,     4.743130,     6.636920,     8.530720,       8.308810,     8.086890,     7.864970,     7.643050,     2.423620,     -2.795810,    -8.015240,    -13.234700,    8.407190,     30.049100,    51.690900,    73.332800,    54.999600,    36.666400,    18.333200,       0.000000,

    0.000000, 0.861750,     1.723500,     2.585250,     3.447000,     3.131850,     2.816700,     2.501550,       2.186400,     2.093310,     2.000220,     1.907130,     1.814040,     2.769300,     3.724560,       4.679820,     5.635080,     6.419880,     7.204670,     7.989470,     8.774260,     6.481610,       4.188950,     1.896290,     -0.396370,    2.007790,     4.411950,     6.816110,     9.220270,       8.670820,     8.121360,     7.571900,     7.022450,     2.123030,     -2.776390,    -7.675800,    -12.575200,    8.888190,     30.351600,    51.815000,    73.278400,    54.958800,    36.639200,    18.319600,       0.000000,

    0.000000, 1.111920,     2.223850,     3.335770,     4.447700,     3.659620,     2.871550,     2.083470,       1.295400,     1.869420,     2.443450,     3.017470,     3.591500,     3.747720,     3.903940,       4.060170,     4.216390,     5.651030,     7.085680,     8.520320,     9.954960,     7.050470,       4.145980,     1.241490,     -1.663000,    1.262420,     4.187850,     7.113270,     10.038700,       9.096120,     8.153540,     7.210960,     6.268390,     1.848520,     -2.571350,    -6.991210,    -11.411100,    9.250290,     29.911700,    50.573000,    71.234400,    53.425800,    35.617200,    17.808600,       0.000000,

    0.000000, 1.374830,     2.749660,     4.124480,     5.499310,     4.229750,     2.960180,     1.690620,       0.421057,     1.717920,     3.014790,     4.311650,     5.608520,     4.837310,     4.066100,       3.294890,     2.523680,     4.650850,     6.778020,     8.905190,     11.032400,    7.566750,       4.101150,     0.635542,     -2.830060,    0.617976,     4.066020,     7.514060,     10.962100,       9.586200,     8.210300,     6.834400,     5.458500,     1.571870,     -2.314750,    -6.201380,    -10.088000,    9.535590,     29.159200,    48.782800,    68.406400,    51.304800,    34.203200,    17.101600,       0.000000,

    0.000000, 1.574430,     3.148860,     4.723280,     6.297710,     4.691290,     3.084870,     1.478460,       -0.127963,    1.791800,     3.711560,     5.631320,     7.551080,     5.872680,     4.194270,       2.515860,     0.837459,     3.594160,     6.350860,     9.107560,     11.864300,    7.918740,       3.973230,     0.027715,     -3.917800,    0.052698,     4.023200,     7.993690,     11.964200,       10.146000,    8.327750,     6.509530,     4.691320,     1.256700,     -2.177910,    -5.612530,    -9.047140,     9.816660,     28.680500,    47.544300,    66.408100,    49.806000,    33.204000,    16.602000,       0.000000,

    0.000000, 1.634690,     3.269380,     4.904060,     6.538750,     4.893310,     3.247870,     1.602440,       -0.043004,    2.244040,     4.531080,     6.818130,     9.105170,     6.688440,     4.271700,       1.854960,     -0.561771,    2.655790,     5.873340,     9.090900,     12.308500,    7.994730,       3.681010,     -0.632715,    -4.946440,    -0.455162,    4.036120,     8.527390,     13.018700,       10.780400,    8.542030,     6.303710,     4.065390,     0.866624,     -2.332140,    -5.530900,    -8.729660,     10.166100,    29.061800,    47.957500,    66.853300,    50.139900,    33.426600,    16.713300,       0.000000,

    0.000000, 1.517420,     3.034840,     4.552260,     6.069680,     4.793880,     3.518090,     2.242290,       0.966495,     3.203280,     5.440070,     7.676860,     9.913650,     7.082750,     4.251860,       1.420970,     -1.409920,    2.009760,     5.429430,     8.849110,     12.268800,    7.733820,       3.198860,     -1.336110,    -5.871070,    -0.897135,    4.076800,     9.050740,     14.024700,       11.442900,    8.861170,     6.279410,     3.697650,     0.416972,     -2.863710,    -6.144380,    -9.425060,     10.632900,    30.690900,    50.748900,    70.806800,    53.105100,    35.403400,    17.701700,       0.000000,

    0.000000, 1.335820,     2.671640,     4.007460,     5.343280,     4.787130,     4.230980,     3.674840,       3.118690,     4.700720,     6.282760,     7.864790,     9.446820,     6.707980,     3.969140,       1.230290,     -1.508550,    1.826880,     5.162310,     8.497740,     11.833200,    7.278330,       2.723490,     -1.831360,    -6.386200,    -1.143880,    4.098440,     9.340760,     14.583100,       11.882000,    9.180860,     6.479750,     3.778640,     0.129947,     -3.518750,    -7.167450,    -10.816100,    11.172500,    33.161100,    55.149700,    77.138300,    57.853700,    38.569100,    19.284600,       0.000000,

    0.000000, 1.240930,     2.481850,     3.722780,     4.963710,     5.376190,     5.788680,     6.201160,       6.613640,     6.743200,     6.872760,     7.002320,     7.131890,     5.180010,     3.228120,       1.276240,     -0.675637,    2.277160,     5.229950,     8.182750,     11.135500,    6.821370,       2.507200,     -1.806980,    -6.121150,    -1.035820,    4.049500,     9.134830,     14.220200,       11.794400,    9.368740,     6.943020,     4.517310,     0.279468,     -3.958380,    -8.196220,    -12.434100,    11.716500,    35.867100,    60.017600,    84.168200,    63.126200,    42.084100,    21.042100,       0.000000,

    0.000000, 1.383790,     2.767570,     4.151360,     5.535150,     7.064210,     8.593280,     10.122300,       11.651400,    9.337560,     7.023720,     4.709870,     2.396030,     2.114730,     1.833430,       1.552140,     1.270840,     3.530590,     5.790350,     8.050100,     10.309900,    6.556070,       2.802290,     -0.951496,    -4.705280,    -0.413408,    3.878460,     8.170340,     12.462200,       10.877300,    9.292400,     7.707490,     6.122590,     1.139450,     -3.843680,    -8.826820,    -13.810000,    12.196800,    38.203600,    64.210400,    90.217200,    67.662900,    45.108600,    22.554300,       0.000000,

    0.000000, 1.837020,     3.674030,     5.511050,     7.348070,     10.013000,    12.678000,    15.343000,       18.007900,    12.300100,    6.592180,     0.884315,     -4.823550,    -2.543430,    -0.263303,       2.016820,     4.296950,     5.590880,     6.884800,     8.178730,     9.472660,     6.617050,       3.761440,     0.905823,     -1.949790,    0.803099,     3.555990,     6.308880,     9.061760,       8.978980,     8.896190,     8.813410,     8.730620,     2.878710,     -2.973210,    -8.825120,    -14.677000,    12.586900,    39.850800,    67.114700,    94.378600,    70.784000,    47.189300,    23.594700,       0.000000,

    0.000000, 2.359550,     4.719090,     7.078640,     9.438180,     13.019300,    16.600400,    20.181500,       23.762600,    14.684500,    5.606500,     -3.471560,    -12.549600,   -7.535980,    -2.522360,       2.491270,     7.504900,     7.796410,     8.087920,     8.379440,     8.670950,     6.904890,       5.138820,     3.372750,     1.606690,     2.374100,     3.141510,     3.908920,     4.676330,       6.553840,     8.431350,     10.308900,    12.186400,    5.245560,     -1.695260,    -8.636080,    -15.576900,    13.026400,    41.629800,    70.233100,    98.836500,    74.127300,    49.418200,    24.709100,       0.000000,

    0.000000, 2.631880,     5.263760,     7.895640,     10.527500,    14.538300,    18.549100,    22.559900,       26.570700,    15.354300,    4.137930,     -7.078470,    -18.294900,   -11.275900,   -4.256970,       2.761990,     9.780940,     9.319290,     8.857640,     8.395980,     7.934330,     7.261660,       6.588990,     5.916320,     5.243650,     3.980150,     2.716650,     1.453150,     0.189652,       4.207750,     8.225840,     12.243900,    16.262000,    7.883230,     -0.495569,    -8.874370,    -17.253200,    13.697000,    44.647100,    75.597300,    106.547000,   79.910600,    53.273700,    26.636900,       0.000000,

    0.000000, 2.334520,     4.669050,     7.003570,     9.338090,     13.025500,    16.712900,    20.400300,       24.087600,    13.172700,    2.257770,     -8.657170,    -19.572100,   -12.176200,   -4.780380,       2.615470,     10.011300,    9.331600,     8.651860,     7.972130,     7.292400,     7.529450,       7.766510,     8.003560,     8.240610,     5.301830,     2.363050,     -0.575736,    -3.514520,       2.546550,     8.607620,     14.668700,    20.729800,    10.434900,    0.140117,     -10.154700,   -20.449500,    14.779900,    50.009400,    85.238900,    120.468000,   90.351300,    60.234200,    30.117100,       0.000000,

    0.000000, 1.238280,     2.476560,     3.714840,     4.953120,     7.384230,     9.815330,     12.246400,       14.677500,    7.358560,     0.039587,     -7.279390,    -14.598400,   -9.098870,    -3.599390,       1.900100,     7.399590,     7.240850,     7.082120,     6.923380,     6.764640,     7.586940,       8.409250,     9.231550,     10.053800,    6.098100,     2.142350,     -1.813410,    -5.769160,       2.021350,     9.811850,     17.602400,    25.392900,    12.621800,    -0.149243,    -12.920300,   -25.691400,    16.401800,    58.494900,    100.588000,   142.681000,   107.011000,   71.340600,    35.670300,       0.000000,

    0.000000, -0.524849,    -1.049700,    -1.574550,    -2.099400,    -1.689430,    -1.279460,    -0.869493,       -0.459526,    -1.446680,    -2.433840,    -3.420990,    -4.408150,    -2.701530,    -0.994921,       0.711691,     2.418300,     3.396250,     4.374210,     5.352160,     6.330110,     7.459260,       8.588410,     9.717560,     10.846700,    6.441510,     2.036300,     -2.368900,    -6.774110,       2.464210,     11.702500,    20.940800,    30.179200,    14.476600,    -1.225910,    -16.928400,   -32.631000,    18.469000,    69.568900,    120.669000,   171.769000,   128.827000,   85.884400,    42.942200,       0.000000,

    0.000000, -2.732570,    -5.465140,    -8.197710,    -10.930300,   -13.051400,   -15.172400,   -17.293500,    -19.414600,   -12.246000,   -5.077420,    2.091150,     9.259730,     5.909120,     2.558520,       -0.792082,    -4.142680,    -1.617580,    0.907516,     3.432620,     5.957720,     7.208120,       8.458520,     9.708920,     10.959300,    6.483010,     2.006700,     -2.469620,    -6.945930,       3.552450,     14.050800,    24.549200,    35.047600,    16.110100,    -2.827450,    -21.765000,   -40.702500,    20.833100,    82.368700,    143.904000,   205.440000,   154.080000,   102.720000,   51.360000,       0.000000,

    0.000000, -5.162590,    -10.325200,   -15.487800,   -20.650400,   -25.557400,   -30.464500,   -35.371500,    -40.278600,   -24.042300,   -7.806080,    8.430170,     24.666400,    15.626400,    6.586450,       -2.453530,    -11.493500,   -7.216040,    -2.938560,    1.338920,     5.616400,     6.895250,       8.174100,     9.452960,     10.731800,    6.373560,     2.015310,     -2.342940,    -6.701190,       4.963370,     16.627900,    28.292500,    39.957000,    17.632800,    -4.691420,    -27.015700,   -49.339900,    23.345700,    96.031300,    168.717000,   241.402000,   181.052000,   120.701000,   60.350600,       0.000000;

 

    return 1;

}

 

 

 

 

 

Tel./Tél. : (613) 990-5019 | Fax/Télécopieur : (613) 952-7136

yan.bombardier@...

 

 

Structures and Materials Performance Laboratory

Laboratoire de performance des structures et des matériaux

Institute for Aerospace Research

Institut de recherche aérospatiale

National Research Council Canada

Conseil national de recherches Canada

M-14, 1200 Montreal Road, Ottawa, ON K1A 0R6

M-14, 1200 chemin Montréal, Ottawa (Ont.) K1A 0R6

Government of Canada

Gouvernement du Canada

 http://www.nrcaerospace.com

 http://www.cnrcaerospatiale.com

Using STL in conjunction with MTL is the solution I've used.  Below is my
code for Blitz vectors - not elegant, but it works.

/**********************************************************************
Sort.h

Used for sorting elements of Blitz++ vector, computing the sortIndex
of a Blitz++ vector, and computing the unique elements of a Blitz++ vector.
*********************************************************************** */

#ifndef BLITZ_SORT_H
#define BLITZ_SORT_H


#include "blitz/array.h"
using namespace blitz;

#include <vector>
#include "mtl/mtl.h"
#include "mtl/dense1D.h"
using namespace mtl;

#include "nr.h"


namespace bul {

template<typename T_numtype>
Array<int,1> sortIndex(const Array<T_numtype,1>& x);

template<typename T_numtype>
Array<T_numtype,1> sort(const Array<T_numtype,1>& x);

template<typename T_numtype>
Array<T_numtype,1> unique(const Array<T_numtype,1>& x);



// Function definitions *************************************************
template<typename T_numtype>
Array<int,1> sortIndex(const Array<T_numtype,1>& x) {

  // copy x to a MTL vector
  mtl::dense1D<double> vec_MTL(x.size());
  int cIndex = 0;
  for (int i=x.lbound(firstDim); i<=x.ubound(firstDim); i++) {
    vec_MTL[cIndex] = x(i);
    cIndex++;
  }

  // call Numerical Recipes routine
  mtl::dense1D<int> index_MTL(x.size());
  NR::indexx(vec_MTL, index_MTL);

  // copy sort index to a Blitz++ vector
  Array<int,1> xIndex(index_MTL.data(), shape(index_MTL.size()),
  neverDeleteData);

  return xIndex.copy();
}


template<typename T_numtype>
Array<T_numtype,1> sort(const Array<T_numtype,1>& x) {

  // copy x to a STL vector
  std::vector<T_numtype> vec_STL(x.size());
  int i, cIndex = 0;
  for (i=x.lbound(firstDim); i<=x.ubound(firstDim); i++) {
    vec_STL[cIndex] = x(i);
    cIndex++;
  }

  // sort the vector
  sort(vec_STL.begin(), vec_STL.end());

  // copy STL vector to a Blitz++ vector
  Array<T_numtype,1> xSorted(vec_STL.size());
  cIndex = 0;
  for (i=xSorted.lbound(firstDim); i<=xSorted.ubound(firstDim); i++) {
    xSorted(i) = vec_STL[cIndex];
    cIndex++;
  }

  return xSorted.copy();
}


template<typename T_numtype>
Array<T_numtype,1> unique(const Array<T_numtype,1>& x) {

  // sort x and copy it to a STL vector
  // (sorting is required by std::unique())
  Array<T_numtype,1> xcopy(x.size());
  xcopy = sort(x);
  std::vector<T_numtype> vec_STL(x.size());
  int i, cIndex = 0;
  for (i=x.lbound(firstDim); i<=x.ubound(firstDim); i++) {
    vec_STL[cIndex] = xcopy(i);
    cIndex++;
  }

  // compute unique elements of STL vector
  std::vector<T_numtype>::iterator new_end;
  new_end = unique(vec_STL.begin(), vec_STL.end());
  vec_STL.erase(new_end, vec_STL.end());


  // copy STL vector to a Blitz++ vector
  Array<T_numtype,1> xUnique(vec_STL.size());
  cIndex = 0;
  for (i=xUnique.lbound(firstDim); i<=xUnique.ubound(firstDim); i++) {
    xUnique(i) = vec_STL[cIndex];
    cIndex++;
  }

  return xUnique.copy();
}

} /* namespace bul */

#endif //  BLITZ_SORT_H




Bob Dorazio
Department of Statistics
University of Florida
P.O. Box 110339
Gainesville, FL  32611-0339
Phone:  (352) 392-3081
Web site:  http://ifasstat.ifas.ufl.edu/BDorazio.html



-----blitz-support-bounces@... wrote: -----


To: blitz-support@...
From: Christopher Subich <csubich@...>
Sent by: blitz-support-bounces@...
Date: 04/03/2008 02:50AM
Subject: Re: [Blitz-support] (no subject)

Ying Wu wrote:

>     GeneralArrayStorage<2> storage;
>     storage.ordering() = firstRank;
>     storage.base() = 0;
>     storage.ascendingFlag() = true;

This is a storage object, determining the array's memory layout; it has
nothing to do with the data itself.

The ascendingFlag(), for example, just means that array(1) is stored at
a higher memory address than array(0).

The suggestion to use the STL sort is probably the correct one.  If
Blitz's iterators aren't up to the task, then a 'lower-level' sorting is
possible by treating array.data() directly as a C-array.  This is often
required, most notably when interfacing with outside libraries.

If you go that route, then pay careful attention to array.stride() and
the lower bounds by dimension (if not all 0.)

-------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services for
just about anything Open Source.
http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marketplace

_______________________________________________
Blitz-support mailing list
Blitz-support@...
https://lists.sourceforge.net/lists/listinfo/blitz-support

Hi,

My area is Data Mining and recently, I am trying to move from matlab to
c++ platform, so I find blitz++. Of course, I am a new one using this
excellent c++ library. 

I have a problem. I try to sort a 1*5 vector like

    GeneralArrayStorage<2> storage;
    storage.ordering() = firstRank;
    storage.base() = 0;
    storage.ascendingFlag() = true;
    Array<int,2> D(1,5,storage);
    D=4,6,8,2,3;
    cout<< "D = " <<D<<endl;

The output should be 2,3,4,6,8
But why doesn't it work?

Could I ask whether there is some function like sort() in matlab?

Thanks very much,

Best wishes,

Ying Wu

***************************************************************************************************************

Following the successful merger between the University of Paisley and Bell College, the University has been renamed as University of the West of Scotland.

Please note that as part of the next stage in the University's development, we have changed name to University of the West of Scotland. 

As a result, our web and email addresses have changed and I would be obliged if you would update your records accordingly. 

New University Web address: 

http://www.uws.ac.uk 

New University e-mail address format: 

firstname.surname@... 

The University of the West of Scotland is a registered Scottish charity.  Charity number SC002520.

***************************************************************************************************************

Legal disclaimer
--------------------------

The information transmitted is the property of the University of the West of Scotland and is intended only for the person or entity 
to which it is addressed and may contain confidential and/or privileged material.  Statements and opinions expressed in this e-mail may not represent those of the company.  Any review, retransmission, dissemination and other use of, or taking of any action in reliance upon, this information by persons or entities other than the intended recipient is prohibited.
If you received this in error, please contact the sender immediately and delete the material from any computer. 

--------------------------