| 1 | //                               -*- Mode: C -*- 
 | 
|---|
| 2 | // 
 | 
|---|
| 3 | // Cforall Version 1.0.0 Copyright (C) 2017 University of Waterloo
 | 
|---|
| 4 | //
 | 
|---|
| 5 | // The contents of this file are covered under the licence agreement in the
 | 
|---|
| 6 | // file "LICENCE" distributed with Cforall.
 | 
|---|
| 7 | // 
 | 
|---|
| 8 | // matrixSum.c -- 
 | 
|---|
| 9 | // 
 | 
|---|
| 10 | // Author           : Peter A. Buhr
 | 
|---|
| 11 | // Created On       : Mon Oct  9 08:29:28 2017
 | 
|---|
| 12 | // Last Modified By : Peter A. Buhr
 | 
|---|
| 13 | // Last Modified On : Fri May 25 09:34:27 2018
 | 
|---|
| 14 | // Update Count     : 10
 | 
|---|
| 15 | // 
 | 
|---|
| 16 | 
 | 
|---|
| 17 | #include <fstream>
 | 
|---|
| 18 | #include <kernel>
 | 
|---|
| 19 | #include <thread>
 | 
|---|
| 20 | 
 | 
|---|
| 21 | thread Adder {
 | 
|---|
| 22 |         int * row, cols, & subtotal;                                            // communication
 | 
|---|
| 23 | };
 | 
|---|
| 24 | 
 | 
|---|
| 25 | void ?{}( Adder & adder, int row[], int cols, int & subtotal ) {
 | 
|---|
| 26 |         adder.[ row, cols ] = [ row, cols ];                            // expression disallowed in multi-member access
 | 
|---|
| 27 |         &adder.subtotal = &subtotal;
 | 
|---|
| 28 | }
 | 
|---|
| 29 | 
 | 
|---|
| 30 | void main( Adder & adder ) with( adder ) {                              // thread starts here
 | 
|---|
| 31 |         subtotal = 0;
 | 
|---|
| 32 |         for ( int c = 0; c < cols; c += 1 ) {
 | 
|---|
| 33 |                 subtotal += row[c];
 | 
|---|
| 34 |         } // for
 | 
|---|
| 35 | }
 | 
|---|
| 36 | 
 | 
|---|
| 37 | int main() {
 | 
|---|
| 38 |         const int rows = 10, cols = 1000;
 | 
|---|
| 39 |         int matrix[rows][cols], subtotals[rows], total = 0;
 | 
|---|
| 40 |         processor p;                                                                            // add kernel thread
 | 
|---|
| 41 | 
 | 
|---|
| 42 |         for ( int r = 0; r < rows; r += 1 ) {
 | 
|---|
| 43 |                 for ( int c = 0; c < cols; c += 1 ) {
 | 
|---|
| 44 |                         matrix[r][c] = 1;
 | 
|---|
| 45 |                 } // for
 | 
|---|
| 46 |         } // for
 | 
|---|
| 47 |         Adder * adders[rows];
 | 
|---|
| 48 |         for ( int r = 0; r < rows; r += 1 ) {                           // start threads to sum rows
 | 
|---|
| 49 |                 adders[r] = &(*malloc()){ matrix[r], cols, subtotals[r] };
 | 
|---|
| 50 | //              adders[r] = new( matrix[r], cols, &subtotals[r] );
 | 
|---|
| 51 |         } // for
 | 
|---|
| 52 |         for ( int r = 0; r < rows; r += 1 ) {                           // wait for threads to finish
 | 
|---|
| 53 |                 delete( adders[r] );
 | 
|---|
| 54 |                 total += subtotals[r];                                                  // total subtotals
 | 
|---|
| 55 |         } // for
 | 
|---|
| 56 |         sout | total | endl;
 | 
|---|
| 57 | }
 | 
|---|
| 58 | 
 | 
|---|
| 59 | // Local Variables: //
 | 
|---|
| 60 | // tab-width: 4 //
 | 
|---|
| 61 | // compile-command: "cfa matrixSum.c" //
 | 
|---|
| 62 | // End: //
 | 
|---|