← Back to branch summary

~ubuntu-branches/ubuntu/raring/rheolef/raring-proposed

« back to all changes in this revision

Viewing changes to nfem/basis/geo_element_v2.cc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito
  • Date: 2012-04-06 09:12:21 UTC
  • mfrom: (1.1.5)
  • Revision ID: package-import@ubuntu.com-20120406091221-m58me99p1nxqui49
Tags: 6.0-1
http://bugs.debian.org/667356
http://bugs.debian.org/661689
* New upstream release 6.0 (major changes):
  - massively distributed and parallel support
  - full FEM characteristic method (Lagrange-Gakerkin method) support
  - enhanced users documentation 
  - source code supports g++-4.7 (closes: #667356)
* debian/control: dependencies for MPI distributed solvers added
* debian/rules: build commands simplified
* debian/librheolef-dev.install: man1/* to man9/* added
* debian/changelog: package description rewritted (closes: #661689)

Show diffs side-by-side

added added

removed removed

Lines of Context:
nfem/basis/geo_element_v2.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
 
//
22
 
// the geo_element hierarchy of classes
23
 
//
24
 
// author: Pierre.Saramito@imag.fr
25
 
//
26
 
// date: 19 december 2010
27
 
//
28
 
#include "rheolef/geo_element_v2.h"
29
 
 
30
 
namespace rheolef {
31
 
 
32
 
char
33
 
skip_blancs_and_tabs (std::istream& is)
34
 
{
35
 
    if (!is.good()) return 0;
36
 
    char c = is.peek();
37
 
    while (is.good() && (c == ' ' || c == '\t')) {
38
 
        is >> c;
39
 
    } while (is.good() && (c == ' ' || c == '\t'));
40
 
    return c;
41
 
}
42
 
/// @brief return orientation and shift between *this element and S
43
 
/** assume that vertices of *this and S match
44
 
 * this routine is used for domain elements identification
45
 
 *  shift     : such that vertex (*this)[0] == S[shift] matches
46
 
 *  orient=+1 : when (*this)[1] == S[shift+1]
47
 
 * return true when all is ok (does not check fully that elements match)
48
 
 */
49
 
bool
50
 
geo_element::get_orientation_and_shift (
51
 
    const geo_element& S,
52
 
    orientation_type&  orient,
53
 
    size_type&         shift) const
54
 
{
55
 
    orient = 0;
56
 
    shift  = std::numeric_limits<size_type>::max();
57
 
    if (S.dimension() != dimension() || S.size() != size()) return false;
58
 
    size_type n = size();
59
 
    for (shift = 0; shift < n; shift++) {
60
 
      if (operator[](0) == S[shift]) break;
61
 
    }
62
 
    if (shift == n) {
63
 
      shift = std::numeric_limits<size_type>::max();
64
 
      return false;
65
 
    }
66
 
    if ((*this)[(shift+1)%n] == S[(shift+1)%n]) {
67
 
      orient = 1;
68
 
    } else {
69
 
      orient = -1;
70
 
    }
71
 
    return true;
72
 
}
73
 
 
74
 
} // namespace rheolef