2
/// This file is part of Rheolef.
4
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
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.
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.
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
20
/// =========================================================================
22
// draft-1 version of the geo_element class with arbitrarily order
23
// solution with std::vector inside each geo_element
26
// - souplesse : gere la memoire en tas dans geo (efficace)
27
// mais aussi peut creer des elements temporaires, si besoin
28
// - inconvenient : tableau de geo_element, qui contiennent des tableaux
29
// alloues dynamiquement => la memoire n'est pas contigue globalement
31
#include "rheolef/hack_array.h"
32
#include "rheolef/geo_element.h"
33
using namespace rheolef;
36
int main(int argc, char**argv)
38
environment env (argc, argv);
39
check_macro (communicator().size() == 2, "expect nproc = 2");
41
distributor ownership (dis_size, communicator(), distributor::decide);
42
check_macro (sizeof(geo_element_hack) == sizeof(geo_element::size_type), "portability problem");
44
geo_element::parameter_type param (reference_element::e, order); // edge, order 1
45
hack_array<geo_element_hack> ge_e (ownership,param);
46
geo_element_auto<> E (reference_element::e, order);
47
size_t jproc = (communicator().rank() == 0) ? 1 : 0;
48
if (communicator().rank() == 0) {
55
for (size_t dis_i = ownership.first_index(jproc), dis_n = ownership.last_index(jproc); dis_i < dis_n; dis_i++) {
56
ge_e.dis_entry(dis_i) = E;
58
ge_e.dis_entry_assembly();
59
ge_e.put_values (dout);