← Back to branch summary

~ubuntu-branches/ubuntu/saucy/rheolef/saucy

« back to all changes in this revision

Viewing changes to skit/ptst2/array_of_idxset_scatter_tst.cc

  • Committer: Bazaar Package Importer
  • Author(s): Pierre Saramito
  • Date: 2011-03-23 11:14:26 UTC
  • mfrom: (1.1.4 upstream)
  • Revision ID: james.westby@ubuntu.com-20110323111426-cjvhey7lxt6077ty
Tags: 5.93-1
* New upstream release (minor changes):
  - some extra warning message deleted in heap_allocator
  - graphic output with mayavi2 fixed
  - add doc refman in .info and .pdf format

Show diffs side-by-side

added added

removed removed

Lines of Context:
skit/ptst2/array_of_idxset_scatter_tst.cc
 
1
///
 
2
/// This file is part of Rheolef.
 
3
///
 
4
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
 
5
///
 
6
/// Rheolef is free software; you can redistribute it and/or modify
 
7
/// it under the terms of the GNU General Public License as published by
 
8
/// the Free Software Foundation; either version 2 of the License, or
 
9
/// (at your option) any later version.
 
10
///
 
11
/// Rheolef is distributed in the hope that it will be useful,
 
12
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
/// GNU General Public License for more details.
 
15
///
 
16
/// You should have received a copy of the GNU General Public License
 
17
/// along with Rheolef; if not, write to the Free Software
 
18
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
19
///
 
20
/// =========================================================================
 
21
#include "rheolef/config.h"
 
22
#ifndef _RHEOLEF_HAVE_MPI
 
23
int main() { return 0; }
 
24
#else // _RHEOLEF_HAVE_MPI
 
25
#include "rheolef/distributed.h"
 
26
#include "rheolef/array.h"
 
27
#include "rheolef/index_set.h"
 
28
using namespace rheolef;
 
29
using namespace std;
 
30
 
 
31
int main(int argc, char**argv) {
 
32
    environment distributed(argc, argv);
 
33
    communicator comm;
 
34
    check_macro (comm.size() == 2 || comm.size() == 3, "expect np=2 or 3");
 
35
    size_t dis_n;
 
36
    dcin >> dis_n;
 
37
    array<index_set> x (dis_n);
 
38
    x.get_values (dcin);
 
39
    set<size_t> ext_idx_set;
 
40
    switch (comm.rank()) {
 
41
       case 0:
 
42
                ext_idx_set.insert (3);
 
43
                ext_idx_set.insert (5);
 
44
                break;
 
45
       case 1:
 
46
                ext_idx_set.insert (0);
 
47
                ext_idx_set.insert (1);
 
48
                break;
 
49
       case 2:
 
50
                ext_idx_set.insert (0);
 
51
                ext_idx_set.insert (3);
 
52
                break;
 
53
       default: error_macro ("unsupported np >= 4");
 
54
    }
 
55
    map <size_t,index_set> ext_idx_map;
 
56
    x.get_dis_entry (ext_idx_set, ext_idx_map);
 
57
    for (map <size_t,index_set>::const_iterator iter = ext_idx_map.begin(), last = ext_idx_map.end();
 
58
        iter != last; iter++) {
 
59
        cout << "proc " << comm.rank() << ": x[" << (*iter).first << "] = " << (*iter).second << endl;
 
60
    }
 
61
}
 
62
#endif // _RHEOLEF_HAVE_MPI