| 1 | #include <iostream>
|
|---|
| 2 | using namespace std;
|
|---|
| 3 | #include <chrono>
|
|---|
| 4 | using namespace chrono;
|
|---|
| 5 |
|
|---|
| 6 | #include "caf/actor_ostream.hpp"
|
|---|
| 7 | #include "caf/caf_main.hpp"
|
|---|
| 8 | #include "caf/event_based_actor.hpp"
|
|---|
| 9 | #include "caf/all.hpp"
|
|---|
| 10 | using namespace caf;
|
|---|
| 11 |
|
|---|
| 12 | struct WorkMsg;
|
|---|
| 13 |
|
|---|
| 14 | CAF_BEGIN_TYPE_ID_BLOCK(custom_types_1, first_custom_type_id)
|
|---|
| 15 | CAF_ADD_TYPE_ID(custom_types_1, (WorkMsg))
|
|---|
| 16 | CAF_END_TYPE_ID_BLOCK(custom_types_1)
|
|---|
| 17 |
|
|---|
| 18 | // --(rst-foo-begin)--
|
|---|
| 19 | struct WorkMsg {
|
|---|
| 20 | unsigned int b;
|
|---|
| 21 | };
|
|---|
| 22 |
|
|---|
| 23 | template <class Inspector>
|
|---|
| 24 | bool inspect(Inspector& f, WorkMsg& x) {
|
|---|
| 25 | return f.object(x).fields(f.field("b", x.b));
|
|---|
| 26 | }
|
|---|
| 27 |
|
|---|
| 28 | unsigned int xr = 3'072, xc = 3'072, yc = 3'072, Processors = 1, MaxProcs = 48; // default values
|
|---|
| 29 |
|
|---|
| 30 | time_point<steady_clock> starttime;
|
|---|
| 31 | actor * actors;
|
|---|
| 32 | WorkMsg ** work_msgs;
|
|---|
| 33 | unsigned int actorCnt = 0;
|
|---|
| 34 |
|
|---|
| 35 | int wCount = 0;
|
|---|
| 36 |
|
|---|
| 37 | int ** w_Z, ** w_X, ** w_Y;
|
|---|
| 38 |
|
|---|
| 39 | class MatrixMult : public event_based_actor {
|
|---|
| 40 | behavior make_behavior() override {
|
|---|
| 41 | return {
|
|---|
| 42 | [=]( const WorkMsg & val ) -> void {
|
|---|
| 43 | unsigned int w = val.b;
|
|---|
| 44 | for ( unsigned int i = 0; i < yc; i += 1 ) { // multiply X_row by Y_col and sum products
|
|---|
| 45 | w_Z[w][i] = 0;
|
|---|
| 46 | for ( unsigned int j = 0; j < xc; j += 1 ) {
|
|---|
| 47 | w_Z[w][i] += w_X[w][j] * w_Y[j][i];
|
|---|
| 48 | } // for
|
|---|
| 49 | } // for
|
|---|
| 50 |
|
|---|
| 51 | if ( __atomic_add_fetch( &actorCnt, 1, __ATOMIC_SEQ_CST ) == xr ) {
|
|---|
| 52 | aout(this) << (steady_clock::now() - starttime).count() * ( MaxProcs / Processors ) / 1'000'000'000.0 << endl;
|
|---|
| 53 | } // if
|
|---|
| 54 | this->quit();
|
|---|
| 55 | return;
|
|---|
| 56 | }
|
|---|
| 57 | };
|
|---|
| 58 | }
|
|---|
| 59 | public:
|
|---|
| 60 | MatrixMult( caf::actor_config & cfg ) : event_based_actor( cfg ) {}
|
|---|
| 61 | }; // MatrixMult
|
|---|
| 62 |
|
|---|
| 63 | void caf_main( actor_system & sys ) {
|
|---|
| 64 | starttime = steady_clock::now();
|
|---|
| 65 |
|
|---|
| 66 | for ( unsigned int i = 0; i < xr; i += 1 ) { // create actors
|
|---|
| 67 | actors[i] = sys.spawn<MatrixMult>();
|
|---|
| 68 | } // for
|
|---|
| 69 |
|
|---|
| 70 | caf::scoped_actor self{sys};
|
|---|
| 71 | for ( unsigned int i = 0; i < xr; i += 1 ) { // start actors
|
|---|
| 72 | self->send( actors[i], WorkMsg{i} );
|
|---|
| 73 | } // for
|
|---|
| 74 | } // caf_main
|
|---|
| 75 |
|
|---|
| 76 | int main( int argc, char * argv[] ) {
|
|---|
| 77 | switch ( argc ) {
|
|---|
| 78 | case 6:
|
|---|
| 79 | if ( strcmp( argv[5], "d" ) != 0 ) { // default ?
|
|---|
| 80 | MaxProcs = stoi( argv[5] );
|
|---|
| 81 | if ( MaxProcs < 1 ) goto Usage;
|
|---|
| 82 | } // if
|
|---|
| 83 | case 5:
|
|---|
| 84 | if ( strcmp( argv[4], "d" ) != 0 ) { // default ?
|
|---|
| 85 | Processors = stoi( argv[4] );
|
|---|
| 86 | if ( Processors < 1 ) goto Usage;
|
|---|
| 87 | } // if
|
|---|
| 88 | case 4:
|
|---|
| 89 | if ( strcmp( argv[3], "d" ) != 0 ) { // default ?
|
|---|
| 90 | xr = stoi( argv[3] );
|
|---|
| 91 | if ( xr < 1 ) goto Usage;
|
|---|
| 92 | } // if
|
|---|
| 93 | case 3:
|
|---|
| 94 | if ( strcmp( argv[2], "d" ) != 0 ) { // default ?
|
|---|
| 95 | xc = stoi( argv[2] );
|
|---|
| 96 | if ( xc < 1 ) goto Usage;
|
|---|
| 97 | } // if
|
|---|
| 98 | case 2:
|
|---|
| 99 | if ( strcmp( argv[1], "d" ) != 0 ) { // default ?
|
|---|
| 100 | yc = stoi( argv[1] );
|
|---|
| 101 | if ( yc < 1 ) goto Usage;
|
|---|
| 102 | } // if
|
|---|
| 103 | case 1: // use defaults
|
|---|
| 104 | break;
|
|---|
| 105 | default:
|
|---|
| 106 | Usage:
|
|---|
| 107 | cerr << "Usage: " << argv[0]
|
|---|
| 108 | << " [ yc (> 0) | 'd' (default " << yc
|
|---|
| 109 | << ") ] [ xc (> 0) | 'd' (default " << xc
|
|---|
| 110 | << ") ] [ xr (> 0) | 'd' (default " << xr
|
|---|
| 111 | << ") ] [ processors (> 0) | 'd' (default " << Processors
|
|---|
| 112 | << ") ]" << endl;
|
|---|
| 113 | exit( EXIT_FAILURE );
|
|---|
| 114 | } // switch
|
|---|
| 115 |
|
|---|
| 116 | xr = xr / ( MaxProcs / Processors );
|
|---|
| 117 |
|
|---|
| 118 | //cout << Actors << " " << Set << " " << Rounds << " " << Processors << endl;
|
|---|
| 119 | actors = new actor[xr];
|
|---|
| 120 |
|
|---|
| 121 | caf::core::init_global_meta_objects();
|
|---|
| 122 | caf::exec_main_init_meta_objects<id_block::custom_types_1>();
|
|---|
| 123 |
|
|---|
| 124 | caf::actor_system_config cfg;
|
|---|
| 125 | cfg.set( "caf.scheduler.max-threads", Processors );
|
|---|
| 126 | caf::actor_system system { cfg };
|
|---|
| 127 |
|
|---|
| 128 | unsigned int r, c;
|
|---|
| 129 | int * Z[xr], * X[xr], * Y[xc];
|
|---|
| 130 |
|
|---|
| 131 | w_Z = Z;
|
|---|
| 132 | w_Y = Y;
|
|---|
| 133 | w_X = X;
|
|---|
| 134 |
|
|---|
| 135 | for ( r = 0; r < xr; r += 1 ) { // create/initialize X matrix
|
|---|
| 136 | X[r] = new int[xc];
|
|---|
| 137 | for ( c = 0; c < xc; c += 1 ) {
|
|---|
| 138 | X[r][c] = r * c % 37; // for timing
|
|---|
| 139 | } // for
|
|---|
| 140 | } // for
|
|---|
| 141 | for ( r = 0; r < xc; r += 1 ) { // create/initialize Y matrix
|
|---|
| 142 | Y[r] = new int[yc];
|
|---|
| 143 | for ( c = 0; c < yc; c += 1 ) {
|
|---|
| 144 | Y[r][c] = r * c % 37; // for timing
|
|---|
| 145 | } // for
|
|---|
| 146 | } // for
|
|---|
| 147 | for ( r = 0; r < xr; r += 1 ) { // create Z matrix
|
|---|
| 148 | Z[r] = new int[yc];
|
|---|
| 149 | } // for
|
|---|
| 150 |
|
|---|
| 151 | caf::exec_main<>(caf_main, argc, argv);
|
|---|
| 152 |
|
|---|
| 153 | for ( r = 0; r < xr; r += 1 ) { // deallocate X and Z matrices
|
|---|
| 154 | delete [] X[r];
|
|---|
| 155 | delete [] Z[r];
|
|---|
| 156 | } // for
|
|---|
| 157 | for ( r = 0; r < xc; r += 1 ) { // deallocate Y matrix
|
|---|
| 158 | delete [] Y[r];
|
|---|
| 159 | } // for
|
|---|
| 160 | } // main
|
|---|
| 161 | //CAF_MAIN()
|
|---|
| 162 |
|
|---|
| 163 | // /usr/bin/time -f "%Uu %Ss %Er %Mkb" a.out
|
|---|
| 164 |
|
|---|
| 165 | // Local Variables: //
|
|---|
| 166 | // compile-command: "g++-10 -Wall -O3 -std=c++17 -ICAF/actor-framework/libcaf_core -ICAF/actor-framework/libcaf_core/caf -ICAF/actor-framework/build/libcaf_core -LCAF/actor-framework/build/libcaf_core -LCAF/actor-framework/build/libcaf_io CAFMatrix.cpp -lcaf_io -lcaf_core -Wl,-rpath=CAF/actor-framework/build/libcaf_core" //
|
|---|
| 167 | // End: //
|
|---|
