root/trunk/blade/BladeDemo.d

//  Written in the D programming language 1.0.
/**
* BLADE Alpha -- Basic Linear Algebra D Expressions
*
*/

module BladeDemo;
import blade.Blade;
import std.stdio;

// Local arrays in D aren't currently aligned to 128-bit boundaries.
// In such cases, the library generates an 'SSE misalignment' assert error,
// to avoid segfaults.
// Use heap-allocated arrays, or static arrays (DMD 1.023 or later)
// cdouble[] always remains aligned, even when sliced.

void main()
{
    static z = [3.4, 565, 31.3, 0];
    double [] a = new double[4];
    double [] d = [0.5, 2.8, 17.0, 28.25, 1, 56.2, 3.4];
    a[0..$] = [3.4, 565, 31.3, 41.8];
    auto w2=z.ptr;
    float[] q = new float[4];
    q[0..$]= [17.0f, 28.25, 1, 0];
    float [4] r;
    real [] p = [2.3, 254, 0.1, 1.2];
    for(int i=0; i<r.length;++i) {
        r[i]= q[i]*2213.3L;
    }
    double [4][] another = [[33.1, 4543, 43, 878.7],
                            [5.14, 455, 554, 2.43]];
    real k=3.4;
    mixin(vectorize(` a += (d[2..$-1]*2.01*a[2]-another[][1])["abc".length-3..$]`));
    mixin(vectorize(" a-= 2.01*(        3.04+k)*r"));

    mixin(vectorize("q+= q*2.01"));
    mixin(vectorize("another[0]+=r-=another[0]+another[0]"));

    // All of the next four are equivalent
    mixin(vectorize("a+=6*another[1,0..$]"));
    mixin(vectorize("a+=6*(another[1,0..$]+another[1,0..$])"));


    mixin(vectorize("a+=6*another[1][0..$]"));
    mixin(vectorize("a+=6*another[1]"));
    mixin(vectorize("a+=6*another[1][]"));
    // I don't think I'll support this syntax long-term.
    mixin(vectorize("a+=6*another[1,[0,$]]"));

    // Strided vector
    mixin(vectorize("another[0..$,1]=6*a[0..2]"));

// Simplifies to q*= 2*dot(q,q)*dot(q,q).
   mixin(vectorize("q *=dot(q,q*dot(2*q,q))"));
   double u;
   mixin(vectorize("u = dot(q,q*dot(q,q))"));
   mixin(vectorize("u = dot(a, q)"));
   mixin(vectorize("u = sum(abs(-p))"));
   mixin(vectorize("a = -a"));
   mixin(vectorize("u = sum(sqrt(abs(p))) + sum(sqrt(abs(q)))"));
   mixin(vectorize("u = prod(q)"));
    writefln("a=", a);

}