Threading Example

Example of Threaded Code
in GAUSS™ 9.0

Below is some sample code written first as a single thread, then divided into four threads.

// Procedure

proc mloop(x,N,its,pos);
  local m;

  m = 0;
  for i(1,its,1);
    m = m + moment(x[pos:pos+N-1,.],0);
    pos = pos + N;
  endfor;

  retp(m);
endp;

N = 4000;
x = rndn(N*N,4);

// Single-Threaded Calculation

mlp = mloop( x, N, N, 1 );

// Multi-Threaded Calculation (split into four threads for a quad-core machine)

R = rows(x)/4;
ThreadStat m1 = mloop(x, N, N/4, 1    );
ThreadStat m2 = mloop(x, N, N/4, 1+R  );
ThreadStat m3 = mloop(x, N, N/4, 1+R*2);
ThreadStat m4 = mloop(x, N, N/4, 1+R*3);
ThreadJoin;
mtmlp = m1 + m2 + m3 + m4;

Running time trials with these calculations on a quad-core machine, the multi-threaded calculation averaged 3.3 times faster than the single-threaded call!

Additional examples of how to incorporate threading into programs will be added periodically, so check back here for new examples.

GAUSS 9.0

GAUSS 9.0 Update Flyer [112k]